How to: GDAL with Martian Data

I’ve seen a lot of tutorial on working with GDAL with Earth data, but not a lot for planetary data, Mars included. As I use GDAL a lot (and in need of reminder myself), why not write a how-to with Martian data.

MOLA, but in different projections (https://core2.gsfc.nasa.gov/PGDA/images/MarsTopography_GMM3.jpeg, https://astrogeology.usgs.gov/search/map/Mars/GlobalSurveyor/MOLA/Mars_MGS_MOLA_ClrShade_merge_global_463m)

There is a reason why there is not a lot of discussion about GDAL in planetary community. It is because a lot of data processing, conversion, projection, are done in ISIS3. However, ISIS3 doesn’t work that well with HRSC data (my main Martian data), VICAR-formatted data, and GeoTIFF you produce yourself. Because of that, I suggest to have gdal installed as well.

Preparation: install GDAL. You can download GDAL here. To work with VICAR data, gdal version has to be higher than 2.0.0

Checking gdal version

gdal –version

Preparation: set projection information.

Compared to Earth where proj4 for projection are easily available as well as definable with just a few characters like WGS84 and EPSG:4926, you need to define yours for Mars. Instead of typing a new CRS/proj4 every time you convert a strip, why not save them in some files

Generally, we can have MOLA Sphere (a,b,c=3396000) or areoid (a,b=3396190 c=3376200). For the projection, we can have latlon, equirectangular, sinusoidal, or stereographic

Create a new textfile (in Notepad, gedit, vi, or Atom), write the proj4, save it under name you easily remember (remember, no space!) such as MarsSouthStereo, MarsEQC, MarsMOLAEQC, etc.

Some examples

Equirectangular areoid projection, centered in Lat=0 and Lon=0

+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +a=3396190 +b=3396190 +units=m +no_defs

LatLon projection, MOLA sphere

+proj=longlat +a=3396000 +b=3396000 +no_defs

gdal_translate

Is useful to convert IMG to GeoTIFF, set no data value, change bit depth

Some examples

Converting a HRSC vicar file to GeoTIFF with world file (TFW), world file is defined by projection information

Gdal_translate H2288_0000_DT4.vic -co TFW=YES h2288_000_dt4.tif

Changing bit depth to 8bit and defining no data value as -32768

Gdal_translate H2288_0000_DT4.tif –a_nodata –32768 –ot Int H2288_0000_8bit_DT4.tif

Removing georeference information

Gdal_translate H2288_0000_DT4.tif –co PROFILE=BASELINE H2288_0000_DT4_base.tif

gdalinfo

Obtain projection information, size, bit depth, no data

gdalinfo -mm

The same with gdalinfo, with minimum and maximum values

gdalinfo –stats

Example

gdalinfo -stats H2288_0000_DT4.tif

Gdalwarp

Useful to change projection and fitting a strip to another strip

Example (with MArsMOLAEQC and MarsSouthStereo as projection files in the same folder with the image strip)

gdalwarp -s_srs MarsMOLAEQC -t_srs MarsSouthStereo H2288_0000_DT4.tif H2288_0000_D4_south_DT4.tif

(Will add more when I’m working with gdal)

About alfiahrizky

Currently Indonesian university student trying to submit more useful contents to share one's limited knowledge. Rather has wide range of interest with quite limited knowledge, the writings will be as sporadic as it can be. Expect more on technical related writings especially softwares and apps, cultures, music, and random things.

Posted on March 10, 2018, in Learning, Plugging, Spacing and tagged , , , , , . Bookmark the permalink. Leave a comment.

Something to be curious about?

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

w

Connecting to %s

%d bloggers like this: