[gdal-dev] Problem reprojecting Clark 1866 to UTM 34 South

[Frank Warmerdam]

Rainer M Krug wrote:
> Hi
> 
> I have a ESRI asc raster in +proj=longlat +ellps=D_Clarke_1866 +no_defs
...
> I want to reproject / warp it into UTM 34 S. I use the following command:
> 
> gdalwarp -s_srs '+proj=longlat +ellps=D_Clarke_1866 +no_defs' -t_srs
> EPSG:32734 hak_ser2.asc hak_ser2_utm34.asc
> 
> where EPSG:32734 is the definition for UTM 34 S WGS 84
> 
> But I get:
> 
> rkrug at ecolmod1:~/Documents/Projects/AlienSpread/GIS.local/Rainer_subset>
> gdalinfo hak_ser2_utm34.asc
> Driver: GTiff/GeoTIFF
> Files: hak_ser2_utm34.asc
> Size is 1028, 880
> Coordinate System is:
> PROJCS["WGS 84 / UTM zone 34N",
>    GEOGCS["WGS 84",

Rainer,

I suspect there is a problem with finding some of the support libgeotiff
or GDAL coordinate system files.

Do you have files in /usr/local/share/epsg_csv?

Do you have files in /usr/local/share/gdal?

I tried a similar command and the South aspect was properly preserved, so
I'm getting that when the support files are not found (either by libgeotiff
or GDAL) the fallback logic mishandles the south modifier.

PS. +ellps=D_Clark_1866 isn't likely to be recognised.  Use +ellps=clrk66.

Best regards,
-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | President OSGeo, http://osgeo.org