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

[Rainer M Krug]

Hi

I have a ESRI asc raster in +proj=longlat +ellps=D_Clarke_1866 +no_defs

rkrug at ecolmod1:~/> gdalinfo hak_ser2.asc
Driver: AAIGrid/Arc/Info ASCII Grid
Files: hak_ser2.asc
Size is 1021, 781
Coordinate System is `'
Origin = (15.991666666664999,-21.991833330730000)
Pixel Size = (0.016666666670000,-0.016666666670000)
Corner Coordinates:
Upper Left  (  15.9916667, -21.9918333)
Lower Left  (  15.9916667, -35.0085000)
Upper Right (  33.0083333, -21.9918333)
Lower Right (  33.0083333, -35.0085000)
Center      (  24.5000000, -28.5001667)
Band 1 Block=1021x1 Type=Float32, ColorInterp=Undefined
 NoData Value=-9999

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",
       DATUM["WGS_1984",
           SPHEROID["WGS 84",6378137,298.2572235629972,
               AUTHORITY["EPSG","7030"]],
           AUTHORITY["EPSG","6326"]],
       PRIMEM["Greenwich",0],
       UNIT["degree",0.0174532925199433],
       AUTHORITY["EPSG","4326"]],
   PROJECTION["Transverse_Mercator"],
   PARAMETER["latitude_of_origin",0],
   PARAMETER["central_meridian",21],
   PARAMETER["scale_factor",0.9996],
   PARAMETER["false_easting",500000],
   PARAMETER["false_northing",0],
   UNIT["metre",1,
       AUTHORITY["EPSG","9001"]],
   AUTHORITY["EPSG","32634"]]
Origin = (-17471.391705359546904,-2431926.174459759611636)
Pixel Size = (1715.270741799227608,-1715.270741799227608)
Metadata:
 AREA_OR_POINT=Area
Image Structure Metadata:
 INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  -17471.392,-2431926.174) ( 15d59'39.65"E, 21d54'55.81"S)
Lower Left  (  -17471.392,-3941364.427) ( 15d17'55.27"E, 35d28'51.77"S)
Upper Right ( 1745826.931,-2431926.174) ( 32d58'5.17"E, 21d33'19.99"S)
Lower Right ( 1745826.931,-3941364.427) ( 34d34'59.95"E, 34d50'47.92"S)
Center      (  864177.770,-3186645.301) ( 24d43'43.82"E, 28d45'21.15"S)
Band 1 Block=1028x1 Type=Float32, ColorInterp=Gray

i.e. UTM 34 N  and it does not overlap with other layers in UTM 34 S.

Same result for

/> gdalwarp -s_srs '+proj=longlat +ellps=D_Clarke_1866 +no_defs'
-t_srs '+proj=utm +zone=34 +south +ellps=WGS84 +datum=WGS84 +
units=m +no_defs' hak_ser2.asc hak_ser2_utm34.asc

Any ideas what the problem is?

Thanks

Rainer

-- 
Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation
Biology, UCT), Dipl. Phys. (Germany)

Plant Conservation Unit
Department of Botany
University of Cape Town
Rondebosch 7701
South Africa