Sorry Micha<br><br><div>My original Landsat image:</div><div><div>Driver: GTiff/GeoTIFF</div><div>Files: L71203033_03320000601_B10.TIF</div><div> L71203033_03320000601_B10.aux</div><div> L71203033_03320000601_B10.rrd</div>
<div>Size is 8161, 7141</div><div>Coordinate System is:</div><div>PROJCS["WGS 84 / UTM zone 29N",</div><div> GEOGCS["WGS 84",</div><div> DATUM["WGS_1984",</div><div> SPHEROID["WGS 84",6378137,298.257223563,</div>
<div> AUTHORITY["EPSG","7030"]],</div><div> AUTHORITY["EPSG","6326"]],</div><div> PRIMEM["Greenwich",0],</div><div> UNIT["degree",0.0174532925199433],</div>
<div> AUTHORITY["EPSG","4326"]],</div><div> PROJECTION["Transverse_Mercator"],</div><div> PARAMETER["latitude_of_origin",0],</div><div> PARAMETER["central_meridian",-9],</div>
<div> PARAMETER["scale_factor",0.9996],</div><div> PARAMETER["false_easting",500000],</div><div> PARAMETER["false_northing",0],</div><div> UNIT["metre",1,</div><div> AUTHORITY["EPSG","9001"]],</div>
<div> AUTHORITY["EPSG","32629"]]</div><div>Origin = (521999.999999999941792,4414200.000000000000000)</div><div>Pixel Size = (30.000000000000000,-30.000000000000000)</div><div>Metadata:</div><div> AREA_OR_POINT=Point</div>
<div>Image Structure Metadata:</div><div> INTERLEAVE=BAND</div><div>Corner Coordinates:</div><div>Upper Left ( 522000.000, 4414200.000) ( 8d44'33.81"W, 39d52'39.24"N)</div><div>Lower Left ( 522000.000, 4199970.000) ( 8d44'58.57"W, 37d56'49.39"N)</div>
<div>Upper Right ( 766830.000, 4414200.000) ( 5d52'54.42"W, 39d50'9.31"N)</div><div>Lower Right ( 766830.000, 4199970.000) ( 5d57'53.91"W, 37d54'29.39"N)</div><div>Center ( 644415.000, 4307085.000) ( 7d20'4.58"W, 38d54'2.83"N)</div>
<div>Band 1 Block=8161x1 Type=Byte, ColorInterp=Gray</div><div> Min=0.000 Max=255.000 </div><div> Minimum=0.000, Maximum=255.000, Mean=60.374, StdDev=43.530</div><div> Overviews: 2041x1786, 1021x893, 511x447, 256x224, 128x112, 64x56</div>
<div> Metadata:</div><div> STATISTICS_MINIMUM=0</div><div> STATISTICS_MAXIMUM=255</div><div> STATISTICS_MEAN=60.374122874202</div><div> STATISTICS_MEDIAN=8.5986839406512e-287</div><div> STATISTICS_MODE=1.7801025473311e-306</div>
<div> STATISTICS_STDDEV=43.529818998446</div><div> LAYER_TYPE=athematic</div><div><br></div><div>And the obtained file:</div><div><div>Driver: GTiff/GeoTIFF</div><div>Files: destination.TIF</div><div>Size is 8229, 7218</div>
<div>Coordinate System is:</div><div>PROJCS["ETRS89 / Portugal TM06",</div><div> GEOGCS["ETRS89",</div><div> DATUM["European_Terrestrial_Reference_System_1989",</div><div> SPHEROID["GRS 1980",6378137,298.2572221010042,</div>
<div> AUTHORITY["EPSG","7019"]],</div><div> AUTHORITY["EPSG","6258"]],</div><div> PRIMEM["Greenwich",0],</div><div> UNIT["degree",0.0174532925199433],</div>
<div> AUTHORITY["EPSG","4258"]],</div><div> PROJECTION["Transverse_Mercator"],</div><div> PARAMETER["latitude_of_origin",39.66825833333333],</div><div> PARAMETER["central_meridian",-8.133108333333334],</div>
<div> PARAMETER["scale_factor",1],</div><div> PARAMETER["false_easting",0],</div><div> PARAMETER["false_northing",0],</div><div> UNIT["metre",1,</div><div> AUTHORITY["EPSG","9001"]],</div>
<div> AUTHORITY["EPSG","3763"]]</div><div>Origin = (-54187.411212536622770,23417.314857349840167)</div><div>Pixel Size = (30.006036956002887,-30.006036956002887)</div><div>Metadata:</div><div> AREA_OR_POINT=Area</div>
<div>Image Structure Metadata:</div><div> INTERLEAVE=BAND</div><div>Corner Coordinates:</div><div>Upper Left ( -54187.411, 23417.315) ( 8d45'59.52"W, 39d52'38.78"N)</div><div>Lower Left ( -54187.411, -193166.260) ( 8d44'57.96"W, 37d55'35.73"N)</div>
<div>Upper Right ( 192732.267, 23417.315) ( 5d52'51.29"W, 39d51'26.27"N)</div><div>Lower Right ( 192732.267, -193166.260) ( 5d56'29.99"W, 37d54'28.07"N)</div><div>Center ( 69272.428, -84874.473) ( 7d20'4.40"W, 38d54'3.73"N)</div>
<div>Band 1 Block=8229x1 Type=Byte, ColorInterp=Gray</div><div><br></div><div>Hope this helps you (and mee :D )</div><div><br></div><div>Thank you</div><div>Jenny</div></div><div><br><div class="gmail_quote">2010/5/25 Micha Silver <span dir="ltr"><<a href="mailto:micha@arava.co.il">micha@arava.co.il</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div text="#000000" bgcolor="#ffffff"><div class="im">
On 05/25/2010 11:17 AM, Jenny Turner wrote:
</div><div class="im"><blockquote type="cite">Yes Sure :)
<div><br>
</div>
<div>Attached to this email you can find the log</div>
</blockquote>
<br></div><div><div></div><div class="h5">
I was thinking of just the results of gdalinfo. This way we can see if
the tiff is being recognized, and if its CRS is is available in the
tiff header.<br>
<br>
Here's an example:<br>
<br>
~>gdalinfo geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif <br>
Driver: GTiff/GeoTIFF<br>
Files: geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif<br>
Size is 8577, 7573<br>
Coordinate System is:<br>
PROJCS["WGS 84 / UTM zone 36N",<br>
GEOGCS["WGS 84",<br>
DATUM["WGS_1984",<br>
SPHEROID["WGS 84",6378137,298.257223563,<br>
AUTHORITY["EPSG","7030"]],<br>
AUTHORITY["EPSG","6326"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4326"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",0],<br>
PARAMETER["central_meridian",33],<br>
PARAMETER["scale_factor",0.9996],<br>
PARAMETER["false_easting",500000],<br>
PARAMETER["false_northing",0],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]],<br>
AUTHORITY["EPSG","32636"]]<br>
Origin = (603487.500000000000000,3623062.500000000000000)<br>
Pixel Size = (28.500000000000000,-28.500000000000000)<br>
Metadata:<br>
TIFFTAG_XRESOLUTION=72<br>
TIFFTAG_YRESOLUTION=72<br>
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)<br>
AREA_OR_POINT=Point<br>
Image Structure Metadata:<br>
INTERLEAVE=BAND<br>
Corner Coordinates:<br>
Upper Left ( 603487.500, 3623062.500) ( 34d 6'16.43"E, 32d44'25.92"N)<br>
Lower Left ( 603487.500, 3407232.000) ( 34d 4'54.01"E, 30d47'37.09"N)<br>
Upper Right ( 847932.000, 3623062.500) ( 36d42'38.00"E, 32d41'25.66"N)<br>
Lower Right ( 847932.000, 3407232.000) ( 36d38'1.81"E, 30d44'49.94"N)<br>
Center ( 725709.750, 3515147.250) ( 35d22'58.18"E, 31d44'58.41"N)<br>
Band 1 Block=8577x1 Type=Byte, ColorInterp=Gray<br>
<br>
<br>
<br>
</div></div><blockquote type="cite">
<div><br>
</div><div><div></div><div class="h5">
<div>Thanks for your help<br>
<br>
<div class="gmail_quote">On Mon, May 24, 2010 at 8:46 PM, Micha
Silver <span dir="ltr"><<a href="mailto:micha@arava.co.il" target="_blank">micha@arava.co.il</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left:1px solid rgb(204, 204, 204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex">
<div bgcolor="#ffffff" text="#000000">
<div>On 05/24/2010 06:26 PM, Jenny Turner wrote:
</div>
<div>
<blockquote type="cite">Hi
<div><br>
</div>
<div>Following this topic I have tested this methodology by
converting from UTM 29N to a specific local coordinate system
(EPSG:3763) (available at PROJ4)</div>
<div>gdalwarp -t_srs "EPSG:3763" band_landsat.TIF destination.tif</div>
<div><br>
</div>
<div>With that I got a segmentation error, one <span style="font-family:monospace;white-space:pre-wrap">TIFFReadDirectory<span style="font-family:arial;white-space:normal"> a lot of <span style="font-family:monospace;white-space:pre-wrap">TIFFWriteEncodedStrip<br>
and the destination.tif has a different size.</span></span></span></div>
<div><span style="font-family:monospace;white-space:pre-wrap"><span style="font-family:arial;white-space:normal"><span style="font-family:monospace;white-space:pre-wrap"><br>
<br>
<br>
<br>
</span></span></span></div>
</blockquote>
</div>
<div>Can you send the output of: gdalinfo
band_landsat.TIF<br>
<br>
</div>
<div>
<div>
<blockquote type="cite">
<div><span style="font-family:monospace;white-space:pre-wrap"><span style="font-family:arial;white-space:normal"><span style="font-family:monospace;white-space:pre-wrap"></span></span></span></div>
<div><span style="font-family:monospace;white-space:pre-wrap"><span style="font-family:arial;white-space:normal"><span style="font-family:monospace;white-space:pre-wrap">Any tips of<br>
what migh be wrong?</span></span></span></div>
<div><span style="font-family:monospace;white-space:pre-wrap"><span style="font-family:arial;white-space:normal"><span style="font-family:monospace;white-space:pre-wrap"><br>
<br>
Thank you</span></span></span></div>
<div><span style="font-family:monospace;white-space:pre-wrap"><span style="font-family:arial;white-space:normal"><span style="font-family:monospace;white-space:pre-wrap">Jenny</span></span></span></div>
<div><br>
<div class="gmail_quote">On Tue, May 18, 2010 at 8:42 PM, Micha
Silver <span dir="ltr"><<a href="mailto:micha@arava.co.il" target="_blank">micha@arava.co.il</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="border-left:1px solid rgb(204, 204, 204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex">
<div>Jenny Turner wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left:1px solid rgb(204, 204, 204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex">Greetings<br>
Firman, but imagine that I have rasters in UTM for differente zones
(e.g. huge countries like USA, etc). How can I produce an automatic
importing tool that uses gdalwarp to project an image to WGS84 and
import to GRASS?<br>
</blockquote>
</div>
Using gdalwarp you don't have to specify the source coordinate system,
*if* the originals have the correct projection defined in their
headers. For Landsat images, this is most likely the case, so you
should be able to get away with just setting the -t_srs option of
gdalwarp and each tile will be re-projected to the WGS84 geographic CRS.<br>
<br>
<blockquote class="gmail_quote" style="border-left:1px solid rgb(204, 204, 204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex">
<div>Thanks for your help<br>
Best regards,<br>
Jenny<br>
<br>
</div>
<div>On Mon, May 17, 2010 at 11:10 PM, Firman Hadi <<a href="mailto:jalmiburung@gmail.com" target="_blank">jalmiburung@gmail.com</a> <mailto:<a href="mailto:jalmiburung@gmail.com" target="_blank">jalmiburung@gmail.com</a>>> wrote:<br>
<br>
Hi,<br>
<br>
As far as I know you must define the source coordinate system.<br>
I think it is pretty straightforward when you are using gdalwarp.<br>
If there are so many images, maybe you can create a shell script<br>
to run the process.<br>
<br>
Cheers,<br>
<br>
><br>
<br>
Firman Hadi<br>
Center for Remote Sensing - ITB<br>
Jl. Ganesha No. 10,<br>
Labtek IX-C, 3rd Floor<br>
Bandung - 40132<br>
INDONESIA<br>
Phone: +62-22-2530701<br>
Fax: +62-22-2530702<br>
</div>
Website : <a href="http://crs.itb.ac.id" target="_blank">http://crs.itb.ac.id</a> ; <a href="http://www.sigro.org" target="_blank">www.sigro.org</a> <<a href="http://www.sigro.org" target="_blank">http://www.sigro.org</a>>
<div><br>
Blog : <a href="http://jalmiburung.wordpress.com" target="_blank">http://jalmiburung.wordpress.com</a><br>
JepratJepret: <a href="http://jalmiburung.fotografer.net" target="_blank">http://jalmiburung.fotografer.net</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
</div>
This mail was received via Mail-SeCure System.<br>
------------------------------------------------------------------------<br>
<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
<br>
This mail was received via Mail-SeCure System.<br>
<br>
<br>
<br>
</blockquote>
<br>
<br>
-- <br>
Micha Silver<br>
Arava Development Co. +972-52-3665918<br>
<a href="http://www.surfaces.co.il" target="_blank">http://www.surfaces.co.il</a><br>
<br>
<br>
</blockquote>
</div>
<br>
</div>
<br>
This mail was received via Mail-SeCure System.<br>
</blockquote>
<br>
<br>
</div>
</div>
<pre cols="72"><div><div>--
Micha Silver
Arava Development Co. +972-52-3665918
</div></div><a href="http://surfaces.co.il" target="_blank">http://surfaces.co.il</a>
</pre>
</div>
</blockquote>
</div>
<br>
</div>
<br>
This mail was received via Mail-SeCure System.<br>
<pre><fieldset></fieldset>
_______________________________________________
grass-user mailing list
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a>
This mail was received via Mail-SeCure System.
</pre>
</div></div></blockquote>
<br>
<br>
<pre cols="72">--
Micha Silver
<a href="http://www.surfaces.co.il/" target="_blank">http://www.surfaces.co.il/</a><div class="im">
Arava Development Co. +972-52-3665918</div></pre>
</div>
</blockquote></div><br></div></div>