No subject
Wed Apr 10 02:58:04 EST 2013
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:
gdal_translate -a_srs epsg:26915 -t_srs epsg:<your code> ...
would transform from UTM Zone 15 to your coord system.
You could also see if you have MrSID support installed in your GDAL setup
with:
gdal_translate --formats
which would allow you to skip the ArcGIS export step and maybe some
confusion.
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.
rich
On Wed, Aug 27, 2014 at 3:19 PM, Mark Volz <MarkVolz at co.lyon.mn.us> wrote:
> Hello,
>
>
>
> 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).
>
>
>
> My original steps include
>
> =C2=B7 Use ArcGIS to export the MrSID in UTM to a GeoTiff in Coun=
ty
> Coordinates (feet).
>
> o This takes care of converting the coordinate system, however at about
> 4GB, the file size is awfully huge.
>
> =C2=B7 Use gdal_translate to tile the image. This also seems to =
make
> the file size much more reasonable=E2=80=A6 about 280 MB.
>
> o gdal_translate -CO TILED=3DYES -CO PHOTOMETRIC=3DYCBCR -CO TFW=3DYES =
-CO
> COMPRESS=3DJPEG -CO BIGTIFF=3DYES LincolnFSA2013.tif, Lin13.tif
>
> =C2=B7 Use gdaladdo to add internal JPEG overviews. This adds a =
small
> amount of data to the file, but it also makes it MUCH quicker.
>
>
>
> 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 t=
o
> use only GDAL to create a GeoTiff in County Coordinates from a Mr Sid fil=
e
> in UTM.
>
>
>
> 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 boundar=
y
> around the edge of the air photos as well.
>
>
>
> Thanks
>
>
>
>
>
> Sincerely,
>
> *Mark Volz*
>
>
>
> _______________________________________________
> FWTools mailing list
> FWTools at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/fwtools
> http://fwtools.maptools.org/
>
>
--=20
Richard W. Greenwood, PLS
www.greenwoodmap.com
--001a11c288f2af9b3d0501a7cc94
Content-Type: text/html; charset=UTF-8
Content-Transfer-Encoding: quoted-printable
<div dir=3D"ltr"><div><div><div><div><div><div><div>From what you show, I d=
on't see that gdal_translate is preforming any coordinate system transf=
ormation. 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 a=
nd (optionally) the target coordinate systems. For example:<br>
</div>=C2=A0=C2=A0 gdal_translate -a_srs=C2=A0epsg:26915 -t_srs epsg:<yo=
ur code> ...<br></div>would transform from UTM Zone 15 to your coord sys=
tem.<br><br></div>You could also see if you have MrSID support installed in=
your GDAL setup with:=C2=A0=C2=A0 <br>
</div>=C2=A0=C2=A0 gdal_translate --formats<br></div>which would allow you =
to skip the ArcGIS export step and maybe some confusion.<br><br></div>There=
is a "nearblack" command line tool in the GDAL suite that can he=
lp 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.<br=
>
<br></div>rich<br><div><div><div><div><div><div><br></div></div></div></div=
></div></div></div><div class=3D"gmail_extra"><br><br><div class=3D"gmail_q=
uote">On Wed, Aug 27, 2014 at 3:19 PM, Mark Volz <span dir=3D"ltr"><<a h=
ref=3D"mailto:MarkVolz at co.lyon.mn.us" target=3D"_blank">MarkVolz at co.lyon.mn=
.us</a>></span> wrote:<br>
<blockquote class=3D"gmail_quote" style=3D"margin:0 0 0 .8ex;border-left:1p=
x #ccc solid;padding-left:1ex">
<div link=3D"blue" vlink=3D"purple" lang=3D"EN-US">
<div>
<p class=3D"MsoNormal">Hello,<u></u><u></u></p>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
<p class=3D"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=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
<p class=3D"MsoNormal">My original steps include<u></u><u></u></p>
<p><u></u><span style=3D"font-family:Symbol"><span>=C2=B7<span style=3D"fon=
t:7.0pt "Times New Roman"">=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=
=C2=A0=C2=A0
</span></span></span><u></u>Use ArcGIS to export the MrSID in UTM to a GeoT=
iff in County Coordinates (feet).<u></u><u></u></p>
<p style=3D"margin-left:1.0in">
<u></u><span style=3D"font-family:"Courier New""><span>o<span sty=
le=3D"font:7.0pt "Times New Roman"">=C2=A0=C2=A0
</span></span></span><u></u>This takes care of converting the coordinate sy=
stem, however at about 4GB, the file size is awfully huge.<u></u><u></u></p=
>
<p><u></u><span style=3D"font-family:Symbol"><span>=C2=B7<span style=3D"fon=
t:7.0pt "Times New Roman"">=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=
=C2=A0=C2=A0
</span></span></span><u></u>Use gdal_translate to tile the image.=C2=A0 Thi=
s also seems to make the file size much more reasonable=E2=80=A6=C2=A0 abou=
t 280 MB.<u></u><u></u></p>
<p style=3D"margin-left:1.0in">
<u></u><span style=3D"font-family:"Courier New""><span>o<span sty=
le=3D"font:7.0pt "Times New Roman"">=C2=A0=C2=A0
</span></span></span><u></u>gdal_translate -CO TILED=3DYES -CO PHOTOMETRIC=
=3DYCBCR -CO TFW=3DYES -CO COMPRESS=3DJPEG -CO BIGTIFF=3DYES LincolnFSA2013=
.tif, =C2=A0Lin13.tif<u></u><u></u></p>
<p><u></u><span style=3D"font-family:Symbol"><span>=C2=B7<span style=3D"fon=
t:7.0pt "Times New Roman"">=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=
=C2=A0=C2=A0
</span></span></span><u></u>Use gdaladdo to add internal JPEG overviews.=C2=
=A0 This adds a small amount of data to the file, but it also makes it MUCH=
quicker.<u></u><u></u></p>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
<p class=3D"MsoNormal">Unfortunately, these steps are NOT working.=C2=A0 It=
appears the coordinates are getting shifted from feet back to meters.=C2=
=A0 I think it may be possible that gdal_translate is not reading or ignori=
ng =C2=A0the .tfw file created by ESRI.=C2=A0 Therefore
I need to either figure out how to force gdal_translate to read the ESRI w=
orld file, or I need to figure out how to figure out how to use only GDAL t=
o create a GeoTiff in County Coordinates from a Mr Sid file in UTM.<u></u><=
u></u></p>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
<p class=3D"MsoNormal">If anyone has any tips on how I can use GDAL to conv=
ert an air photo in a Mr SID UTM file to a GeoTiff in County Coordinates le=
t me know.=C2=A0 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=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
<p class=3D"MsoNormal">Thanks<u></u><u></u></p>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
<p class=3D"MsoNormal">Sincerely,<u></u><u></u></p>
<p class=3D"MsoNormal"><b><span style=3D"font-size:12.0pt">Mark Volz</span>=
</b><u></u><u></u></p>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
</div>
</div>
<br>_______________________________________________<br>
FWTools mailing list<br>
<a href=3D"mailto:FWTools at lists.maptools.org">FWTools at lists.maptools.org</a=
><br>
<a href=3D"http://lists.maptools.org/mailman/listinfo/fwtools" target=3D"_b=
lank">http://lists.maptools.org/mailman/listinfo/fwtools</a><br>
<a href=3D"http://fwtools.maptools.org/" target=3D"_blank">http://fwtools.m=
aptools.org/</a><br>
<br></blockquote></div><br><br clear=3D"all"><br>-- <br><div dir=3D"ltr">Ri=
chard W. Greenwood, PLS<br><a href=3D"http://www.greenwoodmap.com" target=
=3D"_blank">www.greenwoodmap.com</a></div>
</div>
--001a11c288f2af9b3d0501a7cc94--
More information about the FWTools
mailing list