[GRASS-user] Importing Landsat images with different coordinates
Jenny Turner
jennyturner1980 at gmail.com
Wed May 26 06:25:08 EDT 2010
Micha
You are absolutely right. I was running it in a shared folder, and for some
reason (not space beucase it's not full), i was getting that error. Now, I
tried in a Desktop folder and it worked.
Thanks
2010/5/26 Micha Silver <micha at arava.co.il>
> 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>
>
>> 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
>>
>>
>
> 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/20100526/ef9a90a2/attachment-0001.html
More information about the grass-user
mailing list