[Proj] Map projection used for CADRG polar areas?

Even Rouault even.rouault at spatialys.com
Wed Oct 22 02:54:23 PDT 2014


Mikael,

Based on my code in GDAL, the answer would be Azimuthal Equidistant :
https://github.com/OSGeo/gdal/blob/trunk/gdal/frmts/nitf/nitfdataset.cpp#L1005

But it is perhaps wrong. However, my experience with data showed that the 4 
corner coordinates reprojected from long/lat to aeqd form a rectangle. But 
that might also be true of polar stereographic...

Did you check if the maths given in 30.5.1 "North polar region" are consistant 
with an aeqd method ?

Another clue is that the ASRP specification 
(https://www.dgiwg.org/digest/documents/asrp_1_2.pdf), quite close to CADRG 
regarding the ARC system, mentions : "3.2.3 Polar zones are square arrays of 
tiles overlaid mathematically on a polar azimuthal equidistant projection". 
And its maths at "B.2.2 North Polar Case" look very much the same as the one 
of CADRG.

Even

Le mercredi 22 octobre 2014 11:34:52, Mikael Rittri a écrit :
> Hello,
> I am reading the specification for CADRG,
> http://earth-info.nga.mil/publications/specs/printed/89038/89038_CADRG.pdf
> . But I am confused about the polar areas.
> 
> Section 3.5.5 says: "CADRG shall use a polar stereographic projection".
> 
> And Section 3.5.6 says "Polar frames are defined in polar stereographic
> projection and non-polar frames are in the ARC projection."
> 
> But Section 3.6.2 says: "Pixel coordinates in the polar zones are
> proportional to rectangular coordinates of the Azimuthal
> Equidistant projection, polar aspect, spherical form (as
> defined in Map Projections-A Working Manual, page 191)."
> 
> And section 30.1 says: "Polar zone equations are based on
> the Azimuthal Equidistant projection, polar aspect,
> spherical form."
> 
> Since Polar Stereographic is not the same thing as Azimuthal Equidistant
> in polar aspect, I assume that one of the two alternatives is just a typo
> (unless I am missing something). But which one is the typo of
> "Stereographic" and "Azimuthal Equidistant"? Does anyone know?
> 
> Best regards,
> 
> Mikael Rittri
> Carmenta
> Sweden
> http://www.carmenta.com

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com



More information about the Proj mailing list