[gdal-dev] gdal_merge.py on large Tiff's. gdalwarp not working properly.

aphunter a1hunter at ucsd.edu
Sun Jul 8 10:48:14 PDT 2012


Definitely, thanks for the help. I am sort of new to this stuff. Below is the
imagery from the small images I'm testing with.

Here is the info for the reprojected image (EPSG26910 -> EPSG26911)

Driver: GTiff/GeoTIFF
Files: 26910_rs.tif
Size is 4577, 11883
Coordinate System is:
PROJCS["NAD83 / UTM zone 11N",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-117],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","26911"]]
Origin = (-177789.631458136136644,4942643.042761978693306)
Pixel Size = (100.253523439120713,-100.236459348383733)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( -177789.631, 4942643.043) (125d29'56.62"W, 44d19'11.15"N)
Lower Left  ( -177789.631, 3751533.196) (124d18'18.68"W, 33d41'14.38"N)
Upper Right (  281070.745, 4942643.043) (119d45'31.64"W, 44d36'13.72"N)
Lower Right (  281070.745, 3751533.196) (119d22' 1.78"W, 33d52'53.21"N)
Center      (   51640.557, 4347088.120) (122d11'16.55"W, 39d 9'26.88"N)
Band 1 Block=4577x1 Type=Byte, ColorInterp=Red
Band 2 Block=4577x1 Type=Byte, ColorInterp=Green
Band 3 Block=4577x1 Type=Byte, ColorInterp=Blue


Here is the info for the native EPSG26911 image.


PROJCS["NAD83 / UTM zone 11N",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-117],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","26911"]]
Origin = (222400.000000000000000,4654450.000000000000000)
Pixel Size = (100.012896094325725,-100.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  222400.000, 4654450.000) (120d21' 4.97"W, 41d59'35.32"N)
Lower Left  (  222400.000, 3595550.000) (119d57'11.96"W, 32d27'45.18"N)
Upper Right (  765270.000, 4654450.000) (113d47'50.08"W, 41d59'50.61"N)
Lower Right (  765270.000, 3595550.000) (114d10'39.77"W, 32d27'56.01"N)
Center      (  493835.000, 4125000.000) (117d 4'10.33"W, 37d16'17.60"N)
Band 1 Block=5428x1 Type=Byte, ColorInterp=Red
Band 2 Block=5428x1 Type=Byte, ColorInterp=Green
Band 3 Block=5428x1 Type=Byte, ColorInterp=Blue


--
View this message in context: http://osgeo-org.1560.n6.nabble.com/gdal-merge-py-on-large-Tiff-s-gdalwarp-not-working-properly-tp4986729p4986770.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.


More information about the gdal-dev mailing list