Skip to content

ReferencePeaks

Jump to bottom
Eric E Palmer edited this page Oct 23, 2026 · 1 revision

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

Calibration-Test11I

Spacecraft Azimuth 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

Sun Variations

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

Spacecraft Zenith Variations

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...)

Clone this wiki locally