Pds2isis errors/ mro gravity data advice

A transfer from the original Astrodiscuss forum (from sa-walton)

Hi all, I’m very new to using ISIS and trying to view MRO gravity data in ArcMap with a .img and .lbl file. I’m doing pds2isis and then trying to convert to a tiff (isis2std). I’m encountering an error during pds2isis:

PROGRAMMER ERROR Unable to find a translation value for [CoreBitsPerPixel, 64] in file [/home/sawalton/anaconda3/envs/isis/appdata/translations/pdsImage.trn].

I’m not quite sure how to fix this, or if anyone has any advice on using the gravity models in ArcMap I would really appreciate the help.

Thanks!

Unfortunately, using a 64bit floating point type in PDS3 is a little odd to see. Looks like ISIS doesn’t support it for ingestion. I think the max for ISIS3 is 32bit floating point. There is another issue we have seen with these labels by not setting the map_scale. So give this a try.

download both the label: https://pds-geosciences.wustl.edu/mro/mro-m-rss-5-sdp-v1/mrors_1xxx/data/rsdmap/ggmro_120_anom_100.lbl
and pixels: https://pds-geosciences.wustl.edu/mro/mro-m-rss-5-sdp-v1/mrors_1xxx/data/rsdmap/ggmro_120_anom_100.img

Edit the *.lbl, double check the math (the number of lines used below) but redefine MAP_SCALE from:
MAP_SCALE = “N/A”
to be:
radius (km) * PI / 2880 -->> 3396.0 * PI / 2880 -->>
MAP_SCALE = 3.70446133736 <km/px>

Now just load the updated ggmro_120_anom_100.lbl into ArcMap.

BTW, this file uses a Longitude range from 0 to 360, centered on 180. ArcMap can handle that, but it might have more issues or odd quirks as compared to using -180 to 180 longitudes.

To convert to GeoTIFF you can also use gdal_translate on the updated label too.

gdal_translate -of GTIFF ggmro_120_anom_100.lbl ggmro_120_anom_100.tif