[gdal-dev] Warping Issues
Even Rouault
even.rouault at spatialys.com
Wed Aug 20 03:43:21 PDT 2014
Martin,
The difference is likely due to a different version of the MrSID SDK used to
compile the GDAL MrSID plugin/driver. According to the version of the SDK, a
different method is used to decode the projection information. I guess that the
one that leads to buggy result is the one that borrows code from the GeoTIFF
driver, that is also known to have an issue with non-meter linear units (as
discussed just yesterday in
http://lists.osgeo.org/pipermail/gdal-dev/2014-August/039827.html).
Even
> Thanks for the quick reply Evan.
>
> I have done what you asked and it gets a bit strange.
> I ran gdalinfo from the QGIS bin directory and then on the alternate machine
> ran gdalinfo from the 'pure gdal' bin directory on the same sid file.
> There are some differences in the output particularly around the definition
> if the PROJCS and in particular the definition of the 'false_easting'.
> The decoded Corner Coordinates are also different probably due to this
> 'false_easting' discrepency.
> Still at a loss to understand the difference between the QGIS GDAL version
> and the built gdal on the other machine.
> Is it a bad build?
> Wrong Proj4 datasets?
> Or just an issue with GDAL/MrSID integration that QGIS seems to have
> overcome? (If they have overcome it then it would be nice to have an
> indication of how)
>
> Output from QGIS gdal info
> -----------------------------------------------------------------------
> Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)
> Files: E:\Stuff\Maps\georgiaville_se.sid
> Size is 5189, 6945
> Coordinate System is:
> *PROJCS["NAD83 / Rhode Island",*
> 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",41.08333333333334],
> PARAMETER["central_meridian",-71.5],
> PARAMETER["scale_factor",0.99999375],
> * PARAMETER["false_easting",328083.3333333333],*
> PARAMETER["false_northing",0],
> * UNIT["US survey foot",0.3048006096012192,
> AUTHORITY["EPSG","9003"]]]*
> Origin = (311060.109583333310000,311248.359216666660000)
> Pixel Size = (3.280833333333334,-3.280833333333334)
> Metadata:
> GEOTIFF_CHAR__GTModelTypeGeoKey=ModelTypeProjected
> GEOTIFF_CHAR__GTRasterTypeGeoKey=RasterPixelIsArea
> GEOTIFF_CHAR__ProjectedCSTypeGeoKey=PCS_NAD83_Rhode_Island
> GEOTIFF_CHAR__ProjLinearUnitsGeoKey=Linear_Foot_US_Survey
> GEOTIFF_NUM__1024__GTModelTypeGeoKey=1
> GEOTIFF_NUM__1025__GTRasterTypeGeoKey=1
> GEOTIFF_NUM__1026__GTCitationGeoKey=IMAGINE GeoTIFF Support
> Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
> @(#)$RCSfile: egtf.c $ $Revision: 1.10.1.6 $ $Date: 2003/10/17 16:36:51EDT $
> Projection Name = State Plane
> Units = feet
> GeoTIFF Units = feet
> GEOTIFF_NUM__3072__ProjectedCSTypeGeoKey=32130
> GEOTIFF_NUM__3073__PCSCitationGeoKey=IMAGINE GeoTIFF Support
> Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
> @(#)$RCSfile: egtf.c $ $Revision: 1.10.1.6 $ $Date: 2003/10/17 16:36:51EDT $
> State Plane Zone -3800
> NAD = 83
> GEOTIFF_NUM__3076__ProjLinearUnitsGeoKey=9003
> IMAGE__COMPRESSION_BLOCK_SIZE=512
> IMAGE__COMPRESSION_GAMMA=2.000000
> IMAGE__COMPRESSION_NLEV=7
> IMAGE__COMPRESSION_VERSION=2,0,0
> IMAGE__COMPRESSION_WEIGHT=2.000000
> IMAGE__CREATION_DATE=Mon Feb 07 13:53:36 2005
>
> IMAGE__INPUT_FILE_SIZE=108169367.000000
> IMAGE__INPUT_FORMAT=GeoTIFF
> IMAGE__INPUT_NAME=E:\Aerials\USDA\TIFS\4107104_se.tif
> IMAGE__TARGET_COMPRESSION_RATIO=25.000000
> VERSION=MG2
> Image Structure Metadata:
> INTERLEAVE=PIXEL
> Corner Coordinates:
> *Upper Left ( 311060.110, 311248.359) ( 71d33'45.24"W, 41d56'15.01"N)
> Lower Left ( 311060.110, 288462.972) ( 71d33'45.02"W, 41d52'29.91"N)
> Upper Right ( 328084.354, 311248.359) ( 71d29'59.99"W, 41d56'15.07"N)
> Lower Right ( 328084.354, 288462.972) ( 71d29'59.99"W, 41d52'29.97"N)
> Center ( 319572.232, 299855.665) ( 71d31'52.56"W, 41d54'22.51"N)*
> Band 1 Block=1024x128 Type=Byte, ColorInterp=Red
> Minimum=23.000, Maximum=253.000, Mean=76.340, StdDev=27.247
> Overviews: 2595x3473, 1298x1737, 649x869, 325x435, 163x218, 82x109, 41x55
> Band 2 Block=1024x128 Type=Byte, ColorInterp=Green
> Minimum=36.000, Maximum=254.000, Mean=92.388, StdDev=23.622
> Overviews: 2595x3473, 1298x1737, 649x869, 325x435, 163x218, 82x109, 41x55
> Band 3 Block=1024x128 Type=Byte, ColorInterp=Blue
> Minimum=46.000, Maximum=254.000, Mean=82.477, StdDev=21.573
> Overviews: 2595x3473, 1298x1737, 649x869, 325x435, 163x218, 82x109, 41x55
>
>
>
>
> Output from 'pure gdal' gdalinfo
> ------------------------------------------------------------------------
> Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)
> Files: C:\Stuff\MrSid\Georgiaville_SE\georgiaville_se.sid
> Size is 5189, 6945
> Coordinate System is:
> *PROJCS["State Plane",*
> 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",41.08333333333334],
> PARAMETER["central_meridian",-71.5],
> PARAMETER["scale_factor",0.99999375],
> * PARAMETER["false_easting",1076386.736111111],*
> PARAMETER["false_northing",0],
> * UNIT["feet",0.3048006096012192]]*
> Origin = (311060.109583333310000,311248.359216666660000)
> Pixel Size = (3.280833333333334,-3.280833333333334)
> Metadata:
> GEOTIFF_CHAR__GTModelTypeGeoKey=ModelTypeProjected
> GEOTIFF_CHAR__GTRasterTypeGeoKey=RasterPixelIsArea
> GEOTIFF_CHAR__ProjectedCSTypeGeoKey=PCS_NAD83_Rhode_Island
> GEOTIFF_CHAR__ProjLinearUnitsGeoKey=Linear_Foot_US_Survey
> GEOTIFF_NUM__1024__GTModelTypeGeoKey=1
> GEOTIFF_NUM__1025__GTRasterTypeGeoKey=1
> GEOTIFF_NUM__1026__GTCitationGeoKey=IMAGINE GeoTIFF Support
> Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
> @(#)$RCSfile: egtf.c $ $Revision: 1.10.1.6 $ $Date: 2003/10/17 16:36:51EDT $
> Projection Name = State Plane
> Units = feet
> GeoTIFF Units = feet
> GEOTIFF_NUM__3072__ProjectedCSTypeGeoKey=32130
> GEOTIFF_NUM__3073__PCSCitationGeoKey=IMAGINE GeoTIFF Support
> Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
> @(#)$RCSfile: egtf.c $ $Revision: 1.10.1.6 $ $Date: 2003/10/17 16:36:51EDT $
> State Plane Zone -3800
> NAD = 83
> GEOTIFF_NUM__3076__ProjLinearUnitsGeoKey=9003
> IMAGE__COMPRESSION_BLOCK_SIZE=512
> IMAGE__COMPRESSION_GAMMA=2.000000
> IMAGE__COMPRESSION_NLEV=7
> IMAGE__COMPRESSION_VERSION=2,0,0
> IMAGE__COMPRESSION_WEIGHT=2.000000
> IMAGE__CREATION_DATE=Mon Feb 07 13:53:36 2005
>
> IMAGE__INPUT_FILE_SIZE=108169367.000000
> IMAGE__INPUT_FORMAT=GeoTIFF
> IMAGE__INPUT_NAME=E:\Aerials\USDA\TIFS\4107104_se.tif
> IMAGE__TARGET_COMPRESSION_RATIO=25.000000
> VERSION=MG2
> Image Structure Metadata:
> INTERLEAVE=PIXEL
> Corner Coordinates:
> *Upper Left ( 311060.110, 311248.359) ( 74d18'40.33"W, 41d54'11.14"N)
> Lower Left ( 311060.110, 288462.972) ( 74d18'30.47"W, 41d50'26.31"N)
> Upper Right ( 328084.354, 311248.359) ( 74d14'55.46"W, 41d54'16.59"N)
> Lower Right ( 328084.354, 288462.972) ( 74d14'45.82"W, 41d50'31.75"N)
> Center ( 319572.232, 299855.665) ( 74d16'43.02"W, 41d52'21.46"N)*
> Band 1 Block=1024x128 Type=Byte, ColorInterp=Red
> Minimum=23.000, Maximum=253.000, Mean=76.340, StdDev=27.247
> Overviews: 2595x3473, 1298x1737, 649x869, 325x435, 163x218, 82x109, 41x55
> Band 2 Block=1024x128 Type=Byte, ColorInterp=Green
> Minimum=36.000, Maximum=254.000, Mean=92.388, StdDev=23.622
> Overviews: 2595x3473, 1298x1737, 649x869, 325x435, 163x218, 82x109, 41x55
> Band 3 Block=1024x128 Type=Byte, ColorInterp=Blue
> Minimum=46.000, Maximum=254.000, Mean=82.477, StdDev=21.573
> Overviews: 2595x3473, 1298x1737, 649x869, 325x435, 163x218, 82x109, 41x55
>
>
>
>
> --
> View this message in context:
> http://osgeo-org.1560.x6.nabble.com/Re-Warping-Issues-tp5157231p5157277.html
> Sent from the GDAL - Dev mailing list archive at Nabble.com.
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
--
Spatialys - Geospatial professional services
http://www.spatialys.com
More information about the gdal-dev
mailing list