Invalid geometry of feature cause exception in geotrellis.spark.clip.ClipToGrid

EmileSonneveld commented 10 months ago

Reported by Bert De Roo

error. org.locationtech.jts.geom.TopologyException: found non-noded intersection between LINESTRING ( 3125327.8980538994 3571729.6168265673, 3125328.649952571 3571732.4740763605 ) and LINESTRING ( 3125328.627648378 3571732.3893194045, 3126011.1691451804 3571351.0673861904 ) [ (3125328.6276483806, 3571732.389319403, NaN) ]. The workflow is run on different years, however the error only occurs for 2021. Can somebody have a look? job-id:j-231116047d824ec0aa46639f5bbb28d4 backend: openeo3-1.openeo-vlcc-prod.vgt.vito.be

Minimal reproduction graph:

{
  "process_graph": {
    "loadcollection1": {
      "arguments": {
        "bands": [
          "B02",
          "SCL"
        ],
        "id": "SENTINEL2_L2A",
        "properties": {
          "tileId": {
            "process_graph": {
              "eq1": {
                "arguments": {
                  "x": {
                    "from_parameter": "value"
                  },
                  "y": "29UNV"
                },
                "process_id": "eq",
                "result": true
              }
            }
          }
        },
        "spatial_extent": {
          "crs": "EPSG:3035",
          "east": 3140000,
          "north": 3580000,
          "south": 3560000,
          "west": 3120000
        },
        "temporal_extent": [
          "2021-10-01",
          "2021-10-03"
        ]
      },
      "process_id": "load_collection"
    },
    "resamplespatial1": {
      "arguments": {
        "align": "upper-left",
        "data": {
          "from_node": "loadcollection1"
        },
        "method": "near",
        "projection": 3035,
        "resolution": 10
      },
      "process_id": "resample_spatial"
    },
    "saveresult1": {
      "arguments": {
        "data": {
          "from_node": "resamplespatial1"
        },
        "format": "GTiff",
        "options": {
          "filename_prefix": "CROP_2021_E312N356-03035-010m_V100",
          "overview_method": "mode"
        }
      },
      "process_id": "save_result",
      "result": true
    }
  }
}

EmileSonneveld commented 10 months ago

This product contained invalid geometry: /eodata/Sentinel-2/MSI/L2A_N0500/2021/10/02/S2B_MSIL2A_20211002T114349_N0500_R123_T29UNV_20230103T031448.SAFE

jdries commented 9 months ago

@EmileSonneveld did you mention that this was now solved?

EmileSonneveld commented 9 months ago

The MR just needs to get trough

EmileSonneveld commented 9 months ago

@jdries This is merged and went through tests, can i deploy it on vlcc-3-1?

jdries commented 9 months ago

fine for me, thomas also asked for a deploy there

EmileSonneveld commented 9 months ago

This is deployed and works on https://openeo-3-1.openeo-vlcc-prod.vgt.vito.be/

Open-EO / openeo-opensearch-client

Invalid geometry of feature cause exception in geotrellis.spark.clip.ClipToGrid #29