I was wondering if this used a Gnomonic projection but the AEQD makes way more sense here (especially if defined in polar coordinates, as I imagine it must be? Then you only need to project the points you actually use?).

Any chance of writing a QGIS plugin with the algorithm?