[gdal-dev] Cannot get spatial reference system from ozi file
Nik Sands
nixanz at nixanz.com
Sun Nov 11 16:04:36 PST 2012
Unfortunately, I did not get any reply to my message about GDAL being unable to georeference an .ozf2 image using an ozi .map file a user sent me (see below). Is there anybody who could take a look and help me out with this?
Please let me know if there is more information that I can provide in order to progress this issue.
Cheers,
Nik.
On 07/11/2012, at 11:11 AM, Nik Sands <nixanz at nixanz.com> wrote:
> Hi all,
>
> My application is reading in ozi file pairs (.ozf2/.map), warping them and saving out as GeoTIFF or other file formats using GDAL. It works fine for my own ozf files, but not for the files sent to me by one of my users (who apparently has a large number of files they want to use).
>
> It has been unable to get a useful spatial reference system from the ozi file pair (.ozf2/.map) sent to me (.map content below). The following code:
>
> ------------------
> OGRSpatialReferenceH srcSRS = OSRNewSpatialReference(pszSrcWKT);
> char *srcWKT;
>
> OSRExportToWkt(srcSRS, &srcWKT);
> printf("\n%s\n\n", srcWKT);
> ------------------
>
> produces the following output:
>
> ------------------
> LOCAL_CS["unnamed",UNIT["unknown",1]]
> ------------------
>
> Can somebody please help me to use GDAL to obtain a valid spatial reference system from this .map file (content below)?
>
> Note that I'm using a fairly recent build of GDAL trunk that includes all the ozi datums, etc, including "Australian Geodetic 1966" and it is able to get a valid SRS from my own files that use the same ozi datum.
>
> The complete content of the ozi .map file in question is below.
>
> Thanks in anticipation,
> Nik.
>
>
> OziExplorer Map Data File Version 2.2
> r2256.jgw
> D:\Documents and Settings\nelsonp\My Documents\General Admin\Personal\Tasmania 1 to 25,000\North West OZF2\r2256.ozf2
> 1 ,Map Code,
> Australian Geodetic 1966,WGS 84, 0.0000, 0.0000,WGS 84
> Reserved 1
> Reserved 2
> Magnetic Variation,,,E
> Map Projection,(UTM) Universal Transverse Mercator,PolyCal,No,AutoCalOnly,No,BSBUseWPX,No
> Point01,xy, 0, 0,in, deg, , ,, , ,, grid, 55,219999.09067205,5570000.7412442,S
> Point02,xy, 1000, 1000,in, deg, , ,, , ,, grid, 55,223259.16350298,5566740.6684133,S
> Point03,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point04,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point05,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point06,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point07,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point08,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point09,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point10,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point11,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point12,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point13,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point14,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point15,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point16,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point17,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point18,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point19,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point20,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point21,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point22,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point23,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point24,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point25,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point26,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point27,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point28,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point29,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Point30,xy, , ,in, deg, , ,, , ,, grid, , , ,
> Projection Setup, 0.000000000, -120.000000000, 0.999600000, 0.00, -4000000.00, 34.000000000, 40.500000000,,,
> Map Feature = MF ; Map Comment = MC These follow if they exist
> Track File = TF These follow if they exist
> Moving Map Parameters = MM? These follow if they exist
> MM0,Yes
> MMPNUM,4
> MMPXY,1,0,0
> MMPXY,2,6136,0
> MMPXY,3,6136,3069
> MMPXY,4,0,3069
> MMPLL,1, 143.721340, -39.973660
> MMPLL,2, 143.955261, -39.980047
> MMPLL,3, 143.951253, -40.070068
> MMPLL,4, 143.717025, -40.063661
> MM1B,3.251882
> MOP,Map Open Position,0,0
> IWH,Map Image Width/Height,6136,3069
>
> _______________________________________________
> 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