# [Proj] Map projection used for CADRG polar areas?

Even Rouault even.rouault at spatialys.com
Wed Oct 22 04:54:23 EST 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
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

Even

Le mercredi 22 octobre 2014 11:34:52, Mikael Rittri a écrit :
> Hello,
> . 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
```