[gdal-dev] Map warping

[gp.saliola at gmail.com]

Hello,

I'm trying to warp a GeoTiff image from its datum/projection to EPSG:32633
(WGS84 UTM33 Nord).

These are the info about the source tiff file:

 

Driver: GTiff/GeoTIFF

Files: C:\Users\gpsal\Desktop\venezia.tif

Size is 10600, 3904

Coordinate System is:

PROJCS["WGS 84 / Pseudo-Mercator",

    GEOGCS["WGS 84",

        DATUM["WGS_1984",

            SPHEROID["WGS 84",6378137,298.257223563,

                AUTHORITY["EPSG","7030"]],

            AUTHORITY["EPSG","6326"]],

        PRIMEM["Greenwich",0,

            AUTHORITY["EPSG","8901"]],

        UNIT["degree",0.0174532925199433,

            AUTHORITY["EPSG","9122"]],

        AUTHORITY["EPSG","4326"]],

    PROJECTION["Mercator_1SP"],

    PARAMETER["central_meridian",0],

    PARAMETER["scale_factor",1],

    PARAMETER["false_easting",0],

    PARAMETER["false_northing",0],

    UNIT["metre",1,

       AUTHORITY["EPSG","9001"]],

    AXIS["Easting",EAST],

    AXIS["Northing",NORTH],

    EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0
+x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs"],

    AUTHORITY["EPSG","3857"]]

Data axis to CRS axis mapping: 1,2

PROJ.4 string is:

'+proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1
+units=m +nadgrids=@null +wktext +no_defs'

Origin = (1335698.850769308628514,5704032.021438719704747)

Pixel Size = (3.356372516969551,-3.345286391264631)

Metadata:

  AREA_OR_POINT=Area

  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)

  TIFFTAG_XRESOLUTION=72

  TIFFTAG_YRESOLUTION=72

Image Structure Metadata:

  COMPRESSION=LZW

  INTERLEAVE=PIXEL

Corner Coordinates:

Upper Left  ( 1335698.851, 5704032.021) ( 11d59'55.63"E, 45d31'18.17"N)

Lower Left  ( 1335698.851, 5690972.023) ( 11d59'55.63"E, 45d26'22.04"N)

Upper Right ( 1371276.399, 5704032.021) ( 12d19' 6.19"E, 45d31'18.17"N)

Lower Right ( 1371276.399, 5690972.023) ( 12d19' 6.19"E, 45d26'22.04"N)

Center      ( 1353487.625, 5697502.022) ( 12d 9'30.91"E, 45d28'50.16"N)

Band 1 Block=10600x1 Type=Byte, ColorInterp=Red

Band 2 Block=10600x1 Type=Byte, ColorInterp=Green

Band 3 Block=10600x1 Type=Byte, ColorInterp=Blue

 

The result is a warped image but georeferencing is bad (difference of some
kilometers).

The problem is in the input map or what?

Thank you for your support.

Best regards,

Gian Paolo Saliola

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20210105/08b57d95/attachment.html>