Mapping package query

[Lester Anderson]

Hello,

I can extract the data and geotiff info from a file by the following:

>> [Image,CS] =  rasterread('Test.tif')

with the data stored in "Image" and the projection details in CS.
Similar results can be obtained from gdalread('Test.tif')

Does anyone know if it is possible to convert the EPSG Projection
string generated by the above to PROJ4 format (see below) within
Octave or not?

PROJ4 string: +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84
+units=m +no_defs

The octproj package works with the PROJ4 syntax.

Lester