<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=us-ascii">
<META content="MSHTML 6.00.5730.11" name=GENERATOR></HEAD>
<BODY>
<DIV dir=ltr align=left><SPAN class=432522002-09022007><FONT face=Arial
color=#0000ff size=2>John -</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=432522002-09022007><FONT face=Arial
color=#0000ff size=2></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=432522002-09022007><FONT face=Arial
color=#0000ff size=2>I think I answered that question below. Remember that
an "EPSG code" really means "a projection definition with the specified ID
number installed in the text file named epsg on your computer". If the
projection definition I described below for the 102629 projection is not already
present in your epsg file (a text file that's part of your PROJ installation)
you can edit that file by inserting the lines I provided below. The 1xxxxx
series EPSG codes are ESRI-defined extensions to the standard EPSG set, so they
are often not found in user installations that don't use ESRI
software.</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=432522002-09022007><FONT face=Arial
color=#0000ff size=2></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=432522002-09022007> <FONT
face=Arial color=#0000ff size=2>- Ed</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=432522002-09022007><!-- Converted from text/plain format -->
<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: +1 (978) 251-4242<BR>Fax: +1 (978)
251-1396<BR>ed@topozone.com</FONT> </P></SPAN></DIV><BR>
<DIV class=OutlookMessageHeader lang=en-us dir=ltr align=left>
<HR tabIndex=-1>
<FONT face=Tahoma size=2><B>From:</B> John Mitchell
[mailto:mitchelljj98@gmail.com] <BR><B>Sent:</B> Thursday, February 08, 2007
8:52 PM<BR><B>To:</B> Ed McNierney<BR><B>Cc:</B> MAPSERVER-USERS@lists.umn.edu;
Gdal-Dev; mitchelljj98@gmail.com<BR><B>Subject:</B> Re: [Gdal-dev] Alabama
geotiff tiles are showing up in Virginia withinmapserver<BR></FONT><BR></DIV>
<DIV></DIV>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="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">
<DIV>
<DIV dir=ltr align=left><SPAN><FONT face=Arial color=#0000ff size=2>John
-</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN><FONT face=Arial color=#0000ff
size=2></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN><FONT face=Arial color=#0000ff 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 face=Arial color=#0000ff
size=2></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN><FONT face=Arial color=#0000ff 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 face=Arial color=#0000ff
size=2></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN><FONT face=Arial color=#0000ff 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 face=Arial color=#0000ff
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 face=Arial color=#0000ff 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 face=Arial color=#0000ff
size=2></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN><FONT face=Arial color=#0000ff
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)"
class=skype_tb_injection onmousedown="SetCallButtonPressed(this, 1,1)"
id=__skype_highlight_id
onmouseover="SetCallButton(this, 1,1);skype_active=CheckCallButton(this);"
onmouseout="SetCallButton(this, 0,1);HideSkypeMenu();"><SPAN
class=skype_tb_injection_left id=__skype_highlight_id_left
onmouseover="SetCallButtonPart(this, 1);"
title="This is a United States of America phone number. The country code cannot be changed."
style="MARGIN-RIGHT: 0px"
onclick="javascript:if(0){doRunCMD(event, 'chdial','0');}else{doRunCMD(event, 'call','+19782514242');}event.preventBubble();return false;"
onmouseout="SetCallButtonPart(this, 0);"><SPAN
class=skype_tb_injection_left_img id=__skype_highlight_id_left_adge
style="BACKGROUND-IMAGE: url(chrome://skype_ff_toolbar_win/content/cb_normal_l.gif)"><IMG
class=skype_tb_img_adge style="WIDTH: 7px; HEIGHT: 11px" height=11
src="chrome://skype_ff_toolbar_win/content/cb_transparent_l.gif" NOSEND="1">
</SPAN><SPAN class=skype_tb_injection_left_img
id=__skype_highlight_id_left_img style="PADDING-RIGHT: 1px"><IMG
class=skype_tb_img_flag title="" style="WIDTH: 16px"
src="chrome://skype_ff_toolbar_win/content/famfamfam/us.gif"
name=skype_tb_img_f0 NOSEND="1"> </SPAN></SPAN><SPAN
class=skype_tb_injection_right id=__skype_highlight_id_right
onmouseover="SetCallButtonPart(this, 1)"
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)"><SPAN class=skype_tb_innerText
id=__skype_highlight_id_innerText><IMG class=skype_tb_img_space
style="PADDING-RIGHT: 0px; PADDING-LEFT: 0px; PADDING-BOTTOM: 0px; MARGIN: 0px; WIDTH: 1px; PADDING-TOP: 0px; HEIGHT: 1px"
height=1 src="chrome://skype_ff_toolbar_win/content/space.gif" width=1
NOSEND="1"> <IMG class=skype_tb_img_space
style="PADDING-RIGHT: 0px; PADDING-LEFT: 0px; PADDING-BOTTOM: 0px; MARGIN: 0px; WIDTH: 1px; PADDING-TOP: 0px; HEIGHT: 1px"
height=1 src="chrome://skype_ff_toolbar_win/content/space.gif" width=1
NOSEND="1"><IMG class=skype_tb_img_space
style="PADDING-RIGHT: 0px; PADDING-LEFT: 0px; PADDING-BOTTOM: 0px; MARGIN: 0px; WIDTH: 1px; PADDING-TOP: 0px; HEIGHT: 1px"
height=1 src="chrome://skype_ff_toolbar_win/content/space.gif" width=1
NOSEND="1"> <IMG class=skype_tb_img_space
style="PADDING-RIGHT: 0px; PADDING-LEFT: 0px; PADDING-BOTTOM: 0px; MARGIN: 0px; WIDTH: 1px; PADDING-TOP: 0px; HEIGHT: 1px"
height=1 src="chrome://skype_ff_toolbar_win/content/space.gif" width=1
NOSEND="1">+1 (978) 251-4242</SPAN><SPAN class=skype_tb_injection_left_img
id=__skype_highlight_id_right_adge
style="BACKGROUND-IMAGE: url(chrome://skype_ff_toolbar_win/content/cb_normal_r.gif)">
<IMG class=skype_tb_img_adge style="WIDTH: 19px; HEIGHT: 11px" height=11
src="chrome://skype_ff_toolbar_win/content/cb_transparent_r.gif"
NOSEND="1"></SPAN></SPAN></SPAN><BR>Fax: +1 (978) 251-1396<BR><A
onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:ed@topozone.com" target=_blank>ed@topozone.com</A></FONT>
</P></SPAN></DIV><BR>
<DIV lang=en-us dir=ltr align=left>
<HR>
<FONT face=Tahoma size=2><B>From:</B> <A
onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:gdal-dev-bounces@lists.maptools.org"
target=_blank>gdal-dev-bounces@lists.maptools.org</A> [mailto:<A
onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:gdal-dev-bounces@lists.maptools.org"
target=_blank>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 onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:MAPSERVER-USERS@lists.umn.edu"
target=_blank>MAPSERVER-USERS@lists.umn.edu</A>; Gdal-Dev<BR><B>Cc:</B> <A
onclick="return top.js.OpenExtLink(window,event,this)"
href="mailto:mitchelljj98@gmail.com"
target=_blank>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 </BODY></HTML>