[GRASS-user] Importing Landsat images with different coordinates
Micha Silver
micha at arava.co.il
Wed May 26 04:43:33 EDT 2010
This all looks fine.
And you said that the destination tiff was smaller than the original?
Can you display the tiff files?
Maybe it's some disk problem (in the log you sent last time, it seems
there was an error reading the file...)?
On 05/25/2010 04:10 PM, Jenny Turner wrote:
> 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 <micha at arava.co.il <mailto:micha at 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 <micha at arava.co.il
>> <mailto:micha at 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
>>> <micha at arava.co.il <mailto:micha at 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
>>> <jalmiburung at gmail.com
>>> <mailto:jalmiburung at gmail.com>
>>> <mailto:jalmiburung at gmail.com
>>> <mailto:jalmiburung at 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 at lists.osgeo.org
>>> <mailto:grass-user at 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 at lists.osgeo.org <mailto:grass-user at 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.
--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co. +972-52-3665918
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20100526/a413a28d/attachment-0001.html
More information about the grass-user
mailing list