<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 12 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
{mso-style-priority:34;
margin-top:0in;
margin-right:0in;
margin-bottom:0in;
margin-left:.5in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:2124302824;
mso-list-type:hybrid;
mso-list-template-ids:-2134229714 67698703 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l0:level1
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-.25in;}
ol
{margin-bottom:0in;}
ul
{margin-bottom:0in;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal>Hello again,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Been having good success reprojecting TIFs from/to various different projections (UTMs, State Plane, Mercator etc.) but I have run into a problem with gdalwarp generating a particular Mercator GeoTIFF. At the bottom of the mail is the console output from three operations:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='mso-list:Ignore'>1.<span style='font:7.0pt "Times New Roman"'> </span></span><![endif]>Gdalinfo output from the input, a standard UTM Zone 35N Geotiff<o:p></o:p></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='mso-list:Ignore'>2.<span style='font:7.0pt "Times New Roman"'> </span></span><![endif]>Gdalwarp output warping to a Mercator Projection.<o:p></o:p></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='mso-list:Ignore'>3.<span style='font:7.0pt "Times New Roman"'> </span></span><![endif]>Gdalinfo output from the resulting Tif.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Although the TIF image looks OK and is rotated about as you would expect, as you can see there are some invalid fields in the output from step #3,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>In case it helps the three Proj4 strings are:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>“+proj=utm +zone=35 +datum=WGS84 +units=m +no_defs” (from the original TIF)<o:p></o:p></p><p class=MsoNormal>"+proj=merc +datum=WGS84 +lon_0=0 +lat_ts=41.5 +x_0=0 +y_0=2000000" (input to gdawarp)<o:p></o:p></p><p class=MsoNormal>“+proj=merc +lon_0=0 +k=1 +x_0=-1.#IND +y_0=0 +datum=WGS84 +to_meter=0 +no_defs” (the proj4 string from the final TIF)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>If it helps I can provide links to the TIF files.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Can anyone spot what I am doing wrong?<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Many thanks,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Geoff<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>STEP#1:<o:p></o:p></p><p class=MsoNormal>------------------------------------------------------------------------------<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>gdalinfo -proj4 G:\Tifftest\Marmara\testTrans_0_TRANS.tif<o:p></o:p></p><p class=MsoNormal>Driver: GTiff/GeoTIFF<o:p></o:p></p><p class=MsoNormal>Files: G:\Tifftest\Marmara\testTrans_0_TRANS.tif<o:p></o:p></p><p class=MsoNormal> G:\Tifftest\Marmara\testTrans_0_TRANS.tfw<o:p></o:p></p><p class=MsoNormal>Size is 1030, 740<o:p></o:p></p><p class=MsoNormal>Coordinate System is:<o:p></o:p></p><p class=MsoNormal>PROJCS["WGS 84 / UTM zone 35N",<o:p></o:p></p><p class=MsoNormal> GEOGCS["WGS 84",<o:p></o:p></p><p class=MsoNormal> DATUM["WGS_1984",<o:p></o:p></p><p class=MsoNormal> SPHEROID["WGS 84",6378137,298.257223563,<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","7030"]],<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","6326"]],<o:p></o:p></p><p class=MsoNormal> PRIMEM["Greenwich",0],<o:p></o:p></p><p class=MsoNormal> UNIT["degree",0.0174532925199433],<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","4326"]],<o:p></o:p></p><p class=MsoNormal> PROJECTION["Transverse_Mercator"],<o:p></o:p></p><p class=MsoNormal> PARAMETER["latitude_of_origin",0],<o:p></o:p></p><p class=MsoNormal> PARAMETER["central_meridian",27],<o:p></o:p></p><p class=MsoNormal> PARAMETER["scale_factor",0.9996],<o:p></o:p></p><p class=MsoNormal> PARAMETER["false_easting",500000],<o:p></o:p></p><p class=MsoNormal> PARAMETER["false_northing",0],<o:p></o:p></p><p class=MsoNormal> UNIT["metre",1,<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","9001"]],<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","32635"]]<o:p></o:p></p><p class=MsoNormal>PROJ.4 string is:<o:p></o:p></p><p class=MsoNormal>'+proj=utm +zone=35 +datum=WGS84 +units=m +no_defs '<o:p></o:p></p><p class=MsoNormal>Origin = (640314.119512796050000,4530851.144891739800000)<o:p></o:p></p><p class=MsoNormal>Pixel Size = (0.999893837330017,-0.999390728012785)<o:p></o:p></p><p class=MsoNormal>Metadata:<o:p></o:p></p><p class=MsoNormal> AREA_OR_POINT=Area<o:p></o:p></p><p class=MsoNormal> TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)<o:p></o:p></p><p class=MsoNormal> TIFFTAG_XRESOLUTION=96<o:p></o:p></p><p class=MsoNormal> TIFFTAG_YRESOLUTION=96<o:p></o:p></p><p class=MsoNormal>Image Structure Metadata:<o:p></o:p></p><p class=MsoNormal> INTERLEAVE=PIXEL<o:p></o:p></p><p class=MsoNormal>Corner Coordinates:<o:p></o:p></p><p class=MsoNormal>Upper Left ( 640314.120, 4530851.145) ( 28d39'58.54"E, 40d55' 0.27"N)<o:p></o:p></p><p class=MsoNormal>Lower Left ( 640314.120, 4530111.596) ( 28d39'57.94"E, 40d54'36.30"N)<o:p></o:p></p><p class=MsoNormal>Upper Right ( 641344.010, 4530851.145) ( 28d40'42.55"E, 40d54'59.63"N)<o:p></o:p></p><p class=MsoNormal>Lower Right ( 641344.010, 4530111.596) ( 28d40'41.94"E, 40d54'35.66"N)<o:p></o:p></p><p class=MsoNormal>Center ( 640829.065, 4530481.370) ( 28d40'20.24"E, 40d54'47.96"N)<o:p></o:p></p><p class=MsoNormal>Band 1 Block=1030x1 Type=Byte, ColorInterp=Red<o:p></o:p></p><p class=MsoNormal> Mask Flags: PER_DATASET ALPHA<o:p></o:p></p><p class=MsoNormal>Band 2 Block=1030x1 Type=Byte, ColorInterp=Green<o:p></o:p></p><p class=MsoNormal> Mask Flags: PER_DATASET ALPHA<o:p></o:p></p><p class=MsoNormal>Band 3 Block=1030x1 Type=Byte, ColorInterp=Blue<o:p></o:p></p><p class=MsoNormal> Mask Flags: PER_DATASET ALPHA<o:p></o:p></p><p class=MsoNormal>Band 4 Block=1030x1 Type=Byte, ColorInterp=Alpha<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>STEP#2<o:p></o:p></p><p class=MsoNormal>------------------------------------------------------------------------------<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Gdalwarp -co TFW=YES -t_srs "+proj=merc +datum=WGS84 +lon_0=0 +lat_ts=41.5 +x_0=0 +y_0=20<o:p></o:p></p><p class=MsoNormal>00000" G:\Tifftest\Marmara\testTrans_0_TRANS.tif G:\Tifftest\Marmara\test.tif<o:p></o:p></p><p class=MsoNormal>Creating output file that is 1044P x 759L.<o:p></o:p></p><p class=MsoNormal>Processing input file G:\Tifftest\Marmara\testTrans_0_TRANS.tif.<o:p></o:p></p><p class=MsoNormal>Using band 4 of source image as alpha.<o:p></o:p></p><p class=MsoNormal>Using band 4 of destination image as alpha.<o:p></o:p></p><p class=MsoNormal>0...10...20...30...40...50...60...70...80...90...100 - done.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>STEP#3<o:p></o:p></p><p class=MsoNormal>-------------------------------------------------------------------------------<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>gdalinfo -proj4 G:\Tifftest\Marmara\test.tif<o:p></o:p></p><p class=MsoNormal>Driver: GTiff/GeoTIFF<o:p></o:p></p><p class=MsoNormal>Files: G:\Tifftest\Marmara\test.tif<o:p></o:p></p><p class=MsoNormal> G:\Tifftest\Marmara\test.tfw<o:p></o:p></p><p class=MsoNormal>Size is 1044, 759<o:p></o:p></p><p class=MsoNormal>Coordinate System is:<o:p></o:p></p><p class=MsoNormal>PROJCS["unnamed",<o:p></o:p></p><p class=MsoNormal> GEOGCS["WGS 84",<o:p></o:p></p><p class=MsoNormal> DATUM["WGS_1984",<o:p></o:p></p><p class=MsoNormal> SPHEROID["WGS 84",6378137,298.257223563,<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","7030"]],<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","6326"]],<o:p></o:p></p><p class=MsoNormal> PRIMEM["Greenwich",0],<o:p></o:p></p><p class=MsoNormal> UNIT["degree",0.0174532925199433],<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","4326"]],<o:p></o:p></p><p class=MsoNormal> PROJECTION["Mercator_1SP"],<o:p></o:p></p><p class=MsoNormal> PARAMETER["central_meridian",0],<o:p></o:p></p><p class=MsoNormal> PARAMETER["scale_factor",1],<o:p></o:p></p><p class=MsoNormal> PARAMETER["false_easting","-1.#IND"],<o:p></o:p></p><p class=MsoNormal> PARAMETER["false_northing","1.#INF"],<o:p></o:p></p><p class=MsoNormal> UNIT["degree minute second hemisphere",0,<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","9108"]]]<o:p></o:p></p><p class=MsoNormal>PROJ.4 string is:<o:p></o:p></p><p class=MsoNormal>'+proj=merc +lon_0=0 +k=1 +x_0=-1.#IND +y_0=0 +datum=WGS84 +to_meter=0 +no_defs<o:p></o:p></p><p class=MsoNormal>'<o:p></o:p></p><p class=MsoNormal>Origin = (2393508.975926194800000,5729350.050399892000000)<o:p></o:p></p><p class=MsoNormal>Pixel Size = (0.990987191737565,-0.990987191737565)<o:p></o:p></p><p class=MsoNormal>Metadata:<o:p></o:p></p><p class=MsoNormal> AREA_OR_POINT=Area<o:p></o:p></p><p class=MsoNormal>Image Structure Metadata:<o:p></o:p></p><p class=MsoNormal> INTERLEAVE=PIXEL<o:p></o:p></p><p class=MsoNormal>Corner Coordinates:<o:p></o:p></p><p class=MsoNormal>Upper Left ( 2393508.976, 5729350.050) ( 0d 0' 0.03"E, 0d 0' 0.01"N)<o:p></o:p></p><p class=MsoNormal>Lower Left ( 2393508.976, 5728597.891) ( 0d 0' 0.03"E, 0d 0' 0.01"N)<o:p></o:p></p><p class=MsoNormal>Upper Right ( 2394543.567, 5729350.050) ( 0d 0' 0.03"E, 0d 0' 0.01"N)<o:p></o:p></p><p class=MsoNormal>Lower Right ( 2394543.567, 5728597.891) ( 0d 0' 0.03"E, 0d 0' 0.01"N)<o:p></o:p></p><p class=MsoNormal>Center ( 2394026.271, 5728973.971) ( 0d 0' 0.03"E, 0d 0' 0.01"N)<o:p></o:p></p><p class=MsoNormal>Band 1 Block=1044x1 Type=Byte, ColorInterp=Red<o:p></o:p></p><p class=MsoNormal> Mask Flags: PER_DATASET ALPHA<o:p></o:p></p><p class=MsoNormal>Band 2 Block=1044x1 Type=Byte, ColorInterp=Green<o:p></o:p></p><p class=MsoNormal> Mask Flags: PER_DATASET ALPHA<o:p></o:p></p><p class=MsoNormal>Band 3 Block=1044x1 Type=Byte, ColorInterp=Blue<o:p></o:p></p><p class=MsoNormal> Mask Flags: PER_DATASET ALPHA<o:p></o:p></p><p class=MsoNormal>Band 4 Block=1044x1 Type=Byte, ColorInterp=Alpha<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:12.0pt;font-family:"Times New Roman","serif"'>Geoff Shipton,<br>Triton Imaging Inc.<br>Customer Support.<br>Capitola CA.<br><a href="mailto:geoff@tritonimaginginc.com"><span style='color:blue'>geoff@tritonimaginginc.com</span></a><br><br></span><span style='font-size:8.0pt;font-family:"Times New Roman","serif"'>The information contained in this message may be privileged or confidential information intended solely for the use of the addressee(s). Any views or opinions presented in this communication do not necessarily represent those of Triton Imaging, Inc. If the reader of this message is not the intended recipient, you are hereby notified that any dissemination, distribution or copying of this communication is strictly prohibited. If you have received this communication in error, please immediately notify me and you are hereby instructed to delete all electronic copies and destroy all printed copies.</span><span style='font-size:12.0pt;font-family:"Times New Roman","serif"'><o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>