[GRASS-user] Off Topic: Importing GRASS tiffs to GMT

Patton, Eric Eric.Patton at NRCan-RNCan.gc.ca
Mon Apr 6 15:27:05 EDT 2009


I'll post here before checking on the GMT list; hopefully I can get some answers from a GRASS/GMT guru...

I've exported a tiff from GRASS using r.out.tiff -t, producing a worldfile in the process. I then converted the tiff to an 8-bit Sun raster image using Imagemagick's convert utility (GMT will only import imagery in this format). So far so good, the image looks fine. My question is how do plot the Sun raster so that it is correctly georeferenced on my GMT basemap? The GMT program psimage doesn't have a -R flag with which I may feed it geographical coordinates! And grd2image is only used for gridding xyz data, not for plotting georeferenced images. 

I already wrote a bash script to extract the upper left coordinates from the worldfile and then convert these positions into GMT-style -R flag like this:

Worldfile:
                  10.000000000000000 
                   0.000000000000000 
                   0.000000000000000 
                 -10.000000000000000 
              710175.000000000000000 
             5418295.000000000000000

Output from my script:

-54.160355/48.399269/-53.359791/48.860100

But can this information actually be used to imoprt the Sun raster? 

Sorry for the off-topic post. I know there are a few GMT users in the GRASS community, and maybe they have encountered this problem before.

Thanks,

~ Eric.


More information about the grass-user mailing list