Downloaded results: Coordinate system WGS84, pixel-shift between images

[TS84]

From an user (Austrian) perspective

1. image-to-image registration is preferred for further processing in ArcGIS

2. national or European projected coordinate system is preferred to match with national / regional data and to optimize data visualisation due to on-the-fly projection. If data is reprojected, categorical data (e.g. SIAM categories) requires nearest neighbor interpolation.

@grmpfhmbl

[whisperingpixel]

Maybe as clarification: The processing is done in the data cube's projection, the EPSG:4326 output is generated for compatibility reasons with different clients. The output in the original projection can be downloaded by replacing the 4326 EPSG code in the filename.

[grmpfhmbl]

We should check the output array in the inference result. The toolbox simply downloads everything it finds, if I interpret that correctly.

https://github.com/ZGIS/arcpy_toolbox_sen2cube/blob/202b3f722a0a1784e8e296a9dd8c636a9439497e/src/dynamic/execution.py#L406-L410

[grmpfhmbl]

Assuming we don't change the format of the output URL at the moment line 408 could be replaced with something like that

regex = r'(?P<filename>.*\/\w+_)4326\.tiff'
if (output["file_type"] == 'geotiff') and re.match(regex, output["data"]):
  url_epsg3035 = re.sub(regex, r'\g<filename>3035.tiff', output["data"])
  output_url.append(url_epsg3035)
else:
  output_url.append(output["data"])

[whisperingpixel]

Thanks @grmpfhmbl !
in addition to switch between the data cube, the epsg code (here: 3035) can be extracted from the WKT in the srs attribute of the factbase endpoint. Using osgeo, we're using the following code but I don't know whether there are methods in arcpy for that:

from osgeo import osr

def get_crs_from_wkt(wkt):
  srs = osr.SpatialReference()
  srs.ImportFromWkt(wkt)
  srs.AutoIdentifyEPSG()
  crs_code = srs.GetAuthorityCode(None)
  #crs_name = srs.GetAttrValue("PROJCS", 0)
  if crs_code is None:
    crs_code = "0000"
  return crs_code
``´

[grmpfhmbl]

Apparently ArcPy has SpatialReference with loadFromText(string) and can then give you the name.