Instructions to convert a MAP file to other formats
Convert to FITS
- Use the AltWG toolset to get the proper format
Maplet2FITS --configFile fitsConfig_orex_spc.ini RAD201.MAP RAD201-fullv1.fits
An example ini file can be found here fitsConfig_spc.ini.
- Notes from Ray Espiritu on 5 May 2020 about the ini file "It may look familiar as it has the same format as the nft config files you send me for generating MLNs. Format is KEYWORD#value. Notice that “//” denotes comment lines. This does not allow you to add custom keywords; the keyword has to exactly match a keyword that is normally populated in the fits header."
There will still be some header information you will want to change. For this use ChangeFitsHeader, also part of the AltWG distribution toolset.
More notes from Ray Espiritu on 5 May 2020 for "This tool allows you to add,delete, or modify keywords in the header. Run “ChangeFitsHeader –examples” to see examples for each operation (add, delete, modify, change). There is a limitation inherent in the fits format in that the new fits header cannot have a size difference of more than 2880 bytes from the old fits header. You shouldn’t encounter that if you are adding or deleting a handful of keywords."
Convert to GeoTIFF
- First, convert to FITS format (and update headers as necessary).
- Use GDAL to convert from FITS to GeoTIFF
gdal_translate -b 3 -a_srs RAD201-v5.prj -a_ullr -8510 8510 8510 -8510 RAD201-fullv1.fits RAD201-fullv1.tif
An example of the prj file can be found here, and can be viewed with a text editor such as "vi". RAD201-v5.prj
The -a_ullr is giving the bounds of the DTM in meters. The math is "Upperleft X = (Q*2+1) * GSD (in meters) / 2 * -1" and "UpperLeft Y = (Q*2+1) * GSD (in meters) / 2". So if the GSD is 20 m with a Q of 425, the math is (425*2+1) * 20 / 2 * -1 = -8510
- The -a_scale option will tell the geoTiff how to convert the DN value to physical units. I am still working out how to determine the value needed.
If you are converting a fits file with multiple bands, you can specify a particular band with the -b <band #>
gdal_translate -a_scale .0023 -a_srs INGV02_jrw.prj -a_ullr -4552.5 4552.5 4552.5 -4552.5 INGV02_DTM.pgm INGV02.tif