[gdal-dev] Cannot read UTM Ozi .map file
Nik Sands
nixanz at nixanz.com
Thu Nov 29 14:32:13 PST 2012
Thanks for the reply. Firstly, it did resolve my problem, so thanks very much for this information.
However, I don't get that warning with gdalinfo (see below). But I set the environment variable in my code, and it did resolve the problem for me there which is all that matters, I guess. :-)
I'm curious though, how does the quality of an image affect this? I had assumed that transforms would work on on pixel x,y coordinates irrespective of the colour values of those pixels.
Anyhow, thanks again for your help. Much appreciated.
Cheers,
Nik.
----------------------------------
$ env | grep OZI
$ gdalinfo 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"]],
...
...
...
----------------------------------
On 29/11/2012, at 10:08 PM, Jean-Claude Repetto <jrepetto at free.fr> wrote:
> 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
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
More information about the gdal-dev
mailing list