I have a grid file in a certain coordinate system. I'd like to change it so I have the same grid in geographic lat/long based on WGS84. This seems like something that gdal can tackle quite easily (maybe with gdalwarp or gdal_translate) but for some reason I can't find the proper way to do it.
I'm trying to run the following:
gdalwarp -s_srs '+proj=cass +lat_0=31.73409694444445 +lon_0=35.21208055555556 +x_0=170251.555 +y_0=1126867.909 +a=6378300.789 +b=6356566.435 +towgs84=-275.722,94.7824,340.894,-8.001,-4.42,-11.821,1 +units=m +no_defs' -t_srs EPSG:4326 thick.grd thick2.grd
The proj4 definition taken from spatialreference.org. I get the following error message:
Warning: The target file has a 'grd' extension, which is normally used by the GSAG, GSBG, GS7BG, NWT_GRD drivers, but the requested output driver is GTiff. Is it really what you want ? ERROR 6: Unable to load PROJ.4 library (libproj.so), creation of OGRCoordinateTransformation failed.
So I have two problems here. First, how do I let gdalwarp know what is the input format (in this case, Surfer grid) and the output format? Also, why can't it load the PROJ.4 library if I have it installed? Is there an alternative to PROJ.4 that will work on my system?