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 > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev > -- Spatialys - Geospatial professional services http://www.spatialys.com _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
