Test 11: ReferencePeaks
Aim and Objectives
Purpose of test: To evaluate SPC's ability to represent a simplified 'boulder' (a peak) which has a constant 20 degree slope.
Additional Objectives
- To investigate the effect of varying the combinations of viewing conditions which comprise the image suite, specifically:
- Variation of spacecraft azimuth;
- Variation of sun zenith (incidence);
- Variation of spacecraft zenith (emission angle);
- To investigate the effect of the number of images and viewing conditions which comprise the image suite.
Methodology
The 'boulder' was simulated by generating a simple axi-symmetric, flat-topped cone (peak) with a constant 20 degree slope. Two boulder simulations were generated:
- A single peak in an otherwise flat landscape of width 5m - Modeled using a single maplet (no overlaps).
- Multiple peaks and (simplified) craters in an otherwise flat landscape of width approx. 22.5m.
Data
- Data was generated by ??.
- Image resolutions: 2cm/pixel.
Planned Tests
Please follow embedded links for test results.
Take PIT020 as the only maplet in MAPLIST.TXT. Generate 5cm images of the topography (which is based at TAG Site 1). Then we build a MAPLET, saving the MAPLETs flatMapVect at several states.
Simple Variations
Groups of images - (adjusted for 5cm)
- Test11A - Detailed Survey. Eric (Has 11 images, very low in sample space)
- Test11B - Min Orbit B (by hand, calculate the minimum requirement)
- Test11C - Average Orbit B (use a random Orbit B solution with 3 unique observation/illumination conditions)
- Test11D - Best Orbit B (as we did before)
- Test11E - Optimum Collection (as we did before)
Name |
Sun Az |
Sun Zenith |
S/C Az |
S/C Zenith |
Optimum Collection1 |
0, 90, 180, 270 |
45 |
every 20 |
45 |
Optimum Collection2 |
135 |
30 |
every 20 |
45 |
- Test11X - Criss Cross (combine both of these sets)
Name |
Sun |
S/C Az |
S/C Zenith |
Criss Cross 1 |
80, 60, 40, 20, 0, -20, -40, -60, -80 |
45 |
45 |
Criss Cross 2 |
45 |
0, 180 |
80, 60, 40, 20, 0 |
Variations
Name |
Sun |
S/C Az |
S/C Zenith |
Total Number of Images |
Minimum CompareOBJ RMS (cm) |
Minimum Formal Uncertainty (cm) |
Test11F |
45, 0, -45 |
every 20 |
45 |
72 |
0.704 |
2.653 |
Test11G |
45, 0, -45 |
every 30 |
45 |
48 |
0.820 |
1.532 |
Test11H |
45, 0, -45 |
every 45 |
45 |
32 |
0.878 |
1.394 |
Test11I |
45, 0, -45 |
every 60 |
45 |
24 |
0.951 |
1.370 |
Test11J |
45, 0, -45 |
every 90 |
45 |
16 |
1.069 |
1.392 |
Name |
Sun |
S/C Az |
S/C Zenith |
Total Number of Images |
Minimum CompareOBJ RMS (cm) |
Minimum Formal Uncertainty (cm) |
Test11K |
30 |
every 45 |
45 |
8 |
|
|
Test11L |
30 -30 |
every 45 |
45 |
16 |
|
|
Test11M |
-45, -45 |
every 45 |
45 |
16 |
|
|
Test11N |
30, 0, -30 |
every 45 |
45 |
32 |
|
|
Test11P |
45, 30, 0, -30, -45 |
every 45 |
45 |
48 |
|
|
Test11Q |
60, 45, 30, 0, -30, -45, -60 |
every 45 |
45 |
64 |
|
|
Test11R |
60, 30, -30, -60 |
every 45 |
45 |
32 |
|
|
Name |
Sun |
S/C Az |
S/C Zenith |
Total Number of Images |
Minimum CompareOBJ RMS (cm) |
Minimum Formal Uncertainty (cm) |
Test11S |
45, 0, -45 |
every 45 |
0, 30, 60 |
80 |
1.035 |
1.413 |
Test11T |
45, 0, -45 |
every 45 |
0, 20, 40, 60 |
112 |
0.874 |
1.522 |
Test11U |
45, 0, -45 |
every 45 |
0, 15, 30, 45, 60 |
144 |
0.905 |
1.598 |
Test Set-Up
Hill Data
I have created the following hills for testing
- REGION - Includes lots of craters and a few 20 deg peaks
Truth Maplets
- -8.027 262.768
- T11TRU.MAP - already built. It has the truth heights loaded from PIT020.
- T11TRU.TXT -- Text file of the maplet. Use for plotting transits
- t11-truth.obj -- OBJ file for testing using CompareOBJ
Testing Maplets
- -8.027 262.768
- Build a model MAPLET using TAG Site 1 center position (5cm resolution and Q=49). Do not load a base map for topography
- Create MOD20#, which # is the test number (letter)
- Align the images with the key location at 50/50.
- Build template with 100
- (Skip until close to perfect) -- Do an auto align
- Build height using the iteration template
- Do flatMapVec and rename it to mod-20#-1
- Do compMapVec between MOD20# and T11TRU
- Repeat as necessary
Iteration Template
0 0 44 2 8 fill no data with shape slope 7 .01 6 y y 1 condition with overlapping maps 2 condition with existing topo .01 6 y y 0 begin integration .025 100 0 v 1 u 1 o overlaps RECENT n 1 i r n n
Useful copy/paste stuff
echo CENTER | /usr/local/localBin/flatMapVec /usr/local/localBin/compMapVec MAPFILES/CENTER MAPFILES/T11TRU.MAP ~/eval/bin/Maplet2FITS MAPFILES/CENTER fit ~/eval/bin/FITS2OBJ --local fit center.obj ~/eval/bin/CompareOBJ --compute-optimal-translation-and-rotation center.obj newtru.obj convert LMRK_DISPLAY1.pgm display.jpg mkdir step02 cp CENTER.TXT MAPFILES/CENTER.MAP step02 gnuplot set term x11 set style data lines plot 'NEW020.TXT' u 0:50 lw 3 lc 0, 'CENTER.TXT' u 0:50 set output 'out.png' set term png replot unset output set term x11 exit mv out.png mod20f.obj MOD20F.TXT comp.txt display.jpg step2
Save States
The goal is to see the evolution of the model over time. You can stop iterations once the peak and base are stable. The most likely tool to evaluate if the model has reached it best fit will be transit.
Save the flatMapVect statistics at the following number of iterations
- 1st
- 2nd
- 3rd
- 5th
- 10th
- 15th
- 20th
- 50th (if needed)
- 100th (if needed)
- 200th (if needed)
Examples
Example of a 14 degree peak, measured at difference locations 
Example of a 20 degree peak, measured over time 
ReferencePeaks_dep Deprecated Reference Peaks (If at 1st you dont succeed...)