<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#ffffff">
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:
<blockquote
cite="mid:42CFFC29FBA4FC48A39DBAFA5A2238D906E0B377@ATLEX03.lggm.llc"
type="cite">
<meta http-equiv="Content-Type"
content="text/html; charset=ISO-8859-1">
<meta name="Generator" content="Microsoft Word 11 (filtered medium)">
<!--[if !mso]>
<style>
v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style>
<![endif]-->
<title>GeoTIFF geographic grid alignment</title>
<o:smarttagtype
namespaceuri="urn:schemas-microsoft-com:office:smarttags"
name="country-region">
<o:smarttagtype
namespaceuri="urn:schemas-microsoft-com:office:smarttags" name="place"><!--[if !mso]>
<style>
st1\:*{behavior:url(#default#ieooui) }
</style>
<![endif]-->
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman";}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{color:blue;
text-decoration:underline;}
p
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman";}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:Arial;
color:navy;}
@page Section1
{size:8.5in 11.0in;
margin:1.0in 1.25in 1.0in 1.25in;}
div.Section1
{page:Section1;}
-->
</style>
<div class="Section1">
<p class="MsoNormal"><font face="Arial" size="2" color="navy"><span
style="font-size: 10pt; font-family: Arial; color: navy;">Tim,<o:p></o:p></span></font></p>
<p class="MsoNormal"><font face="Arial" size="2" color="navy"><span
style="font-size: 10pt; font-family: Arial; color: navy;"><o:p> </o:p></span></font></p>
<p class="MsoNormal"><font face="Arial" size="2" color="navy"><span
style="font-size: 10pt; font-family: Arial; color: navy;">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.<o:p></o:p></span></font></p>
<p class="MsoNormal"><font face="Arial" size="2" color="navy"><span
style="font-size: 10pt; font-family: Arial; color: navy;"><o:p> </o:p></span></font></p>
<p class="MsoNormal"><font face="Arial" size="2" color="navy"><span
style="font-size: 10pt; font-family: Arial; color: navy;">We use a
-0.5, -0.5 translation to
translate pixel is area to pixel is point just has you have done below.<o:p></o:p></span></font></p>
<p class="MsoNormal"><font face="Arial" size="2" color="navy"><span
style="font-size: 10pt; font-family: Arial; color: navy;"><o:p> </o:p></span></font></p>
<p class="MsoNormal"><font face="Arial" size="2" color="navy"><span
style="font-size: 10pt; font-family: Arial; color: navy;">Max<o:p></o:p></span></font></p>
<div>
<div class="MsoNormal" style="text-align: center;" align="center"><font
face="Times New Roman" size="3"><span style="font-size: 12pt;">
<hr tabindex="-1" size="2" width="100%" align="center"></span></font></div>
<p class="MsoNormal"><b><font face="Tahoma" size="2"><span
style="font-size: 10pt; font-family: Tahoma; font-weight: bold;">From:</span></font></b><font
face="Tahoma" size="2"><span
style="font-size: 10pt; font-family: Tahoma;">
<a class="moz-txt-link-abbreviated" href="mailto:geotiff-bounces@lists.maptools.org">geotiff-bounces@lists.maptools.org</a>
[<a class="moz-txt-link-freetext" href="mailto:geotiff-bounces@lists.maptools.org">mailto:geotiff-bounces@lists.maptools.org</a>] <b><span
style="font-weight: bold;">On Behalf Of </span></b><a class="moz-txt-link-abbreviated" href="mailto:dburken@comcast.net">dburken@comcast.net</a><br>
<b><span style="font-weight: bold;">Sent:</span></b> Monday, May 10,
2010 1:48 PM<br>
<b><span style="font-weight: bold;">To:</span></b> tjn98<br>
<b><span style="font-weight: bold;">Cc:</span></b>
<a class="moz-txt-link-abbreviated" href="mailto:Geotiff@lists.maptools.org">Geotiff@lists.maptools.org</a><br>
<b><span style="font-weight: bold;">Subject:</span></b> Re: [Geotiff]
GeoTIFF
geographic grid alignment</span></font><o:p></o:p></p>
</div>
<p class="MsoNormal"><font face="Times New Roman" size="3"><span
style="font-size: 12pt;"><o:p> </o:p></span></font></p>
<div>
<p class="MsoNormal" style="margin-bottom: 12pt;"><font face="Arial"
size="3" color="black"><span
style="font-size: 12pt; font-family: Arial; color: black;">Hi,<br>
<br>
If </span></font><font face="Consolas" size="2" color="black"><span
style="font-size: 10pt; font-family: Consolas; color: black;">RasterPixelIsPoint
I think it would be
(0.5, -0.5).<br>
<br>
Dave<br>
<br>
</span></font><font face="Arial" color="black"><span
style="font-family: Arial; color: black;">----- Original Message -----<br>
From: "tjn98" <a class="moz-txt-link-rfc2396E" href="mailto:tim.nightingale@stfc.ac.uk"><tim.nightingale@stfc.ac.uk></a><br>
To: <a class="moz-txt-link-abbreviated" href="mailto: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>
<br>
</span></font><font face="Consolas" size="2" color="black"><span
style="font-size: 10pt; font-family: Consolas; color: black;">Dear All,<br>
<br>
I’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’ve found, this is a common source of confusion.<br>
<br>
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 “PixelIsArea” raster described by:<br>
<br>
ModelTiepointTag = (0.0, 0.0, 0.0, -180.0, 90.0, 0.0)<br>
ModelPixelScaleTag = (0.002777777778, 0.002777777778,
0.0) <br>
GeoKeyDirectoryTag:<br>
GTModelTypeGeoKey
= 2 (ModelTypeGeographic)<br>
GTRasterTypeGeoKey
= 1 (RasterPixelIsArea)<br>
GeographicTypeGeoKey = 4326
(GCS_WGS_84)<br>
<br>
exactly fits into a –180 -> 180 by –90 -> 90 degree box that
touches<br>
the edges of the extreme cell boundaries, and that a “PixelIsPoint”<br>
raster, whose entries fall exactly in the centre of the
“PixelIsArea”<br>
cells in the first example is described by:<br>
<br>
ModelTiepointTag = (-0.5, -0.5, 0.0, -180.0, 90.0, 0.0)<br>
ModelPixelScaleTag = (0.002777777778, 0.002777777778,
0.0) <br>
GeoKeyDirectoryTag:<br>
GTModelTypeGeoKey
= 2 (ModelTypeGeographic)<br>
GTRasterTypeGeoKey
= 2 (RasterPixelIsPoint)<br>
GeographicTypeGeoKey = 4326
(GCS_WGS_84)<br>
<br>
Variants I have seen include both +0.5 and –0.5 pixel offsets
for<br>
the PixelIsArea case, negative ModelPixelScaleTag values for the<br>
“J” direction and 0.0 pixel offsets for the PixelIsPoint case.
QGIS,<br>
for example, appears to assume the last.<br>
<br>
Can anyone help?<br>
<br>
Many thanks,<br>
<br>
Tim.<br>
<br>
------------------------------------------------------------------------<br>
Dr T.J. Nightingale<br>
Space Science and Technology Department<br>
Rutherford Appleton Laboratory<br>
Chilton, Didcot
Phone:
+44/0 1235 445914<br>
Oxon OX11 0QX
Fax:
+44/0 1235 445848<br>
<st1:country-region w:st="on"><st1:place w:st="on">United Kingdom</st1:place></st1:country-region>
Email:
<a moz-do-not-send="true" href="tim.nightingale@stfc.ac.uk"
target="_blank">tim.nightingale@stfc.ac.uk</a><br>
------------------------------------------------------------------------
<br>
<br>
<br>
</span></font><font face="Arial" color="black"><span
style="font-family: Arial; color: black;"><o:p></o:p></span></font></p>
<p><font face="Arial" size="3" color="black"><span
style="font-size: 12pt; font-family: Arial; color: black;">-- <br>
Scanned by iCritical. <o:p></o:p></span></font></p>
<p class="MsoNormal"><font face="Arial" size="3" color="black"><span
style="font-size: 12pt; font-family: Arial; color: black;"><br>
<br>
_______________________________________________ Geotiff mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Geotiff@lists.maptools.org">Geotiff@lists.maptools.org</a>
<a class="moz-txt-link-freetext" href="http://lists.maptools.org/mailman/listinfo/geotiff">http://lists.maptools.org/mailman/listinfo/geotiff</a><o:p></o:p></span></font></p>
</div>
</div>
</o:smarttagtype></o:smarttagtype></blockquote>
<br>
</body>
</html>