We want to:

• Be able to detect the AOI from the DTM.
• Be able to set the resolution from the GSD rather than the DTM resolution
• Be able to export a sequence of local images for memory efficiency (automatically)

I compared 3 methods to convert geodetic coordinates to ECEF coordinates and noticed some precision issues with convertLatLonToECEF.

Input code snippet:

double lat = 47.37315016988094;
double lon = 7.361895523802324;
double alt = 646.5192265328203;

convertLatLonToECEF(lat, lon, alt, WGS84_Ellipsoid);

Results are:

4.27644e+06
552521
4.68484e+06

Where as pymap3d and pyproj both give the same results:

4292018.66822882 554533.5989806354 4670434.035306832

I will investigate this further and see where the discrepancy comes from.

The best options to not slow everything down too much would be to compute a "normal" vector for each pixel indicating the directions the model extend from the level 0, as most of the assumptions we made for the different algorithms would still be valid like that.