Differences between revisions 36 and 37
 ⇤ ← Revision 36 as of 2016-07-26 09:43:02 → Size: 8982 Editor: BMittan Comment: ← Revision 37 as of 2022-03-14 13:37:31 → ⇥ Size: 8983 Editor: JohnWeirich Comment: Deletions are marked like this. Additions are marked like this. Line 68: Line 68: '''8. Fill with shape''' - This option uses the shape model to fill any "holes" in the maplet, if they exit. '''8. Fill with shape''' - This option uses the shape model to fill any "holes" in the maplet, if they exist.

# Find heights

(lithos Main Menu Option 2. Find heights)

## Description

This Main Menu option solves for the brightness and slope of each pixel. It includes a variety of options, including stereo and photoclinometry.

### Suboptions

When you choose Option 2. Find heights from the Main Menu, you will see this set of Suboptions:

```Main Menu: Q 0123O IUCRSG NVAMBX DEPL FZH? 2
>2

TOOLS
0. Begin integration
1. Auto include neighboring maps
2. Include existing heights
3. Include shape vectors
4. Include limb vectors
5. Include external map
6. Include differential stereo
8. Fill with shape
9. Include control points
Q. Exit```

0. Begin integration - Provide these inputs:

1. Weight for the conditioning heights - Always use .025.

Whether to use weight 0.025 or not depends on a number of factors. Often in the first slope to height integration you may want to use a smaller number, 0.01. Also if you have a large number of overlapping maps and the seeds.pgm plot looks very dense, you may want to decrease the weight so as not to over constrain the solution with too many seed heights. Also, if some of the conditioning heights are derived from maps that have poorly determined heights, you should decrease the value of the weight to even less than 0.01. If this does not improve the heights integration, you may skip options 1 and 2 and use option 6 before you perform the heights integration.

1. Number of iterations - Typically use 30.

2. After the iterations are finished, enter one of these choices:
• '0' to go back to the Main Menu
• '1' to input more iterations
• '2' to change the weight

1. Auto include neighboring maps - Includes data from surrounding maplets that have been added to the "overlap" list.

• If you are building a maplet for the first time, there won't be any overlaps assigned, so this option won't be needed until later.

2. Include existing heights - Use this option most times. You dn't have to use this option if you include an external map.

• This option uses heights from the existing topography as conditioning heights, with the fraction input by the user (typically 1%). Therefore, this option weights data from the current maplet as well, and that way the changes will occur slowly and will remain stable.

3. Include shape vectors - This option is seldom used, perhaps only in the beginning stage of creating the first set of landmarks.

4. Include limb vectors - Constrains the solution by using limb data.

• This option prompts you to enter expansion, res/scale limit, d_hgt limit. Entries can be separated by 'Enter', ' ', or ','. "Expansion" is included only for backwards compatibility, "res/scale limit" refers to the resolution of the pictures that it searches for limbs, and "d_hgt limit" refers to the distance away (in number of pixels) it will search for limbs. d_hgt limit is similar to the Enter spacing prompt from Align landmarks (i.e., you need to enter a number between 1 and 5). Typically, use the values 1, 2.5, 3 for these entries.

5. Include external map - You can use an external map, often a big map, to provide seed heights for the slope to heights integration.

• You have to input a map name and a fraction of points just like when using editing heights.

If the input to the "fraction of points" is < 0.0001 then instead of seed heights use slopes directly derived from the external map. Doing so means that you can skip options 1 and 2 and go directly to option 0 integration.

6. Include differential stereo - This option calculates stereo points at random locations on the maplet. Testing has shown that this option is a good way to keep maplets from tilting at the edges.

• After you enter this option, respond to these prompts:
• use gradients - enter 'y'

• autoheights - enter 'y'

7. Include shadowing - This is still an experimental procedure. It will be tried thoroughly when there's significant shadowing in the images used for the map in question.

8. Fill with shape - This option uses the shape model to fill any "holes" in the maplet, if they exist.

9. Include control points - This option is rarely used. It will use landmarks as seed heights.

Q. Exit - Exit to the Main Menu.

## Typical Example

The following sample shows how to use 2. Find heights when starting at the Main Menu.

• Sometimes options 6 or 7 (or both) will be removed, but the others are typically present.

```2         <- Find heights option
8         <- Fill w/ shape
2         <- Include existing heights
.01       <- weighting
1         <- Include overlaps
4         <- Include limbs
1, 2.5, 3 <- inputs for limbs
6         <- include stereo
y         <- use auto heights
0         <- Begin integration
.025      <- weighting
30        <- number of iterations
0         <- End iterations```

The following sample output begins with the '0' input for Begin integration:

```0
0
gc seeds.pgm

input weight
.025
Enter number of iterations
30
-0.33759     0.34356    -0.34050     0.41873

gc slope.pgm

-0.28767     0.27729    -0.25837     0.36825

gc slope.pgm

...
-0.22599     0.18857    -0.16774     0.26969

gc slope.pgm

-0.22600     0.18859    -0.16777     0.26958

gc slope.pgm

0. end iterations
1. more iterations
2. change weight
0
-8.18594598952327801E-006

Current landmark = EE0001
SCALE =  0.000100 QSZ =   49
Region = EE

gc LMRK_DISPLAY0.pgm
gc LMRK_DISPLAY1.pgm

check display
Main Menu: Q 0123O IUCRSG NVAMBX DEPL FZH?```

1. Auto include neighboring maps

• [From LITHOS.f] - This provides conditioning heights from overlapping maplets. Points on the overlapping maplet are chosen randomly according to a fixed formula that provides more points (up to 5%) if the overlapping maplet is higher resolution than the target maplet, 2.5% if resolutions are the same and dropping to zero at low resolution overlap. The projection is done by the subroutine GET_MAP (see abstract).

3. Include shape vectors

• [From LITHOS.f] - This option provides conditioning heights from the existing shape model, with a fraction set by the user and heights from the subroutine GET_HEIGHTS. Option 3 is almost never used.

4. Include limb vectors

• [From LIMB_HEIGHTS.f] - There are three user input search parameters. ZEXP is an expansion parameter that is always set to 1 but is kept for backward compatibility with old scripts. ZRES restricts the image RESOLUTION (in km/px) relative to the maplet SCALE: RESOLUTION < ZRES*SCALE. ZHGT is the maximum allowed difference between predicted limb point height on a maplet and the measured height in units of maplet SCALE.

5. Include external map

• [From LITHOS.f] - This option gets conditioning heights from an externally produced map, with a fraction set by the user and heights from the subroutine GET_MAP (see abstract). If the fraction is set to zero, the option fills missing data areas with slopes from the external map, in much the same way as option 8, but with less granularity. If used in this context, option 5 should always be used before option 8. Note that this option may be used twice, once to provide conditioning heights and once to fill in missing slopes. Option 5 is the natural interface between SPC and other sources of topographic data such as SPG or scanning lidar.

6. Include differential stereo

• After you enter this option, respond to these prompts:
• use gradients - enter 'y'

• autoheights - enter 'y'

• [From LITHOS.f] - This option determines conditioning heights from differential stereo. Features in a maplet move on trajectories as their heights change, owing to the obliqueness of its appearance in an image. The subroutine STEREO (see abstract) uses differences in these trajectories from image to image to determine conditioning heights.