Sorry Micha
My original Landsat image:
Driver: GTiff/GeoTIFF
Files: L71203033_03320000601_B10.TIF
L71203033_03320000601_B10.aux
L71203033_03320000601_B10.rrd
Size is 8161, 7141
Coordinate System is:
PROJCS["WGS 84 / UTM zone 29N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-9],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32629"]]
Origin = (521999.999999999941792,4414200.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
AREA_OR_POINT=Point
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 522000.000, 4414200.000) ( 8d44'33.81"W, 39d52'39.24"N)
Lower Left ( 522000.000, 4199970.000) ( 8d44'58.57"W, 37d56'49.39"N)
Upper Right ( 766830.000, 4414200.000) ( 5d52'54.42"W, 39d50'9.31"N)
Lower Right ( 766830.000, 4199970.000) ( 5d57'53.91"W, 37d54'29.39"N)
Center ( 644415.000, 4307085.000) ( 7d20'4.58"W, 38d54'2.83"N)
Band 1 Block=8161x1 Type=Byte, ColorInterp=Gray
Min=0.000 Max=255.000
Minimum=0.000, Maximum=255.000, Mean=60.374, StdDev=43.530
Overviews: 2041x1786, 1021x893, 511x447, 256x224, 128x112, 64x56
Metadata:
STATISTICS_MINIMUM=0
STATISTICS_MAXIMUM=255
STATISTICS_MEAN=60.374122874202
STATISTICS_MEDIAN=8.5986839406512e-287
STATISTICS_MODE=1.7801025473311e-306
STATISTICS_STDDEV=43.529818998446
LAYER_TYPE=athematic
And the obtained file:
Driver: GTiff/GeoTIFF
Files: destination.TIF
Size is 8229, 7218
Coordinate System is:
PROJCS["ETRS89 / Portugal TM06",
GEOGCS["ETRS89",
DATUM["European_Terrestrial_Reference_System_1989",
SPHEROID["GRS 1980",6378137,298.2572221010042,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6258"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4258"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",39.66825833333333],
PARAMETER["central_meridian",-8.133108333333334],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3763"]]
Origin = (-54187.411212536622770,23417.314857349840167)
Pixel Size = (30.006036956002887,-30.006036956002887)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -54187.411, 23417.315) ( 8d45'59.52"W, 39d52'38.78"N)
Lower Left ( -54187.411, -193166.260) ( 8d44'57.96"W, 37d55'35.73"N)
Upper Right ( 192732.267, 23417.315) ( 5d52'51.29"W, 39d51'26.27"N)
Lower Right ( 192732.267, -193166.260) ( 5d56'29.99"W, 37d54'28.07"N)
Center ( 69272.428, -84874.473) ( 7d20'4.40"W, 38d54'3.73"N)
Band 1 Block=8229x1 Type=Byte, ColorInterp=Gray
Hope this helps you (and mee :D )
Thank you
Jenny
2010/5/25 Micha Silver <mi...@arava.co.il <mailto:mi...@arava.co.il>>
On 05/25/2010 11:17 AM, Jenny Turner wrote:
Yes Sure :)
Attached to this email you can find the log
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.
Here's an example:
~>gdalinfo geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif
Driver: GTiff/GeoTIFF
Files: geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif
Size is 8577, 7573
Coordinate System is:
PROJCS["WGS 84 / UTM zone 36N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",33],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32636"]]
Origin = (603487.500000000000000,3623062.500000000000000)
Pixel Size = (28.500000000000000,-28.500000000000000)
Metadata:
TIFFTAG_XRESOLUTION=72
TIFFTAG_YRESOLUTION=72
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
AREA_OR_POINT=Point
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 603487.500, 3623062.500) ( 34d 6'16.43"E,
32d44'25.92"N)
Lower Left ( 603487.500, 3407232.000) ( 34d 4'54.01"E,
30d47'37.09"N)
Upper Right ( 847932.000, 3623062.500) ( 36d42'38.00"E,
32d41'25.66"N)
Lower Right ( 847932.000, 3407232.000) ( 36d38'1.81"E, 30d44'49.94"N)
Center ( 725709.750, 3515147.250) ( 35d22'58.18"E,
31d44'58.41"N)
Band 1 Block=8577x1 Type=Byte, ColorInterp=Gray
Thanks for your help
On Mon, May 24, 2010 at 8:46 PM, Micha Silver <mi...@arava.co.il
<mailto:mi...@arava.co.il>> wrote:
On 05/24/2010 06:26 PM, Jenny Turner wrote:
Hi
Following this topic I have tested this methodology by
converting from UTM 29N to a specific local coordinate
system (EPSG:3763) (available at PROJ4)
gdalwarp -t_srs "EPSG:3763" band_landsat.TIF destination.tif
With that I got a segmentation error, one
TIFFReadDirectory a lot of TIFFWriteEncodedStrip
and the destination.tif has a different size.
Can you send the output of: gdalinfo band_landsat.TIF
Any tips of
what migh be wrong?
Thank you
Jenny
On Tue, May 18, 2010 at 8:42 PM, Micha Silver
<mi...@arava.co.il <mailto:mi...@arava.co.il>> wrote:
Jenny Turner wrote:
Greetings
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?
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.
Thanks for your help
Best regards,
Jenny
On Mon, May 17, 2010 at 11:10 PM, Firman Hadi
<jalmibur...@gmail.com
<mailto:jalmibur...@gmail.com>
<mailto:jalmibur...@gmail.com
<mailto:jalmibur...@gmail.com>>> wrote:
Hi,
As far as I know you must define the source
coordinate system.
I think it is pretty straightforward when you are
using gdalwarp.
If there are so many images, maybe you can create
a shell script
to run the process.
Cheers,
>
Firman Hadi
Center for Remote Sensing - ITB
Jl. Ganesha No. 10,
Labtek IX-C, 3rd Floor
Bandung - 40132
INDONESIA
Phone: +62-22-2530701
Fax: +62-22-2530702
Website : http://crs.itb.ac.id ; www.sigro.org
<http://www.sigro.org> <http://www.sigro.org>
Blog : http://jalmiburung.wordpress.com
JepratJepret: http://jalmiburung.fotografer.net
This mail was received via Mail-SeCure System.
------------------------------------------------------------------------
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
<mailto:grass-user@lists.osgeo.org>
http://lists.osgeo.org/mailman/listinfo/grass-user
This mail was received via Mail-SeCure System.
--
Micha Silver
Arava Development Co. +972-52-3665918
http://www.surfaces.co.il
This mail was received via Mail-SeCure System.
-- Micha Silver Arava Development Co. +972-52-3665918
http://surfaces.co.il
This mail was received via Mail-SeCure System.
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org <mailto:grass-user@lists.osgeo.org>
http://lists.osgeo.org/mailman/listinfo/grass-user
This mail was received via Mail-SeCure System.
--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co. +972-52-3665918
This mail was received via Mail-SeCure System.