ogr2ogr: speed-up reprojection in Arrow code path

autotest/cpp/test_ogr.cpp

@@ -4201,4 +4201,73 @@ TEST_F(test_ogr, OGRCurve_reversePoints)
     }
 }
 
+// Test OGRGeometryFactory::transformWithOptions()
+TEST_F(test_ogr, transformWithOptions)
+{
+    // Projected CRS to national geographic CRS (not including poles or antimeridian)
+    {
+        OGRGeometry *poGeom = nullptr;
+        OGRGeometryFactory::createFromWkt(
+            "LINESTRING(700000 6600000, 700001 6600001)", nullptr, &poGeom);
+        ASSERT_NE(poGeom, nullptr);
+
+        OGRSpatialReference oEPSG_2154;
+        oEPSG_2154.importFromEPSG(2154);  // "RGF93 v1 / Lambert-93"
+        OGRSpatialReference oEPSG_4171;
+        oEPSG_4171.importFromEPSG(4171);  // "RGF93 v1"
+        oEPSG_4171.SetAxisMappingStrategy(OAMS_TRADITIONAL_GIS_ORDER);
+        auto poCT = std::unique_ptr<OGRCoordinateTransformation>(
+            OGRCreateCoordinateTransformation(&oEPSG_2154, &oEPSG_4171));
+        OGRGeometryFactory::TransformWithOptionsCache oCache;
+        poGeom = OGRGeometryFactory::transformWithOptions(poGeom, poCT.get(),
+                                                          nullptr, oCache);
+        EXPECT_NEAR(poGeom->toLineString()->getX(0), 3, 1e-8);
+        EXPECT_NEAR(poGeom->toLineString()->getY(0), 46.5, 1e-8);
+
+        delete poGeom;
+    }
+
+#ifdef HAVE_GEOS
+    // Projected CRS to national geographic CRS including antimeridian
+    {
+        OGRGeometry *poGeom = nullptr;
+        OGRGeometryFactory::createFromWkt(
+            "LINESTRING(657630.64 4984896.17,815261.43 4990738.26)", nullptr,
+            &poGeom);
+        ASSERT_NE(poGeom, nullptr);
+
+        OGRSpatialReference oEPSG_6329;
+        oEPSG_6329.importFromEPSG(6329);  // "NAD83(2011) / UTM zone 60N"
+        OGRSpatialReference oEPSG_6318;
+        oEPSG_6318.importFromEPSG(6318);  // "NAD83(2011)"
+        oEPSG_6318.SetAxisMappingStrategy(OAMS_TRADITIONAL_GIS_ORDER);
+        auto poCT = std::unique_ptr<OGRCoordinateTransformation>(
+            OGRCreateCoordinateTransformation(&oEPSG_6329, &oEPSG_6318));
+        OGRGeometryFactory::TransformWithOptionsCache oCache;
+        poGeom = OGRGeometryFactory::transformWithOptions(poGeom, poCT.get(),
+                                                          nullptr, oCache);
+        EXPECT_EQ(poGeom->getGeometryType(), wkbMultiLineString);
+        if (poGeom->getGeometryType() == wkbMultiLineString)
+        {
+            const auto poMLS = poGeom->toMultiLineString();
+            EXPECT_EQ(poMLS->getNumGeometries(), 2);
+            if (poMLS->getNumGeometries() == 2)
+            {
+                const auto poLS = poMLS->getGeometryRef(0);
+                EXPECT_EQ(poLS->getNumPoints(), 2);
+                if (poLS->getNumPoints() == 2)
+                {
+                    EXPECT_NEAR(poLS->getX(0), 179, 1e-6);
+                    EXPECT_NEAR(poLS->getY(0), 45, 1e-6);
+                    EXPECT_NEAR(poLS->getX(1), 180, 1e-6);
+                    EXPECT_NEAR(poLS->getY(1), 45.004384301691303, 1e-6);
+                }
+            }
+        }
+
+        delete poGeom;
+    }
+#endif
+}
+
 }  // namespace