[gdal-dev] "Too many points" for GDALSuggestedWarpOutput with Ozi .map

[Jean-Claude Repetto]

On 05/12/2012 22:19, Nik Sands wrote:
> Thanks for the reply.  Just to clarify though, can you tell me what it is that is wrong about it?
> 
The GCP coordinates look like UTM coordinates (UTM zone 36) :
Point01,xy,3051,328,in,deg,,0,N,,0,E, grid, 36, 673750, 3572910,N

So the projection is UTM, not Lat/Lon.
gdalinfo reports an error (Invalid angle) :

$ OZI_APPROX_GEOTRANSFORM=YES gdalinfo -noct IL50K-0715_ozf.map
Driver: MAP/OziExplorer .MAP
Files: IL50K-0715_ozf.map
       IL50K-0715.ozf2
Size is 5860, 5860
Coordinate System is:
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"]]
GeoTransform =
  647690.1672605723, 8.519190252396545, 0.1802216446798053
  3575200.856193653, 0.2040104927636261, -8.558172983414893
Corner Coordinates:
Upper Left  (  647690.167, 3575200.856) (Invalid angle,Invalid angle)
Lower Left  (  648746.266, 3525049.963) (Invalid angle,Invalid angle)
Upper Right (  697612.622, 3576396.358) (Invalid angle,Invalid angle)
Lower Right (  698668.721, 3526245.464) (Invalid angle,Invalid angle)
Center      (  673179.444, 3550723.160) (Invalid angle,Invalid angle)
Band 1 Block=64x64 Type=Byte, ColorInterp=Palette
  Overviews: 1465x1465, 586x586, 293x293, 146x146, 59x59
  Color Table (RGB with 256 entries)


If you put the right projection in the .map file, it works :

$ OZI_APPROX_GEOTRANSFORM=YES gdalinfo -noct IL50K-0715_ozf.map
Driver: MAP/OziExplorer .MAP
Files: IL50K-0715_ozf.map
Size is 5860, 5860
Coordinate System is:
PROJCS["UTM Zone 36, Northern Hemisphere",
    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["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",33],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]
GeoTransform =
  647690.1672605723, 8.519190252396545, 0.1802216446798053
  3575200.856193653, 0.2040104927636261, -8.558172983414893
Corner Coordinates:
Upper Left  (  647690.167, 3575200.856) ( 34d34' 7.35"E, 32d18'14.00"N)
Lower Left  (  648746.266, 3525049.963) ( 34d34'19.78"E, 31d51' 5.30"N)
Upper Right (  697612.622, 3576396.358) ( 35d 5'56.21"E, 32d18'25.07"N)
Lower Right (  698668.721, 3526245.464) ( 35d 5'59.25"E, 31d51'16.70"N)
Center      (  673179.444, 3550723.160) ( 34d50' 5.64"E, 32d 4'46.28"N)
Band 1 Block=64x64 Type=Byte, ColorInterp=Palette
  Overviews: 1465x1465, 586x586, 293x293, 146x146, 59x59
  Color Table (RGB with 256 entries)