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 <[email protected] <mailto:[email protected]>> 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 <[email protected]
    <mailto:[email protected]>> 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
            <[email protected] <mailto:[email protected]>
            <mailto:[email protected]
            <mailto:[email protected]>>> 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
            [email protected]
            <mailto:[email protected]>
            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
[email protected]
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

_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to