fix: add support and test case for spherical to projection transform

Fix issue pdssp#10

Author: jsorel

build.gradle.kts

@@ -31,7 +31,7 @@ dependencies {
     compileOnly("org.jspecify:jspecify:1.0.0")
 
     // Referencing engine
-    implementation("org.apache.sis.core:sis-referencing:1.5.0-ALPHA-1")
+    implementation(libs.apache.sis.referencing)
     implementation("org.apache.sis.non-free:sis-embedded-data:1.3")
     implementation("org.apache.sis.non-free:sis-epsg:1.3")
     implementation("org.apache.derby:derby")
@@ -41,7 +41,7 @@ dependencies {
     // For client
     testFixturesImplementation(platform(libs.geomatys.backend.bom))
     testFixturesImplementation("org.springframework:spring-core")
-    testFixturesApi("org.apache.sis.core:sis-referencing:1.5.0-ALPHA-1")
+    testFixturesApi(libs.apache.sis.referencing)
     testFixturesImplementation("org.graalvm.js:js:24.1.1")
     testFixturesImplementation("org.graalvm.js:js-scriptengine:24.1.1")
     testFixturesImplementation("org.python:jython-slim:2.7.4")

gradle/libs.versions.toml

@@ -2,8 +2,10 @@
 geomatys-backend = "2025.1"
 springdoc = "2.8.9"
 therapi = "0.15.0"
+apache-sis = "1.5.0-ALPHA-2"
 
 [libraries]
+apache-sis-referencing = { module = "org.apache.sis.core:sis-referencing", version.ref = "apache-sis" }
 geomatys-backend-bom = { module = "com.geomatys.backend:downstream-bom", version.ref = "geomatys-backend" }
 springdoc-starter = { module = "org.springdoc:springdoc-openapi-starter-webmvc-ui", version.ref = "springdoc" }
 therapi-runtime = { module = "com.github.therapi:therapi-runtime-javadoc", version.ref = "therapi" }

src/main/java/com/geomatys/crsservice/service/DefaultCrsOperationService.java

@@ -478,6 +478,8 @@ private static String toPythonClass(MathTransform trs) {
         pythonCode = pythonCode.replaceAll("β", "beta");
         //replace ρ
         pythonCode = pythonCode.replaceAll("ρ", "rho");
+        //replace Ω
+        pythonCode = pythonCode.replaceAll("Ω", "omega");
 
         {//replace NaN
             final Pattern objPattern = Pattern.compile("(Number.math.NaN|Number.NaN|NaN)");

src/test/java/com/geomatys/crsservice/client/ClientTest.java

@@ -17,7 +17,7 @@
 
 /**
  *
- * @author jsorel
+ * @author Johann Sorel (Geomatys)
  */
 public class ClientTest extends AbstractIntegrationTest {
 
@@ -43,6 +43,41 @@ public void operationTest() throws Exception {
         testTransform(factory, "EPSG:4326", "EPSG:2154", new double[]{48, 2}, true, cartTolerance, geogTolerance);
         //Geocentric
         testTransform(factory, "EPSG:4326", "EPSG:4978", new double[]{48, 2}, true, cartTolerance, geogTolerance);
+        //Spherical to projection
+        final String mars1 =
+                "GEODCRS[\"Mars (2015) / Ocentric\",\n" +
+                "  DATUM[\"Mars (2015)\",\n" +
+                "    ELLIPSOID[\"Mars (2015)\", 3396190, 169.8944472236118,\n" +
+                "      LENGTHUNIT[\"metre\", 1, ID[\"EPSG\", 9001]]],\n" +
+                "    ANCHOR[\"Viking 1 lander : 47.95137 W\"]],\n" +
+                "    PRIMEM[\"Reference Meridian\", 0,\n" +
+                "      ANGLEUNIT[\"degree\", 0.0174532925199433, ID[\"EPSG\", 9122]]],\n" +
+                "  CS[spherical, 2],\n" +
+                "    AXIS[\"planetocentric latitude (U)\", north,\n" +
+                "      ANGLEUNIT[\"degree\", 0.0174532925199433]],\n" +
+                "    AXIS[\"planetocentric longitude (V)\", east,\n" +
+                "      ANGLEUNIT[\"degree\", 0.0174532925199433]],\n" +
+                "  ID[\"IAU\", 49902, 2015],\n" +
+                "  REMARK[\"Source of IAU Coordinate systems: doi:10.1007/s10569-017-9805-5\"]]";
+        final String mars2 =
+                "PROJCRS[\"Mars (2015) / Ocentric / Equirectangular, clon = 0\",\n" +
+                "  BASEGEODCRS[\"Mars (2015) / Ocentric\",\n" +
+                "    DATUM[\"Mars (2015)\",\n" +
+                "      ELLIPSOID[\"Mars (2015)\", 3396190, 169.8944472236118,\n" +
+                "        LENGTHUNIT[\"metre\", 1, ID[\"EPSG\", 9001]]],\n" +
+                "      ANCHOR[\"Viking 1 lander : 47.95137 W\"]],\n" +
+                "      PRIMEM[\"Reference Meridian\", 0,\n" +
+                "        ANGLEUNIT[\"degree\", 0.0174532925199433, ID[\"EPSG\", 9122]]],\n" +
+                "    ID[\"IAU\", 49902, 2015]],\n" +
+                "  CONVERSION[\"Equirectangular, clon = 0\",\n" +
+                "    METHOD[\"Equidistant Cylindrical\", ID[\"EPSG\", 1028]]],\n" +
+                "  CS[Cartesian, 2],\n" +
+                "    AXIS[\"Easting (E)\", east,\n" +
+                "      LENGTHUNIT[\"metre\", 1]],\n" +
+                "    AXIS[\"Northing (N)\", north,\n" +
+                "      LENGTHUNIT[\"metre\", 1]],\n" +
+                "  ID[\"IAU\", 49912, 2015]]";
+        testTransform(factory, mars1, mars2, new double[]{40, 120}, true, cartTolerance, geogTolerance);
 
         //test PassThroughTransform : Not working yet
         CoordinateReferenceSystem crs1 = CRS.compound(CRS.forCode("CRS:84"));