Skip to content

Integrate new LIA workflow directly done in S2 geometry

Generating XYZ coordinate positions in S1 geometry (through SARDEMProjection + SARCartesianMeanEstimation) doesn't result in normals good enough. Artefacts are observed in rugged areas. See normlim_sigma0#3 (closed)

Another approach seems to yield much better results:

  1. gdalwrap dem + geiod on S2 geometry
  2. BandMath sum them
  3. Generate XYZ ground positions + satellite position on the S2 geometry with new version of SARDEMProjection (normlim_sigma0#9 (closed))
  4. Compute LIA + sin(LIA) from these informations directly in S2 geometry