Lesson 4

Objectives

• AUTOREGISTER, Create Bigmap from 0.4 pm/px maps, set up batch run tiling (3-4 hours)
  • CONVERT last set of images
  • AUTOREGISTER
  • Create bigmap at 0.4 km/px resolution
  • Set up batch run
    • Map_coverage
    • Make_tilefile
    • Make_scriptT
    • Find_nofitT
  • Iterate using GEOMETRY, RESIDUALS, and LITHOS

Process last set of images

• cd ../DAWN/NEW_IMAGES_2011_246

• ls -1 > $myvesta/newpix.txt

• cd $myvesta
• ../DAWN/UTILITIES/PROCESS_IMG.e

Make Sumfiles

• ../DAWN/UTILITIES/make_sumfiles.e

• y - use B-frame? (y/n)

  • 0. Create SUMFILES & NOMINALS and ..

  • 1. Update old NOMINALS and ..
  • 2. Update old SUMFILES

• 1

DYNAMICS

• emacs DYNAMICS.TXT

• emacs PICTLIST.TXT

• FC21A0006742 through the end from PICTLIST.TXT into ‘SURVEY to HAMO’ in DYNAMICS.TXT

• DYAMICS.e

ETLM=  1800, 1.D-6      APPROACH
VSIG= 0.200, 0.200, 0.200
PSIG= 0.0001, 0.0001, 0.0001
 FC21A0002099
 FC21A0002101
…
FC21A0003895
 FC21A0003910
ETLM=  1800, 1.D-6      SURVEY
VSIG= 0.100, 0.100, 0.040
PSIG= 0.0001, 0.0001, 0.0001
 FC21A0003931
 FC21A0003932
 …
FC21A0006731
 FC21A0006732
ETLM=  1800, 1.D-6      SURVEY TO HAMO
VSIG= 0.200, 0.200, 0.200
PSIG= 0.0001, 0.0001, 0.0001
FC21A0006742
...
FC21A0006781
END

AUTOREGISTER

• AUTOREGISTER.e

• FC21A0006781 - input 12-character picture name. 1 to quit.

• y - add landmarks? (y/n)

• 0.5 - enter fractional width (0 = center)

• n - reject invisibles? (y/n)

• display LMRK_DISPLAY1.pgm

• 1 - auto align

• 1 - enter spacing

• n - new spacing

• 0 - continue

• y - update landmark pixel locations?

• 1 - auto align

• 1 - enter spacing

• n - new spacing?

• 0 - continue

• y - update landmark pixel locations?

• 4 - change flags

• b - use all

• 0 - exit

• q - quit

AUTOREGISTER.e

BIGMAP

• BIGMAP

• m - pixlinpic (p) or ltding (l) or map p/l (m)

• LMK000 - input 6 character map name

• 50,50 - input patch center (map p/l)

• .5,200,1.234 - scale, qsz, seed

• BM0001 - input 6 character map name

• 1 - integrate slopes

• .005 - input fraction

• .025 - input weight

• 1 - more iteration

• 1 - more iteration

• 1 - more iteration

• 1 - more iteration

• 1 - more iteration

• 1 - more iteration

• 1 - more iteration

• 0 - end iteration

• 0 - no template

• SHOWMAP.e

• BM0001 - input map name

• display BM0001.pgm

=== Batch Run Set-Up ===

1. run map_coverage to find blank spots in BM0001

   • map_coverage.e

   • BM0001

   • .4,.4

   • display coverage_m.pgm &

2. Then run make_tilefile.e which takes coverage_m.pgm and determines what you have left to tile

   • make_tilefile.e

   • make_tilefile.e>make_scriptT.in

Batch Run Set-up (cont.)

• Open make_scriptT.in in emacs and add the following to the top of the file
  • BM0001
  • BM0001.SEED

• cp USED_MAPS.TXT LMRKLISTO.TXT

• emacs BM0001.SEED

• ./rem_script.b

• rm -f ./TESTFILES/*

• rm -f ./TESTFILES1/*

• make_scriptT.e

• ls -lrt

• emacs run_script.b

Check batch run before letting it go

• Cut and paste the first 4 lines in terminal
• When it’s finished, look at:

  • cat 001.OOT

  • OR

  • find_nofitT.e

• If it looks pretty good, deleted first 4 lines for 001 (since we already executed it) and

• ./run_script.b

• When it’s all done, run find_nofitT.e

Geometry

• Now let’s solve for the geometry and residuals to check for outliers

• geometry

  • 1: landmark vectors.
  • 2: camera pointing, scobj
  • 3: RA/DEC/TWIST
  • 0: end.

• 120 - input operation list

• 3 - enter number of iterations

• n - change position and pointing aprioris? (y/n)

• n - use limbs for pointing? (y/n)

• n - continue? (y/n)

Residuals

• After running GEOMETRY, we can solve for the residuals
• Check for landmarks with poor residuals

  • possible misregistration => manually adjust

• RESIDUALS.e

• 2,5,10 - enter plim (px,km,km)

• Then look at MAPINFO.TXT, PICINFO.TXT, RESIDUALS.TXT

MAPINFO.TXT