Unfortunately, GDAL does not support a particular projection I am trying to achieve (the ob_tran rotated pole projection), so I am trying to use PROJ4 to convert a geotiff file to the projection I want.
However, I am not sure what the correct syntax for such an operation is. Using the command:
invproj -m 57.2957795130823 +proj=ob_tran +o_proj=latlon +o_lon_p=20 +o_lat_p=20 +lon_0=180 input.tif
I get only a lot of gibberish in the console, such as:
133d13'9.044"E 62d0'32.793"N <Item name="missing_value" sample="0">-32767</Item> 133d13'9.044"E 62d0'32.793"N <Item name="NETCDF_VARNAME" sample="0">Band1</Item> 133d13'9.044"E 62d0'32.793"N</GDALMetadata> 133d13'9.044"E 62d0'32.793"N 133d13'9.044"E 62d0'32.793"N� �D�k�����������k�R�f�����P�w��B�����|�[�=�$
This is not due to the particular projection I am trying to use. Trying a normal Mercator projection with the command:
proj +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +units=m +k=1.0 +nadgrids=@null +no_defs input.tif
I also get a lot of gibberish.
0.00 -0.00�4�L�}�z���`����)�,�$�������� 0.00 -0.00��6�������G�R�f�@� 0.00 -0.00 0.00 -0.00� 0.00 -0.00���������3���<�����; 0.00 -0.00�����
What am I doing wrong? From looking at the manual, I thought I need to put the name of the file(s) at the end of the command line, but apparently, that is wrong?
How do I convert a file to another projection with PROJ4?