make_script.in
Description
This input file lists the set landmarks or images intended for batch processing.
make_script.in is used by these script-makers:
make_scriptA - a script-maker for batch autoregistering new images
make_scriptAP - a script-maker for batch autoregistering new images using parallel processing
make_scriptF - a script-maker for batch lithos processing
make_scriptP - a script-maker for batch standard iteration of landmarks using parallel processing
make_scriptR - a script-maker for batch registration of new images
Here is a sample make_script.in file used by make_scriptA, make_scriptAP, and make_scriptR:
P00045000450 P00045000451 P00045000452 END
![/!\](/wiki/classic/img/alert.png "/!\")
You must precede each image filename with a space when creating make_script.in files for these script-makers.
Here is a sample make_script.in file used by make_scriptF and make_scriptP:
EE0001 EE0002 EE0003 END
- You must not precede each landmark filename with a space when creating make_script.in files for these script-makers.
(Compiled by DL)
make_scriptA.seed
Description
This seed file is used by make_scriptA, the script-maker for batch autoregistering new images.
Here is a sample make_scriptA.seed file annotated with explanations of the commands:
n <- auto-eliminate new landmarks 0,50,.25,0,5 <- INVLIM, SLIM, CLIM, RSMN, RSMX 1 <- auto-align 2 <- spacing 2 n <- no spacing change 0 <- shift by full amount y <- update landmark pixel locations 1 <- auto-align 1 <- spacing 1 n <- no spacing change 0 <- shift by full amount y <- update landmark pixel locations o <- remove low-correlation landmarks .5 <- correlation limit 0.5 1 <- auto-align 1 <- spacing 1 n <- no spacing change 0 <- shift by full amount y <- update landmark pixel locations 4 <- change flags b <- use all 0 <- exit (image) q <- quit autoregister END
(Compiled by DL)
make_scriptF.seed
Description
This seed file is used by make_scriptF, the script-maker for batch lithos processing. make_scriptR.seed contains the lithos commands for processing each landmark.
make_scriptF is a generic script maker for batch lithos processing tasks. You specify the set of landmarks to process in make_script.in and the sequence of commands to process in make_scriptF.seed.
make_scriptF generates the .INN files and run scripts to batch process the landmarks, and generates .OOT files and saves display files during processing for user inspection.
- You should test the command sequence in lithos before running the script generated by make_scriptF.
Writing the make_scriptF.seed
lithos expects a sequence of start and end commands which must be contained in the seed file before and after the processing commands.
Here is a sample make_scriptF.seed file. The required start and end command sequences are explained below.
# Seed file to attach map to overlapping maps #start commands n <- don't check for more images n <- don't include a single image #processing commands o RECENT y 1 o RECENT n 3 y 1, 3, 5 w i RECENT n n v 1 u 1 o RECENT n 1 #end command q <- quit #end-of-file identifier END <- end-of-file
Start Commands
The .INN file starts with the command i: input landmark, and the landmark name. lithos then asks you whether to check for more images. The following samples show the prompts and responses for both 'y' and 'n' responses to this prompt:
Check for more images? y[n] n Include a single image? y[n] n
Check for more images? y[n] y Enter fractional width (0=center). .5 Reject invisibles? y[n] n
When you create the make_scriptF.seed file, you must start with a sequence that deals with this set of options. In the sample, the Start Commands lines are annotated to indicate these responses.
End Command
make_scriptF expects you to quit lithos on completion of each landmark. The make_scriptF.seed file must contain the quit command q on the second-to-last-line.
End-Of-File Identifier
make_scriptF expects the seed file to end with the end-of-file identifier, END. Any commands following this will not be appended to the .INN file.
Comments
Use # at the beginning of any comment line you want to include. make_scriptF does not append these to the .INN file.
(Compiled by DL)
make_scriptP.seed
This seed file is used by make_scriptP for standard iteration of landmarks. make_scriptP.seed contains the lithosP commands for processing each maplet.
Here is an annotated sample make_scriptP.seed file:
# This script performs two maplet iterations in parallel mode. # It reinitializes slope and albedo. # It recomputes maplet overlaps. # It recomputes limb fits. # It recomputes central maplet vector. # Topography determination from slopes is conditioned by: # Initial topography, # Overlapping topography, # Shadowing, # Differential stereo, # Limb heghts. # Dataless area slopes are set from the shape model. # # The script hides most of the screen output. In order to obtain the # full display replace H* by #H*. # # Modified 2011_12_31. Header and comments added. # Modified 2014_02_26. Hide screen output added. # H* <- Toggle ON: suppress standard output x <- turn on extract filter .025 <- extract filter tolerance, 0.025 # # Initialize slope and albedo. # a <- reset albedo or slopes y <- reset slopes y <- reset albedos # # Begin first iteration. # 0 <- find template (build template #1) 40 <- number of iterations, 40 1 <- align landmarks 0 <- auto-align 1 <- auto-align spacing, 1 n <- don't try new spacing 0 <- continue (full shift) y <- update landmark pixel locations u <- update landmark files 0 <- find template, as above (build template #2) 40 # # Begin topo determination. # 2 <- find heights (build topography #1) 8 <- fill no data with shape slope 2 <- include existing heights .01 <- existing heights, fraction of points, 0.01 7 <- include shadowing .025 <- shadowing, fraction of points, 0.025 1 <- auto-include neighboring maps 4 <- include limb vectors 1, 2.5, 3 <- expansion (1), res/scale (2.5), d_hgt limit (3) 6 <- include differential stereo y <- use gradients / auto heights 0 <- begin integration .025 <- weight, 0.025 15 <- number of iterations 0 <- end iterations u <- update landmark files 1 <- align landmarks (as above) 0 2 <- auto-align spacing, 2 n 0 y 1 <- align landmarks (as above) 0 1 n 0 y u <- update landmark files # # Begin second iteration. # 0 <- find template (as above) (build template #3) 40 1 <- align landmarks (as above) 0 1 n 0 y u <- update landmark files 0 <- find template (as above) (build template #4) 40 # # Begin topo determination. # 2 <- find heights, as above (build topography #2) 8 2 .01 7 .025 1 4 1, 2.5, 3 6 y 0 .025 15 0 u <- update landmark files 1 <- align landmarks (as above) 0 2 n 0 y H* <- Toggle OFF: suppress standard output 1 <- align landmarks (as above) 0 1 n 0 y H* <- Toggle ON: suppress standard output u <- update landmark files # # Begin overlaps and limb fits. # o <- attach map to maps or limbs RECENT <- re-load most recent maplet y <- reset all 1 <- automatic o <- attach map to maps or limbs RECENT <- re-load most recent maplet n <- don't reset all 3 <- attach map to limbs y <- clear 1, 3, 5 <- expansion (1), res/scale limit (3), d_hgt limit (5) # # Central vector determination. # v <- find V from Z and pointing (mini-geometry) RECENT <- re-load most recent maplet q <- quit lithosP END
(Compiled by DL)
make_scriptR.seed
Description
This seed file is used by make_scriptR, the script-maker for batch registering new images.
Using make_scriptR.seed
The following annotated make_scriptR.seed file was used for registering 450 Vesta images to Zmaps:
m <- reference = map
0 <- Auto-choose Zmap
1.0 <- Scale = 1 km
a <- Turn off background
3 <- Shift image
y <- Autocorrelate
XSTOP <- Stop on no correlation
1 <- Change scale
0.25 <- New scale = 250 m
3 <- Shift image
y <- Autocorrelate
XSTOP <- Stop on no correlation
0 <- Quit
y <- Save new .SUM data
n <- Don't save rotation history
n <- Don't save nominal
q <- Quit procedure
ENDXSTOP - This is a provision for the autocorrelation procedure to end gracefully in the event that the correlation fails.
- If you enter
3 y XSTOP
and then the correlation fails, the program will reset the .SUM data to the input values and then stop.
XSKIP# - Used to to activate the autocorrelate sequence.
- If you enter
3 y XSKIP#
the program will skip # lines of the script, reset the .SUM data and ask you to input the next image.
The last two entries of 'n' are always so. The rotation history file was introduced to keep track of pointing errors in Clementine data during Lunar orbits in an attempt to quantify systematic shifts and it is rare that we want to set the nominal file equal to the .SUM solution.
(Compiled by DL)
Based on SPOC v3.02A PDF/LITHOSPHERE/REGISTER.f File Reference
make_scriptT.in
Description
This support file is used by make_scriptT, the script-maker for tiling a bigmap or shape.
Basic Tiling
Basic tiling involves tiling a bigmap which is located on the shape. For more details, refer to Basic_tiling.
Here is an annotated sample of a make_scriptT.in file used for basic tiling:
XXXXXX <- For maplet locations, reference to temporary file containing the tiling bigmap
scripts/t11-10.seed <- Reference to the make_scriptT.seed file
50 50 <- bigmap pixel/line location for placement of the center of each new maplet
100 50
150 50
200 50
50 100
100 100
150 100
200 100
50 150
100 150
150 150
200 150
50 200
100 200
150 200
200 200
ENDXXXXXX - You must save the tiling bigmap to a temporary file for tiling using a command line like this:
cp MAPFILES/BGMAP1.MAP MAPFILES/XXXXXX.MAP
- You then include the temporary filename as the first line of the file.
make_scriptT seed file - File containing the lithos commands for generating each new maplet. Refer to make_scriptT.seed for more information.
bigmap pixel/line locations - These are generated using the programs map_coverage and make_tilefile. Refer to those entries for details.
Lat/Long Tiling
Lat/Long tiling involves directly tiling a shape model. For more details, refer to LatLon_tiling.
Here is an annotated sample make_scriptT.in file used for lat/long tiling:
LATLON <- For maplet locations, reference to shape model SCRIPTS/XXX150.SEED <- Reference to the make_scriptT.seed file 00, 000 <- shape model lat/wlong for placement of the center of each new maplet 00, 012 00, 024 00, 036 00, 048 00, 060 00, 072 00, 084 .... -84, 315 90, 000 -90, 000 END
make_scriptT seed file - File containing the lithos commands for generating each new maplet. Refer to make_scriptT.seed for more information.
Image Tiling
- The 1st line must be "IMAGES".
- The 2nd line is the script name
- Following lines are image name, sample and pixel
- Image name must be char 1-12. Sample and pixel are after that with any fortran format (space or comma delimited)
- Ends with the statement "END"
- # symbol will cause that entry to be skipped
Here is an annotated sample make_scriptT.in file used for tiling an image:
IMAGES <- For maplet locations, reference to shape model TIL_eekingLimbsInImage.seed <- Reference to the make_scriptT.seed file D7175061248G 358 511 <- image name followed up pixel/line location D7175061210G 372 515 D7175061210G 387 513 D7175061210G 403 507 .... D7175061210G 566 887 END
(Compiled by DL)
make_scriptT.seed
Description
This seed file is used by make_scriptT, the script-maker for tiling a bigmap or shape with a suite of new maplets. make_scriptT.seed contains the lithos commands for processing each new maplet.
Here is an annotated sample make_scriptT.seed file:
0.00010,49 <- maplet ground sample distance and Q-size g <- turn on auto-rename (based on its shape model lat/wlon region) i <- input landmark a <- new landmark temporary name 'a', will immediately be auto re-named y <- check for more images .5 <- fractional width, 0.5 n <- don't reject invisibles x <- turn on extract filter .025 <- extract filter tolerance, 0.025 0 <- find template a <- use all images b <- get heights from surrounding map n <- don't list all possibilities XXXXXX <- reference map name (temporary file which contains tiling bigmap) n <- find normal 2 <- find normal from heights b <- get heights from surrounding map (as above) n XXXXXX u <- update landmark files 1 <- save map v <- find V, Z, or PTG 2 <- find Z from V and pointing e <- eliminate images from landmark a <- auto-remove 0,60,.25,.25,0,3 <- INVLIM, SLIM, CLIM, RSMN, RSMX 1 <- align landmarks 0 <- auto-align 3 <- auto-align spacing, 3 n <- don't try new spacing 0 <- continue (full shift) y <- update landmark pixel locations 1 <- align landmarks (as above) 0 1 <- auto-align spacing, 1 n 0 y 0 <- find template c <- use correlated images .5 <- correlation limit, 0.5 0 <- find template (build template #1) 0 <- continue 40 <- number of iterations, 40 1 <- align landmarks (as above) 0 3 n 0 y 1 <- align landmarks (as above) 0 1 n 0 y e <- eliminate images from landmark o <- remove low correlating images .5 <- correlation limit, 0.5 0 <- find template a <- use all images u <- update landmark files (as above) 1 0 <- find template (as above) (build template #2) 0 40 <- number of iterations, 40 2 <- find heights (build topography #1) 8 <- fill with shape 2 <- include existing heights .01 <- existing heights, fraction of points, 0.01 1 <- auto-include neighboring maps 4 <- include limb vectors 1, 2.5, 3 <- expansion (1), res/scale (2.5), d_hgt limit (3) 6 <- include differential stereo y <- use gradients y <- use auto heights 0 <- begin iteration .025 <- weight, 0.025 30 <- number of iterations 0 <- end iterations u <- update landmark files (as above) 1 v <- find V, Z, or PTG 1 <- find V from Z and pointing 1 <- align landmarks (as above) 0 1 n 0 y u <- update landmark files (as above) 1 0 <- find template (as above) (build template #3) 0 40 1 <- align landmarks (as above) 0 2 n 0 y 1 <- align landmarks (as above) 0 1 n 0 y e <- eliminate images from landmark (as above) o .5 u <- update landmark files (as above) 1 0 <- find template (as above) (build template #4) 0 40 2 <- find heights (build topography #2) 8 2 .01 1 4 1, 2.5, 3 6 y y 0 .025 30 0 u <- update landmark files (as above) 1 v <- find V, Z, or PTG (as above) 1 1 <- align landmarks (as above) 0 1 n 0 y u <- update landmark files (as above) 1 o <- attach map to maps or limbs RECENT <- re-load most recent maplet y <- reset all 1 <- automatic o <- attach map to maps or limbs RECENT <- re-load most recent maplet n <- don't reset all 3 <- attach map to limbs y <- clear 1, 3, 5 <- expansion (1), res/scale limit (3), d_hgt limit (5) i <- input landmark RECENT <- re-load most recent maplet n <- don't check for more images n <- don't include single image v <- find V, Z, or PTG (as above) 1 u <- update landmark files (as above) 1 v <- find V, Z, or PTG 4 <- find PTG/SCOBJ from Z and V (current images) - update S/C position and pointing o <- attach map to maps or limbs RECENT <- re-load most recent maplet n <- don't reset all 1 <- automatic q <- quit lithos END
(Compiled by DL)
Bibliography
Gaskell R.W., Landmark navigation and target characterization in a simulated Itokawa encounter, AAS paper 05-289, AAS/AIAA Astrodynamics Specialists Conf., Lake Tahoe, CA, 2005. [GASKELL2005]
Gaskell, R.W., O. Barnouin-Jha, and D. Scheeres, Modeling Eros with Stereophotoclinometry, Abstract #1333, 38th LPSC, Houston, TX, 2007. [GASKELLETAL2007]
Gaskell, R.W., O S Barnouin-Jha, D J Scheeres, A S Konopliv, T Mukai, S Abe, J Saito, M Ishiguro, T Kubota, T Hashimoto, J Kawaguchi, M Yoshikawa, K Shirakawa, T Kominato, N Hirata, H Demura (2008). Characterizing and navigating small bodies with imaging data. Meteoritics and Planetary Science 43, 1049-1062.
Kirk, R. L., Barrett, J. M., and Soderblom, L. A., 2003, Photoclinometry made simple…?, ISPRS Working Group IV/9 Workshop "Advances in Planetary Mapping 2003", Houston, March 2003, online at http://astrogeology.usgs.gov/Projects/ISPRS/Meetings/Houston2003/abstracts/Kirk_isprs_mar03.pdf.