Le 28/11/2012 22:54, Nik Sands a écrit :
After applying the suggested patch, it no longer errors out in the same way,
however, I now get another error just a little further on in my code.
---------------------------------
Failed to get geo-transform (error 3) for file: Cradle Mountain_ozf.map
---------------------------------
I'm afraid this is a bit beyond me, so if you can help, I'd be very grateful.
Did you notice that gdalinfo reports the same error ?
$ gdalinfo --debug on -noct Cradle\ Mountain_ozf.map
GDAL: GDALLoadOziMapFile(Cradle Mountain_ozf.map) found file, wasn't
able to derive a first order geotransform. Using points as GCPs.
GDAL: GDALOpen(Cradle Mountain_ozf.map, this=0x155a4a0) succeeds as OZI.
Driver: OZI/OziExplorer Image File
GDAL: GDALDefaultOverviews::OverviewScan()
Files: Cradle Mountain_ozf.map
Size is 6535, 11880
Coordinate System is `'
GCP Projection =
PROJCS["UTM Zone 55, Southern Hemisphere",
GEOGCS["GDA94",
DATUM["Geocentric_Datum_of_Australia_1994",
SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6283"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4283"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",147],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",10000000],
UNIT["Meter",1]]
GCP[ 0]: Id=, Info=
(749,50) -> (410000,5397000,0)
GCP[ 1]: Id=, Info=
(4697,38) -> (415000,5397000,0)
GCP[ 2]: Id=, Info=
(4668,11128) -> (415000,5383000,0)
GCP[ 3]: Id=, Info=
(757,11164) -> (410000,5383000,0)
GCP[ 4]: Id=, Info=
(1532,5600) -> (411000,5390000,0)
GCP[ 5]: Id=, Info=
(3112,1611) -> (413000,5395000,0)
GCP[ 6]: Id=, Info=
(3894,7962) -> (414000,5387000,0)
GCP[ 7]: Id=, Info=
(3890,4794) -> (414000,5391000,0)
GCP[ 8]: Id=, Info=
(2318,9563) -> (412000,5385000,0)
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0,11880.0)
Upper Right ( 6535.0, 0.0)
Lower Right ( 6535.0,11880.0)
Center ( 3267.5, 5940.0)
Band 1 Block=64x64 Type=Byte, ColorInterp=Palette
Overviews: 1634x2970, 654x1188, 327x594, 163x297, 65x119
Color Table (RGB with 256 entries)
GDAL: GDALClose(Cradle Mountain_ozf.map, this=0x155a4a0)
This is because the quality of the scan is very bad. Look at the text
"No Camping Allowed", that appears almost twice. You can ask GDAL to
ignore the problem by setting the environment variable
OZI_APPROX_GEOTRANSFORM to the value YES (In this case, you can drop the
patch I sent you yesterday) :
$ OZI_APPROX_GEOTRANSFORM=YES gdalinfo -noct Cradle\ Mountain_ozf.map
Driver: OZI/OziExplorer Image File
Files: Cradle Mountain_ozf.map
Size is 6535, 11880
Coordinate System is:
PROJCS["UTM Zone 55, Southern Hemisphere",
GEOGCS["GDA94",
DATUM["Geocentric_Datum_of_Australia_1994",
SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6283"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4283"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",147],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",10000000],
UNIT["Meter",1]]
GeoTransform =
409038.6971553778, 1.27227951543291, 0.00115086216077703
5397071.561063163, -0.008135856054426157, -1.260559301556319
Corner Coordinates:
Upper Left ( 409038.697, 5397071.561) (145d54'32.19"E, 41d34'22.31"S)
Lower Left ( 409052.369, 5382096.117) (145d54'24.58"E, 41d42'27.83"S)
Upper Right ( 417353.044, 5397018.393) (146d 0'31.13"E, 41d34'27.29"S)
Lower Right ( 417366.716, 5382042.949) (146d 0'24.27"E, 41d42'32.82"S)
Center ( 413202.707, 5389557.255) (145d57'28.05"E, 41d38'27.60"S)
Band 1 Block=64x64 Type=Byte, ColorInterp=Palette
Overviews: 1634x2970, 654x1188, 327x594, 163x297, 65x119
Color Table (RGB with 256 entries)
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
