<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=utf-8">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","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
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></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"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Richard and others,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">I think I am close to projecting using gdal, but I am having issues with the to srs parameter. I first tried to supply to srs code, which I thought should
work because the code 103749 is defined in c:\ms4w\proj\nad\epsg, and works in my mapserver applications. Next, I tried to supply the Proj.4 definition. That did not work either. In both cases I received the following error: Option -t_srs incomplete, or
not recognised. Please let me know if there is a type-o in the following commands.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">gdal_translate -a_srs epsg:26915 -t_srs espg:103749 -CO TILED=YES -CO PHOTOMETRIC=YCBCR -CO TFW=YES -CO COMPRESS=JPEG -CO BIGTIFF=YES ortho_1-1_1n_s_mn081_2013_1.sid
lin13.tif<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">gdal_translate -a_srs epsg:26915 -t_srs "+proj=lcc +lat_1=44.25 +lat_2=44.58333333333334 +lat_0=44.19555555555555 +lon_0=-95.85000000000001 +x_0=152400.3048006096
+y_0=30480.06096012192 +a=6378559.758 + b=6357173.65471281 +units=ft +no_defs" -CO TILED=YES -CO PHOTOMETRIC=YCBCR -CO TFW=YES -CO COMPRESS=JPEG -CO BIGTIFF=YES ortho_1-1_1n_s_mn081_2013_1.sid lin13.tif<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">p.s. I noticed that “recognised” is either misspelled or possibly another English variant of the “recognized”<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Sincerely,<o:p></o:p></span></p>
<p class="MsoNormal"><b><span style="font-family:"Calibri","sans-serif";color:#1F497D">Mark Volz, GISP<o:p></o:p></span></b></p>
<p class="MsoNormal"><b><span style="font-family:"Calibri","sans-serif";color:#1F497D">Lyon County GIS Coordinator<o:p></o:p></span></b></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><a href="https://www.google.com/maps/?daddr=504%20Fairgrounds%20Rd,%20Marshall,%20MN%2056258?z=6"><span style="color:blue">504 Fairgrounds Rd</span></a></span><u><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:blue"><o:p></o:p></span></u></p>
<p class="MsoNormal"><u><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:blue"><a href="https://www.google.com/maps/?daddr=504%20Fairgrounds%20Rd,%20Marshall,%20MN%2056258?z=6"><span style="color:blue">Marshall, MN 56258</span></a></span></u><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Ph: (507) 532-8218<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Fax: (507) 532-8217<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><a href="http://lyonco.org/"><span style="color:blue">http://lyonco.org/</span></a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><a href="http://geomoose.lyonco.org/"><span style="color:blue">http://geomoose.lyonco.org/</span></a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div style="border:none;border-left:solid blue 1.5pt;padding:0in 0in 0in 4.0pt">
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> Richard Greenwood [mailto:richard.greenwood@gmail.com]
<br>
<b>Sent:</b> Wednesday, August 27, 2014 10:05 PM<br>
<b>To:</b> Mark Volz<br>
<b>Cc:</b> fwtools@lists.maptools.org<br>
<b>Subject:</b> Re: [FWTools] converting mr sid in UTM to GeoTiff in County Coordinates<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">From what you show, I don't see that gdal_translate is preforming any coordinate system transformation. However, the world file doesn't specify the coordinate system so I would suggest using -a_srs and -t_srs which will specify the source
and (optionally) the target coordinate systems. For example:<o:p></o:p></p>
</div>
<p class="MsoNormal"> gdal_translate -a_srs epsg:26915 -t_srs epsg:<your code> ...<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">would transform from UTM Zone 15 to your coord system.<o:p></o:p></p>
</div>
<p class="MsoNormal">You could also see if you have MrSID support installed in your GDAL setup with:
<o:p></o:p></p>
</div>
<p class="MsoNormal"> gdal_translate --formats<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">which would allow you to skip the ArcGIS export step and maybe some confusion.<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">There is a "nearblack" command line tool in the GDAL suite that can help with removing the black pixels around the image but I haven't found it too useful. There are a few other options, but I'd suggest you
work thru your coordinate system issues before dealing with the black pixels.<o:p></o:p></p>
</div>
<p class="MsoNormal">rich<o:p></o:p></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p>
<div>
<p class="MsoNormal">On Wed, Aug 27, 2014 at 3:19 PM, Mark Volz <<a href="mailto:MarkVolz@co.lyon.mn.us" target="_blank">MarkVolz@co.lyon.mn.us</a>> wrote:<o:p></o:p></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Hello,<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I would like to know how to use gdal to convert a Mr Sid air photo in UTM 15 to a GeoTiff in County Coordinates (feet).<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">My original steps include<o:p></o:p></p>
<p><span style="font-family:Symbol">·</span><span style="font-size:7.0pt">
</span>Use ArcGIS to export the MrSID in UTM to a GeoTiff in County Coordinates (feet).<o:p></o:p></p>
<p style="margin-left:1.0in"><span style="font-family:"Courier New"">o</span><span style="font-size:7.0pt">
</span>This takes care of converting the coordinate system, however at about 4GB, the file size is awfully huge.<o:p></o:p></p>
<p><span style="font-family:Symbol">·</span><span style="font-size:7.0pt">
</span>Use gdal_translate to tile the image. This also seems to make the file size much more reasonable… about 280 MB.<o:p></o:p></p>
<p style="margin-left:1.0in"><span style="font-family:"Courier New"">o</span><span style="font-size:7.0pt">
</span>gdal_translate -CO TILED=YES -CO PHOTOMETRIC=YCBCR -CO TFW=YES -CO COMPRESS=JPEG -CO BIGTIFF=YES LincolnFSA2013.tif, Lin13.tif<o:p></o:p></p>
<p><span style="font-family:Symbol">·</span><span style="font-size:7.0pt">
</span>Use gdaladdo to add internal JPEG overviews. This adds a small amount of data to the file, but it also makes it MUCH quicker.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Unfortunately, these steps are NOT working. It appears the coordinates are getting shifted from feet back to meters. I think it may be possible that gdal_translate is not reading
or ignoring the .tfw file created by ESRI. Therefore I need to either figure out how to force gdal_translate to read the ESRI world file, or I need to figure out how to figure out how to use only GDAL to create a GeoTiff in County Coordinates from a Mr Sid
file in UTM.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">If anyone has any tips on how I can use GDAL to convert an air photo in a Mr SID UTM file to a GeoTiff in County Coordinates let me know. In addition, I would like to also figure
out how to remove the black boundary around the edge of the air photos as well.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Thanks<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Sincerely,<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b>Mark Volz</b><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
_______________________________________________<br>
FWTools mailing list<br>
<a href="mailto:FWTools@lists.maptools.org">FWTools@lists.maptools.org</a><br>
<a href="http://lists.maptools.org/mailman/listinfo/fwtools" target="_blank">http://lists.maptools.org/mailman/listinfo/fwtools</a><br>
<a href="http://fwtools.maptools.org/" target="_blank">http://fwtools.maptools.org/</a><o:p></o:p></p>
</div>
<p class="MsoNormal"><br>
<br clear="all">
<br>
-- <o:p></o:p></p>
<div>
<p class="MsoNormal">Richard W. Greenwood, PLS<br>
<a href="http://www.greenwoodmap.com" target="_blank">www.greenwoodmap.com</a><o:p></o:p></p>
</div>
</div>
</div>
</div>
</body>
</html>