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:323734 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.
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