[GRASS-user] Importing Landsat images with different coordinates
Jenny Turner
jennyturner1980 at gmail.com
Tue May 25 09:10:32 EDT 2010
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>
> 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> 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> 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>> 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>
>>>>
>>>> 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
>>>> 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 listgrass-user at lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user
>
> This mail was received via Mail-SeCure System.
>
>
>
>
>
> --
> Micha Silverhttp://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/20100525/28fd1927/attachment-0001.html
More information about the grass-user
mailing list