Hello,<br />
<br />
It would appear that GDAL14 does not understand arc-second units. I created<br />
an elevation file (GeoTIFF) using Global Mapper. When I look at it with<br />
gdalinfo I get:<br />
<br />
Driver: GTiff/GeoTIFF<br />
Size is 777, 657<br />
Coordinate System is:<br />
GEOGCS["WGS 84",<br />
DATUM["WGS_1984",<br />
SPHEROID["WGS 84",6378137,298.2572235630016,<br />
AUTHORITY["EPSG","7030"]],<br />
AUTHORITY["EPSG","6326"]],<br />
PRIMEM["Greenwich",0],<br />
UNIT[,0.0174532925199433],<br />
AUTHORITY["EPSG","4326"]]<br />
Origin = (-399424.610220731760000,116211.336444500940000)<br />
Pixel Size = (0.647917999999991,-0.647917999999996)<br />
Metadata:<br />
AREA_OR_POINT=Point<br />
Corner Coordinates:<br />
Upper Left ( -399424.610, 116211.336)<br />
(399424d36'36.79"W,116211d20'11.20"N)<br />
Lower Left ( -399424.610, 115785.654)<br />
(399424d36'36.79"W,115785d39'15.55"N)<br />
Upper Right ( -398921.178, 116211.336)<br />
(398921d10'40.57"W,116211d20'11.20"N)<br />
Lower Right ( -398921.178, 115785.654)<br />
(398921d10'40.57"W,115785d39'15.55"N)<br />
Center ( -399172.894, 115998.495)<br />
(399172d53'38.68"W,115998d29'43.37"N)<br />
Band 1 Block=777x1 Type=Float32, ColorInterp=Gray<br />
<br />
<script><!--
D(["mb","The UNIT is empty, but it clearly thinks the unit is DEGREES as the corner<br />coordinates are obviously completely wrong.<br /><br />When looking at the same file with a different reader (a GeoTIFF tag<br />reader), the information returned is:<br /><br />Elevation file:<br />Geotiff_Information:<br /> &nbsp; &nbsp;Version: 1<br /> &nbsp; &nbsp;Key_Revision: 1.0<br /> &nbsp; &nbsp;Tagged_Information:<br /> &nbsp; &nbsp; &nbsp; ModelTiepointTag (2,3):<br /> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0<br /> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;-399424.61 &nbsp; &nbsp; &nbsp; 116211.336 &nbsp; &nbsp; &nbsp; 0<br /> &nbsp; &nbsp; &nbsp; ModelPixelScaleTag (1,3):<br /> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0.647918 &nbsp; &nbsp; &nbsp; &nbsp; 0.647918 &nbsp; &nbsp; &nbsp; &nbsp; 1<br /> &nbsp; &nbsp; &nbsp; End_Of_Tags.<br /> &nbsp; &nbsp;Keyed_Information:<br /> &nbsp; &nbsp; &nbsp; GTModelTypeGeoKey (Short,1): ModelTypeGeographic<br /> &nbsp; &nbsp; &nbsp; GTRasterTypeGeoKey (Short,1): RasterPixelIsPoint<br /> &nbsp; &nbsp; &nbsp; GeographicTypeGeoKey (Short,1): GCS_WGS_84<br /> &nbsp; &nbsp; &nbsp; GeogAngularUnitsGeoKey (Short,1): Angular_Arc_Second<br /> &nbsp; &nbsp; &nbsp; VerticalUnitsGeoKey (Short,1): Unknown-9001<br /> &nbsp; &nbsp; &nbsp; End_Of_Keys.<br /> &nbsp; &nbsp;End_Of_Geotiff.<br /><br />The units are clearly set to Angular_Arc_Second. &nbsp;If the coordinates above<br />are interpreted as
arc-seconds then the corner coordinates are correct.<br /><br />First, am I understanding things correctly, and is my interpretation that<br />GDAL14 is not handling this correctly, correct? &nbsp;Or am I missing something?<br />I\'m new to GDAL and certainly could easily be overlooking something very<br />basic.<br /><br />If GDAL does in fact not handle this type of coordinate correctly, can<br />someone point me to where in the sources I would need to make a patch? &nbsp;I\'m<br />more than willing to go in and extend things as necessary.<br /><br />Thanks in advance,<br /></div>",1]
);
//--></script>The UNIT is empty, but it clearly thinks the unit is DEGREES as the corner<br />
coordinates are obviously completely wrong.<br />
<br />
When looking at the same file with a different reader (a GeoTIFF tag<br />
reader), the information returned is:<br />
<br />
Elevation file:<br />
Geotiff_Information:<br />
Version: 1<br />
Key_Revision: 1.0<br />
Tagged_Information:<br />
ModelTiepointTag (2,3):<br />
0 0 0<br />
-399424.61 116211.336 0<br />
ModelPixelScaleTag (1,3):<br />
0.647918 0.647918 1<br />
End_Of_Tags.<br />
Keyed_Information:<br />
GTModelTypeGeoKey (Short,1): ModelTypeGeographic<br />
GTRasterTypeGeoKey (Short,1): RasterPixelIsPoint<br />
GeographicTypeGeoKey (Short,1): GCS_WGS_84<br />
GeogAngularUnitsGeoKey (Short,1): Angular_Arc_Second<br />
VerticalUnitsGeoKey (Short,1): Unknown-9001<br />
End_Of_Keys.<br />
End_Of_Geotiff.<br />
<br />
The units are clearly set to Angular_Arc_Second. If the coordinates above<br />
are interpreted as arc-seconds then the corner coordinates are correct.<br />
<br />
First, am I understanding things correctly, and is my interpretation that<br />
GDAL14 is not handling this correctly, correct? Or am I missing something?<br />
I'm new to GDAL and certainly could easily be overlooking something very<br />
basic.<br />
<br />
If GDAL does in fact not handle this type of coordinate correctly, can<br />
someone point me to where in the sources I would need to make a patch? I'm<br />
more than willing to go in and extend things as necessary.<br />
<br />
Thanks in advance,<br />
<br />
Nick<br />