[GRASS-user] Problem reprojecting to UTM 34 S

Rainer M Krug r.m.krug at gmail.com
Thu May 15 10:01:38 EDT 2008


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


More information about the grass-user mailing list