= 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 ===

{{attachment:lesson4_autoreseduals.jpg||width=600}}

=== 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 ===
{{attachment:lesson4_mapinfo.jpg||width=600}}

=== PICINFO.TXT ===
{{attachment:lesson4_PICINFO.jpg||width=750}}

=== RESIDUALS.TXT ===
{{attachment:lesson4_RESIDUALS.jpg||width=600}}

=== A typical list of procedures when a new set of images comes in (DAWN Ops) ===
{{attachment:lesson4_pic5.jpg||width=600}}