<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Mark,<br>
<br>
If you want to reproject the data, I think you need to use
gdalwarp.<br>
<pre class="moz-signature" cols="72">Best Regards,
Brent Fraser</pre>
On 8/29/2014 7:54 AM, Mark Volz wrote:<br>
</div>
<blockquote
cite="mid:36CA828A36E29F45B7CF0A1766E5DFA33AA7A6CE@swmail01.r8nssis.local"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<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]-->
<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
moz-do-not-send="true"
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
moz-do-not-send="true"
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
moz-do-not-send="true" 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
moz-do-not-send="true" 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 [<a class="moz-txt-link-freetext" href="mailto:richard.greenwood@gmail.com">mailto:richard.greenwood@gmail.com</a>]
<br>
<b>Sent:</b> Wednesday, August 27, 2014 10:05 PM<br>
<b>To:</b> Mark Volz<br>
<b>Cc:</b> <a class="moz-txt-link-abbreviated" href="mailto:fwtools@lists.maptools.org">fwtools@lists.maptools.org</a><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 moz-do-not-send="true"
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 moz-do-not-send="true"
href="mailto:FWTools@lists.maptools.org">FWTools@lists.maptools.org</a><br>
<a moz-do-not-send="true"
href="http://lists.maptools.org/mailman/listinfo/fwtools"
target="_blank">http://lists.maptools.org/mailman/listinfo/fwtools</a><br>
<a moz-do-not-send="true"
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 moz-do-not-send="true"
href="http://www.greenwoodmap.com" target="_blank">www.greenwoodmap.com</a><o:p></o:p></p>
</div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
FWTools mailing list
<a class="moz-txt-link-abbreviated" href="mailto:FWTools@lists.maptools.org">FWTools@lists.maptools.org</a>
<a class="moz-txt-link-freetext" href="http://lists.maptools.org/mailman/listinfo/fwtools">http://lists.maptools.org/mailman/listinfo/fwtools</a>
<a class="moz-txt-link-freetext" href="http://fwtools.maptools.org/">http://fwtools.maptools.org/</a>
</pre>
</blockquote>
<br>
</body>
</html>