Assert sensible extends in UTM zones. #279

openeo-geotrellis/src/main/scala/org/openeo/geotrellis/layers/FileLayerProvider.scala

@@ -886,6 +886,12 @@ class FileLayerProvider private(openSearch: OpenSearchClient, openSearchCollecti
     val reprojectedBoundingBox: ProjectedExtent = DatacubeSupport.targetBoundingBox(fullBBox, layoutScheme)
     val alignedExtent = worldLayout.createAlignedRasterExtent(reprojectedBoundingBox.extent)
 
+    val polygonIsUTM = polygons_crs.proj4jCrs.getProjection.getName == "utm"
+    if (polygonIsUTM) {
+      // This is an extend that has the highest sensible values for northern and/or southern hemisphere UTM zones
+      val utmProjectedBounds = Extent(166021.44, 0000000.00, 833978.56, 10000000)
+      assert(alignedExtent.extent.intersects(utmProjectedBounds), "Extend should be in valid values of UTM zone to avoid distortion when reprojecting.")
+    }
 
     logger.info(s"Loading ${openSearchCollectionId} with params ${datacubeParams.getOrElse(new DataCubeParameters)} and bands ${openSearchLinkTitles.toList.mkString(";")} initial layout: ${worldLayout}")
 

openeo-geotrellis/src/test/scala/org/openeo/geotrellis/layers/FileLayerProviderTest.scala

@@ -1165,6 +1165,11 @@ class FileLayerProviderTest extends RasterMatchers{
     cubeSpatial.writeGeoTiff(outDir + "/testMissingS2DateLine_" + crsName.replace(":", "_") + ".tiff")
   }
 
+  @Test
+  def testMissingS2DateLineOutside(): Unit = {
+    assertThrows[AssertionError](testMissingS2DateLine("EPSG:32631"))
+  }
+
   private def keysForLargeArea(useBBox:Boolean=false) = {
     val date = LocalDate.of(2022, 2, 11).atStartOfDay(UTC)
     val crs = CRS.fromEpsgCode(32630)