GCP optimization - UTM zone ValueError

[MatteaE]

I am trying to run the full orthorectification workflow on a single SPOT 5 HRG scene, using a Pléiades orthorectified scene + DEM as reference.
When I call cGCPOptimization(), I get the attached error message. For some reason, there is also a warning "UserWarning: Reproject DEM from 32642-->32641", but the entire domain lies well within the EPSG:32642 region.

gcpOptimization_errorMessage.txt

[MatteaE]

Updated findings: the error happens only if the reference DEM has gaps and/or does not cover the full extent of the SPOT scene. Using a complete global DEM with a large extent, the workflow runs without issues.

[SaifAati]

Great! thank you @MatteaE for the update. I will add a TODO comment to handle NaNs or non full cover DEMs.