[GRASS-user] Problem reprojecting to UTM 34 S

[Michael Perdue]

 From what I see, gdalinfo hak_ser2.asc is listing all the latitudes  
as being positive which means that they are in the northern  
hemisphere. It looks to me like gdal is doing the right thing and  
there is something wrong with the ASCII grid.

Mike


On 15-May-08, at 8:01 AM, Rainer M Krug wrote:

> 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
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user