[GRASS-user] Problem reprojecting to UTM 34 S

[Rainer M Krug]

On Thu, May 15, 2008 at 4:46 PM, Michael Perdue <michael_perdue at yahoo.ca> wrote:
> 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.

A friend of mine reprojected the grid successfully using Arc GIS - so
there should be a way of doing it with gdal as well.



>
> 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
>
>
>



-- 
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