<HTML>
<HEAD>
<TITLE>Re: [Geotiff] GeoTIFF geographic grid alignment</TITLE>
</HEAD>
<BODY>
<FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'>Dear Max, David,<BR>
<BR>
&nbsp;&nbsp;Thanks for getting back to me. I see I&#8217;m not alone in my confusion.<BR>
I think we all agree that ModelTiePointTag 0.5 offsets (of whatever<BR>
sign) ARE required for RasterPixelIsPoint and ARE NOT required for<BR>
RasterPixelIsArea. This is a big step forward for me!<BR>
<BR>
&nbsp;&nbsp;We also agree that the ModelPixelScaleTag entries should be positive<BR>
(though I think this is inconsistent with the left-handed I/J system<BR>
and right-handed lon/lat system).<BR>
<BR>
&nbsp;&nbsp;I would have thought that the offset sign would have to be the same for<BR>
both the I and J directions (and K for that matter) as the handedness<BR>
of the IJK axis system doesn&#8217;t have any effect on this. Even though J<BR>
is positive down, the edge of the cell centred on the first point still<BR>
has to start &#8220;before&#8221; the beginning of the axis, as would be the case<BR>
for I also. Consequently, I&#8217;d think both must be negative.<BR>
<BR>
&nbsp;&nbsp;If the offsets were applied instead to the geographic tie point<BR>
information, then one adjustment would have to be positive (Sx) and<BR>
the other negative (Sy).<BR>
<BR>
&nbsp;&nbsp;How does this sound?<BR>
<BR>
&nbsp;&nbsp;&nbsp;&nbsp;Tim.<BR>
<BR>
------------------------------------------------------------------------<BR>
&nbsp;Tim 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>
&nbsp;<BR>
<BR>
<BR>
On 14/05/2010 09:03, &quot;David Burken&quot; &lt;<a href="dburken@comcast.net">dburken@comcast.net</a>&gt; wrote:<BR>
<BR>
</SPAN></FONT><BLOCKQUOTE><FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'>Max, Tim,<BR>
<BR>
If your tie upper left corner of the pixel (AREA) and relative to an Easting Northing of (10,000, 10,000) and your pixel scale is (100, 100) meters then to convert to point (center) it would be:<BR>
<BR>
Move the Easting 50 meters to the east or right ( +.5 pixel)<BR>
<BR>
Move the Northing 50 meters to the South or down ( -.5 pixel) <BR>
<BR>
So doing ( -.5, -5 ) is wrong in a map projection.<BR>
<BR>
We actually did a drawing of this on ossim.org to show the error when ignored but I cannot find it. Sorry.<BR>
<BR>
Hope that helps,<BR>
Dave<BR>
<BR>
<BR>
On 05/13/2010 05:15 PM, Max Martinez wrote: <BR>
</SPAN></FONT><BLOCKQUOTE><FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'> &nbsp;&nbsp;GeoTIFF geographic grid alignment &nbsp;&nbsp;&nbsp;<BR>
&nbsp;<BR>
<BR>
</SPAN></FONT><SPAN STYLE='font-size:10pt'><FONT COLOR="#000080"><FONT FACE="Arial">Tim,<BR>
</FONT></FONT><FONT FACE="Consolas, Courier New, Courier"> <BR>
</FONT><FONT COLOR="#000080"><FONT FACE="Arial"> <BR>
</FONT></FONT><FONT FACE="Consolas, Courier New, Courier"> <BR>
</FONT><FONT COLOR="#000080"><FONT FACE="Arial">I think the GeoTIFF spec tries to show through its examples that although the orientation of the positive J raster space axis is usually in the opposite direction of the positive Y model space axis, it is expected that both pixel scales will be positive under these conditions. (They probably could have been clearer if they described full map extents of the images in the examples). Negation should be used appropriately as the current condition differs from that usual condition. This is also supported by the example on page 27 where the ModelTransformationTag content equivalents of an image with tiepoint (I,J,K,X,Y,Z) and scale (Sx, Sy, Sz) is provided. So I think you have that right in your example.<BR>
</FONT></FONT><FONT FACE="Consolas, Courier New, Courier"> <BR>
</FONT><FONT COLOR="#000080"><FONT FACE="Arial"> <BR>
</FONT></FONT><FONT FACE="Consolas, Courier New, Courier"> <BR>
</FONT><FONT COLOR="#000080"><FONT FACE="Arial">We use a -0.5, -0.5 translation to translate pixel is area to pixel is point just has you have done below.<BR>
</FONT></FONT><FONT FACE="Consolas, Courier New, Courier"> <BR>
</FONT><FONT COLOR="#000080"><FONT FACE="Arial"> <BR>
</FONT></FONT><FONT FACE="Consolas, Courier New, Courier"> <BR>
</FONT><FONT COLOR="#000080"><FONT FACE="Arial">Max<BR>
</FONT></FONT><FONT FACE="Consolas, Courier New, Courier"> <BR>
&nbsp;
</FONT></SPAN>
<P ALIGN=CENTER>
<FONT SIZE="4"><FONT FACE="Times New Roman"><SPAN STYLE='font-size:12pt'> <BR>
<HR ALIGN=CENTER SIZE="2" WIDTH="100%"></SPAN></FONT></FONT>
<P>
<FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'> <BR>
<BR>
</SPAN></FONT><SPAN STYLE='font-size:10pt'><FONT FACE="Tahoma, Verdana, Helvetica, Arial"><B>From:</B> <a href="geotiff-bounces@lists.maptools.org">geotiff-bounces@lists.maptools.org</a> [<a href="mailto:geotiff-bounces@lists.maptools.org">mailto:geotiff-bounces@lists.maptools.org</a>] <B>On Behalf Of </B><a href="dburken@comcast.net">dburken@comcast.net</a><BR>
&nbsp;<B>Sent:</B> Monday, May 10, 2010 1:48 PM<BR>
&nbsp;<B>To:</B> tjn98<BR>
&nbsp;<B>Cc:</B> <a href="Geotiff@lists.maptools.org">Geotiff@lists.maptools.org</a><BR>
&nbsp;<B>Subject:</B> Re: [Geotiff] GeoTIFF geographic grid alignment<BR>
</FONT><FONT FACE="Consolas, Courier New, Courier"> <BR>
&nbsp;<BR>
</FONT></SPAN><FONT SIZE="4"><FONT FACE="Times New Roman"><SPAN STYLE='font-size:12pt'> <BR>
</SPAN></FONT></FONT><FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'> <BR>
&nbsp;<BR>
<BR>
</SPAN></FONT><FONT SIZE="4"><FONT FACE="Arial"><SPAN STYLE='font-size:12pt'>Hi,<BR>
&nbsp;<BR>
If </SPAN></FONT></FONT><FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'>RasterPixelIsPoint I think it would be (0.5, -0.5).<BR>
&nbsp;<BR>
Dave<BR>
&nbsp;<BR>
&nbsp;</SPAN></FONT><FONT SIZE="4"><FONT FACE="Arial"><SPAN STYLE='font-size:12pt'>----- Original Message -----<BR>
From: &quot;tjn98&quot; &lt;<a href="tim.nightingale@stfc.ac.uk">tim.nightingale@stfc.ac.uk</a>&gt; &lt;<a href="mailto:tim.nightingale@stfc.ac.uk">mailto:tim.nightingale@stfc.ac.uk</a>&gt; <BR>
To: <a href="Geotiff@lists.maptools.org">Geotiff@lists.maptools.org</a><BR>
Sent: Monday, May 10, 2010 10:28:07 AM GMT -05:00 US/Canada Eastern<BR>
Subject: [Geotiff] GeoTIFF geographic grid alignment<BR>
&nbsp;<BR>
&nbsp;</SPAN></FONT></FONT><FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'>Dear All,<BR>
&nbsp;<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>
&nbsp;<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>
&nbsp;<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>
&nbsp;<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>
&nbsp;<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>
&nbsp;<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>
&nbsp;<BR>
&nbsp;&nbsp;Can anyone help?<BR>
&nbsp;<BR>
&nbsp;&nbsp;&nbsp;&nbsp;Many thanks,<BR>
&nbsp;<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Tim.<BR>
&nbsp;<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>
&nbsp;<BR>
&nbsp;<BR>
&nbsp;<BR>
&nbsp;<BR>
</SPAN></FONT></BLOCKQUOTE></BLOCKQUOTE>

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