Help needed: absolute triangulation HiRISE SOCET GXP

[marlens123]

Help needed: absolute triangulation HiRISE SOCET GXP

[thareUSGS]

For GXP questions, the best methods for questions is via our email "planetaryphotogrammetry at usgs.gov" or direct emails. There is also the Open Planetary forum: https://www.openplanetary.org/

Anyway, MOLA topography can be extracted from the file above for GXP. While map projected GeoTIFFs seem to work, we still often import the clipped area using an Arc ASCII grid format (created using GDAL) in degrees. Likely the easiest method to clip from the MOLA is visually in QGIS or using gdal "-projwin" parameters. GDAL can be a little tricky if you want to clip a file in a projected meters using Long/Lat bounds. But you can do that using gdal_translate (projwin and projwin_srs) or gdalwarp (te and te_srs).

e.g.;

gdal_translate -projwin 120 45 115 40 -projwin_srs IAU_2015:49900 MOLA_in_meters_projection.tif out_clip.tif

• Here you are clipping by a lon/lat range defined by the projwin_srs set to Mars in degrees. IAU_2015:49900.

or (exporting to ASCII for GXP ASCII terrain import using degrees)

gdal_warp -te 115 40 120 45 -te_srs IAU_2015:49900 -of AAIGrid MOLA_in_meters_projection.tif out_clip_degrees.asc

• note the order of "-te" is different from gdal_translate's "-projwin"

To go from elevations in "topographic" heights to "geoid" heights (optional but recommended) you can run ASP's dem_geoid.