Hi,
[yet another "funny" image problem. - yay]
Using Ubuntu 9.10 Karmic with gdal 1.5.4-4 (I think - how can I check
for sure? --help doesn't help)
I've tried various -t_srs permutations trying to convert a tiff ortho
from LATLONG to UTM34S but I keep ending up with PROJCS in the Northern
hemisphere, even though the Corner Coordinates show a latitude in the
South. The UTM coords are not in the expected value range, either.
Note I am not using -s_srs as I am expecting that from the TIFF headers.
gdalwarp -t_srs 'epsg:32734' 3319bc23.tif bc23.tif
Notice the UTM zone, Central Meridian and False Northing in the target
bc23.tif file.
Here are gdalinfo outputs for the following 3 images: (Third being a
reference image that is slightly East of the bc23 image)
-rw-rw-r-- 1 1001 geograph 376 2010-02-20 17:43 3319bc23.IF
-rw-rw-r-- 1 1001 geograph 118 2010-02-19 09:57 3319bc23.tfw
-rw-rw-r-- 1 1001 geograph 158072320 2010-02-19 09:57 3319bc23.tif
-rw-r--r-- 1 1001 geograph 295016802 2010-02-25 11:02 bc23.tif
-rw-r--r-- 1 1001 geograph 355 2010-02-24 09:19 worcester.IF
-rw-r--r-- 1 1001 geograph 207421516 2010-02-24 09:22 worcester.MR
-rw-r--r-- 1 1001 geograph 585550792 2010-02-24 09:11 worcester.tif
z...@gl1:/home/gisortho/t4a/breederivier$
gdalinfo 3319bc23.tif
Driver: GTiff/GeoTIFF
Files: 3319bc23.tif
Size is 7259, 7256
Coordinate System is:
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"]]
Origin = (19.598375571136735,-33.448385541900883)
Pixel Size = (0.000007335933852,-0.000007335933852)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 19.5983756, -33.4483855) ( 19d35'54.15"E, 33d26'54.19"S)
Lower Left ( 19.5983756, -33.5016151) ( 19d35'54.15"E, 33d30'5.81"S)
Upper Right ( 19.6516271, -33.4483855) ( 19d39'5.86"E, 33d26'54.19"S)
Lower Right ( 19.6516271, -33.5016151) ( 19d39'5.86"E, 33d30'5.81"S)
Center ( 19.6250013, -33.4750003) ( 19d37'30.00"E, 33d28'30.00"S)
Band 1 Block=7259x1 Type=Byte, ColorInterp=Red
Band 2 Block=7259x1 Type=Byte, ColorInterp=Green
Band 3 Block=7259x1 Type=Byte, ColorInterp=Blue
z...@gl1:/home/gisortho/t4a/breederivier$
z...@gl1:/home/gisortho/t4a/breederivier$
gdalinfo bc23.tif
Driver: GTiff/GeoTIFF
Files: bc23.tif
Size is 9688, 10117
Coordinate System is:
PROJCS["UTM Zone -34, Northern Hemisphere",
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",-387],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (4985652.260930398479104,-4864755.980756304226816)
Pixel Size = (0.945609074139101,-0.945609074139101)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 4985652.261,-4864755.981) ( 18d25'30.82"E, 33d43'5.31"S)
Lower Left ( 4985652.261,-4874322.708) ( 18d26'54.45"E, 33d46'51.98"S)
Upper Right ( 4994813.322,-4864755.981) ( 18d28'22.86"E, 33d41'23.95"S)
Lower Right ( 4994813.322,-4874322.708) ( 18d29'45.37"E, 33d45'10.62"S)
Center ( 4990232.791,-4869539.344) ( 18d27'38.64"E, 33d44'7.92"S)
Band 1 Block=9688x1 Type=Byte, ColorInterp=Red
Band 2 Block=9688x1 Type=Byte, ColorInterp=Green
Band 3 Block=9688x1 Type=Byte, ColorInterp=Blue
z...@gl1:/home/gisortho/t4a/breederivier$
Reference image:
z...@gl1:/home/gisortho/t4a/breederivier$
gdalinfo worcester.tif
Driver: GTiff/GeoTIFF
Files: worcester.tif
Size is 15556, 12545
Coordinate System is:
PROJCS["WGS 84 / UTM zone 34S",
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",10000000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32734"]]
Origin = (352199.486889008665457,6280000.033236353658140)
Pixel Size = (0.597834417917855,-0.597834417917858)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 352199.487, 6280000.033) ( 19d24'24.68"E, 33d36'34.29"S)
Lower Left ( 352199.487, 6272500.200) ( 19d24'20.20"E, 33d40'37.72"S)
Upper Right ( 361499.399, 6280000.033) ( 19d30'25.45"E, 33d36'38.79"S)
Lower Right ( 361499.399, 6272500.200) ( 19d30'21.25"E, 33d40'42.23"S)
Center ( 356849.443, 6276250.117) ( 19d27'22.89"E, 33d38'38.30"S)
Band 1 Block=15556x1 Type=Byte, ColorInterp=Red
Band 2 Block=15556x1 Type=Byte, ColorInterp=Green
Band 3 Block=15556x1 Type=Byte, ColorInterp=Blue
z...@gl1:/home/gisortho/t4a/breederivier$
TIA,
Zoltan
--
===========================================
Zoltan Szecsei
Director
Geograph (Pty) Ltd.
P.O. Box 7, Muizenberg 7950, South Africa.
65 Main Road, Muizenberg 7945
Western Cape, South Africa.
34° 6'16.35"S 18°28'5.62"E
Tel: +27-21-7884897 Mobile: +27-83-6004028
Fax: +27-86-6115323 www.geograph.co.za
===========================================
|
