get80 will determine the best 80% landmarks. The output from this program is used to calculate FormU.
Must have run RESIDUALS!
usage: get80.sh [x] - (x may be used to limit the landmarks used to a specific GSD, ex. 0.0003 will only use the 30 cm to 39cm landmarks. Currently the 5th decimal place isn't used. If you give it something with 5 decimals it will hang.)
This script will produce files in your ~/ directory. Produces
- tmp80a - list of landmarks
- tmp80b - list of landmarks in order, along with residual values
- tmp80c - list of LMK files in order; best 80%
- tmp80d - list of landmark files in order; best 80%. Same as tmp80c, but without the .LMK at the end of each maplet file. (JRW - at least I think that's the only difference)
The output to the screen will be for the landmark at the 80% level (same line for that landmark in tmp80b). Per the example below, Col 5 is the valuable data you want from get80.sh. It indicates that 80% of the landmarks are better than that residual. If you want a more holistic number (i.e. FormU from Weirich et al. (2022)), take the RMS of the smallest GSD maplet values in Col 5 from tmp80b.
Example of tmp80b: Col 1: Landmark name; Col 2-4: Col 2 should be pointing, Col 3 should be limb, and Col 4 should be overlap, but these should be confirmed via the source code; Col 5: combination of pointing and limb errors that includes all maplets above and including the current line (use this value for internal consistency!); Col 6-7: Col 6 is a label for Col7, and Col 7 should be the combination of all 3 error, and should the same value found in RESIDUALS.TXT at the end of each landmark, but this should be confirmed via the source code; Col 8: Maplet resolution (stops at 4th decimal)
AR0001 0.10387 0.36778 0.10081 0.33276 all: 0.32265 Res:0.0003
GJ0004 0.57802 0.12566 0.07926 0.34191 all: 0.30685 Res:0.0003
AP0001 0.10221 0.39447 0.15224 0.34759 all: 0.33704 Res:0.0003
AG0002 0.09340 0.45815 0.16250 0.35798 all: 0.34324 Res:0.0003
FK0005 0.44616 0.13738 0.06904 0.36796 all: 0.34842 Res:0.0003
AC0001 0.09889 0.44376 0.16125 0.38225 all: 0.37034 Res:0.0003
AE0002 0.06820 0.48915 0.15824 0.39687 all: 0.38315 Res:0.0003
HR0001 0.60162 0.11211 0.07913 0.42346 all: 0.39629 Res:0.0003
AM0003 0.07107 0.51015 0.14339 0.42508 all: 0.41137 Res:0.0003
GK0002 1.21825 0.09294 0.09338 0.75481 all: 0.66459 Res:0.0003Note: The last landmark in RESIDUALS.TXT isn't used, so there must be an off by 1 error in one of the loops. Shouldn't be a big deal as long as you have many landmarks.