Ed,<br><br>Thanks for the information. <br><br>What should I use for the EPSG code when using mapserver or gdalwarp?<br><br>John<br><br><div><span class="gmail_quote">On 2/8/07, <b class="gmail_sendername">Ed McNierney</b>
<<a href="mailto:ed@topozone.com">ed@topozone.com</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">John -</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">Alabama East SPCS as EPSG:26929 is in units of meters, but
your GeoTIFF seems to think it's in feet. This will cause the reprojection
to move the file 3 times further from the CS origin than you
think.</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">It seems like you have a "hybrid" GeoTIFF with one
projection definition and a different EPSG code.</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">The ESRI-defined code of 102629 is an Alabama East NAD83
Feet projection. The definition is:</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span># NAD 1983 StatePlane
Alabama East FIPS 0101 Feet<br><102629> +proj=tmerc +lat_0=30.5
+lon_0=-85.83333333333333 +k=0.999960 +x_0=200000 +y_0=0 +ellps=GRS80
+datum=NAD83 +to_meter=0.3048006096012192 no_defs <></span></div>
<div dir="ltr" align="left"><span></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">Which is identical to the 26929 definition except for the
"to_meter" parameter.</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"> - Ed</font></span></div>
<div dir="ltr" align="left"><span>
<p><font size="2">Ed McNierney<br>President and Chief Mapmaker<br>TopoZone.com /
Maps a la carte, Inc.<br>73 Princeton Street, Suite 305<br>North Chelmsford,
MA 01863<br>Phone: <span onmouseup="SetCallButtonPressed(this, 0,1)" onmousedown="SetCallButtonPressed(this, 1,1)" onmouseover="SetCallButton(this, 1,1);skype_active=CheckCallButton(this);" onmouseout="SetCallButton(this, 0,1);HideSkypeMenu();" class="skype_tb_injection" id="__skype_highlight_id">
<span style="margin-right: 0px;" title="This is a United States of America phone number. The country code cannot be changed." onclick="javascript:if(0){doRunCMD(event, 'chdial','0');}else{doRunCMD(event, 'call','+19782514242');}event.preventBubble();return false;" onmouseout="SetCallButtonPart(this, 0);" onmouseover="SetCallButtonPart(this, 1);" class="skype_tb_injection_left" id="__skype_highlight_id_left">
<span style="background-image: url(chrome://skype_ff_toolbar_win/content/cb_normal_l.gif);" class="skype_tb_injection_left_img" id="__skype_highlight_id_left_adge"><img src="chrome://skype_ff_toolbar_win/content/cb_transparent_l.gif" style="height: 11px; width: 7px;" class="skype_tb_img_adge" height="11">
</span><span style="padding-right: 1px;" class="skype_tb_injection_left_img" id="__skype_highlight_id_left_img"><img style="width: 16px;" src="chrome://skype_ff_toolbar_win/content/famfamfam/us.gif" title="" class="skype_tb_img_flag" name="skype_tb_img_f0">
</span></span><span title="Call this phone number in United States of America with Skype: +19782514242" onclick="javascript:doRunCMD(event, 'call','+19782514242');event.preventBubble();return false;" onmouseout="SetCallButtonPart(this, 0)" onmouseover="SetCallButtonPart(this, 1)" class="skype_tb_injection_right" id="__skype_highlight_id_right">
<span class="skype_tb_innerText" id="__skype_highlight_id_innerText"><img src="chrome://skype_ff_toolbar_win/content/space.gif" style="margin: 0px; padding: 0px; height: 1px; width: 1px;" class="skype_tb_img_space" height="1" width="1">
<img src="chrome://skype_ff_toolbar_win/content/space.gif" style="margin: 0px; padding: 0px; height: 1px; width: 1px;" class="skype_tb_img_space" height="1" width="1"><img src="chrome://skype_ff_toolbar_win/content/space.gif" style="margin: 0px; padding: 0px; height: 1px; width: 1px;" class="skype_tb_img_space" height="1" width="1">
<img src="chrome://skype_ff_toolbar_win/content/space.gif" style="margin: 0px; padding: 0px; height: 1px; width: 1px;" class="skype_tb_img_space" height="1" width="1">+1 (978) 251-4242</span><span style="background-image: url(chrome://skype_ff_toolbar_win/content/cb_normal_r.gif);" class="skype_tb_injection_left_img" id="__skype_highlight_id_right_adge">
<img src="chrome://skype_ff_toolbar_win/content/cb_transparent_r.gif" style="height: 11px; width: 19px;" class="skype_tb_img_adge" height="11"></span></span></span><br>Fax: +1 (978)
251-1396<br><a href="mailto:ed@topozone.com" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">ed@topozone.com</a></font> </p></span></div><br>
<div dir="ltr" align="left" lang="en-us">
<hr>
<font face="Tahoma" size="2"><b>From:</b> <a href="mailto:gdal-dev-bounces@lists.maptools.org" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">gdal-dev-bounces@lists.maptools.org</a>
[mailto:<a href="mailto:gdal-dev-bounces@lists.maptools.org" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">gdal-dev-bounces@lists.maptools.org</a>] <b>On Behalf Of </b>John
Mitchell<br><b>Sent:</b> Thursday, February 08, 2007 6:38 PM<br><b>To:</b>
<a href="mailto:MAPSERVER-USERS@lists.umn.edu" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">MAPSERVER-USERS@lists.umn.edu</a>; Gdal-Dev<br><b>Cc:</b>
<a href="mailto:mitchelljj98@gmail.com" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">mitchelljj98@gmail.com</a><br><b>Subject:</b> [Gdal-dev] Alabama geotiff tiles are
showing up in Virginia withinmapserver<br></font><br></div><div><span class="e" id="q_110a3d3a1497d006_1">
<div></div>Hi,<br><br>I have Alabama geotiff tiles which are showing up in
Virginia within mapserver. If you look at the original gdalinfo against
one of the tiles you will see that the corner coordinates are within Alabama,
but within MapServer if you do on the fly projection(EPSG:26929 to EPSG:4326) it
instead shows up within the state of Virginia. I also used gdalwarp to
change the projection from EPSG:26929 to EPSG:4326 and then displayed within
mapserver without having to reproject on the fly and it also shows up within
Virginia. Note: when running gdalinfo on the gdalwarp file which was reprojected
to EPSG:4326 it also has the coordinates within the state of Virginia.
<br><br>Below are the gdalinfo dumps first the original file and then the file
that was converted from EPSG:26929 to EPSG:4326<br><br>Driver:
GTiff/GeoTIFF<br>Size is 10000, 10000<br>Coordinate System is:<br>PROJCS["NAD83
/ Alabama East", <br>
GEOGCS["NAD83",<br>
DATUM["North_American_Datum_1983",<br>
SPHEROID["GRS
1980",6378137,298.2572221010002,<br>
AUTHORITY["EPSG","7019"]],
<br>
AUTHORITY["EPSG","6269"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4269"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",30.5],<br>
PARAMETER["central_meridian",-85.83333333333333],<br>
PARAMETER["scale_factor",0.99996], <br>
PARAMETER["false_easting",656166.6666666665],<br>
PARAMETER["false_northing",0],<br> UNIT["US survey
foot",0.3048006096012192,<br>
AUTHORITY["EPSG","9003"]], <br>
AUTHORITY["EPSG","26929"]]<br>Origin =
(723240.000000000000000,820340.000000000000000)<br>Pixel Size =
(0.500000000000000,-0.500000000000000)<br>Metadata:<br>
AREA_OR_POINT=Area<br> TIFFTAG_XRESOLUTION=72 <br>
TIFFTAG_YRESOLUTION=72<br> TIFFTAG_RESOLUTIONUNIT=2
(pixels/inch)<br>Corner Coordinates:<br>Upper Left (
723240.000, 820340.000) ( 85d36'54.59"W, 32d45'17.80"N)<br>Lower
Left ( 723240.000, 815340.000 ) ( 85d36'54.71"W,
32d44'28.33"N)<br>Upper Right ( 728240.000, 820340.000) (
85d35'56.04"W, 32d45'17.69"N)<br>Lower Right ( 728240.000,
815340.000) ( 85d35'56.17"W, 32d44'
28.22"N)<br>Center ( 725740.000,
817840.000) ( 85d36'25.38"W, 32d44'53.01"N)<br>Band 1 Block=10000x1 Type=Byte,
ColorInterp=Red<br> Overviews: 1250x1250<br>Band 2 Block=10000x1
Type=Byte, ColorInterp=Green <br> Overviews: 1250x1250<br>Band 3
Block=10000x1 Type=Byte, ColorInterp=Blue<br> Overviews:
1250x1250<br><br>-------------------------------------------------------------------------------------------------------------------
<br><br>Driver: GTiff/GeoTIFF<br>Size is 11930, 9476<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["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4326"]]<br>Origin = (-79.902533251890048
,37.745881181682833)<br>Pixel Size =
(0.000005030208747,-0.000005030208747)<br>Metadata:<br>
AREA_OR_POINT=Area<br>Corner Coordinates:<br>Upper Left (
-79.9025333, 37.7458812) ( 79d54'9.12"W, 37d44'45.17"N) <br>Lower
Left ( -79.9025333, 37.6982149) ( 79d54'9.12"W,
37d41'53.57"N)<br>Upper Right ( -79.8425229, 37.7458812) ( 79d50'33.08"W,
37d44'45.17"N)<br>Lower Right ( -79.8425229, 37.6982149 ) ( 79d50'33.08"W,
37d41'53.57"N)<br>Center ( -79.8725281,
37.7220481) ( 79d52'21.10"W, 37d43'19.37"N)<br>Band 1 Block=256x256 Type=Byte,
ColorInterp=Red<br>Band 2 Block=256x256 Type=Byte, ColorInterp=Green <br>Band 3
Block=256x256 Type=Byte, ColorInterp=Blue<br><br clear="all"><br>-- <br>John J.
Mitchell </span></div></div>
</blockquote></div><br><br clear="all"><br>-- <br>John J. Mitchell