<div dir="ltr"><div>Ah yes. Sorry I didn't even notice that in the original post. As Brent notes, gdalwarp does coordinate transformations (reprojects), gdal_translate does not.<br><br></div>Rich<br></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Fri, Aug 29, 2014 at 5:05 PM, Brent Fraser <span dir="ltr"><<a href="mailto:bfraser@geoanalytic.com" target="_blank">bfraser@geoanalytic.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<div>Mark,<br>
<br>
If you want to reproject the data, I think you need to use
gdalwarp.<br>
<pre cols="72">Best Regards,
Brent Fraser</pre><div><div class="h5">
On 8/29/2014 7:54 AM, Mark Volz wrote:<br>
</div></div></div><div><div class="h5">
<blockquote type="cite">
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Richard
and others,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></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.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></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<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></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<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></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”<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Sincerely,<u></u><u></u></span></p>
<p class="MsoNormal"><b><span style="font-family:"Calibri","sans-serif";color:#1f497d">Mark
Volz, GISP<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span style="font-family:"Calibri","sans-serif";color:#1f497d">Lyon
County GIS Coordinator<u></u><u></u></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" target="_blank"><span style="color:blue">504 Fairgrounds Rd</span></a></span><u><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:blue"><u></u><u></u></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" target="_blank"><span style="color:blue">Marshall, MN 56258</span></a></span></u><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Ph:
<a href="tel:%28507%29%20532-8218" value="+15075328218" target="_blank">(507) 532-8218</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Fax:
<a href="tel:%28507%29%20532-8217" value="+15075328217" target="_blank">(507) 532-8217</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><a href="http://lyonco.org/" target="_blank"><span style="color:blue">http://lyonco.org/</span></a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><a href="http://geomoose.lyonco.org/" target="_blank"><span style="color:blue">http://geomoose.lyonco.org/</span></a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></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 href="mailto:richard.greenwood@gmail.com" target="_blank">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 href="mailto:fwtools@lists.maptools.org" target="_blank">fwtools@lists.maptools.org</a><br>
<b>Subject:</b> Re: [FWTools] converting mr sid in UTM
to GeoTiff in County Coordinates<u></u><u></u></span></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></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:<u></u><u></u></p>
</div>
<p class="MsoNormal"> gdal_translate
-a_srs epsg:26915 -t_srs epsg:<your
code> ...<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">would
transform from UTM Zone 15 to your coord system.<u></u><u></u></p>
</div>
<p class="MsoNormal">You could also see if you have
MrSID support installed in your GDAL setup with:
<u></u><u></u></p>
</div>
<p class="MsoNormal"> gdal_translate --formats<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">which
would allow you to skip the ArcGIS export step and
maybe some confusion.<u></u><u></u></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.<u></u><u></u></p>
</div>
<p class="MsoNormal">rich<u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><u></u> <u></u></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:<u></u><u></u></p>
<div>
<div>
<p class="MsoNormal">Hello,<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">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).<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">My
original steps include<u></u><u></u></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).<u></u><u></u></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.<u></u><u></u></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.<u></u><u></u></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<u></u><u></u></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.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">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.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">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.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Thanks<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">Sincerely,<u></u><u></u></p>
<p class="MsoNormal"><b>Mark
Volz</b><u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
_______________________________________________<br>
FWTools mailing list<br>
<a href="mailto:FWTools@lists.maptools.org" target="_blank">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><u></u><u></u></p>
</div>
<p class="MsoNormal"><br>
<br clear="all">
<br>
-- <u></u><u></u></p>
<div>
<p class="MsoNormal">Richard W. Greenwood, PLS<br>
<a href="http://www.greenwoodmap.com" target="_blank">www.greenwoodmap.com</a><u></u><u></u></p>
</div>
</div>
</div>
</div>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
FWTools mailing list
<a href="mailto:FWTools@lists.maptools.org" target="_blank">FWTools@lists.maptools.org</a>
<a href="http://lists.maptools.org/mailman/listinfo/fwtools" target="_blank">http://lists.maptools.org/mailman/listinfo/fwtools</a>
<a href="http://fwtools.maptools.org/" target="_blank">http://fwtools.maptools.org/</a>
</pre>
</blockquote>
<br>
</div></div></div>
<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><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr">Richard W. Greenwood, PLS<br><a href="http://www.greenwoodmap.com" target="_blank">www.greenwoodmap.com</a></div>
</div>