<HTML>
<HEAD>
<TITLE>GeoTIFF geographic grid alignment</TITLE>
</HEAD>
<BODY>
<FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'>Dear All,<BR>
<BR>
&nbsp;&nbsp;I&#8217;m looking for some guidance on the correct interpretation of<BR>
GeoTIFF raster spaces, in particular their correct alignments at the<BR>
half-pixel level. Judging by the number of mutually contradictory<BR>
examples I&#8217;ve found, this is a common source of confusion.<BR>
<BR>
&nbsp;&nbsp;My interpretation of the rather terse text in the GeoTIFF Format<BR>
Specification (Revision 1.0) is that, for the example of a global<BR>
map measuring 129,600 longitude pixels by 64,800 latitude pixels,<BR>
a &#8220;PixelIsArea&#8221; raster described by:<BR>
<BR>
&nbsp;&nbsp;&nbsp;&nbsp;ModelTiepointTag = (0.0, 0.0, 0.0, -180.0, 90.0, 0.0)<BR>
&nbsp;&nbsp;&nbsp;&nbsp;ModelPixelScaleTag = (0.002777777778, 0.002777777778, 0.0) <BR>
&nbsp;&nbsp;&nbsp;&nbsp;GeoKeyDirectoryTag:<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;GTModelTypeGeoKey &nbsp;&nbsp;&nbsp;= 2 &nbsp;&nbsp;&nbsp;(ModelTypeGeographic)<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;GTRasterTypeGeoKey &nbsp;&nbsp;= 1 &nbsp;&nbsp;&nbsp;(RasterPixelIsArea)<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;GeographicTypeGeoKey = 4326 (GCS_WGS_84)<BR>
<BR>
exactly fits into a &#8211;180 -&gt; 180 by &#8211;90 -&gt; 90 degree box that touches<BR>
the edges of the extreme cell boundaries, and that a &#8220;PixelIsPoint&#8221;<BR>
raster, whose entries fall exactly in the centre of the &#8220;PixelIsArea&#8221;<BR>
cells in the first example is described by:<BR>
<BR>
&nbsp;&nbsp;&nbsp;&nbsp;ModelTiepointTag = (-0.5, -0.5, 0.0, -180.0, 90.0, 0.0)<BR>
&nbsp;&nbsp;&nbsp;&nbsp;ModelPixelScaleTag = (0.002777777778, 0.002777777778, 0.0) <BR>
&nbsp;&nbsp;&nbsp;&nbsp;GeoKeyDirectoryTag:<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;GTModelTypeGeoKey &nbsp;&nbsp;&nbsp;= 2 &nbsp;&nbsp;&nbsp;(ModelTypeGeographic)<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;GTRasterTypeGeoKey &nbsp;&nbsp;= 2 &nbsp;&nbsp;&nbsp;(RasterPixelIsPoint)<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;GeographicTypeGeoKey = 4326 (GCS_WGS_84)<BR>
<BR>
&nbsp;&nbsp;Variants I have seen include both +0.5 and &#8211;0.5 pixel offsets for<BR>
the PixelIsArea case, negative ModelPixelScaleTag values for the<BR>
&#8220;J&#8221; direction and 0.0 pixel offsets for the PixelIsPoint case. QGIS,<BR>
for example, appears to assume the last.<BR>
<BR>
&nbsp;&nbsp;Can anyone help?<BR>
<BR>
&nbsp;&nbsp;&nbsp;&nbsp;Many thanks,<BR>
<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Tim.<BR>
<BR>
------------------------------------------------------------------------<BR>
&nbsp;Dr T.J. Nightingale<BR>
&nbsp;&nbsp;Space Science and Technology Department<BR>
&nbsp;&nbsp;&nbsp;Rutherford Appleton Laboratory<BR>
&nbsp;&nbsp;&nbsp;&nbsp;Chilton, Didcot &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Phone: +44/0 1235 445914<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Oxon OX11 0QX &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Fax: &nbsp;&nbsp;+44/0 1235 445848<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;United Kingdom &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Email: <a href="tim.nightingale@stfc.ac.uk">tim.nightingale@stfc.ac.uk</a><BR>
------------------------------------------------------------------------ <BR>
<BR>
<BR>
</SPAN></FONT>

<br><p>-- 
<BR>Scanned by iCritical.
</p>
<br></BODY>
</HTML>