[GRASS-user] Importing Landsat images with different coordinates
Micha Silver
micha at arava.co.il
Tue May 25 08:50:29 EDT 2010
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
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20100525/2e0e42f4/attachment-0001.html
More information about the grass-user
mailing list