LCOV - code coverage report
Current view: top level - autotest/cpp - test_ogr.cpp (source / functions) Hit Total Coverage
Test: gdal_filtered.info Lines: 2765 2782 99.4 %
Date: 2026-09-14 13:07:48 Functions: 381 381 100.0 %

          Line data    Source code
       1             : ///////////////////////////////////////////////////////////////////////////////
       2             : //
       3             : // Project:  C++ Test Suite for GDAL/OGR
       4             : // Purpose:  Test general OGR features.
       5             : // Author:   Mateusz Loskot <mateusz@loskot.net>
       6             : //
       7             : ///////////////////////////////////////////////////////////////////////////////
       8             : // Copyright (c) 2006, Mateusz Loskot <mateusz@loskot.net>
       9             : /*
      10             :  * SPDX-License-Identifier: MIT
      11             :  ****************************************************************************/
      12             : 
      13             : #include "gdal_unit_test.h"
      14             : 
      15             : #include "ogr_p.h"
      16             : #include "ogrsf_frmts.h"
      17             : #include "../../ogr/ogrsf_frmts/osm/gpb.h"
      18             : #include "ogr_recordbatch.h"
      19             : #include "ogrlayerarrow.h"
      20             : #include "ogr_geos.h"
      21             : 
      22             : #include <string>
      23             : #include <algorithm>
      24             : #include <cmath>
      25             : #include <fstream>
      26             : #include <limits>
      27             : 
      28             : #ifdef HAVE_SQLITE3
      29             : #include <sqlite3.h>
      30             : #endif
      31             : 
      32             : #include "gtest_include.h"
      33             : 
      34             : namespace
      35             : {
      36             : 
      37             : // Common fixture with test data
      38             : struct test_ogr : public ::testing::Test
      39             : {
      40             :     std::string data_{tut::common::data_basedir};
      41             :     std::string data_tmp_{tut::common::tmp_basedir};
      42             : };
      43             : 
      44             : // Test OGR driver registrar access
      45           4 : TEST_F(test_ogr, GetGDALDriverManager)
      46             : {
      47           1 :     ASSERT_TRUE(nullptr != GetGDALDriverManager());
      48             : }
      49             : 
      50             : template <class T>
      51          16 : void testSpatialReferenceLeakOnCopy(OGRSpatialReference *poSRS)
      52             : {
      53          16 :     ASSERT_EQ(1, poSRS->GetReferenceCount());
      54             :     {
      55             :         int nCurCount;
      56          16 :         int nLastCount = 1;
      57          16 :         T value;
      58          16 :         value.assignSpatialReference(poSRS);
      59          16 :         nCurCount = poSRS->GetReferenceCount();
      60          16 :         ASSERT_GT(nCurCount, nLastCount);
      61          16 :         nLastCount = nCurCount;
      62             : 
      63          16 :         T value2(value);
      64          16 :         nCurCount = poSRS->GetReferenceCount();
      65          16 :         ASSERT_GT(nCurCount, nLastCount);
      66          16 :         nLastCount = nCurCount;
      67             : 
      68          16 :         T value3;
      69          16 :         value3 = value;
      70          16 :         nCurCount = poSRS->GetReferenceCount();
      71          16 :         ASSERT_GT(nCurCount, nLastCount);
      72          16 :         nLastCount = nCurCount;
      73             : 
      74          16 :         value3 = value;
      75             :         // avoid Coverity Scan warning about above assignment being better
      76             :         // replaced with a move.
      77          16 :         EXPECT_NE(value.getSpatialReference(), nullptr);
      78          16 :         ASSERT_EQ(nLastCount, poSRS->GetReferenceCount());
      79             :     }
      80          16 :     ASSERT_EQ(1, poSRS->GetReferenceCount());
      81             : }
      82             : 
      83             : // Test if copy does not leak or double delete the spatial reference
      84           4 : TEST_F(test_ogr, SpatialReference_leak)
      85             : {
      86           1 :     OGRSpatialReference *poSRS = new OGRSpatialReference();
      87           1 :     ASSERT_TRUE(nullptr != poSRS);
      88             : 
      89           1 :     testSpatialReferenceLeakOnCopy<OGRPoint>(poSRS);
      90           1 :     testSpatialReferenceLeakOnCopy<OGRLineString>(poSRS);
      91           1 :     testSpatialReferenceLeakOnCopy<OGRLinearRing>(poSRS);
      92           1 :     testSpatialReferenceLeakOnCopy<OGRCircularString>(poSRS);
      93           1 :     testSpatialReferenceLeakOnCopy<OGRCompoundCurve>(poSRS);
      94           1 :     testSpatialReferenceLeakOnCopy<OGRCurvePolygon>(poSRS);
      95           1 :     testSpatialReferenceLeakOnCopy<OGRPolygon>(poSRS);
      96           1 :     testSpatialReferenceLeakOnCopy<OGRGeometryCollection>(poSRS);
      97           1 :     testSpatialReferenceLeakOnCopy<OGRMultiSurface>(poSRS);
      98           1 :     testSpatialReferenceLeakOnCopy<OGRMultiPolygon>(poSRS);
      99           1 :     testSpatialReferenceLeakOnCopy<OGRMultiPoint>(poSRS);
     100           1 :     testSpatialReferenceLeakOnCopy<OGRMultiCurve>(poSRS);
     101           1 :     testSpatialReferenceLeakOnCopy<OGRMultiLineString>(poSRS);
     102           1 :     testSpatialReferenceLeakOnCopy<OGRTriangle>(poSRS);
     103           1 :     testSpatialReferenceLeakOnCopy<OGRPolyhedralSurface>(poSRS);
     104           1 :     testSpatialReferenceLeakOnCopy<OGRTriangulatedSurface>(poSRS);
     105             : 
     106           1 :     delete poSRS;
     107             : 
     108             :     // Check that assignSpatialReference() works when passed the SRS
     109             :     // object it already owns and whose has a single reference.
     110           1 :     poSRS = new OGRSpatialReference();
     111           2 :     OGRPoint oPoint;
     112           1 :     oPoint.assignSpatialReference(poSRS);
     113           1 :     poSRS->Release();
     114           1 :     oPoint.assignSpatialReference(oPoint.getSpatialReference());
     115             : }
     116             : 
     117             : template <class T> T *make();
     118             : 
     119           6 : template <> OGRPoint *make()
     120             : {
     121           6 :     return new OGRPoint(1.0, 2.0, 3.0);
     122             : }
     123             : 
     124          18 : template <> OGRLineString *make()
     125             : {
     126          18 :     OGRLineString *poLineString = new OGRLineString();
     127             : 
     128          18 :     poLineString->addPoint(1.0, 2.0, 3.0);
     129          18 :     poLineString->addPoint(1.1, 2.1, 3.1);
     130          18 :     poLineString->addPoint(1.2, 2.2, 3.2);
     131             : 
     132          18 :     return poLineString;
     133             : }
     134             : 
     135          22 : template <> OGRLinearRing *make()
     136             : {
     137          22 :     OGRLinearRing *poLinearRing = new OGRLinearRing();
     138             : 
     139          22 :     poLinearRing->addPoint(1.0, 2.0, 3.0);
     140          22 :     poLinearRing->addPoint(1.1, 2.1, 3.1);
     141          22 :     poLinearRing->addPoint(1.2, 2.2, 3.2);
     142          22 :     poLinearRing->addPoint(1.0, 2.0, 3.0);
     143             : 
     144          22 :     return poLinearRing;
     145             : }
     146             : 
     147          14 : template <> OGRCircularString *make()
     148             : {
     149          14 :     OGRCircularString *poCircularString = new OGRCircularString();
     150             : 
     151          14 :     poCircularString->addPoint(1.0, 2.0, 3.0);
     152          14 :     poCircularString->addPoint(1.1, 2.1, 3.1);
     153          14 :     poCircularString->addPoint(1.2, 2.2, 3.2);
     154             : 
     155          14 :     return poCircularString;
     156             : }
     157             : 
     158          12 : template <> OGRCompoundCurve *make()
     159             : {
     160          12 :     OGRCompoundCurve *poCompoundCurve = new OGRCompoundCurve();
     161             : 
     162          12 :     poCompoundCurve->addCurveDirectly(make<OGRLineString>());
     163          12 :     OGRCircularString *poCircularString = make<OGRCircularString>();
     164          12 :     poCircularString->reversePoints();
     165          12 :     poCompoundCurve->addCurveDirectly(poCircularString);
     166             : 
     167          12 :     return poCompoundCurve;
     168             : }
     169             : 
     170           4 : template <> OGRCurvePolygon *make()
     171             : {
     172           4 :     OGRCurvePolygon *poCurvePolygon = new OGRCurvePolygon();
     173             : 
     174           4 :     poCurvePolygon->addRingDirectly(make<OGRCompoundCurve>());
     175           4 :     poCurvePolygon->addRingDirectly(make<OGRCompoundCurve>());
     176             : 
     177           4 :     return poCurvePolygon;
     178             : }
     179             : 
     180           8 : template <> OGRPolygon *make()
     181             : {
     182           8 :     OGRPolygon *poPolygon = new OGRPolygon();
     183             : 
     184           8 :     poPolygon->addRingDirectly(make<OGRLinearRing>());
     185           8 :     poPolygon->addRingDirectly(make<OGRLinearRing>());
     186             : 
     187           8 :     return poPolygon;
     188             : }
     189             : 
     190           2 : template <> OGRGeometryCollection *make()
     191             : {
     192           2 :     OGRGeometryCollection *poCollection = new OGRGeometryCollection();
     193             : 
     194           2 :     poCollection->addGeometryDirectly(make<OGRPoint>());
     195           2 :     poCollection->addGeometryDirectly(make<OGRLinearRing>());
     196             : 
     197           2 :     return poCollection;
     198             : }
     199             : 
     200           2 : template <> OGRMultiSurface *make()
     201             : {
     202           2 :     OGRMultiSurface *poCollection = new OGRMultiSurface();
     203             : 
     204           2 :     poCollection->addGeometryDirectly(make<OGRPolygon>());
     205           2 :     poCollection->addGeometryDirectly(make<OGRCurvePolygon>());
     206             : 
     207           2 :     return poCollection;
     208             : }
     209             : 
     210           2 : template <> OGRMultiPolygon *make()
     211             : {
     212           2 :     OGRMultiPolygon *poCollection = new OGRMultiPolygon();
     213             : 
     214           2 :     poCollection->addGeometryDirectly(make<OGRPolygon>());
     215             : 
     216           2 :     return poCollection;
     217             : }
     218             : 
     219           2 : template <> OGRMultiPoint *make()
     220             : {
     221           2 :     OGRMultiPoint *poCollection = new OGRMultiPoint();
     222             : 
     223           2 :     poCollection->addGeometryDirectly(make<OGRPoint>());
     224             : 
     225           2 :     return poCollection;
     226             : }
     227             : 
     228           2 : template <> OGRMultiCurve *make()
     229             : {
     230           2 :     OGRMultiCurve *poCollection = new OGRMultiCurve();
     231             : 
     232           2 :     poCollection->addGeometryDirectly(make<OGRLineString>());
     233           2 :     poCollection->addGeometryDirectly(make<OGRCompoundCurve>());
     234             : 
     235           2 :     return poCollection;
     236             : }
     237             : 
     238           2 : template <> OGRMultiLineString *make()
     239             : {
     240           2 :     OGRMultiLineString *poCollection = new OGRMultiLineString();
     241             : 
     242           2 :     poCollection->addGeometryDirectly(make<OGRLineString>());
     243           2 :     poCollection->addGeometryDirectly(make<OGRLinearRing>());
     244             : 
     245           2 :     return poCollection;
     246             : }
     247             : 
     248           4 : template <> OGRTriangle *make()
     249             : {
     250           8 :     OGRPoint p1(0, 0), p2(0, 1), p3(1, 1);
     251           8 :     return new OGRTriangle(p1, p2, p3);
     252             : }
     253             : 
     254           2 : template <> OGRTriangulatedSurface *make()
     255             : {
     256           2 :     OGRTriangulatedSurface *poTS = new OGRTriangulatedSurface();
     257           2 :     poTS->addGeometryDirectly(make<OGRTriangle>());
     258           2 :     return poTS;
     259             : }
     260             : 
     261           2 : template <> OGRPolyhedralSurface *make()
     262             : {
     263           2 :     OGRPolyhedralSurface *poPS = new OGRPolyhedralSurface();
     264           2 :     poPS->addGeometryDirectly(make<OGRPolygon>());
     265           2 :     return poPS;
     266             : }
     267             : 
     268          16 : template <class T> void testCopyEquals()
     269             : {
     270          16 :     auto poOrigin = std::unique_ptr<T>(make<T>());
     271          16 :     ASSERT_TRUE(nullptr != poOrigin);
     272             : 
     273          16 :     T value2(*poOrigin);
     274             : 
     275          16 :     ASSERT_TRUE(poOrigin->Equals(&value2))
     276           0 :         << poOrigin->getGeometryName() << ": copy constructor changed a value";
     277             : 
     278          16 :     T value3;
     279          16 :     value3 = *poOrigin;
     280          16 :     value3 = *poOrigin;
     281          16 :     auto &value3Ref(value3);
     282          16 :     value3 = value3Ref;
     283             : 
     284             : #ifdef DEBUG_VERBOSE
     285             :     char *wkt1 = NULL, *wkt2 = NULL;
     286             :     poOrigin->exportToWkt(&wkt1);
     287             :     value3.exportToWkt(&wkt2);
     288             :     printf("%s %s\n", wkt1, wkt2);
     289             :     CPLFree(wkt1);
     290             :     CPLFree(wkt2);
     291             : #endif
     292          16 :     ASSERT_TRUE(poOrigin->Equals(&value3))
     293           0 :         << poOrigin->getGeometryName()
     294           0 :         << ": assignment operator changed a value";
     295             : 
     296          16 :     value3 = T();
     297          16 :     ASSERT_TRUE(value3.IsEmpty());
     298             : }
     299             : 
     300             : // Test if copy constructor and assignment operators succeeds on copying the
     301             : // geometry data
     302           4 : TEST_F(test_ogr, SpatialReference_leak_copy_constructor)
     303             : {
     304           1 :     testCopyEquals<OGRPoint>();
     305           1 :     testCopyEquals<OGRLineString>();
     306           1 :     testCopyEquals<OGRLinearRing>();
     307           1 :     testCopyEquals<OGRCircularString>();
     308           1 :     testCopyEquals<OGRCompoundCurve>();
     309           1 :     testCopyEquals<OGRCurvePolygon>();
     310           1 :     testCopyEquals<OGRPolygon>();
     311           1 :     testCopyEquals<OGRGeometryCollection>();
     312           1 :     testCopyEquals<OGRMultiSurface>();
     313           1 :     testCopyEquals<OGRMultiPolygon>();
     314           1 :     testCopyEquals<OGRMultiPoint>();
     315           1 :     testCopyEquals<OGRMultiCurve>();
     316           1 :     testCopyEquals<OGRMultiLineString>();
     317           1 :     testCopyEquals<OGRTriangle>();
     318           1 :     testCopyEquals<OGRPolyhedralSurface>();
     319           1 :     testCopyEquals<OGRTriangulatedSurface>();
     320           1 : }
     321             : 
     322             : // Test crazy usage of OGRGeometryCollection copy constructor
     323           4 : TEST_F(test_ogr, OGRGeometryCollection_copy_constructor_illegal_use)
     324             : {
     325           2 :     OGRGeometryCollection gc;
     326           1 :     gc.addGeometryDirectly(new OGRPoint(1, 2));
     327             : 
     328           2 :     OGRMultiPolygon mp;
     329           1 :     mp.addGeometryDirectly(new OGRPolygon());
     330             : 
     331           1 :     OGRGeometryCollection *mp_as_gc = &mp;
     332           1 :     CPLErrorReset();
     333             :     {
     334           2 :         CPLErrorHandlerPusher oPusher(CPLQuietErrorHandler);
     335           1 :         *mp_as_gc = gc;
     336             :         // avoid Coverity Scan warning
     337           1 :         EXPECT_EQ(gc.getSpatialReference(), nullptr);
     338             :     }
     339           1 :     EXPECT_STREQ(CPLGetLastErrorMsg(),
     340             :                  "Illegal use of OGRGeometryCollection::operator=(): trying to "
     341             :                  "assign an incompatible sub-geometry");
     342           1 :     EXPECT_TRUE(mp.IsEmpty());
     343           1 : }
     344             : 
     345             : // Test crazy usage of OGRCurvePolygon copy constructor
     346           4 : TEST_F(test_ogr, OGRCurvePolygon_copy_constructor_illegal_use)
     347             : {
     348           2 :     OGRCurvePolygon cp;
     349           1 :     auto poCC = new OGRCircularString();
     350           1 :     poCC->addPoint(0, 0);
     351           1 :     poCC->addPoint(1, 1);
     352           1 :     poCC->addPoint(2, 0);
     353           1 :     poCC->addPoint(1, -1);
     354           1 :     poCC->addPoint(0, 0);
     355           1 :     cp.addRingDirectly(poCC);
     356             : 
     357           2 :     OGRPolygon poly;
     358           1 :     auto poLR = new OGRLinearRing();
     359           1 :     poLR->addPoint(0, 0);
     360           1 :     poLR->addPoint(1, 1);
     361           1 :     poLR->addPoint(2, 0);
     362           1 :     poLR->addPoint(1, -1);
     363           1 :     poLR->addPoint(0, 0);
     364           1 :     poly.addRingDirectly(poLR);
     365             : 
     366           1 :     OGRCurvePolygon *poly_as_cp = &poly;
     367           1 :     CPLErrorReset();
     368             :     {
     369           2 :         CPLErrorHandlerPusher oPusher(CPLQuietErrorHandler);
     370           1 :         *poly_as_cp = cp;
     371             :         // avoid Coverity Scan warning
     372           1 :         EXPECT_EQ(cp.getSpatialReference(), nullptr);
     373             :     }
     374           1 :     EXPECT_STREQ(CPLGetLastErrorMsg(),
     375             :                  "Illegal use of OGRCurvePolygon::operator=(): trying to "
     376             :                  "assign an incompatible sub-geometry");
     377           1 :     EXPECT_TRUE(poly.IsEmpty());
     378           1 : }
     379             : 
     380          16 : template <class T> void testMove()
     381             : {
     382          16 :     auto poSRS = new OGRSpatialReference();
     383             :     {
     384          16 :         auto poOrigin = std::unique_ptr<T>(make<T>());
     385          16 :         ASSERT_TRUE(nullptr != poOrigin);
     386          16 :         poOrigin->assignSpatialReference(poSRS);
     387             : 
     388          16 :         T valueCopy(*poOrigin);
     389          16 :         const int refCountBefore = poSRS->GetReferenceCount();
     390          16 :         T fromMoved(std::move(*poOrigin));
     391          16 :         EXPECT_EQ(poSRS->GetReferenceCount(), refCountBefore);
     392             : 
     393          16 :         ASSERT_TRUE(fromMoved.Equals(&valueCopy))
     394           0 :             << valueCopy.getGeometryName()
     395           0 :             << ": move constructor changed a value";
     396          16 :         EXPECT_EQ(fromMoved.getSpatialReference(), poSRS);
     397             : 
     398          16 :         T valueCopy2(valueCopy);
     399          16 :         EXPECT_EQ(valueCopy.getSpatialReference(), poSRS);
     400          16 :         T value3;
     401          16 :         const int refCountBefore2 = poSRS->GetReferenceCount();
     402          16 :         value3 = std::move(valueCopy);
     403          16 :         EXPECT_EQ(poSRS->GetReferenceCount(), refCountBefore2);
     404             : 
     405          16 :         ASSERT_TRUE(value3.Equals(&valueCopy2))
     406           0 :             << valueCopy2.getGeometryName()
     407           0 :             << ": move assignment operator changed a value";
     408          16 :         EXPECT_EQ(value3.getSpatialReference(), poSRS);
     409             :     }
     410          16 :     EXPECT_EQ(poSRS->GetReferenceCount(), 1);
     411          16 :     poSRS->Release();
     412             : }
     413             : 
     414           4 : TEST_F(test_ogr, geometry_move)
     415             : {
     416           1 :     testMove<OGRPoint>();
     417           1 :     testMove<OGRLineString>();
     418           1 :     testMove<OGRLinearRing>();
     419           1 :     testMove<OGRCircularString>();
     420           1 :     testMove<OGRCompoundCurve>();
     421           1 :     testMove<OGRCurvePolygon>();
     422           1 :     testMove<OGRPolygon>();
     423           1 :     testMove<OGRGeometryCollection>();
     424           1 :     testMove<OGRMultiSurface>();
     425           1 :     testMove<OGRMultiPolygon>();
     426           1 :     testMove<OGRMultiPoint>();
     427           1 :     testMove<OGRMultiCurve>();
     428           1 :     testMove<OGRMultiLineString>();
     429           1 :     testMove<OGRTriangle>();
     430           1 :     testMove<OGRPolyhedralSurface>();
     431           1 :     testMove<OGRTriangulatedSurface>();
     432           1 : }
     433             : 
     434           4 : TEST_F(test_ogr, geometry_get_point)
     435             : {
     436             :     {
     437           1 :         OGRPoint p;
     438           1 :         double x = 1, y = 2;
     439           1 :         OGR_G_SetPoints((OGRGeometryH)&p, 1, &x, 0, &y, 0, nullptr, 0);
     440           1 :         ASSERT_EQ(p.getCoordinateDimension(), 2);
     441           1 :         ASSERT_EQ(p.getX(), 1);
     442           1 :         ASSERT_EQ(p.getY(), 2);
     443           1 :         ASSERT_EQ(p.getZ(), 0);
     444             :     }
     445             : 
     446             :     {
     447           1 :         OGRPoint p;
     448           1 :         double x = 1, y = 2, z = 3;
     449           1 :         OGR_G_SetPoints((OGRGeometryH)&p, 1, &x, 0, &y, 0, &z, 0);
     450           1 :         ASSERT_EQ(p.getCoordinateDimension(), 3);
     451           1 :         ASSERT_EQ(p.getX(), 1);
     452           1 :         ASSERT_EQ(p.getY(), 2);
     453           1 :         ASSERT_EQ(p.getZ(), 3);
     454             :     }
     455             : 
     456             :     {
     457           2 :         OGRPoint p;
     458           1 :         CPLPushErrorHandler(CPLQuietErrorHandler);
     459           1 :         OGR_G_SetPoints((OGRGeometryH)&p, 1, nullptr, 0, nullptr, 0, nullptr,
     460             :                         0);
     461           1 :         CPLPopErrorHandler();
     462             :     }
     463             : 
     464             :     {
     465           1 :         OGRLineString ls;
     466           1 :         double x = 1, y = 2;
     467           1 :         OGR_G_SetPoints((OGRGeometryH)&ls, 1, &x, 0, &y, 0, nullptr, 0);
     468           1 :         ASSERT_EQ(ls.getCoordinateDimension(), 2);
     469           1 :         ASSERT_EQ(ls.getX(0), 1);
     470           1 :         ASSERT_EQ(ls.getY(0), 2);
     471           1 :         ASSERT_EQ(ls.getZ(0), 0);
     472             :     }
     473             : 
     474             :     {
     475           1 :         OGRLineString ls;
     476           1 :         double x = 1, y = 2;
     477           1 :         OGR_G_SetPoints((OGRGeometryH)&ls, 1, &x, 0, &y, 0, nullptr, 0);
     478           1 :         ASSERT_EQ(ls.getCoordinateDimension(), 2);
     479           1 :         ASSERT_EQ(ls.getX(0), 1);
     480           1 :         ASSERT_EQ(ls.getY(0), 2);
     481           1 :         ASSERT_EQ(ls.getZ(0), 0);
     482             :     }
     483             : 
     484             :     {
     485           1 :         OGRLineString ls;
     486           1 :         double x = 1, y = 2;
     487           1 :         OGR_G_SetPoints((OGRGeometryH)&ls, 1, &x, 8, &y, 8, nullptr, 0);
     488           1 :         ASSERT_EQ(ls.getCoordinateDimension(), 2);
     489           1 :         ASSERT_EQ(ls.getX(0), 1);
     490           1 :         ASSERT_EQ(ls.getY(0), 2);
     491           1 :         ASSERT_EQ(ls.getZ(0), 0);
     492             :     }
     493             : 
     494             :     {
     495           1 :         OGRLineString ls;
     496           1 :         double x = 1, y = 2, z = 3;
     497           1 :         OGR_G_SetPoints((OGRGeometryH)&ls, 1, &x, 0, &y, 0, &z, 0);
     498           1 :         ASSERT_EQ(ls.getCoordinateDimension(), 3);
     499           1 :         ASSERT_EQ(ls.getX(0), 1);
     500           1 :         ASSERT_EQ(ls.getY(0), 2);
     501           1 :         ASSERT_EQ(ls.getZ(0), 3);
     502             :     }
     503             : 
     504             :     {
     505           1 :         OGRLineString ls;
     506           1 :         double x = 1, y = 2, z = 3;
     507           1 :         OGR_G_SetPoints((OGRGeometryH)&ls, 1, &x, 8, &y, 8, &z, 8);
     508           1 :         ASSERT_EQ(ls.getCoordinateDimension(), 3);
     509           1 :         ASSERT_EQ(ls.getX(0), 1);
     510           1 :         ASSERT_EQ(ls.getY(0), 2);
     511           1 :         ASSERT_EQ(ls.getZ(0), 3);
     512             :     }
     513             : 
     514             :     {
     515           2 :         OGRLineString ls;
     516           1 :         CPLPushErrorHandler(CPLQuietErrorHandler);
     517           1 :         OGR_G_SetPoints((OGRGeometryH)&ls, 1, nullptr, 0, nullptr, 0, nullptr,
     518             :                         0);
     519           1 :         CPLPopErrorHandler();
     520             :     }
     521             : }
     522             : 
     523           4 : TEST_F(test_ogr, OGR_G_CreateGeometry_unknown)
     524             : {
     525           1 :     EXPECT_EQ(OGR_G_CreateGeometry(wkbUnknown), nullptr);
     526           1 : }
     527             : 
     528           4 : TEST_F(test_ogr, style_manager)
     529             : {
     530           1 :     OGRStyleMgrH hSM = OGR_SM_Create(nullptr);
     531           1 :     EXPECT_TRUE(OGR_SM_InitStyleString(
     532             :         hSM, "PEN(w:2px,c:#000000,id:\"mapinfo-pen-2,ogr-pen-0\")"));
     533           1 :     OGRStyleToolH hTool = OGR_SM_GetPart(hSM, 0, nullptr);
     534           1 :     EXPECT_TRUE(hTool != nullptr);
     535           1 :     if (hTool)
     536             :     {
     537             :         int bValueIsNull;
     538             : 
     539           1 :         EXPECT_NEAR(OGR_ST_GetParamDbl(hTool, OGRSTPenWidth, &bValueIsNull),
     540             :                     2.0 * (1.0 / (72.0 * 39.37)) * 1000, 1e-6);
     541           1 :         EXPECT_EQ(OGR_ST_GetUnit(hTool), OGRSTUMM);
     542             : 
     543           1 :         OGR_ST_SetUnit(hTool, OGRSTUPixel, 1.0);
     544           1 :         EXPECT_EQ(OGR_ST_GetParamDbl(hTool, OGRSTPenWidth, &bValueIsNull), 2.0);
     545           1 :         EXPECT_EQ(OGR_ST_GetUnit(hTool), OGRSTUPixel);
     546           1 :         OGR_ST_Destroy(hTool);
     547             :     }
     548             : 
     549           1 :     OGR_SM_Destroy(hSM);
     550           1 : }
     551             : 
     552           4 : TEST_F(test_ogr, OGRParseDate)
     553             : {
     554             :     const auto OGRParseDateWrapper =
     555          64 :         [](const char *str, OGRField *psField, int nFlags)
     556             :     {
     557             :         // Putting inside a std::string helps Valgrind and other checkers to
     558             :         // detect out-of-bounds access
     559          64 :         return OGRParseDate(std::string(str).c_str(), psField, nFlags);
     560             :     };
     561             : 
     562             :     OGRField sField;
     563           1 :     EXPECT_EQ(OGRParseDateWrapper("2017/11/31 12:34:56", &sField, 0), TRUE);
     564           1 :     EXPECT_EQ(sField.Date.Year, 2017);
     565           1 :     EXPECT_EQ(sField.Date.Month, 11);
     566           1 :     EXPECT_EQ(sField.Date.Day, 31);
     567           1 :     EXPECT_EQ(sField.Date.Hour, 12);
     568           1 :     EXPECT_EQ(sField.Date.Minute, 34);
     569           1 :     EXPECT_EQ(sField.Date.Second, 56.0f);
     570           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     571             : 
     572           1 :     EXPECT_EQ(OGRParseDateWrapper("2017/11/31 12:34:56+00", &sField, 0), TRUE);
     573           1 :     EXPECT_EQ(sField.Date.TZFlag, 100);
     574             : 
     575           1 :     EXPECT_EQ(OGRParseDateWrapper("2017/11/31 12:34:56+12:00", &sField, 0),
     576             :               TRUE);
     577           1 :     EXPECT_EQ(sField.Date.TZFlag, 100 + 12 * 4);
     578             : 
     579           1 :     EXPECT_EQ(OGRParseDateWrapper("2017/11/31 12:34:56+1200", &sField, 0),
     580             :               TRUE);
     581           1 :     EXPECT_EQ(sField.Date.TZFlag, 100 + 12 * 4);
     582             : 
     583           1 :     EXPECT_EQ(OGRParseDateWrapper("2017/11/31 12:34:56+815", &sField, 0), TRUE);
     584           1 :     EXPECT_EQ(sField.Date.TZFlag, 100 + 8 * 4 + 1);
     585             : 
     586           1 :     EXPECT_EQ(OGRParseDateWrapper("2017/11/31 12:34:56-12:00", &sField, 0),
     587             :               TRUE);
     588           1 :     EXPECT_EQ(sField.Date.TZFlag, 100 - 12 * 4);
     589             : 
     590           1 :     EXPECT_EQ(OGRParseDateWrapper(" 2017/11/31 12:34:56", &sField, 0), TRUE);
     591           1 :     EXPECT_EQ(sField.Date.Year, 2017);
     592             : 
     593           1 :     EXPECT_EQ(OGRParseDateWrapper("2017/11/31 12:34:56.789", &sField, 0), TRUE);
     594           1 :     EXPECT_EQ(sField.Date.Second, 56.789f);
     595             : 
     596             :     // Leap second
     597           1 :     EXPECT_EQ(OGRParseDateWrapper("2017/11/31 12:34:60", &sField, 0), TRUE);
     598           1 :     EXPECT_EQ(sField.Date.Second, 60.0f);
     599             : 
     600           1 :     EXPECT_EQ(OGRParseDateWrapper("2017-11-31T12:34:56", &sField, 0), TRUE);
     601           1 :     EXPECT_EQ(sField.Date.Year, 2017);
     602           1 :     EXPECT_EQ(sField.Date.Month, 11);
     603           1 :     EXPECT_EQ(sField.Date.Day, 31);
     604           1 :     EXPECT_EQ(sField.Date.Hour, 12);
     605           1 :     EXPECT_EQ(sField.Date.Minute, 34);
     606           1 :     EXPECT_EQ(sField.Date.Second, 56.0f);
     607           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     608             : 
     609           1 :     EXPECT_EQ(OGRParseDateWrapper("2017-11-31T12:34:56Z", &sField, 0), TRUE);
     610           1 :     EXPECT_EQ(sField.Date.Second, 56.0f);
     611           1 :     EXPECT_EQ(sField.Date.TZFlag, 100);
     612             : 
     613           1 :     EXPECT_EQ(OGRParseDateWrapper("2017-11-31T12:34:56.789Z", &sField, 0),
     614             :               TRUE);
     615           1 :     EXPECT_EQ(sField.Date.Second, 56.789f);
     616           1 :     EXPECT_EQ(sField.Date.TZFlag, 100);
     617             : 
     618           1 :     EXPECT_EQ(OGRParseDateWrapper("2017-11-31T23:59:59.999999Z", &sField, 0),
     619             :               TRUE);
     620           1 :     EXPECT_EQ(sField.Date.Hour, 23);
     621           1 :     EXPECT_EQ(sField.Date.Minute, 59);
     622           1 :     EXPECT_EQ(sField.Date.Second, 59.999f);
     623             : 
     624           1 :     EXPECT_EQ(OGRParseDateWrapper("2017-11-31", &sField, 0), TRUE);
     625           1 :     EXPECT_EQ(sField.Date.Year, 2017);
     626           1 :     EXPECT_EQ(sField.Date.Month, 11);
     627           1 :     EXPECT_EQ(sField.Date.Day, 31);
     628           1 :     EXPECT_EQ(sField.Date.Hour, 0);
     629           1 :     EXPECT_EQ(sField.Date.Minute, 0);
     630           1 :     EXPECT_EQ(sField.Date.Second, 0.0f);
     631           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     632             : 
     633           1 :     EXPECT_EQ(OGRParseDateWrapper("2017-11-31Z", &sField, 0), TRUE);
     634           1 :     EXPECT_EQ(sField.Date.Year, 2017);
     635           1 :     EXPECT_EQ(sField.Date.Month, 11);
     636           1 :     EXPECT_EQ(sField.Date.Day, 31);
     637           1 :     EXPECT_EQ(sField.Date.Hour, 0);
     638           1 :     EXPECT_EQ(sField.Date.Minute, 0);
     639           1 :     EXPECT_EQ(sField.Date.Second, 0.0f);
     640           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     641             : 
     642           1 :     EXPECT_EQ(OGRParseDateWrapper("12:34", &sField, 0), TRUE);
     643           1 :     EXPECT_EQ(sField.Date.Year, 0);
     644           1 :     EXPECT_EQ(sField.Date.Month, 0);
     645           1 :     EXPECT_EQ(sField.Date.Day, 0);
     646           1 :     EXPECT_EQ(sField.Date.Hour, 12);
     647           1 :     EXPECT_EQ(sField.Date.Minute, 34);
     648           1 :     EXPECT_EQ(sField.Date.Second, 0.0f);
     649           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     650             : 
     651           1 :     EXPECT_EQ(OGRParseDateWrapper("12:34:56", &sField, 0), TRUE);
     652           1 :     EXPECT_EQ(OGRParseDateWrapper("12:34:56.789", &sField, 0), TRUE);
     653             : 
     654           1 :     EXPECT_EQ(OGRParseDateWrapper("T12:34:56", &sField, 0), TRUE);
     655           1 :     EXPECT_EQ(sField.Date.Year, 0);
     656           1 :     EXPECT_EQ(sField.Date.Month, 0);
     657           1 :     EXPECT_EQ(sField.Date.Day, 0);
     658           1 :     EXPECT_EQ(sField.Date.Hour, 12);
     659           1 :     EXPECT_EQ(sField.Date.Minute, 34);
     660           1 :     EXPECT_EQ(sField.Date.Second, 56.0f);
     661           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     662             : 
     663           1 :     EXPECT_EQ(OGRParseDateWrapper("T123456", &sField, 0), TRUE);
     664           1 :     EXPECT_EQ(sField.Date.Year, 0);
     665           1 :     EXPECT_EQ(sField.Date.Month, 0);
     666           1 :     EXPECT_EQ(sField.Date.Day, 0);
     667           1 :     EXPECT_EQ(sField.Date.Hour, 12);
     668           1 :     EXPECT_EQ(sField.Date.Minute, 34);
     669           1 :     EXPECT_EQ(sField.Date.Second, 56.0f);
     670           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     671             : 
     672           1 :     EXPECT_EQ(OGRParseDateWrapper("T123456.789", &sField, 0), TRUE);
     673           1 :     EXPECT_EQ(sField.Date.Year, 0);
     674           1 :     EXPECT_EQ(sField.Date.Month, 0);
     675           1 :     EXPECT_EQ(sField.Date.Day, 0);
     676           1 :     EXPECT_EQ(sField.Date.Hour, 12);
     677           1 :     EXPECT_EQ(sField.Date.Minute, 34);
     678           1 :     EXPECT_EQ(sField.Date.Second, 56.789f);
     679           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     680             : 
     681           1 :     CPLPushErrorHandler(CPLQuietErrorHandler);
     682           1 :     EXPECT_TRUE(!OGRParseDateWrapper("123456-01-01", &sField, 0));
     683           1 :     CPLPopErrorHandler();
     684           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017", &sField, 0));
     685           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017x-01-01", &sField, 0));
     686           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-1-01", &sField, 0));
     687           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-1", &sField, 0));
     688           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-01x", &sField, 0));
     689           1 :     EXPECT_TRUE(!OGRParseDateWrapper("12:", &sField, 0));
     690           1 :     EXPECT_TRUE(!OGRParseDateWrapper("12:3", &sField, 0));
     691           1 :     EXPECT_TRUE(!OGRParseDateWrapper("1:23", &sField, 0));
     692           1 :     EXPECT_TRUE(!OGRParseDateWrapper("12:34:5", &sField, 0));
     693           1 :     EXPECT_TRUE(!OGRParseDateWrapper("1a:34", &sField, 0));
     694           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-a-31T12:34:56", &sField, 0));
     695           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-00-31T12:34:56", &sField, 0));
     696           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-13-31T12:34:56", &sField, 0));
     697           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-00T12:34:56", &sField, 0));
     698           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-aT12:34:56", &sField, 0));
     699           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-32T12:34:56", &sField, 0));
     700           1 :     EXPECT_TRUE(!OGRParseDateWrapper("a:34:56", &sField, 0));
     701           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-01Ta:34:56", &sField, 0));
     702           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-01T25:34:56", &sField, 0));
     703           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-01T00:a:00", &sField, 0));
     704           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-01T00: 34:56", &sField, 0));
     705           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-01T00:61:00", &sField, 0));
     706           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-01T00:00:61", &sField, 0));
     707           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-01-01T00:00:a", &sField, 0));
     708             : 
     709             :     // Test OGRPARSEDATE_OPTION_LAX
     710           1 :     EXPECT_EQ(OGRParseDateWrapper("2017-1-9", &sField, OGRPARSEDATE_OPTION_LAX),
     711             :               TRUE);
     712           1 :     EXPECT_EQ(sField.Date.Year, 2017);
     713           1 :     EXPECT_EQ(sField.Date.Month, 1);
     714           1 :     EXPECT_EQ(sField.Date.Day, 9);
     715             : 
     716           1 :     EXPECT_EQ(
     717             :         OGRParseDateWrapper("2017-1-31", &sField, OGRPARSEDATE_OPTION_LAX),
     718             :         TRUE);
     719           1 :     EXPECT_EQ(sField.Date.Year, 2017);
     720           1 :     EXPECT_EQ(sField.Date.Month, 1);
     721           1 :     EXPECT_EQ(sField.Date.Day, 31);
     722             : 
     723           1 :     EXPECT_EQ(OGRParseDateWrapper("2017-1-31T1:2:3", &sField,
     724             :                                   OGRPARSEDATE_OPTION_LAX),
     725             :               TRUE);
     726           1 :     EXPECT_EQ(sField.Date.Year, 2017);
     727           1 :     EXPECT_EQ(sField.Date.Month, 1);
     728           1 :     EXPECT_EQ(sField.Date.Day, 31);
     729           1 :     EXPECT_EQ(sField.Date.Hour, 1);
     730           1 :     EXPECT_EQ(sField.Date.Minute, 2);
     731           1 :     EXPECT_EQ(sField.Date.Second, 3.0f);
     732           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     733             : 
     734           1 :     EXPECT_EQ(
     735             :         OGRParseDateWrapper("2017-1-31T1:3", &sField, OGRPARSEDATE_OPTION_LAX),
     736             :         TRUE);
     737           1 :     EXPECT_EQ(sField.Date.Year, 2017);
     738           1 :     EXPECT_EQ(sField.Date.Month, 1);
     739           1 :     EXPECT_EQ(sField.Date.Day, 31);
     740           1 :     EXPECT_EQ(sField.Date.Hour, 1);
     741           1 :     EXPECT_EQ(sField.Date.Minute, 3);
     742           1 :     EXPECT_EQ(sField.Date.Second, 0.0f);
     743           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     744             : 
     745           1 :     EXPECT_EQ(OGRParseDateWrapper("1:3", &sField, OGRPARSEDATE_OPTION_LAX),
     746             :               TRUE);
     747           1 :     EXPECT_EQ(sField.Date.Year, 0);
     748           1 :     EXPECT_EQ(sField.Date.Month, 0);
     749           1 :     EXPECT_EQ(sField.Date.Day, 0);
     750           1 :     EXPECT_EQ(sField.Date.Hour, 1);
     751           1 :     EXPECT_EQ(sField.Date.Minute, 3);
     752           1 :     EXPECT_EQ(sField.Date.Second, 0.0f);
     753           1 :     EXPECT_EQ(sField.Date.TZFlag, 0);
     754             : 
     755           1 :     EXPECT_TRUE(
     756             :         !OGRParseDateWrapper("2017-a-01", &sField, OGRPARSEDATE_OPTION_LAX));
     757           1 :     EXPECT_TRUE(
     758             :         !OGRParseDateWrapper("2017-0-01", &sField, OGRPARSEDATE_OPTION_LAX));
     759           1 :     EXPECT_TRUE(
     760             :         !OGRParseDateWrapper("2017-1", &sField, OGRPARSEDATE_OPTION_LAX));
     761           1 :     EXPECT_TRUE(
     762             :         !OGRParseDateWrapper("2017-1-", &sField, OGRPARSEDATE_OPTION_LAX));
     763           1 :     EXPECT_TRUE(
     764             :         !OGRParseDateWrapper("2017-1-a", &sField, OGRPARSEDATE_OPTION_LAX));
     765           1 :     EXPECT_TRUE(
     766             :         !OGRParseDateWrapper("2017-1-0", &sField, OGRPARSEDATE_OPTION_LAX));
     767           1 :     EXPECT_TRUE(
     768             :         !OGRParseDateWrapper("2017-1-32", &sField, OGRPARSEDATE_OPTION_LAX));
     769           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-1-1Ta:00:00", &sField,
     770             :                                      OGRPARSEDATE_OPTION_LAX));
     771           1 :     EXPECT_TRUE(
     772             :         !OGRParseDateWrapper("2017-1-1T1", &sField, OGRPARSEDATE_OPTION_LAX));
     773           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-1-1T00:a:00", &sField,
     774             :                                      OGRPARSEDATE_OPTION_LAX));
     775           1 :     EXPECT_TRUE(
     776             :         !OGRParseDateWrapper("2017-1-1T1:", &sField, OGRPARSEDATE_OPTION_LAX));
     777           1 :     EXPECT_TRUE(!OGRParseDateWrapper("2017-1-1T00:00:a", &sField,
     778             :                                      OGRPARSEDATE_OPTION_LAX));
     779           1 :     EXPECT_TRUE(!OGRParseDateWrapper("1a:3", &sField, OGRPARSEDATE_OPTION_LAX));
     780           1 : }
     781             : 
     782             : // Test OGRPolygon::IsPointOnSurface()
     783           4 : TEST_F(test_ogr, IsPointOnSurface)
     784             : {
     785           1 :     OGRPolygon oPoly;
     786             : 
     787           1 :     OGRPoint oEmptyPoint;
     788           1 :     ASSERT_TRUE(!oPoly.IsPointOnSurface(&oEmptyPoint));
     789             : 
     790           1 :     OGRPoint oPoint;
     791           1 :     oPoint.setX(1);
     792           1 :     oPoint.setY(1);
     793           1 :     ASSERT_TRUE(!oPoly.IsPointOnSurface(&oPoint));
     794             : 
     795           1 :     const char *pszPoly =
     796             :         "POLYGON((0 0,0 10,10 10,10 0,0 0),(4 4,4 6,6 6,6 4,4 4))";
     797           1 :     oPoly.importFromWkt(&pszPoly);
     798             : 
     799           1 :     ASSERT_TRUE(!oPoly.IsPointOnSurface(&oEmptyPoint));
     800             : 
     801           1 :     ASSERT_EQ(oPoly.IsPointOnSurface(&oPoint), TRUE);
     802             : 
     803           1 :     oPoint.setX(5);
     804           1 :     oPoint.setY(5);
     805           1 :     ASSERT_TRUE(!oPoly.IsPointOnSurface(&oPoint));
     806             : }
     807             : 
     808             : // Test gpb.h
     809           4 : TEST_F(test_ogr, gpb_h)
     810             : {
     811           1 :     ASSERT_EQ(GetVarUIntSize(0), 1);
     812           1 :     ASSERT_EQ(GetVarUIntSize(127), 1);
     813           1 :     ASSERT_EQ(GetVarUIntSize(128), 2);
     814           1 :     ASSERT_EQ(GetVarUIntSize((1 << 14) - 1), 2);
     815           1 :     ASSERT_EQ(GetVarUIntSize(1 << 14), 3);
     816           1 :     ASSERT_EQ(GetVarUIntSize(GUINT64_MAX), 10);
     817             : 
     818           1 :     ASSERT_EQ(GetVarIntSize(0), 1);
     819           1 :     ASSERT_EQ(GetVarIntSize(127), 1);
     820           1 :     ASSERT_EQ(GetVarIntSize(128), 2);
     821           1 :     ASSERT_EQ(GetVarIntSize((1 << 14) - 1), 2);
     822           1 :     ASSERT_EQ(GetVarIntSize(1 << 14), 3);
     823           1 :     ASSERT_EQ(GetVarIntSize(GINT64_MAX), 9);
     824           1 :     ASSERT_EQ(GetVarIntSize(-1), 10);
     825           1 :     ASSERT_EQ(GetVarIntSize(GINT64_MIN), 10);
     826             : 
     827           1 :     ASSERT_EQ(GetVarSIntSize(0), 1);
     828           1 :     ASSERT_EQ(GetVarSIntSize(63), 1);
     829           1 :     ASSERT_EQ(GetVarSIntSize(64), 2);
     830           1 :     ASSERT_EQ(GetVarSIntSize(-1), 1);
     831           1 :     ASSERT_EQ(GetVarSIntSize(-64), 1);
     832           1 :     ASSERT_EQ(GetVarSIntSize(-65), 2);
     833           1 :     ASSERT_EQ(GetVarSIntSize(GINT64_MIN), 10);
     834           1 :     ASSERT_EQ(GetVarSIntSize(GINT64_MAX), 10);
     835             : 
     836           1 :     ASSERT_EQ(GetTextSize(""), 1);
     837           1 :     ASSERT_EQ(GetTextSize(" "), 2);
     838           1 :     ASSERT_EQ(GetTextSize(std::string(" ")), 2);
     839             : 
     840           1 :     GByte abyBuffer[11] = {0};
     841             :     GByte *pabyBuffer;
     842             :     const GByte *pabyBufferRO;
     843             : 
     844           1 :     pabyBuffer = abyBuffer;
     845           1 :     WriteVarUInt(&pabyBuffer, 0);
     846           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 1);
     847           1 :     pabyBufferRO = abyBuffer;
     848           1 :     ASSERT_EQ(ReadVarUInt64(&pabyBufferRO), 0U);
     849             : 
     850           1 :     pabyBuffer = abyBuffer;
     851           1 :     WriteVarUInt(&pabyBuffer, 127);
     852           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 1);
     853           1 :     pabyBufferRO = abyBuffer;
     854           1 :     ASSERT_EQ(ReadVarUInt64(&pabyBufferRO), 127U);
     855             : 
     856           1 :     pabyBuffer = abyBuffer;
     857           1 :     WriteVarUInt(&pabyBuffer, 0xDEADBEEFU);
     858           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 5);
     859           1 :     pabyBufferRO = abyBuffer;
     860           1 :     ASSERT_EQ(ReadVarUInt64(&pabyBufferRO), 0xDEADBEEFU);
     861             : 
     862           1 :     pabyBuffer = abyBuffer;
     863           1 :     WriteVarUInt(&pabyBuffer, GUINT64_MAX);
     864           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 10);
     865           1 :     pabyBufferRO = abyBuffer;
     866           1 :     ASSERT_EQ(ReadVarUInt64(&pabyBufferRO), GUINT64_MAX);
     867             : 
     868           1 :     pabyBuffer = abyBuffer;
     869           1 :     WriteVarInt(&pabyBuffer, GINT64_MAX);
     870           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 9);
     871           1 :     pabyBufferRO = abyBuffer;
     872           1 :     ASSERT_EQ(ReadVarInt64(&pabyBufferRO), GINT64_MAX);
     873             : 
     874           1 :     pabyBuffer = abyBuffer;
     875           1 :     WriteVarInt(&pabyBuffer, -1);
     876           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 10);
     877           1 :     pabyBufferRO = abyBuffer;
     878           1 :     ASSERT_EQ(ReadVarInt64(&pabyBufferRO), -1);
     879             : 
     880           1 :     pabyBuffer = abyBuffer;
     881           1 :     WriteVarInt(&pabyBuffer, GINT64_MIN);
     882           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 10);
     883           1 :     pabyBufferRO = abyBuffer;
     884           1 :     ASSERT_EQ(ReadVarInt64(&pabyBufferRO), GINT64_MIN);
     885             : 
     886           1 :     pabyBuffer = abyBuffer;
     887           1 :     WriteVarSInt(&pabyBuffer, 0);
     888           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 1);
     889             :     {
     890             :         GIntBig nVal;
     891           1 :         pabyBufferRO = abyBuffer;
     892           1 :         READ_VARSINT64(pabyBufferRO, abyBuffer + 10, nVal);
     893           1 :         ASSERT_EQ(nVal, 0);
     894             :     }
     895             : 
     896           1 :     pabyBuffer = abyBuffer;
     897           1 :     WriteVarSInt(&pabyBuffer, 1);
     898           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 1);
     899             :     {
     900             :         GIntBig nVal;
     901           1 :         pabyBufferRO = abyBuffer;
     902           1 :         READ_VARSINT64(pabyBufferRO, abyBuffer + 10, nVal);
     903           1 :         ASSERT_EQ(nVal, 1);
     904             :     }
     905             : 
     906           1 :     pabyBuffer = abyBuffer;
     907           1 :     WriteVarSInt(&pabyBuffer, -1);
     908           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 1);
     909             :     {
     910             :         GIntBig nVal;
     911           1 :         pabyBufferRO = abyBuffer;
     912           1 :         READ_VARSINT64(pabyBufferRO, abyBuffer + 10, nVal);
     913           1 :         ASSERT_EQ(nVal, -1);
     914             :     }
     915             : 
     916           1 :     pabyBuffer = abyBuffer;
     917           1 :     WriteVarSInt(&pabyBuffer, GINT64_MAX);
     918           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 10);
     919             :     {
     920             :         GIntBig nVal;
     921           1 :         pabyBufferRO = abyBuffer;
     922           1 :         READ_VARSINT64(pabyBufferRO, abyBuffer + 10, nVal);
     923           1 :         ASSERT_EQ(nVal, GINT64_MAX);
     924             :     }
     925             : 
     926           1 :     pabyBuffer = abyBuffer;
     927           1 :     WriteVarSInt(&pabyBuffer, GINT64_MIN);
     928           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 10);
     929             :     {
     930             :         GIntBig nVal;
     931           1 :         pabyBufferRO = abyBuffer;
     932           1 :         READ_VARSINT64(pabyBufferRO, abyBuffer + 10, nVal);
     933           1 :         ASSERT_EQ(nVal, GINT64_MIN);
     934             :     }
     935             : 
     936           1 :     pabyBuffer = abyBuffer;
     937           1 :     WriteText(&pabyBuffer, "x");
     938           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 2);
     939           1 :     ASSERT_EQ(abyBuffer[0], 1);
     940           1 :     ASSERT_EQ(abyBuffer[1], 'x');
     941             : 
     942           1 :     pabyBuffer = abyBuffer;
     943           1 :     WriteText(&pabyBuffer, std::string("x"));
     944           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 2);
     945           1 :     ASSERT_EQ(abyBuffer[0], 1);
     946           1 :     ASSERT_EQ(abyBuffer[1], 'x');
     947             : 
     948           1 :     pabyBuffer = abyBuffer;
     949           1 :     WriteFloat32(&pabyBuffer, 1.25f);
     950           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 4);
     951           1 :     pabyBufferRO = abyBuffer;
     952           1 :     ASSERT_EQ(ReadFloat32(&pabyBufferRO, abyBuffer + 4), 1.25f);
     953             : 
     954           1 :     pabyBuffer = abyBuffer;
     955           1 :     WriteFloat64(&pabyBuffer, 1.25);
     956           1 :     ASSERT_EQ(pabyBuffer - abyBuffer, 8);
     957           1 :     pabyBufferRO = abyBuffer;
     958           1 :     ASSERT_EQ(ReadFloat64(&pabyBufferRO, abyBuffer + 8), 1.25);
     959             : }
     960             : 
     961             : // Test OGRGeometry::toXXXXX()
     962           4 : TEST_F(test_ogr, OGRGeometry_toXXXXX)
     963             : {
     964             : #define CONCAT(X, Y) X##Y
     965             : #define TEST_OGRGEOMETRY_TO(X)                                                 \
     966             :     {                                                                          \
     967             :         CONCAT(OGR, X) o;                                                      \
     968             :         OGRGeometry *poGeom = &o;                                              \
     969             :         ASSERT_EQ(poGeom->CONCAT(to, X)(), &o);                                \
     970             :     }
     971             : 
     972           1 :     TEST_OGRGEOMETRY_TO(Point);
     973           1 :     TEST_OGRGEOMETRY_TO(LineString);
     974           1 :     TEST_OGRGEOMETRY_TO(LinearRing);
     975           1 :     TEST_OGRGEOMETRY_TO(CircularString);
     976           1 :     TEST_OGRGEOMETRY_TO(CompoundCurve);
     977           1 :     TEST_OGRGEOMETRY_TO(CurvePolygon);
     978           1 :     TEST_OGRGEOMETRY_TO(Polygon);
     979           1 :     TEST_OGRGEOMETRY_TO(GeometryCollection);
     980           1 :     TEST_OGRGEOMETRY_TO(MultiSurface);
     981           1 :     TEST_OGRGEOMETRY_TO(MultiPolygon);
     982           1 :     TEST_OGRGEOMETRY_TO(MultiPoint);
     983           1 :     TEST_OGRGEOMETRY_TO(MultiCurve);
     984           1 :     TEST_OGRGEOMETRY_TO(MultiLineString);
     985           1 :     TEST_OGRGEOMETRY_TO(Triangle);
     986           1 :     TEST_OGRGEOMETRY_TO(PolyhedralSurface);
     987           1 :     TEST_OGRGEOMETRY_TO(TriangulatedSurface);
     988             :     {
     989           1 :         OGRLineString o;
     990           1 :         OGRGeometry *poGeom = &o;
     991           1 :         ASSERT_EQ(poGeom->toCurve(), &o);
     992             :     }
     993             :     {
     994           1 :         OGRPolygon o;
     995           1 :         OGRGeometry *poGeom = &o;
     996           1 :         ASSERT_EQ(poGeom->toSurface(), &o);
     997             :     }
     998             : 
     999             :     {
    1000           1 :         OGRPoint o;
    1001             :         // ASSERT_EQ(o.toPoint(), &o);
    1002             :     }
    1003             : 
    1004             :     {
    1005           1 :         OGRLineString o;
    1006           1 :         ASSERT_EQ(o.toCurve(), &o);
    1007           1 :         ASSERT_EQ(o.toSimpleCurve(), &o);
    1008             :         // ASSERT_EQ(o.toLineString(), &o);
    1009             : 
    1010             :         {
    1011           1 :             OGRCurve &oRef = o;
    1012           1 :             ASSERT_EQ(oRef.toLineString(), &o);
    1013             :         }
    1014             : 
    1015             :         {
    1016           1 :             OGRSimpleCurve &oRef = o;
    1017           1 :             ASSERT_EQ(oRef.toLineString(), &o);
    1018             :         }
    1019             :     }
    1020             : 
    1021             :     {
    1022           1 :         OGRLinearRing o;
    1023           1 :         ASSERT_EQ(o.toCurve(), &o);
    1024           1 :         ASSERT_EQ(o.toSimpleCurve(), &o);
    1025             :         // ASSERT_EQ(o.toLinearRing(), &o);
    1026             : 
    1027             :         {
    1028           1 :             OGRCurve &oRef = o;
    1029           1 :             ASSERT_EQ(oRef.toLinearRing(), &o);
    1030             :         }
    1031             :         {
    1032           1 :             OGRSimpleCurve &oRef = o;
    1033           1 :             ASSERT_EQ(oRef.toLinearRing(), &o);
    1034             :         }
    1035             :         {
    1036           1 :             OGRLineString &oRef = o;
    1037           1 :             ASSERT_EQ(oRef.toLinearRing(), &o);
    1038             :         }
    1039             :     }
    1040             : 
    1041             :     {
    1042           1 :         OGRCircularString o;
    1043           1 :         ASSERT_EQ(o.toCurve(), &o);
    1044           1 :         ASSERT_EQ(o.toSimpleCurve(), &o);
    1045             :         // ASSERT_EQ(o.toCircularString(), &o);
    1046             : 
    1047             :         {
    1048           1 :             OGRCurve &oRef = o;
    1049           1 :             ASSERT_EQ(oRef.toCircularString(), &o);
    1050             :         }
    1051             : 
    1052             :         {
    1053           1 :             OGRSimpleCurve &oRef = o;
    1054           1 :             ASSERT_EQ(oRef.toCircularString(), &o);
    1055             :         }
    1056             :     }
    1057             : 
    1058             :     {
    1059           1 :         OGRCompoundCurve o;
    1060           1 :         ASSERT_EQ(o.toCurve(), &o);
    1061             :         // ASSERT_EQ(o.toCompoundCurve(), &o);
    1062             : 
    1063             :         {
    1064           1 :             OGRCurve &oRef = o;
    1065           1 :             ASSERT_EQ(oRef.toCompoundCurve(), &o);
    1066             :         }
    1067             :     }
    1068             : 
    1069             :     {
    1070           1 :         OGRCurvePolygon o;
    1071           1 :         ASSERT_EQ(o.toSurface(), &o);
    1072             :         // ASSERT_EQ(o.toCurvePolygon(), &o);
    1073             : 
    1074             :         {
    1075           1 :             OGRSurface &oRef = o;
    1076           1 :             ASSERT_EQ(oRef.toCurvePolygon(), &o);
    1077             :         }
    1078             :     }
    1079             : 
    1080             :     {
    1081           1 :         OGRPolygon o;
    1082           1 :         ASSERT_EQ(o.toSurface(), &o);
    1083           1 :         ASSERT_EQ(o.toCurvePolygon(), &o);
    1084             :         // ASSERT_EQ(o.toPolygon(), &o);
    1085             : 
    1086             :         {
    1087           1 :             OGRSurface &oRef = o;
    1088           1 :             ASSERT_EQ(oRef.toPolygon(), &o);
    1089             :         }
    1090             : 
    1091             :         {
    1092           1 :             OGRCurvePolygon &oRef = o;
    1093           1 :             ASSERT_EQ(oRef.toPolygon(), &o);
    1094             :         }
    1095             :     }
    1096             : 
    1097             :     {
    1098           1 :         OGRTriangle o;
    1099           1 :         ASSERT_EQ(o.toSurface(), &o);
    1100           1 :         ASSERT_EQ(o.toCurvePolygon(), &o);
    1101           1 :         ASSERT_EQ(o.toPolygon(), &o);
    1102             :         // ASSERT_EQ(o.toTriangle(), &o);
    1103             : 
    1104             :         {
    1105           1 :             OGRSurface &oRef = o;
    1106           1 :             ASSERT_EQ(oRef.toTriangle(), &o);
    1107             :         }
    1108             : 
    1109             :         {
    1110           1 :             OGRCurvePolygon &oRef = o;
    1111           1 :             ASSERT_EQ(oRef.toTriangle(), &o);
    1112             :         }
    1113             : 
    1114             :         {
    1115           1 :             OGRPolygon &oRef = o;
    1116           1 :             ASSERT_EQ(oRef.toTriangle(), &o);
    1117             :         }
    1118             :     }
    1119             : 
    1120             :     {
    1121           1 :         OGRMultiPoint o;
    1122           1 :         ASSERT_EQ(o.toGeometryCollection(), &o);
    1123             :         // ASSERT_EQ(o.toMultiPoint(), &o);
    1124             : 
    1125             :         {
    1126           1 :             OGRGeometryCollection &oRef = o;
    1127           1 :             ASSERT_EQ(oRef.toMultiPoint(), &o);
    1128             :         }
    1129             :     }
    1130             : 
    1131             :     {
    1132           1 :         OGRMultiCurve o;
    1133           1 :         ASSERT_EQ(o.toGeometryCollection(), &o);
    1134             :         // ASSERT_EQ(o.toMultiCurve(), &o);
    1135             : 
    1136             :         {
    1137           1 :             OGRGeometryCollection &oRef = o;
    1138           1 :             ASSERT_EQ(oRef.toMultiCurve(), &o);
    1139             :         }
    1140             :     }
    1141             : 
    1142             :     {
    1143           1 :         OGRMultiLineString o;
    1144           1 :         ASSERT_EQ(o.toGeometryCollection(), &o);
    1145           1 :         ASSERT_EQ(o.toMultiCurve(), &o);
    1146             :         // ASSERT_EQ(o.toMultiLineString(), &o);
    1147             : 
    1148             :         {
    1149           1 :             OGRMultiCurve &oRef = o;
    1150           1 :             ASSERT_EQ(oRef.toMultiLineString(), &o);
    1151             :         }
    1152             : 
    1153             :         {
    1154           1 :             OGRGeometryCollection &oRef = o;
    1155           1 :             ASSERT_EQ(oRef.toMultiLineString(), &o);
    1156             :         }
    1157             :     }
    1158             : 
    1159             :     {
    1160           1 :         OGRMultiSurface o;
    1161           1 :         ASSERT_EQ(o.toGeometryCollection(), &o);
    1162             :         // ASSERT_EQ(o.toMultiSurface(), &o);
    1163             : 
    1164             :         {
    1165           1 :             OGRGeometryCollection &oRef = o;
    1166           1 :             ASSERT_EQ(oRef.toMultiSurface(), &o);
    1167             :         }
    1168             :     }
    1169             : 
    1170             :     {
    1171           1 :         OGRMultiPolygon o;
    1172           1 :         ASSERT_EQ(o.toGeometryCollection(), &o);
    1173           1 :         ASSERT_EQ(o.toMultiSurface(), &o);
    1174             :         // ASSERT_EQ(o.toMultiPolygon(), &o);
    1175             : 
    1176             :         {
    1177           1 :             OGRMultiSurface &oRef = o;
    1178           1 :             ASSERT_EQ(oRef.toMultiPolygon(), &o);
    1179             :         }
    1180             : 
    1181             :         {
    1182           1 :             OGRGeometryCollection &oRef = o;
    1183           1 :             ASSERT_EQ(oRef.toMultiPolygon(), &o);
    1184             :         }
    1185             :     }
    1186             : 
    1187             :     {
    1188           1 :         OGRPolyhedralSurface o;
    1189           1 :         ASSERT_EQ(o.toSurface(), &o);
    1190             :         // ASSERT_EQ(o.toPolyhedralSurface(), &o);
    1191             : 
    1192             :         {
    1193           1 :             OGRSurface &oRef = o;
    1194           1 :             ASSERT_EQ(oRef.toPolyhedralSurface(), &o);
    1195             :         }
    1196             :     }
    1197             : 
    1198             :     {
    1199           1 :         OGRTriangulatedSurface o;
    1200           1 :         ASSERT_EQ(o.toSurface(), &o);
    1201           1 :         ASSERT_EQ(o.toPolyhedralSurface(), &o);
    1202             :         // ASSERT_EQ(o.toTriangulatedSurface(), &o);
    1203             : 
    1204             :         {
    1205           1 :             OGRSurface &oRef = o;
    1206           1 :             ASSERT_EQ(oRef.toTriangulatedSurface(), &o);
    1207             :         }
    1208             : 
    1209             :         {
    1210           1 :             OGRPolyhedralSurface &oRef = o;
    1211           1 :             ASSERT_EQ(oRef.toTriangulatedSurface(), &o);
    1212             :         }
    1213             :     }
    1214             : }
    1215             : 
    1216          34 : template <typename T> void TestIterator(T *obj, int nExpectedPointCount)
    1217             : {
    1218          34 :     int nCount = 0;
    1219          71 :     for (auto &elt : obj)
    1220             :     {
    1221          25 :         nCount++;
    1222          25 :         CPL_IGNORE_RET_VAL(elt);
    1223             :     }
    1224          34 :     ASSERT_EQ(nCount, nExpectedPointCount);
    1225             : 
    1226          34 :     nCount = 0;
    1227          34 :     const T *const_obj(obj);
    1228          71 :     for (const auto &elt : const_obj)
    1229             :     {
    1230          25 :         nCount++;
    1231          25 :         CPL_IGNORE_RET_VAL(elt);
    1232             :     }
    1233          34 :     ASSERT_EQ(nCount, nExpectedPointCount);
    1234             : }
    1235             : 
    1236             : template <typename Concrete, typename Abstract = Concrete>
    1237          34 : void TestIterator(const char *pszWKT = nullptr, int nExpectedPointCount = 0)
    1238             : {
    1239          68 :     Concrete obj;
    1240          34 :     if (pszWKT)
    1241             :     {
    1242          18 :         obj.importFromWkt(&pszWKT);
    1243             :     }
    1244          34 :     TestIterator<Abstract>(&obj, nExpectedPointCount);
    1245          34 : }
    1246             : 
    1247             : // Test geometry visitor
    1248           4 : TEST_F(test_ogr, OGRGeometry_visitor)
    1249             : {
    1250             :     static const struct
    1251             :     {
    1252             :         const char *pszWKT;
    1253             :         int nExpectedPointCount;
    1254             :     } asTests[] = {
    1255             :         {"POINT(0 0)", 1},
    1256             :         {"LINESTRING(0 0)", 1},
    1257             :         {"POLYGON((0 0),(0 0))", 2},
    1258             :         {"MULTIPOINT(0 0)", 1},
    1259             :         {"MULTILINESTRING((0 0))", 1},
    1260             :         {"MULTIPOLYGON(((0 0)))", 1},
    1261             :         {"GEOMETRYCOLLECTION(POINT(0 0))", 1},
    1262             :         {"CIRCULARSTRING(0 0,1 1,0 0)", 3},
    1263             :         {"COMPOUNDCURVE((0 0,1 1))", 2},
    1264             :         {"CURVEPOLYGON((0 0,1 1,1 0,0 0))", 4},
    1265             :         {"MULTICURVE((0 0))", 1},
    1266             :         {"MULTISURFACE(((0 0)))", 1},
    1267             :         {"TRIANGLE((0 0,0 1,1 1,0 0))", 4},
    1268             :         {"POLYHEDRALSURFACE(((0 0,0 1,1 1,0 0)))", 4},
    1269             :         {"TIN(((0 0,0 1,1 1,0 0)))", 4},
    1270             :     };
    1271             : 
    1272             :     class PointCounterVisitor : public OGRDefaultGeometryVisitor
    1273             :     {
    1274             :         int m_nPoints = 0;
    1275             : 
    1276             :       public:
    1277          15 :         PointCounterVisitor()
    1278          15 :         {
    1279          15 :         }
    1280             : 
    1281             :         using OGRDefaultGeometryVisitor::visit;
    1282             : 
    1283          31 :         void visit(OGRPoint *) override
    1284             :         {
    1285          31 :             m_nPoints++;
    1286          31 :         }
    1287             : 
    1288          15 :         int getNumPoints() const
    1289             :         {
    1290          15 :             return m_nPoints;
    1291             :         }
    1292             :     };
    1293             : 
    1294             :     class PointCounterConstVisitor : public OGRDefaultConstGeometryVisitor
    1295             :     {
    1296             :         int m_nPoints = 0;
    1297             : 
    1298             :       public:
    1299          15 :         PointCounterConstVisitor()
    1300          15 :         {
    1301          15 :         }
    1302             : 
    1303             :         using OGRDefaultConstGeometryVisitor::visit;
    1304             : 
    1305          31 :         void visit(const OGRPoint *) override
    1306             :         {
    1307          31 :             m_nPoints++;
    1308          31 :         }
    1309             : 
    1310          15 :         int getNumPoints() const
    1311             :         {
    1312          15 :             return m_nPoints;
    1313             :         }
    1314             :     };
    1315             : 
    1316          16 :     for (size_t i = 0; i < CPL_ARRAYSIZE(asTests); i++)
    1317             :     {
    1318          15 :         OGRGeometry *poGeom = nullptr;
    1319          15 :         OGRGeometryFactory::createFromWkt(asTests[i].pszWKT, nullptr, &poGeom);
    1320          15 :         PointCounterVisitor oVisitor;
    1321          15 :         poGeom->accept(&oVisitor);
    1322          15 :         ASSERT_EQ(oVisitor.getNumPoints(), asTests[i].nExpectedPointCount);
    1323          15 :         PointCounterConstVisitor oConstVisitor;
    1324          15 :         poGeom->accept(&oConstVisitor);
    1325          15 :         ASSERT_EQ(oConstVisitor.getNumPoints(), asTests[i].nExpectedPointCount);
    1326          15 :         delete poGeom;
    1327             :     }
    1328             : 
    1329             :     {
    1330           2 :         OGRLineString ls;
    1331           1 :         ls.setNumPoints(2);
    1332           2 :         auto oIter1 = ls.begin();
    1333           1 :         EXPECT_TRUE(oIter1 != ls.end());
    1334           1 :         EXPECT_TRUE(!(oIter1 != ls.begin()));
    1335           2 :         auto oIter2 = ls.begin();
    1336           1 :         EXPECT_TRUE(!(oIter1 != oIter2));
    1337           1 :         ++oIter2;
    1338           1 :         EXPECT_TRUE(oIter1 != oIter2);
    1339           1 :         ++oIter2;
    1340           1 :         EXPECT_TRUE(oIter1 != oIter2);
    1341             :     }
    1342             : 
    1343             :     {
    1344           2 :         OGRLineString ls;
    1345           1 :         EXPECT_TRUE(!(ls.begin() != ls.end()));
    1346             :     }
    1347             : 
    1348           1 :     TestIterator<OGRLineString>();
    1349           1 :     TestIterator<OGRLineString>("LINESTRING(0 0)", 1);
    1350           1 :     TestIterator<OGRLineString, OGRCurve>("LINESTRING(0 0)", 1);
    1351           1 :     TestIterator<OGRLineString, OGRCurve>();
    1352           1 :     TestIterator<OGRLinearRing>();
    1353           1 :     TestIterator<OGRCircularString>();
    1354           1 :     TestIterator<OGRCircularString>("CIRCULARSTRING(0 0,1 1,0 0)", 3);
    1355           1 :     TestIterator<OGRCircularString, OGRCurve>("CIRCULARSTRING(0 0,1 1,0 0)", 3);
    1356           1 :     TestIterator<OGRCompoundCurve>();
    1357           1 :     TestIterator<OGRCompoundCurve>("COMPOUNDCURVE((0 0,1 1))", 1);
    1358           1 :     TestIterator<OGRCompoundCurve, OGRCurve>(
    1359             :         "COMPOUNDCURVE((0 0,1 1),CIRCULARSTRING(1 1,2 2,3 3))", 4);
    1360           1 :     TestIterator<OGRCompoundCurve>("COMPOUNDCURVE(CIRCULARSTRING EMPTY)", 1);
    1361           1 :     TestIterator<OGRCurvePolygon>();
    1362           1 :     TestIterator<OGRCurvePolygon>("CURVEPOLYGON((0 0,1 1,1 0,0 0))", 1);
    1363           1 :     TestIterator<OGRPolygon>();
    1364           1 :     TestIterator<OGRPolygon>("POLYGON((0 0,1 1,1 0,0 0))", 1);
    1365           1 :     TestIterator<OGRGeometryCollection>();
    1366           1 :     TestIterator<OGRGeometryCollection>("GEOMETRYCOLLECTION(POINT(0 0))", 1);
    1367           1 :     TestIterator<OGRMultiSurface>();
    1368           1 :     TestIterator<OGRMultiSurface>("MULTISURFACE(((0 0)))", 1);
    1369           1 :     TestIterator<OGRMultiPolygon>();
    1370           1 :     TestIterator<OGRMultiPolygon>("MULTIPOLYGON(((0 0)))", 1);
    1371           1 :     TestIterator<OGRMultiPoint>();
    1372           1 :     TestIterator<OGRMultiPoint>("MULTIPOINT(0 0)", 1);
    1373           1 :     TestIterator<OGRMultiCurve>();
    1374           1 :     TestIterator<OGRMultiCurve>("MULTICURVE((0 0))", 1);
    1375           1 :     TestIterator<OGRMultiLineString>();
    1376           1 :     TestIterator<OGRMultiLineString>("MULTILINESTRING((0 0))", 1);
    1377           1 :     TestIterator<OGRTriangle>();
    1378           1 :     TestIterator<OGRTriangle>("TRIANGLE((0 0,0 1,1 1,0 0))", 1);
    1379           1 :     TestIterator<OGRPolyhedralSurface>();
    1380           1 :     TestIterator<OGRPolyhedralSurface>("POLYHEDRALSURFACE(((0 0,0 1,1 1,0 0)))",
    1381             :                                        1);
    1382           1 :     TestIterator<OGRTriangulatedSurface>();
    1383           1 :     TestIterator<OGRTriangulatedSurface>("TIN(((0 0,0 1,1 1,0 0)))", 1);
    1384             : 
    1385             :     // Test that the update of the iterated point of a linestring is
    1386             :     // immediately taken into account
    1387             :     // (https://github.com/OSGeo/gdal/issues/6215)
    1388             :     {
    1389           1 :         OGRLineString oLS;
    1390           1 :         oLS.addPoint(1, 2);
    1391           1 :         oLS.addPoint(3, 4);
    1392           1 :         int i = 0;
    1393           3 :         for (auto &&p : oLS)
    1394             :         {
    1395           2 :             p.setX(i * 10);
    1396           2 :             p.setY(i * 10 + 1);
    1397           2 :             p.setZ(i * 10 + 2);
    1398           2 :             p.setM(i * 10 + 3);
    1399           2 :             ASSERT_EQ(oLS.getX(i), p.getX());
    1400           2 :             ASSERT_EQ(oLS.getY(i), p.getY());
    1401           2 :             ASSERT_EQ(oLS.getZ(i), p.getZ());
    1402           2 :             ASSERT_EQ(oLS.getM(i), p.getM());
    1403           2 :             ++i;
    1404             :         }
    1405             :     }
    1406             : 
    1407             :     {
    1408             :         class PointCounterVisitorAndUpdate : public OGRDefaultGeometryVisitor
    1409             :         {
    1410             :           public:
    1411             :             PointCounterVisitorAndUpdate() = default;
    1412             : 
    1413             :             using OGRDefaultGeometryVisitor::visit;
    1414             : 
    1415           2 :             void visit(OGRPoint *poPoint) override
    1416             :             {
    1417           2 :                 poPoint->setZ(100);
    1418           2 :                 poPoint->setM(1000);
    1419           2 :             }
    1420             :         };
    1421             : 
    1422           1 :         OGRLineString oLS;
    1423           1 :         oLS.addPoint(1, 2);
    1424           1 :         oLS.addPoint(3, 4);
    1425           0 :         PointCounterVisitorAndUpdate oVisitor;
    1426           1 :         oLS.accept(&oVisitor);
    1427             : 
    1428           1 :         ASSERT_EQ(oLS.getZ(0), 100.0);
    1429           1 :         ASSERT_EQ(oLS.getZ(1), 100.0);
    1430           1 :         ASSERT_EQ(oLS.getM(0), 1000.0);
    1431           1 :         ASSERT_EQ(oLS.getM(1), 1000.0);
    1432             :     }
    1433             : }
    1434             : 
    1435             : // Test OGRToOGCGeomType()
    1436           4 : TEST_F(test_ogr, OGRToOGCGeomType)
    1437             : {
    1438           1 :     EXPECT_STREQ(OGRToOGCGeomType(wkbPoint), "POINT");
    1439           1 :     EXPECT_STREQ(OGRToOGCGeomType(wkbPointM), "POINT");
    1440           1 :     EXPECT_STREQ(OGRToOGCGeomType(wkbPoint, /*bCamelCase=*/true), "Point");
    1441           1 :     EXPECT_STREQ(
    1442             :         OGRToOGCGeomType(wkbPoint, /*bCamelCase=*/true, /*bAddZM=*/true),
    1443             :         "Point");
    1444           1 :     EXPECT_STREQ(
    1445             :         OGRToOGCGeomType(wkbPoint25D, /*bCamelCase=*/true, /*bAddZM=*/true),
    1446             :         "PointZ");
    1447           1 :     EXPECT_STREQ(
    1448             :         OGRToOGCGeomType(wkbPointM, /*bCamelCase=*/true, /*bAddZM=*/true),
    1449             :         "PointM");
    1450           1 :     EXPECT_STREQ(
    1451             :         OGRToOGCGeomType(wkbPointZM, /*bCamelCase=*/true, /*bAddZM=*/true),
    1452             :         "PointZM");
    1453           1 :     EXPECT_STREQ(OGRToOGCGeomType(wkbPointZM, /*bCamelCase=*/true,
    1454             :                                   /*bAddZM=*/true, /*bAddSpaceBeforeZM=*/true),
    1455             :                  "Point ZM");
    1456           1 : }
    1457             : 
    1458             : // Test layer, dataset-feature and layer-feature iterators
    1459           4 : TEST_F(test_ogr, DatasetFeature_and_LayerFeature_iterators)
    1460             : {
    1461           1 :     if (!GDALGetDriverByName("ESRI Shapefile"))
    1462             :     {
    1463           0 :         GTEST_SKIP() << "ESRI Shapefile driver missing";
    1464             :     }
    1465             :     else
    1466             :     {
    1467           1 :         std::string file(data_ + SEP + "poly.shp");
    1468             :         GDALDatasetUniquePtr poDS(
    1469           1 :             GDALDataset::Open(file.c_str(), GDAL_OF_VECTOR));
    1470           1 :         ASSERT_TRUE(poDS != nullptr);
    1471             : 
    1472             :         {
    1473           1 :             GIntBig nExpectedFID = 0;
    1474          11 :             for (const auto &oFeatureLayerPair : poDS->GetFeatures())
    1475             :             {
    1476          10 :                 ASSERT_EQ(oFeatureLayerPair.feature->GetFID(), nExpectedFID);
    1477          10 :                 nExpectedFID++;
    1478          10 :                 ASSERT_EQ(oFeatureLayerPair.layer, poDS->GetLayer(0));
    1479             :             }
    1480           1 :             ASSERT_EQ(nExpectedFID, 10);
    1481             :         }
    1482             : 
    1483           1 :         ASSERT_EQ(poDS->GetLayers().size(), 1U);
    1484           1 :         ASSERT_EQ(poDS->GetLayers()[0], poDS->GetLayer(0));
    1485           1 :         ASSERT_EQ(poDS->GetLayers()[static_cast<size_t>(0)], poDS->GetLayer(0));
    1486           1 :         ASSERT_EQ(poDS->GetLayers()["poly"], poDS->GetLayer(0));
    1487             : 
    1488           2 :         for (auto poLayer : poDS->GetLayers())
    1489             :         {
    1490           1 :             GIntBig nExpectedFID = 0;
    1491          11 :             for (const auto &poFeature : poLayer)
    1492             :             {
    1493          10 :                 ASSERT_EQ(poFeature->GetFID(), nExpectedFID);
    1494          10 :                 nExpectedFID++;
    1495             :             }
    1496           1 :             ASSERT_EQ(nExpectedFID, 10);
    1497             : 
    1498           1 :             nExpectedFID = 0;
    1499          11 :             for (const auto &oFeatureLayerPair : poDS->GetFeatures())
    1500             :             {
    1501          10 :                 ASSERT_EQ(oFeatureLayerPair.feature->GetFID(), nExpectedFID);
    1502          10 :                 nExpectedFID++;
    1503          10 :                 ASSERT_EQ(oFeatureLayerPair.layer, poLayer);
    1504             :             }
    1505           1 :             ASSERT_EQ(nExpectedFID, 10);
    1506             : 
    1507           1 :             nExpectedFID = 0;
    1508           5 :             OGR_FOR_EACH_FEATURE_BEGIN(hFeat,
    1509             :                                        reinterpret_cast<OGRLayerH>(poLayer))
    1510             :             {
    1511           5 :                 if (nExpectedFID == 0)
    1512             :                 {
    1513           1 :                     nExpectedFID = 1;
    1514           1 :                     continue;
    1515             :                 }
    1516           4 :                 ASSERT_EQ(OGR_F_GetFID(hFeat), nExpectedFID);
    1517           4 :                 nExpectedFID++;
    1518           4 :                 if (nExpectedFID == 5)
    1519           1 :                     break;
    1520             :             }
    1521           5 :             OGR_FOR_EACH_FEATURE_END(hFeat)
    1522           1 :             ASSERT_EQ(nExpectedFID, 5);
    1523             : 
    1524           1 :             auto oIter = poLayer->begin();
    1525           1 :             CPLPushErrorHandler(CPLQuietErrorHandler);
    1526             :             // Only one feature iterator can be active at a time
    1527           1 :             auto oIter2 = poLayer->begin();
    1528           1 :             CPLPopErrorHandler();
    1529           1 :             ASSERT_TRUE(!(oIter2 != poLayer->end()));
    1530           1 :             ASSERT_TRUE(oIter != poLayer->end());
    1531             :         }
    1532             : 
    1533           1 :         poDS.reset(GetGDALDriverManager()->GetDriverByName("MEM")->Create(
    1534             :             "", 0, 0, 0, GDT_Unknown, nullptr));
    1535           1 :         int nCountLayers = 0;
    1536           1 :         for (auto poLayer : poDS->GetLayers())
    1537             :         {
    1538           0 :             CPL_IGNORE_RET_VAL(poLayer);
    1539           0 :             nCountLayers++;
    1540             :         }
    1541           1 :         ASSERT_EQ(nCountLayers, 0);
    1542             : 
    1543           1 :         poDS->CreateLayer("foo");
    1544           1 :         poDS->CreateLayer("bar", nullptr);
    1545           3 :         for (auto poLayer : poDS->GetLayers())
    1546             :         {
    1547           2 :             if (nCountLayers == 0)
    1548             :             {
    1549           1 :                 EXPECT_STREQ(poLayer->GetName(), "foo")
    1550           0 :                     << "layer " << poLayer->GetName();
    1551             :             }
    1552           1 :             else if (nCountLayers == 1)
    1553             :             {
    1554           1 :                 EXPECT_STREQ(poLayer->GetName(), "bar")
    1555           0 :                     << "layer " << poLayer->GetName();
    1556             :             }
    1557           2 :             nCountLayers++;
    1558             :         }
    1559           1 :         ASSERT_EQ(nCountLayers, 2);
    1560             : 
    1561           1 :         auto layers = poDS->GetLayers();
    1562             :         {
    1563             :             // std::copy requires a InputIterator
    1564           1 :             std::vector<OGRLayer *> oTarget;
    1565           1 :             oTarget.resize(2);
    1566           1 :             std::copy(layers.begin(), layers.end(), oTarget.begin());
    1567           1 :             ASSERT_EQ(oTarget[0], layers[0]);
    1568           1 :             ASSERT_EQ(oTarget[1], layers[1]);
    1569             : 
    1570             :             // but in practice not necessarily uses the postincrement iterator.
    1571           1 :             oTarget.clear();
    1572           1 :             oTarget.resize(2);
    1573           1 :             auto input_iterator = layers.begin();
    1574           1 :             auto output_iterator = oTarget.begin();
    1575           3 :             while (input_iterator != layers.end())
    1576             :             {
    1577           2 :                 *output_iterator++ = *input_iterator++;
    1578             :             }
    1579           1 :             ASSERT_EQ(oTarget[0], layers[0]);
    1580           1 :             ASSERT_EQ(oTarget[1], layers[1]);
    1581             :         }
    1582             : 
    1583             :         // Test copy constructor
    1584             :         {
    1585           1 :             GDALDataset::Layers::Iterator srcIter(poDS->GetLayers().begin());
    1586           1 :             ++srcIter;
    1587           1 :             GDALDataset::Layers::Iterator newIter(srcIter);
    1588           1 :             srcIter = layers.begin();  // avoid Coverity Scan warning
    1589           1 :             ASSERT_EQ(*newIter, layers[1]);
    1590             :         }
    1591             : 
    1592             :         // Test assignment operator
    1593             :         {
    1594           1 :             GDALDataset::Layers::Iterator srcIter(poDS->GetLayers().begin());
    1595           1 :             ++srcIter;
    1596           1 :             GDALDataset::Layers::Iterator newIter;
    1597           1 :             newIter = srcIter;
    1598           1 :             srcIter = layers.begin();  // avoid Coverity Scan warning
    1599           1 :             ASSERT_EQ(*newIter, layers[1]);
    1600             :         }
    1601             : 
    1602             :         // Test move constructor
    1603             :         {
    1604           1 :             GDALDataset::Layers::Iterator srcIter(poDS->GetLayers().begin());
    1605           1 :             ++srcIter;
    1606           1 :             GDALDataset::Layers::Iterator newIter(std::move(srcIter));
    1607           1 :             ASSERT_EQ(*newIter, layers[1]);
    1608             :         }
    1609             : 
    1610             :         // Test move assignment operator
    1611             :         {
    1612           1 :             GDALDataset::Layers::Iterator srcIter(poDS->GetLayers().begin());
    1613           1 :             ++srcIter;
    1614           1 :             GDALDataset::Layers::Iterator newIter;
    1615           1 :             newIter = std::move(srcIter);
    1616           1 :             ASSERT_EQ(*newIter, layers[1]);
    1617             :         }
    1618             : 
    1619           1 :         const GDALDataset *poConstDS = poDS.get();
    1620           1 :         ASSERT_EQ(poConstDS->GetLayers().size(), 2U);
    1621           1 :         ASSERT_EQ(poConstDS->GetLayers()[0], poConstDS->GetLayer(0));
    1622           1 :         ASSERT_EQ(poConstDS->GetLayers()[static_cast<size_t>(0)],
    1623             :                   poConstDS->GetLayer(0));
    1624           1 :         ASSERT_EQ(poConstDS->GetLayers()["foo"], poConstDS->GetLayer(0));
    1625           1 :         nCountLayers = 0;
    1626           3 :         for (auto &&poLayer : poDS->GetLayers())
    1627             :         {
    1628           2 :             if (nCountLayers == 0)
    1629             :             {
    1630           1 :                 EXPECT_STREQ(poLayer->GetName(), "foo")
    1631           0 :                     << "layer " << poLayer->GetName();
    1632             :             }
    1633           1 :             else if (nCountLayers == 1)
    1634             :             {
    1635           1 :                 EXPECT_STREQ(poLayer->GetName(), "bar")
    1636           0 :                     << "layer " << poLayer->GetName();
    1637             :             }
    1638           2 :             nCountLayers++;
    1639             :         }
    1640           1 :         ASSERT_EQ(nCountLayers, 2);
    1641             : 
    1642           1 :         auto constLayers = poConstDS->GetLayers();
    1643             :         {
    1644             :             // std::copy requires a InputIterator
    1645           1 :             std::vector<const OGRLayer *> oTarget;
    1646           1 :             oTarget.resize(2);
    1647           1 :             std::copy(constLayers.begin(), constLayers.end(), oTarget.begin());
    1648           1 :             ASSERT_EQ(oTarget[0], constLayers[0]);
    1649           1 :             ASSERT_EQ(oTarget[1], constLayers[1]);
    1650             : 
    1651             :             // but in practice not necessarily uses the postincrement iterator.
    1652           1 :             oTarget.clear();
    1653           1 :             oTarget.resize(2);
    1654           1 :             auto input_iterator = constLayers.begin();
    1655           1 :             auto output_iterator = oTarget.begin();
    1656           3 :             while (input_iterator != constLayers.end())
    1657             :             {
    1658           2 :                 *output_iterator++ = *input_iterator++;
    1659             :             }
    1660           1 :             ASSERT_EQ(oTarget[0], constLayers[0]);
    1661           1 :             ASSERT_EQ(oTarget[1], constLayers[1]);
    1662             :         }
    1663             : 
    1664             :         // Test copy constructor
    1665             :         {
    1666           1 :             auto srcIter(poConstDS->GetLayers().begin());
    1667           1 :             ++srcIter;
    1668           1 :             auto newIter(srcIter);
    1669           1 :             srcIter = constLayers.begin();  // avoid Coverity Scan warning
    1670           1 :             ASSERT_EQ(*newIter, constLayers[1]);
    1671             :         }
    1672             : 
    1673             :         // Test assignment operator
    1674             :         {
    1675           1 :             auto srcIter(poConstDS->GetLayers().begin());
    1676           1 :             ++srcIter;
    1677           1 :             GDALDataset::ConstLayers::Iterator newIter;
    1678           1 :             newIter = srcIter;
    1679           1 :             srcIter = constLayers.begin();  // avoid Coverity Scan warning
    1680           1 :             ASSERT_EQ(*newIter, constLayers[1]);
    1681             :         }
    1682             : 
    1683             :         // Test move constructor
    1684             :         {
    1685           1 :             auto srcIter(poConstDS->GetLayers().begin());
    1686           1 :             ++srcIter;
    1687           1 :             auto newIter(std::move(srcIter));
    1688           1 :             ASSERT_EQ(*newIter, constLayers[1]);
    1689             :         }
    1690             : 
    1691             :         // Test move assignment operator
    1692             :         {
    1693           1 :             auto srcIter(poConstDS->GetLayers().begin());
    1694           1 :             ++srcIter;
    1695           1 :             GDALDataset::ConstLayers::Iterator newIter;
    1696           1 :             newIter = std::move(srcIter);
    1697           1 :             ASSERT_EQ(*newIter, constLayers[1]);
    1698             :         }
    1699             :     }
    1700             : }
    1701             : 
    1702             : // Test field iterator
    1703           4 : TEST_F(test_ogr, field_iterator)
    1704             : {
    1705           1 :     OGRFeatureDefn *poFeatureDefn = new OGRFeatureDefn();
    1706           1 :     poFeatureDefn->Reference();
    1707             :     {
    1708           2 :         OGRFieldDefn oFieldDefn("str_field", OFTString);
    1709           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1710             :     }
    1711             :     {
    1712           2 :         OGRFieldDefn oFieldDefn("int_field", OFTInteger);
    1713           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1714             :     }
    1715             :     {
    1716           2 :         OGRFieldDefn oFieldDefn("int64_field", OFTInteger64);
    1717           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1718             :     }
    1719             :     {
    1720           2 :         OGRFieldDefn oFieldDefn("double_field", OFTReal);
    1721           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1722             :     }
    1723             :     {
    1724           2 :         OGRFieldDefn oFieldDefn("null_field", OFTReal);
    1725           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1726             :     }
    1727             :     {
    1728           2 :         OGRFieldDefn oFieldDefn("unset_field", OFTReal);
    1729           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1730             :     }
    1731             :     {
    1732           2 :         OGRFieldDefn oFieldDefn("dt_field", OFTDateTime);
    1733           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1734             :     }
    1735             :     {
    1736           2 :         OGRFieldDefn oFieldDefn("strlist_field", OFTStringList);
    1737           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1738             :     }
    1739             :     {
    1740           2 :         OGRFieldDefn oFieldDefn("intlist_field", OFTIntegerList);
    1741           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1742             :     }
    1743             :     {
    1744           2 :         OGRFieldDefn oFieldDefn("int64list_field", OFTInteger64List);
    1745           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1746             :     }
    1747             :     {
    1748           2 :         OGRFieldDefn oFieldDefn("doublelist_field", OFTRealList);
    1749           1 :         poFeatureDefn->AddFieldDefn(&oFieldDefn);
    1750             :     }
    1751           1 :     OGRFeature oFeature(poFeatureDefn);
    1752             : 
    1753             :     {
    1754           1 :         OGRFeature oFeatureTmp(poFeatureDefn);
    1755           1 :         oFeatureTmp[0] = "bar";
    1756           2 :         ASSERT_STREQ(oFeatureTmp[0].GetString(), "bar");
    1757             :         {
    1758             :             // Proxy reference
    1759           1 :             auto &&x = oFeatureTmp[0];
    1760           1 :             auto &xRef(x);
    1761           1 :             x = xRef;
    1762           2 :             ASSERT_STREQ(oFeatureTmp[0].GetString(), "bar");
    1763             :         }
    1764             :         {
    1765           1 :             oFeatureTmp[0] = oFeatureTmp[0];
    1766           2 :             ASSERT_STREQ(oFeatureTmp[0].GetString(), "bar");
    1767             :         }
    1768             :         {
    1769             :             // Proxy reference
    1770           1 :             auto &&x = oFeatureTmp[0];
    1771           1 :             x = "baz";
    1772           1 :             ASSERT_STREQ(x.GetString(), "baz");
    1773             :         }
    1774           1 :         oFeatureTmp["str_field"] = std::string("foo");
    1775           1 :         oFeatureTmp["int_field"] = 123;
    1776           1 :         oFeatureTmp["int64_field"] = oFeatureTmp["int_field"];
    1777           1 :         ASSERT_EQ(oFeatureTmp["int64_field"].GetInteger64(), 123);
    1778           1 :         oFeatureTmp["int64_field"] = static_cast<GIntBig>(1234567890123);
    1779           1 :         oFeatureTmp["double_field"] = 123.45;
    1780           1 :         oFeatureTmp["null_field"].SetNull();
    1781           1 :         oFeatureTmp["unset_field"].clear();
    1782           1 :         oFeatureTmp["unset_field"].Unset();
    1783           1 :         oFeatureTmp["dt_field"].SetDateTime(2018, 4, 5, 12, 34, 56.75f, 0);
    1784           1 :         oFeatureTmp["strlist_field"] = CPLStringList().List();
    1785           1 :         oFeatureTmp["strlist_field"] = std::vector<std::string>();
    1786           3 :         oFeatureTmp["strlist_field"] = std::vector<std::string>{"foo", "bar"};
    1787             :         oFeatureTmp["strlist_field"] =
    1788           1 :             static_cast<CSLConstList>(oFeatureTmp["strlist_field"]);
    1789           1 :         ASSERT_EQ(
    1790             :             CSLCount(static_cast<CSLConstList>(oFeatureTmp["strlist_field"])),
    1791             :             2);
    1792           1 :         oFeatureTmp["intlist_field"] = std::vector<int>();
    1793           1 :         oFeatureTmp["intlist_field"] = std::vector<int>{12, 34};
    1794           1 :         oFeatureTmp["int64list_field"] = std::vector<GIntBig>();
    1795             :         oFeatureTmp["int64list_field"] =
    1796           1 :             std::vector<GIntBig>{1234567890123, 34};
    1797           1 :         oFeatureTmp["doublelist_field"] = std::vector<double>();
    1798           1 :         oFeatureTmp["doublelist_field"] = std::vector<double>{12.25, 56.75};
    1799             : 
    1800          12 :         for (const auto &oField : oFeatureTmp)
    1801             :         {
    1802          11 :             oFeature[oField.GetIndex()] = oField;
    1803             :         }
    1804             :     }
    1805             : 
    1806             :     {
    1807           1 :         int x = oFeature[1];
    1808           1 :         ASSERT_EQ(x, 123);
    1809             :     }
    1810             :     {
    1811           1 :         int x = oFeature["int_field"];
    1812           1 :         ASSERT_EQ(x, 123);
    1813             :     }
    1814             :     {
    1815           1 :         GIntBig x = oFeature["int64_field"];
    1816           1 :         ASSERT_EQ(x, static_cast<GIntBig>(1234567890123));
    1817             :     }
    1818             :     {
    1819           1 :         double x = oFeature["double_field"];
    1820           1 :         ASSERT_EQ(x, 123.45);
    1821             :     }
    1822             :     {
    1823           1 :         const char *x = oFeature["str_field"];
    1824           1 :         ASSERT_STREQ(x, "foo");
    1825             :     }
    1826           1 :     bool bExceptionHit = false;
    1827             :     try
    1828             :     {
    1829           1 :         oFeature["inexisting_field"];
    1830             :     }
    1831           1 :     catch (const OGRFeature::FieldNotFoundException &)
    1832             :     {
    1833           1 :         bExceptionHit = true;
    1834             :     }
    1835           1 :     ASSERT_TRUE(bExceptionHit);
    1836             : 
    1837           1 :     int iIter = 0;
    1838           1 :     const OGRFeature *poConstFeature = &oFeature;
    1839          12 :     for (const auto &oField : *poConstFeature)
    1840             :     {
    1841          11 :         ASSERT_EQ(oField.GetIndex(), iIter);
    1842          11 :         ASSERT_EQ(oField.GetDefn(), poFeatureDefn->GetFieldDefn(iIter));
    1843          22 :         ASSERT_EQ(CPLString(oField.GetName()),
    1844             :                   CPLString(oField.GetDefn()->GetNameRef()));
    1845          11 :         ASSERT_EQ(oField.GetType(), oField.GetDefn()->GetType());
    1846          11 :         ASSERT_EQ(oField.GetSubType(), oField.GetDefn()->GetSubType());
    1847          11 :         if (iIter == 0)
    1848             :         {
    1849           1 :             ASSERT_EQ(oField.IsUnset(), false);
    1850           1 :             ASSERT_EQ(oField.IsNull(), false);
    1851           2 :             ASSERT_EQ(CPLString(oField.GetRawValue()->String),
    1852             :                       CPLString("foo"));
    1853           2 :             ASSERT_EQ(CPLString(oField.GetString()), CPLString("foo"));
    1854           2 :             ASSERT_EQ(CPLString(oField.GetAsString()), CPLString("foo"));
    1855             :         }
    1856          10 :         else if (iIter == 1)
    1857             :         {
    1858           1 :             ASSERT_EQ(oField.GetRawValue()->Integer, 123);
    1859           1 :             ASSERT_EQ(oField.GetInteger(), 123);
    1860           1 :             ASSERT_EQ(oField.GetAsInteger(), 123);
    1861           1 :             ASSERT_EQ(oField.GetAsInteger64(), 123);
    1862           1 :             ASSERT_EQ(oField.GetAsDouble(), 123.0);
    1863           2 :             ASSERT_EQ(CPLString(oField.GetAsString()), CPLString("123"));
    1864             :         }
    1865           9 :         else if (iIter == 2)
    1866             :         {
    1867           1 :             ASSERT_EQ(oField.GetRawValue()->Integer64, 1234567890123);
    1868           1 :             ASSERT_EQ(oField.GetInteger64(), 1234567890123);
    1869           1 :             ASSERT_EQ(oField.GetAsInteger(), 2147483647);
    1870           1 :             ASSERT_EQ(oField.GetAsInteger64(), 1234567890123);
    1871           1 :             ASSERT_EQ(oField.GetAsDouble(), 1234567890123.0);
    1872           2 :             ASSERT_EQ(CPLString(oField.GetAsString()),
    1873             :                       CPLString("1234567890123"));
    1874             :         }
    1875           8 :         else if (iIter == 3)
    1876             :         {
    1877           1 :             ASSERT_EQ(oField.GetRawValue()->Real, 123.45);
    1878           1 :             ASSERT_EQ(oField.GetDouble(), 123.45);
    1879           1 :             ASSERT_EQ(oField.GetAsInteger(), 123);
    1880           1 :             ASSERT_EQ(oField.GetAsInteger64(), 123);
    1881           1 :             ASSERT_EQ(oField.GetAsDouble(), 123.45);
    1882           2 :             ASSERT_EQ(CPLString(oField.GetAsString()), CPLString("123.45"));
    1883             :         }
    1884           7 :         else if (iIter == 4)
    1885             :         {
    1886           1 :             ASSERT_EQ(oField.IsUnset(), false);
    1887           1 :             ASSERT_EQ(oField.IsNull(), true);
    1888             :         }
    1889           6 :         else if (iIter == 5)
    1890             :         {
    1891           1 :             ASSERT_EQ(oField.IsUnset(), true);
    1892           1 :             ASSERT_EQ(oField.empty(), true);
    1893           1 :             ASSERT_EQ(oField.IsNull(), false);
    1894             :         }
    1895           5 :         else if (iIter == 6)
    1896             :         {
    1897             :             int nYear, nMonth, nDay, nHour, nMin, nTZFlag;
    1898             :             float fSec;
    1899           1 :             ASSERT_EQ(oField.GetDateTime(&nYear, &nMonth, &nDay, &nHour, &nMin,
    1900             :                                          &fSec, &nTZFlag),
    1901             :                       true);
    1902           1 :             ASSERT_EQ(nYear, 2018);
    1903           1 :             ASSERT_EQ(nMonth, 4);
    1904           1 :             ASSERT_EQ(nDay, 5);
    1905           1 :             ASSERT_EQ(nHour, 12);
    1906           1 :             ASSERT_EQ(nMin, 34);
    1907           1 :             ASSERT_EQ(fSec, 56.75f);
    1908           1 :             ASSERT_EQ(nTZFlag, 0);
    1909             :         }
    1910           4 :         else if (iIter == 7)
    1911             :         {
    1912             :             std::vector<std::string> oExpected{std::string("foo"),
    1913           5 :                                                std::string("bar")};
    1914           1 :             decltype(oExpected) oGot = oField;
    1915           1 :             ASSERT_EQ(oGot.size(), oExpected.size());
    1916           3 :             for (size_t i = 0; i < oExpected.size(); i++)
    1917           2 :                 ASSERT_EQ(oGot[i], oExpected[i]);
    1918             :         }
    1919           3 :         else if (iIter == 8)
    1920             :         {
    1921           1 :             std::vector<int> oExpected{12, 34};
    1922           1 :             decltype(oExpected) oGot = oField;
    1923           1 :             ASSERT_EQ(oGot.size(), oExpected.size());
    1924           3 :             for (size_t i = 0; i < oExpected.size(); i++)
    1925           2 :                 ASSERT_EQ(oGot[i], oExpected[i]);
    1926             :         }
    1927           2 :         else if (iIter == 9)
    1928             :         {
    1929           1 :             std::vector<GIntBig> oExpected{1234567890123, 34};
    1930           1 :             decltype(oExpected) oGot = oField;
    1931           1 :             ASSERT_EQ(oGot.size(), oExpected.size());
    1932           3 :             for (size_t i = 0; i < oExpected.size(); i++)
    1933           2 :                 ASSERT_EQ(oGot[i], oExpected[i]);
    1934             :         }
    1935           1 :         else if (iIter == 10)
    1936             :         {
    1937           1 :             std::vector<double> oExpected{12.25, 56.75};
    1938           1 :             decltype(oExpected) oGot = oField;
    1939           1 :             ASSERT_EQ(oGot.size(), oExpected.size());
    1940           3 :             for (size_t i = 0; i < oExpected.size(); i++)
    1941           2 :                 ASSERT_EQ(oGot[i], oExpected[i]);
    1942             :         }
    1943          11 :         iIter++;
    1944             :     }
    1945           1 :     poFeatureDefn->Release();
    1946             : }
    1947             : 
    1948             : // Test OGRLinearRing::isPointOnRingBoundary()
    1949           4 : TEST_F(test_ogr, isPointOnRingBoundary)
    1950             : {
    1951           1 :     OGRPolygon oPoly;
    1952           1 :     const char *pszPoly = "POLYGON((10 9,11 10,10 11,9 10,10 9))";
    1953           1 :     oPoly.importFromWkt(&pszPoly);
    1954           1 :     auto poRing = oPoly.getExteriorRing();
    1955             : 
    1956             :     // On first vertex
    1957             :     {
    1958           1 :         OGRPoint p(10, 9);
    1959           1 :         ASSERT_TRUE(poRing->isPointOnRingBoundary(&p, false));
    1960             :     }
    1961             : 
    1962             :     // On second vertex
    1963             :     {
    1964           1 :         OGRPoint p(11, 10);
    1965           1 :         ASSERT_TRUE(poRing->isPointOnRingBoundary(&p, false));
    1966             :     }
    1967             : 
    1968             :     // Middle of first segment
    1969             :     {
    1970           1 :         OGRPoint p(10.5, 9.5);
    1971           1 :         ASSERT_TRUE(poRing->isPointOnRingBoundary(&p, false));
    1972             :     }
    1973             : 
    1974             :     // "Before" first segment
    1975             :     {
    1976           1 :         OGRPoint p(10 - 1, 9 - 1);
    1977           1 :         ASSERT_TRUE(!poRing->isPointOnRingBoundary(&p, false));
    1978             :     }
    1979             : 
    1980             :     // "After" first segment
    1981             :     {
    1982           1 :         OGRPoint p(11 + 1, 10 + 1);
    1983           1 :         ASSERT_TRUE(!poRing->isPointOnRingBoundary(&p, false));
    1984             :     }
    1985             : 
    1986             :     // On third vertex
    1987             :     {
    1988           1 :         OGRPoint p(10, 11);
    1989           1 :         ASSERT_TRUE(poRing->isPointOnRingBoundary(&p, false));
    1990             :     }
    1991             : 
    1992             :     // Middle of second segment
    1993             :     {
    1994           1 :         OGRPoint p(10.5, 10.5);
    1995           1 :         ASSERT_TRUE(poRing->isPointOnRingBoundary(&p, false));
    1996             :     }
    1997             : 
    1998             :     // On fourth vertex
    1999             :     {
    2000           1 :         OGRPoint p(9, 10);
    2001           1 :         ASSERT_TRUE(poRing->isPointOnRingBoundary(&p, false));
    2002             :     }
    2003             : 
    2004             :     // Middle of third segment
    2005             :     {
    2006           1 :         OGRPoint p(9.5, 10.5);
    2007           1 :         ASSERT_TRUE(poRing->isPointOnRingBoundary(&p, false));
    2008             :     }
    2009             : 
    2010             :     // Middle of fourth segment
    2011             :     {
    2012           1 :         OGRPoint p(9.5, 9.5);
    2013           1 :         ASSERT_TRUE(poRing->isPointOnRingBoundary(&p, false));
    2014             :     }
    2015             : }
    2016             : 
    2017             : // Test OGRGeometry::exportToWkt()
    2018           4 : TEST_F(test_ogr, OGRGeometry_exportToWkt)
    2019             : {
    2020           1 :     char *pszWKT = nullptr;
    2021           1 :     OGRPoint p(1, 2);
    2022           1 :     p.exportToWkt(&pszWKT);
    2023           1 :     ASSERT_TRUE(pszWKT != nullptr);
    2024           1 :     EXPECT_STREQ(pszWKT, "POINT (1 2)");
    2025           1 :     CPLFree(pszWKT);
    2026             : }
    2027             : 
    2028             : // Test OGRGeometry::clone()
    2029           4 : TEST_F(test_ogr, OGRGeometry_clone)
    2030             : {
    2031           1 :     const char *apszWKT[] = {
    2032             :         "POINT (0 0)",
    2033             :         "POINT ZM EMPTY",
    2034             :         "LINESTRING (0 0)",
    2035             :         "LINESTRING ZM EMPTY",
    2036             :         "POLYGON ((0 0),(0 0))",
    2037             :         "MULTIPOLYGON ZM EMPTY",
    2038             :         "MULTIPOINT ((0 0))",
    2039             :         "MULTIPOINT ZM EMPTY",
    2040             :         "MULTILINESTRING ((0 0))",
    2041             :         "MULTILINESTRING ZM EMPTY",
    2042             :         "MULTIPOLYGON (((0 0)))",
    2043             :         "MULTIPOLYGON ZM EMPTY",
    2044             :         "GEOMETRYCOLLECTION (POINT (0 0))",
    2045             :         "GEOMETRYCOLLECTION ZM EMPTY",
    2046             :         "CIRCULARSTRING (0 0,1 1,0 0)",
    2047             :         "CIRCULARSTRING Z EMPTY",
    2048             :         "CIRCULARSTRING ZM EMPTY",
    2049             :         "COMPOUNDCURVE ((0 0,1 1))",
    2050             :         "COMPOUNDCURVE ZM EMPTY",
    2051             :         "CURVEPOLYGON ((0 0,1 1,1 0,0 0))",
    2052             :         "CURVEPOLYGON ZM EMPTY",
    2053             :         "MULTICURVE ((0 0))",
    2054             :         "MULTICURVE ZM EMPTY",
    2055             :         "MULTISURFACE (((0 0)))",
    2056             :         "MULTISURFACE ZM EMPTY",
    2057             :         "TRIANGLE ((0 0,0 1,1 1,0 0))",
    2058             :         "TRIANGLE ZM EMPTY",
    2059             :         "POLYHEDRALSURFACE (((0 0,0 1,1 1,0 0)))",
    2060             :         "POLYHEDRALSURFACE ZM EMPTY",
    2061             :         "TIN (((0 0,0 1,1 1,0 0)))",
    2062             :         "TIN ZM EMPTY",
    2063             :     };
    2064           1 :     OGRSpatialReference oSRS;
    2065          32 :     for (const char *pszWKT : apszWKT)
    2066             :     {
    2067          31 :         OGRGeometry *poGeom = nullptr;
    2068          31 :         OGRGeometryFactory::createFromWkt(pszWKT, &oSRS, &poGeom);
    2069          31 :         auto poClone = poGeom->clone();
    2070          31 :         ASSERT_TRUE(poClone != nullptr);
    2071          31 :         char *outWKT = nullptr;
    2072          31 :         poClone->exportToWkt(&outWKT, wkbVariantIso);
    2073          31 :         EXPECT_STREQ(pszWKT, outWKT);
    2074          31 :         CPLFree(outWKT);
    2075          31 :         delete poClone;
    2076          31 :         delete poGeom;
    2077             :     }
    2078             : }
    2079             : 
    2080             : // Test OGRLineString::removePoint()
    2081           4 : TEST_F(test_ogr, OGRLineString_removePoint)
    2082             : {
    2083             :     {
    2084           1 :         OGRLineString ls;
    2085           1 :         ls.addPoint(0, 1);
    2086           1 :         ls.addPoint(2, 3);
    2087           1 :         ls.addPoint(4, 5);
    2088           1 :         ASSERT_TRUE(!ls.removePoint(-1));
    2089           1 :         ASSERT_TRUE(!ls.removePoint(3));
    2090           1 :         ASSERT_EQ(ls.getNumPoints(), 3);
    2091           1 :         ASSERT_TRUE(ls.removePoint(1));
    2092           1 :         ASSERT_EQ(ls.getNumPoints(), 2);
    2093           1 :         ASSERT_EQ(ls.getX(0), 0.0);
    2094           1 :         ASSERT_EQ(ls.getY(0), 1.0);
    2095           1 :         ASSERT_EQ(ls.getX(1), 4.0);
    2096           1 :         ASSERT_EQ(ls.getY(1), 5.0);
    2097           1 :         ASSERT_TRUE(ls.removePoint(1));
    2098           1 :         ASSERT_EQ(ls.getNumPoints(), 1);
    2099           1 :         ASSERT_TRUE(ls.removePoint(0));
    2100           1 :         ASSERT_EQ(ls.getNumPoints(), 0);
    2101             :     }
    2102             :     {
    2103             :         // With Z, M
    2104           1 :         OGRLineString ls;
    2105           1 :         ls.addPoint(0, 1, 20, 30);
    2106           1 :         ls.addPoint(2, 3, 40, 50);
    2107           1 :         ls.addPoint(4, 5, 60, 70);
    2108           1 :         ASSERT_TRUE(!ls.removePoint(-1));
    2109           1 :         ASSERT_TRUE(!ls.removePoint(3));
    2110           1 :         ASSERT_EQ(ls.getNumPoints(), 3);
    2111           1 :         ASSERT_TRUE(ls.removePoint(1));
    2112           1 :         ASSERT_EQ(ls.getNumPoints(), 2);
    2113           1 :         ASSERT_EQ(ls.getX(0), 0.0);
    2114           1 :         ASSERT_EQ(ls.getY(0), 1.0);
    2115           1 :         ASSERT_EQ(ls.getZ(0), 20.0);
    2116           1 :         ASSERT_EQ(ls.getM(0), 30.0);
    2117           1 :         ASSERT_EQ(ls.getX(1), 4.0);
    2118           1 :         ASSERT_EQ(ls.getY(1), 5.0);
    2119           1 :         ASSERT_EQ(ls.getZ(1), 60.0);
    2120           1 :         ASSERT_EQ(ls.getM(1), 70.0);
    2121           1 :         ASSERT_TRUE(ls.removePoint(1));
    2122           1 :         ASSERT_EQ(ls.getNumPoints(), 1);
    2123           1 :         ASSERT_TRUE(ls.removePoint(0));
    2124           1 :         ASSERT_EQ(ls.getNumPoints(), 0);
    2125             :     }
    2126             : }
    2127             : 
    2128             : // Test effect of MarkSuppressOnClose() on DXF
    2129           4 : TEST_F(test_ogr, DXF_MarkSuppressOnClose)
    2130             : {
    2131           1 :     CPLString tmpFilename(CPLGenerateTempFilename(nullptr));
    2132           1 :     tmpFilename += ".dxf";
    2133           1 :     auto poDrv = GDALDriver::FromHandle(GDALGetDriverByName("DXF"));
    2134           1 :     if (poDrv)
    2135             :     {
    2136             :         auto poDS(GDALDatasetUniquePtr(
    2137           1 :             poDrv->Create(tmpFilename, 0, 0, 0, GDT_Unknown, nullptr)));
    2138           1 :         ASSERT_TRUE(poDS != nullptr);
    2139             : 
    2140             :         OGRLayer *poLayer =
    2141           1 :             poDS->CreateLayer("test", nullptr, wkbPoint, nullptr);
    2142           1 :         ASSERT_TRUE(poLayer != nullptr);
    2143             : 
    2144         101 :         for (double x = 0; x < 100; x++)
    2145             :         {
    2146             :             OGRFeature *poFeature =
    2147         100 :                 OGRFeature::CreateFeature(poLayer->GetLayerDefn());
    2148         100 :             ASSERT_TRUE(poFeature != nullptr);
    2149         100 :             OGRPoint pt(x, 42);
    2150         100 :             ASSERT_EQ(OGRERR_NONE, poFeature->SetGeometry(&pt));
    2151         100 :             ASSERT_EQ(OGRERR_NONE, poLayer->CreateFeature(poFeature));
    2152         100 :             OGRFeature::DestroyFeature(poFeature);
    2153             :         }
    2154             : 
    2155           1 :         poDS->MarkSuppressOnClose();
    2156             : 
    2157           1 :         poDS.reset();
    2158             :         VSIStatBufL sStat;
    2159           1 :         ASSERT_TRUE(0 != VSIStatL(tmpFilename, &sStat));
    2160             :     }
    2161             : }
    2162             : 
    2163             : // Test OGREnvelope
    2164           4 : TEST_F(test_ogr, OGREnvelope)
    2165             : {
    2166           1 :     OGREnvelope s1;
    2167           1 :     ASSERT_TRUE(!s1.IsInit());
    2168             :     {
    2169           1 :         OGREnvelope s2(s1);
    2170           1 :         ASSERT_TRUE(s1 == s2);
    2171           1 :         ASSERT_TRUE(!(s1 != s2));
    2172             :     }
    2173             : 
    2174           1 :     s1.MinX = 0;
    2175           1 :     s1.MinY = 1;
    2176           1 :     s1.MaxX = 2;
    2177           1 :     s1.MaxY = 3;
    2178           1 :     ASSERT_TRUE(s1.IsInit());
    2179             :     {
    2180           1 :         OGREnvelope s2(s1);
    2181           1 :         ASSERT_TRUE(s1 == s2);
    2182           1 :         ASSERT_TRUE(!(s1 != s2));
    2183           1 :         s2.MinX += 1;
    2184           1 :         ASSERT_TRUE(s1 != s2);
    2185           1 :         ASSERT_TRUE(!(s1 == s2));
    2186             :     }
    2187             : }
    2188             : 
    2189             : // Test OGREnvelope3D
    2190           4 : TEST_F(test_ogr, OGREnvelope3D)
    2191             : {
    2192           1 :     OGREnvelope3D s1;
    2193           1 :     EXPECT_TRUE(!s1.IsInit());
    2194             :     {
    2195           1 :         OGREnvelope3D s2(s1);
    2196           1 :         EXPECT_TRUE(s1 == s2);
    2197           1 :         EXPECT_TRUE(!(s1 != s2));
    2198             :     }
    2199             : 
    2200           1 :     s1.MinX = 0;
    2201           1 :     s1.MinY = 1;
    2202           1 :     s1.MaxX = 2;
    2203           1 :     s1.MaxY = 3;
    2204           1 :     EXPECT_TRUE(s1.IsInit());
    2205           1 :     EXPECT_FALSE(s1.Is3D());
    2206           1 :     s1.MinZ = 4;
    2207           1 :     s1.MaxZ = 5;
    2208           1 :     EXPECT_TRUE(s1.Is3D());
    2209             :     {
    2210           1 :         OGREnvelope3D s2(s1);
    2211           1 :         EXPECT_TRUE(s1 == s2);
    2212           1 :         EXPECT_TRUE(!(s1 != s2));
    2213           1 :         s2.MinX += 1;
    2214           1 :         EXPECT_TRUE(s1 != s2);
    2215           1 :         EXPECT_TRUE(!(s1 == s2));
    2216             :     }
    2217           1 : }
    2218             : 
    2219             : // Test OGRStyleMgr::InitStyleString() with a style name
    2220             : // (https://github.com/OSGeo/gdal/issues/5555)
    2221           4 : TEST_F(test_ogr, InitStyleString_with_style_name)
    2222             : {
    2223           1 :     OGRStyleTableH hStyleTable = OGR_STBL_Create();
    2224           1 :     OGR_STBL_AddStyle(hStyleTable, "@my_style", "PEN(c:#FF0000,w:5px)");
    2225           1 :     OGRStyleMgrH hSM = OGR_SM_Create(hStyleTable);
    2226           1 :     EXPECT_EQ(OGR_SM_GetPartCount(hSM, nullptr), 0);
    2227           1 :     EXPECT_TRUE(OGR_SM_InitStyleString(hSM, "@my_style"));
    2228           1 :     EXPECT_EQ(OGR_SM_GetPartCount(hSM, nullptr), 1);
    2229           1 :     EXPECT_TRUE(!OGR_SM_InitStyleString(hSM, "@i_do_not_exist"));
    2230           1 :     OGR_SM_Destroy(hSM);
    2231           1 :     OGR_STBL_Destroy(hStyleTable);
    2232           1 : }
    2233             : 
    2234             : // Test OGR_L_GetArrowStream
    2235           4 : TEST_F(test_ogr, OGR_L_GetArrowStream)
    2236             : {
    2237             :     auto poDS = std::unique_ptr<GDALDataset>(
    2238             :         GetGDALDriverManager()->GetDriverByName("MEM")->Create(
    2239           1 :             "", 0, 0, 0, GDT_Unknown, nullptr));
    2240           1 :     auto poLayer = poDS->CreateLayer("test");
    2241             :     {
    2242           2 :         OGRFieldDefn oFieldDefn("str", OFTString);
    2243           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2244             :     }
    2245             :     {
    2246           2 :         OGRFieldDefn oFieldDefn("bool", OFTInteger);
    2247           1 :         oFieldDefn.SetSubType(OFSTBoolean);
    2248           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2249             :     }
    2250             :     {
    2251           2 :         OGRFieldDefn oFieldDefn("int16", OFTInteger);
    2252           1 :         oFieldDefn.SetSubType(OFSTInt16);
    2253           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2254             :     }
    2255             :     {
    2256           2 :         OGRFieldDefn oFieldDefn("int32", OFTInteger);
    2257           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2258             :     }
    2259             :     {
    2260           2 :         OGRFieldDefn oFieldDefn("int64", OFTInteger64);
    2261           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2262             :     }
    2263             :     {
    2264           2 :         OGRFieldDefn oFieldDefn("float32", OFTReal);
    2265           1 :         oFieldDefn.SetSubType(OFSTFloat32);
    2266           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2267             :     }
    2268             :     {
    2269           2 :         OGRFieldDefn oFieldDefn("float64", OFTReal);
    2270           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2271             :     }
    2272             :     {
    2273           2 :         OGRFieldDefn oFieldDefn("date", OFTDate);
    2274           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2275             :     }
    2276             :     {
    2277           2 :         OGRFieldDefn oFieldDefn("time", OFTTime);
    2278           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2279             :     }
    2280             :     {
    2281           2 :         OGRFieldDefn oFieldDefn("datetime", OFTDateTime);
    2282           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2283             :     }
    2284             :     {
    2285           2 :         OGRFieldDefn oFieldDefn("binary", OFTBinary);
    2286           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2287             :     }
    2288             :     {
    2289           2 :         OGRFieldDefn oFieldDefn("strlist", OFTStringList);
    2290           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2291             :     }
    2292             :     {
    2293           2 :         OGRFieldDefn oFieldDefn("boollist", OFTIntegerList);
    2294           1 :         oFieldDefn.SetSubType(OFSTBoolean);
    2295           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2296             :     }
    2297             :     {
    2298           2 :         OGRFieldDefn oFieldDefn("int16list", OFTIntegerList);
    2299           1 :         oFieldDefn.SetSubType(OFSTInt16);
    2300           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2301             :     }
    2302             :     {
    2303           2 :         OGRFieldDefn oFieldDefn("int32list", OFTIntegerList);
    2304           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2305             :     }
    2306             :     {
    2307           2 :         OGRFieldDefn oFieldDefn("int64list", OFTInteger64List);
    2308           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2309             :     }
    2310             :     {
    2311           2 :         OGRFieldDefn oFieldDefn("float32list", OFTRealList);
    2312           1 :         oFieldDefn.SetSubType(OFSTFloat32);
    2313           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2314             :     }
    2315             :     {
    2316           2 :         OGRFieldDefn oFieldDefn("float64list", OFTRealList);
    2317           1 :         EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    2318             :     }
    2319           1 :     auto poFDefn = poLayer->GetLayerDefn();
    2320             :     struct ArrowArrayStream stream;
    2321           1 :     ASSERT_TRUE(
    2322             :         OGR_L_GetArrowStream(OGRLayer::ToHandle(poLayer), &stream, nullptr));
    2323             :     {
    2324             :         // Cannot start a new stream while one is active
    2325             :         struct ArrowArrayStream stream2;
    2326           1 :         CPLPushErrorHandler(CPLQuietErrorHandler);
    2327           1 :         ASSERT_TRUE(OGR_L_GetArrowStream(OGRLayer::ToHandle(poLayer), &stream2,
    2328             :                                          nullptr) == false);
    2329           1 :         CPLPopErrorHandler();
    2330             :     }
    2331           1 :     ASSERT_TRUE(stream.release != nullptr);
    2332             : 
    2333             :     struct ArrowSchema schema;
    2334           1 :     CPLErrorReset();
    2335           1 :     ASSERT_TRUE(stream.get_last_error(&stream) == nullptr);
    2336           1 :     ASSERT_EQ(stream.get_schema(&stream, &schema), 0);
    2337           1 :     ASSERT_TRUE(stream.get_last_error(&stream) == nullptr);
    2338           1 :     ASSERT_TRUE(schema.release != nullptr);
    2339           1 :     ASSERT_EQ(schema.n_children,
    2340             :               1 + poFDefn->GetFieldCount() + poFDefn->GetGeomFieldCount());
    2341           1 :     schema.release(&schema);
    2342             : 
    2343             :     struct ArrowArray array;
    2344             :     // Next batch ==> End of stream
    2345           1 :     ASSERT_EQ(stream.get_next(&stream, &array), 0);
    2346           1 :     ASSERT_TRUE(array.release == nullptr);
    2347             : 
    2348             :     // Release stream
    2349           1 :     stream.release(&stream);
    2350             : 
    2351             :     {
    2352           1 :         auto poFeature = std::unique_ptr<OGRFeature>(new OGRFeature(poFDefn));
    2353           1 :         poFeature->SetField("bool", 1);
    2354           1 :         poFeature->SetField("int16", -12345);
    2355           1 :         poFeature->SetField("int32", 12345678);
    2356           1 :         poFeature->SetField("int64", static_cast<GIntBig>(12345678901234));
    2357           1 :         poFeature->SetField("float32", 1.25);
    2358           1 :         poFeature->SetField("float64", 1.250123);
    2359           1 :         poFeature->SetField("str", "abc");
    2360           1 :         poFeature->SetField("date", "2022-05-31");
    2361           1 :         poFeature->SetField("time", "12:34:56.789");
    2362           1 :         poFeature->SetField("datetime", "2022-05-31T12:34:56.789Z");
    2363           1 :         poFeature->SetField("boollist", "[False,True]");
    2364           1 :         poFeature->SetField("int16list", "[-12345,12345]");
    2365           1 :         poFeature->SetField("int32list", "[-12345678,12345678]");
    2366           1 :         poFeature->SetField("int64list", "[-12345678901234,12345678901234]");
    2367           1 :         poFeature->SetField("float32list", "[-1.25,1.25]");
    2368           1 :         poFeature->SetField("float64list", "[-1.250123,1.250123]");
    2369           1 :         poFeature->SetField("strlist", "[\"abc\",\"defghi\"]");
    2370           1 :         poFeature->SetField(poFDefn->GetFieldIndex("binary"), 2, "\xDE\xAD");
    2371           1 :         OGRGeometry *poGeom = nullptr;
    2372           1 :         OGRGeometryFactory::createFromWkt("POINT(1 2)", nullptr, &poGeom);
    2373           1 :         poFeature->SetGeometryDirectly(poGeom);
    2374           1 :         ASSERT_EQ(poLayer->CreateFeature(poFeature.get()), OGRERR_NONE);
    2375             :     }
    2376             : 
    2377             :     // Get a new stream now that we've released it
    2378           1 :     ASSERT_TRUE(
    2379             :         OGR_L_GetArrowStream(OGRLayer::ToHandle(poLayer), &stream, nullptr));
    2380           1 :     ASSERT_TRUE(stream.release != nullptr);
    2381             : 
    2382           1 :     ASSERT_EQ(stream.get_next(&stream, &array), 0);
    2383           1 :     ASSERT_TRUE(array.release != nullptr);
    2384           1 :     ASSERT_EQ(array.n_children,
    2385             :               1 + poFDefn->GetFieldCount() + poFDefn->GetGeomFieldCount());
    2386           1 :     ASSERT_EQ(array.length, poLayer->GetFeatureCount(false));
    2387           1 :     ASSERT_EQ(array.null_count, 0);
    2388           1 :     ASSERT_EQ(array.n_buffers, 1);
    2389           1 :     ASSERT_TRUE(array.buffers[0] == nullptr);  // no bitmap
    2390          21 :     for (int i = 0; i < array.n_children; i++)
    2391             :     {
    2392          20 :         ASSERT_TRUE(array.children[i]->release != nullptr);
    2393          20 :         ASSERT_EQ(array.children[i]->length, array.length);
    2394          20 :         ASSERT_TRUE(array.children[i]->n_buffers >= 2);
    2395          20 :         ASSERT_TRUE(array.children[i]->buffers[0] == nullptr);  // no bitmap
    2396          20 :         ASSERT_EQ(array.children[i]->null_count, 0);
    2397          20 :         ASSERT_TRUE(array.children[i]->buffers[1] != nullptr);
    2398          20 :         if (array.children[i]->n_buffers == 3)
    2399             :         {
    2400           3 :             ASSERT_TRUE(array.children[i]->buffers[2] != nullptr);
    2401             :         }
    2402             :     }
    2403           1 :     array.release(&array);
    2404             : 
    2405             :     // Next batch ==> End of stream
    2406           1 :     ASSERT_EQ(stream.get_next(&stream, &array), 0);
    2407           1 :     ASSERT_TRUE(array.release == nullptr);
    2408             : 
    2409             :     // Release stream
    2410           1 :     stream.release(&stream);
    2411             : 
    2412             :     // Insert 2 empty features
    2413             :     {
    2414           1 :         auto poFeature = std::unique_ptr<OGRFeature>(new OGRFeature(poFDefn));
    2415           1 :         ASSERT_EQ(poLayer->CreateFeature(poFeature.get()), OGRERR_NONE);
    2416             :     }
    2417             : 
    2418             :     {
    2419           1 :         auto poFeature = std::unique_ptr<OGRFeature>(new OGRFeature(poFDefn));
    2420           1 :         ASSERT_EQ(poLayer->CreateFeature(poFeature.get()), OGRERR_NONE);
    2421             :     }
    2422             : 
    2423             :     // Get a new stream now that we've released it
    2424             :     {
    2425             :         char **papszOptions =
    2426           1 :             CSLSetNameValue(nullptr, "MAX_FEATURES_IN_BATCH", "2");
    2427           1 :         ASSERT_TRUE(OGR_L_GetArrowStream(OGRLayer::ToHandle(poLayer), &stream,
    2428             :                                          papszOptions));
    2429           1 :         CSLDestroy(papszOptions);
    2430             :     }
    2431           1 :     ASSERT_TRUE(stream.release != nullptr);
    2432             : 
    2433           1 :     ASSERT_EQ(stream.get_next(&stream, &array), 0);
    2434           1 :     ASSERT_TRUE(array.release != nullptr);
    2435           1 :     ASSERT_EQ(array.n_children,
    2436             :               1 + poFDefn->GetFieldCount() + poFDefn->GetGeomFieldCount());
    2437           1 :     ASSERT_EQ(array.length, 2);
    2438          21 :     for (int i = 0; i < array.n_children; i++)
    2439             :     {
    2440          20 :         ASSERT_TRUE(array.children[i]->release != nullptr);
    2441          20 :         ASSERT_EQ(array.children[i]->length, array.length);
    2442          20 :         ASSERT_TRUE(array.children[i]->n_buffers >= 2);
    2443          20 :         if (i > 0)
    2444             :         {
    2445          19 :             ASSERT_TRUE(array.children[i]->buffers[0] !=
    2446             :                         nullptr);  // we have a bitmap
    2447          19 :             ASSERT_EQ(array.children[i]->null_count, 1);
    2448             :         }
    2449          20 :         ASSERT_TRUE(array.children[i]->buffers[1] != nullptr);
    2450          20 :         if (array.children[i]->n_buffers == 3)
    2451             :         {
    2452           3 :             ASSERT_TRUE(array.children[i]->buffers[2] != nullptr);
    2453             :         }
    2454             :     }
    2455           1 :     array.release(&array);
    2456             : 
    2457             :     // Next batch
    2458           1 :     ASSERT_EQ(stream.get_next(&stream, &array), 0);
    2459           1 :     ASSERT_TRUE(array.release != nullptr);
    2460           1 :     ASSERT_EQ(array.n_children,
    2461             :               1 + poFDefn->GetFieldCount() + poFDefn->GetGeomFieldCount());
    2462           1 :     ASSERT_EQ(array.length, 1);
    2463           1 :     array.release(&array);
    2464             : 
    2465             :     // Next batch ==> End of stream
    2466           1 :     ASSERT_EQ(stream.get_next(&stream, &array), 0);
    2467           1 :     ASSERT_TRUE(array.release == nullptr);
    2468             : 
    2469             :     // Release stream
    2470           1 :     stream.release(&stream);
    2471             : 
    2472             :     // Get a new stream now that we've released it
    2473           1 :     ASSERT_TRUE(
    2474             :         OGR_L_GetArrowStream(OGRLayer::ToHandle(poLayer), &stream, nullptr));
    2475           1 :     ASSERT_TRUE(stream.release != nullptr);
    2476             : 
    2477             :     // Free dataset & layer
    2478           1 :     poDS.reset();
    2479             : 
    2480             :     // Test releasing the stream after the dataset/layer has been closed
    2481           1 :     CPLPushErrorHandler(CPLQuietErrorHandler);
    2482           1 :     CPLErrorReset();
    2483           1 :     ASSERT_TRUE(stream.get_schema(&stream, &schema) != 0);
    2484           1 :     ASSERT_TRUE(stream.get_last_error(&stream) != nullptr);
    2485           1 :     ASSERT_TRUE(stream.get_next(&stream, &array) != 0);
    2486           1 :     CPLPopErrorHandler();
    2487           1 :     stream.release(&stream);
    2488             : }
    2489             : 
    2490             : // Test field domain cloning
    2491           4 : TEST_F(test_ogr, field_domain_cloning)
    2492             : {
    2493             :     // range domain
    2494             :     OGRField min;
    2495           1 :     min.Real = 5.5;
    2496             :     OGRField max;
    2497           1 :     max.Real = 6.5;
    2498             :     OGRRangeFieldDomain oRange("name", "description", OGRFieldType::OFTReal,
    2499             :                                OGRFieldSubType::OFSTBoolean, min, true, max,
    2500           2 :                                true);
    2501           1 :     oRange.SetMergePolicy(OGRFieldDomainMergePolicy::OFDMP_GEOMETRY_WEIGHTED);
    2502           1 :     oRange.SetSplitPolicy(OGRFieldDomainSplitPolicy::OFDSP_GEOMETRY_RATIO);
    2503           1 :     std::unique_ptr<OGRRangeFieldDomain> poClonedRange(oRange.Clone());
    2504           1 :     ASSERT_EQ(poClonedRange->GetName(), oRange.GetName());
    2505           1 :     ASSERT_EQ(poClonedRange->GetDescription(), oRange.GetDescription());
    2506           1 :     bool originalInclusive = false;
    2507           1 :     bool cloneInclusive = false;
    2508           1 :     ASSERT_EQ(poClonedRange->GetMin(originalInclusive).Real,
    2509             :               oRange.GetMin(cloneInclusive).Real);
    2510           1 :     ASSERT_EQ(originalInclusive, cloneInclusive);
    2511           1 :     ASSERT_EQ(poClonedRange->GetMax(originalInclusive).Real,
    2512             :               oRange.GetMax(cloneInclusive).Real);
    2513           1 :     ASSERT_EQ(originalInclusive, cloneInclusive);
    2514           1 :     ASSERT_EQ(poClonedRange->GetFieldType(), oRange.GetFieldType());
    2515           1 :     ASSERT_EQ(poClonedRange->GetFieldSubType(), oRange.GetFieldSubType());
    2516           1 :     ASSERT_EQ(poClonedRange->GetSplitPolicy(), oRange.GetSplitPolicy());
    2517           1 :     ASSERT_EQ(poClonedRange->GetMergePolicy(), oRange.GetMergePolicy());
    2518             : 
    2519             :     // glob domain
    2520             :     OGRGlobFieldDomain oGlob("name", "description", OGRFieldType::OFTString,
    2521           2 :                              OGRFieldSubType::OFSTBoolean, "*a*");
    2522           1 :     oGlob.SetMergePolicy(OGRFieldDomainMergePolicy::OFDMP_GEOMETRY_WEIGHTED);
    2523           1 :     oGlob.SetSplitPolicy(OGRFieldDomainSplitPolicy::OFDSP_GEOMETRY_RATIO);
    2524           1 :     std::unique_ptr<OGRGlobFieldDomain> poClonedGlob(oGlob.Clone());
    2525           1 :     ASSERT_EQ(poClonedGlob->GetName(), oGlob.GetName());
    2526           1 :     ASSERT_EQ(poClonedGlob->GetDescription(), oGlob.GetDescription());
    2527           1 :     ASSERT_EQ(poClonedGlob->GetGlob(), oGlob.GetGlob());
    2528           1 :     ASSERT_EQ(poClonedGlob->GetFieldType(), oGlob.GetFieldType());
    2529           1 :     ASSERT_EQ(poClonedGlob->GetFieldSubType(), oGlob.GetFieldSubType());
    2530           1 :     ASSERT_EQ(poClonedGlob->GetSplitPolicy(), oGlob.GetSplitPolicy());
    2531           1 :     ASSERT_EQ(poClonedGlob->GetMergePolicy(), oGlob.GetMergePolicy());
    2532             : 
    2533             :     // coded value domain
    2534             :     OGRCodedFieldDomain oCoded("name", "description", OGRFieldType::OFTString,
    2535           2 :                                OGRFieldSubType::OFSTBoolean, {OGRCodedValue()});
    2536           1 :     oCoded.SetMergePolicy(OGRFieldDomainMergePolicy::OFDMP_GEOMETRY_WEIGHTED);
    2537           1 :     oCoded.SetSplitPolicy(OGRFieldDomainSplitPolicy::OFDSP_GEOMETRY_RATIO);
    2538           1 :     std::unique_ptr<OGRCodedFieldDomain> poClonedCoded(oCoded.Clone());
    2539           1 :     ASSERT_EQ(poClonedCoded->GetName(), oCoded.GetName());
    2540           1 :     ASSERT_EQ(poClonedCoded->GetDescription(), oCoded.GetDescription());
    2541           1 :     ASSERT_EQ(poClonedCoded->GetFieldType(), oCoded.GetFieldType());
    2542           1 :     ASSERT_EQ(poClonedCoded->GetFieldSubType(), oCoded.GetFieldSubType());
    2543           1 :     ASSERT_EQ(poClonedCoded->GetSplitPolicy(), oCoded.GetSplitPolicy());
    2544           1 :     ASSERT_EQ(poClonedCoded->GetMergePolicy(), oCoded.GetMergePolicy());
    2545             : }
    2546             : 
    2547             : // Test field comments
    2548           4 : TEST_F(test_ogr, field_comments)
    2549             : {
    2550           1 :     OGRFieldDefn oFieldDefn("field1", OFTString);
    2551           1 :     ASSERT_EQ(oFieldDefn.GetComment(), "");
    2552           1 :     oFieldDefn.SetComment("my comment");
    2553           1 :     ASSERT_EQ(oFieldDefn.GetComment(), "my comment");
    2554             : 
    2555           1 :     OGRFieldDefn oFieldDefn2(&oFieldDefn);
    2556           1 :     ASSERT_EQ(oFieldDefn2.GetComment(), "my comment");
    2557           1 :     ASSERT_TRUE(oFieldDefn.IsSame(&oFieldDefn2));
    2558             : 
    2559           1 :     oFieldDefn2.SetComment("my comment 2");
    2560           1 :     ASSERT_FALSE(oFieldDefn.IsSame(&oFieldDefn2));
    2561             : }
    2562             : 
    2563             : // Test OGRFeatureDefn C++ GetFields() iterator
    2564           4 : TEST_F(test_ogr, feature_defn_fields_iterator)
    2565             : {
    2566           2 :     OGRFeatureDefn oFDefn;
    2567             :     {
    2568           2 :         OGRFieldDefn oFieldDefn("field1", OFTString);
    2569           1 :         oFDefn.AddFieldDefn(&oFieldDefn);
    2570             :     }
    2571             :     {
    2572           2 :         OGRFieldDefn oFieldDefn("field2", OFTString);
    2573           1 :         oFDefn.AddFieldDefn(&oFieldDefn);
    2574             :     }
    2575           1 :     EXPECT_EQ(oFDefn.GetFields().size(),
    2576             :               static_cast<size_t>(oFDefn.GetFieldCount()));
    2577           1 :     int i = 0;
    2578           3 :     for (const auto *poFieldDefn : oFDefn.GetFields())
    2579             :     {
    2580           2 :         EXPECT_EQ(oFDefn.GetFields()[i], oFDefn.GetFieldDefn(i));
    2581           2 :         EXPECT_EQ(poFieldDefn, oFDefn.GetFieldDefn(i));
    2582           2 :         ++i;
    2583             :     }
    2584           1 :     EXPECT_EQ(i, oFDefn.GetFieldCount());
    2585           1 : }
    2586             : 
    2587             : // Test OGRFeatureDefn C++ GetGeomFields() iterator
    2588           4 : TEST_F(test_ogr, feature_defn_geomfields_iterator)
    2589             : {
    2590           2 :     OGRFeatureDefn oFDefn;
    2591             :     {
    2592           2 :         OGRGeomFieldDefn oGeomFieldDefn("field1", wkbUnknown);
    2593           1 :         oFDefn.AddGeomFieldDefn(&oGeomFieldDefn);
    2594             :     }
    2595             :     {
    2596           2 :         OGRGeomFieldDefn oGeomFieldDefn("field2", wkbUnknown);
    2597           1 :         oFDefn.AddGeomFieldDefn(&oGeomFieldDefn);
    2598             :     }
    2599           1 :     EXPECT_EQ(oFDefn.GetGeomFields().size(),
    2600             :               static_cast<size_t>(oFDefn.GetGeomFieldCount()));
    2601           1 :     int i = 0;
    2602           4 :     for (const auto *poGeomFieldDefn : oFDefn.GetGeomFields())
    2603             :     {
    2604           3 :         EXPECT_EQ(oFDefn.GetGeomFields()[i], oFDefn.GetGeomFieldDefn(i));
    2605           3 :         EXPECT_EQ(poGeomFieldDefn, oFDefn.GetGeomFieldDefn(i));
    2606           3 :         ++i;
    2607             :     }
    2608           1 :     EXPECT_EQ(i, oFDefn.GetGeomFieldCount());
    2609           1 : }
    2610             : 
    2611             : // Test OGRGeomFieldDefn copy constructor
    2612           4 : TEST_F(test_ogr, geom_field_defn_copy_constructor)
    2613             : {
    2614             :     {
    2615           2 :         OGRGeomFieldDefn oGeomFieldDefn("field1", wkbPoint);
    2616           1 :         oGeomFieldDefn.SetNullable(false);
    2617           2 :         OGRGeomFieldDefn oGeomFieldDefn2("field2", wkbLineString);
    2618           1 :         oGeomFieldDefn2 = oGeomFieldDefn;
    2619           1 :         EXPECT_TRUE(oGeomFieldDefn2.IsSame(&oGeomFieldDefn));
    2620             :     }
    2621             : 
    2622             :     {
    2623           2 :         OGRSpatialReference oSRS;
    2624           1 :         oSRS.SetFromUserInput("WGS84");
    2625           1 :         EXPECT_EQ(oSRS.GetReferenceCount(), 1);
    2626           2 :         OGRGeomFieldDefn oGeomFieldDefn("field1", wkbPoint);
    2627           1 :         oGeomFieldDefn.SetSpatialRef(&oSRS);
    2628           1 :         EXPECT_EQ(oSRS.GetReferenceCount(), 2);
    2629           2 :         OGRGeomFieldDefn oGeomFieldDefn2("field2", wkbLineString);
    2630           1 :         oGeomFieldDefn2 = oGeomFieldDefn;
    2631           1 :         EXPECT_EQ(oSRS.GetReferenceCount(), 3);
    2632           1 :         EXPECT_TRUE(oGeomFieldDefn2.IsSame(&oGeomFieldDefn));
    2633             : 
    2634             :         // oGeomFieldDefn2 already points to oSRS
    2635           1 :         oGeomFieldDefn2 = oGeomFieldDefn;
    2636           1 :         EXPECT_EQ(oSRS.GetReferenceCount(), 3);
    2637           1 :         EXPECT_TRUE(oGeomFieldDefn2.IsSame(&oGeomFieldDefn));
    2638             :     }
    2639           1 : }
    2640             : 
    2641             : // Test GDALDataset QueryLoggerFunc callback
    2642           4 : TEST_F(test_ogr, GDALDatasetSetQueryLoggerFunc)
    2643             : {
    2644           1 :     if (GDALGetDriverByName("GPKG") == nullptr)
    2645             :     {
    2646           0 :         GTEST_SKIP() << "GPKG driver missing";
    2647             :     }
    2648             : 
    2649           1 :     auto tmpGPKG{testing::TempDir() + "/poly-1-feature.gpkg"};
    2650             :     {
    2651           2 :         std::string srcfilename(data_ + SEP + "poly-1-feature.gpkg");
    2652           2 :         std::ifstream src(srcfilename, std::ios::binary);
    2653           2 :         std::ofstream dst(tmpGPKG, std::ios::binary);
    2654           1 :         dst << src.rdbuf();
    2655             :     }
    2656             : 
    2657             :     struct QueryLogEntry
    2658             :     {
    2659             :         std::string sql;
    2660             :         std::string error;
    2661             :         int64_t numRecords;
    2662             :         int64_t executionTimeMilliseconds;
    2663             :     };
    2664             : 
    2665             :     // Note: this must be constructed before poDS or the order
    2666             :     //       of destruction will make the callback call the already
    2667             :     //       destructed vector
    2668           1 :     std::vector<QueryLogEntry> queryLog;
    2669             : 
    2670             :     auto poDS = std::unique_ptr<GDALDataset>(
    2671           1 :         GDALDataset::Open(tmpGPKG.c_str(), GDAL_OF_VECTOR | GDAL_OF_UPDATE));
    2672           1 :     ASSERT_TRUE(poDS);
    2673           1 :     auto hDS = GDALDataset::ToHandle(poDS.get());
    2674           1 :     ASSERT_TRUE(hDS);
    2675             : 
    2676           1 :     const bool retVal = GDALDatasetSetQueryLoggerFunc(
    2677             :         hDS,
    2678          18 :         [](const char *pszSQL, const char *pszError, int64_t lNumRecords,
    2679             :            int64_t lExecutionTimeMilliseconds, void *pQueryLoggerArg)
    2680             :         {
    2681          18 :             std::vector<QueryLogEntry> *queryLogLocal{
    2682             :                 reinterpret_cast<std::vector<QueryLogEntry> *>(
    2683             :                     pQueryLoggerArg)};
    2684          36 :             QueryLogEntry entryLocal;
    2685          18 :             if (pszSQL)
    2686             :             {
    2687          18 :                 entryLocal.sql = pszSQL;
    2688             :             }
    2689          18 :             entryLocal.numRecords = lNumRecords;
    2690          18 :             entryLocal.executionTimeMilliseconds = lExecutionTimeMilliseconds;
    2691          18 :             if (pszError)
    2692             :             {
    2693           1 :                 entryLocal.error = pszError;
    2694             :             }
    2695          18 :             queryLogLocal->push_back(std::move(entryLocal));
    2696          18 :         },
    2697           1 :         &queryLog);
    2698             : 
    2699           1 :     ASSERT_TRUE(retVal);
    2700           1 :     auto hLayer{GDALDatasetGetLayer(hDS, 0)};
    2701           1 :     ASSERT_TRUE(hLayer);
    2702           1 :     ASSERT_STREQ(OGR_L_GetName(hLayer), "poly");
    2703             :     auto poFeature = std::unique_ptr<OGRFeature>(
    2704           1 :         OGRFeature::FromHandle(OGR_L_GetNextFeature(hLayer)));
    2705           1 :     auto hFeature = OGRFeature::ToHandle(poFeature.get());
    2706           1 :     ASSERT_TRUE(hFeature);
    2707           1 :     ASSERT_GT(queryLog.size(), 1U);
    2708             : 
    2709           1 :     QueryLogEntry entry{queryLog.back()};
    2710           1 :     ASSERT_EQ(entry.sql.find("SELECT", 0), 0);
    2711           1 :     ASSERT_TRUE(entry.executionTimeMilliseconds >= 0);
    2712           1 :     ASSERT_EQ(entry.numRecords, -1);
    2713           1 :     ASSERT_TRUE(entry.error.empty());
    2714             : 
    2715             :     // Test erroneous query
    2716           1 :     OGRLayerH queryResultLayerH{GDALDatasetExecuteSQL(
    2717             :         hDS, "SELECT * FROM not_existing_table", nullptr, nullptr)};
    2718           1 :     GDALDatasetReleaseResultSet(hDS, queryResultLayerH);
    2719           1 :     ASSERT_FALSE(queryResultLayerH);
    2720             : 
    2721           1 :     entry = queryLog.back();
    2722           1 :     ASSERT_EQ(entry.sql.find("SELECT * FROM not_existing_table", 0), 0);
    2723           1 :     ASSERT_EQ(entry.executionTimeMilliseconds, -1);
    2724           1 :     ASSERT_EQ(entry.numRecords, -1);
    2725           1 :     ASSERT_FALSE(entry.error.empty());
    2726             : 
    2727             :     // Test prepared arg substitution
    2728           1 :     hFeature = OGR_F_Create(OGR_L_GetLayerDefn(hLayer));
    2729           1 :     poFeature.reset(OGRFeature::FromHandle(hFeature));
    2730           1 :     OGR_F_SetFieldInteger(hFeature, 1, 123);
    2731           1 :     OGRErr err = OGR_L_CreateFeature(hLayer, hFeature);
    2732           1 :     ASSERT_EQ(OGRERR_NONE, err);
    2733             : 
    2734             :     auto insertEntry = std::find_if(
    2735          16 :         queryLog.cbegin(), queryLog.cend(), [](const QueryLogEntry &e)
    2736          17 :         { return e.sql.find(R"sql(INSERT INTO "poly")sql", 0) == 0; });
    2737             : 
    2738           1 :     ASSERT_TRUE(insertEntry != queryLog.end());
    2739           1 :     ASSERT_EQ(
    2740             :         insertEntry->sql.find(
    2741             :             R"sql(INSERT INTO "poly" ( "geom", "AREA", "EAS_ID", "PRFEDEA") VALUES (NULL, NULL, 123, NULL))sql",
    2742             :             0),
    2743             :         0);
    2744             : }
    2745             : 
    2746           4 : TEST_F(test_ogr, OGRParseDateTimeYYYYMMDDTHHMMZ)
    2747             : {
    2748             :     {
    2749           1 :         char szInput[] = "2023-07-11T17:27Z";
    2750             :         OGRField sField;
    2751           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMZ(szInput, &sField), true);
    2752           1 :         EXPECT_EQ(sField.Date.Year, 2023);
    2753           1 :         EXPECT_EQ(sField.Date.Month, 7);
    2754           1 :         EXPECT_EQ(sField.Date.Day, 11);
    2755           1 :         EXPECT_EQ(sField.Date.Hour, 17);
    2756           1 :         EXPECT_EQ(sField.Date.Minute, 27);
    2757           1 :         EXPECT_EQ(sField.Date.Second, 0.0f);
    2758           1 :         EXPECT_EQ(sField.Date.TZFlag, 100);
    2759             :     }
    2760             :     {
    2761           1 :         char szInput[] = "2023-07-11T17:27";
    2762             :         OGRField sField;
    2763           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMZ(szInput, &sField), true);
    2764           1 :         EXPECT_EQ(sField.Date.Year, 2023);
    2765           1 :         EXPECT_EQ(sField.Date.Month, 7);
    2766           1 :         EXPECT_EQ(sField.Date.Day, 11);
    2767           1 :         EXPECT_EQ(sField.Date.Hour, 17);
    2768           1 :         EXPECT_EQ(sField.Date.Minute, 27);
    2769           1 :         EXPECT_EQ(sField.Date.Second, 0.0f);
    2770           1 :         EXPECT_EQ(sField.Date.TZFlag, 0);
    2771             :     }
    2772             :     {
    2773             :         // Invalid
    2774           1 :         char szInput[] = "2023-07-11T17:2";
    2775             :         OGRField sField;
    2776           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMZ(szInput, &sField), false);
    2777             :     }
    2778             :     {
    2779             :         // Invalid
    2780           1 :         char szInput[] = "2023-07-11T17:99";
    2781             :         OGRField sField;
    2782           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMZ(szInput, &sField), false);
    2783             :     }
    2784           1 : }
    2785             : 
    2786           4 : TEST_F(test_ogr, OGRParseDateTimeYYYYMMDDTHHMMSSZ)
    2787             : {
    2788             :     {
    2789           1 :         char szInput[] = "2023-07-11T17:27:34Z";
    2790             :         OGRField sField;
    2791           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMSSZ(szInput, &sField), true);
    2792           1 :         EXPECT_EQ(sField.Date.Year, 2023);
    2793           1 :         EXPECT_EQ(sField.Date.Month, 7);
    2794           1 :         EXPECT_EQ(sField.Date.Day, 11);
    2795           1 :         EXPECT_EQ(sField.Date.Hour, 17);
    2796           1 :         EXPECT_EQ(sField.Date.Minute, 27);
    2797           1 :         EXPECT_EQ(sField.Date.Second, 34.0f);
    2798           1 :         EXPECT_EQ(sField.Date.TZFlag, 100);
    2799             :     }
    2800             :     {
    2801           1 :         char szInput[] = "2023-07-11T17:27:34";
    2802             :         OGRField sField;
    2803           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMSSZ(szInput, &sField), true);
    2804           1 :         EXPECT_EQ(sField.Date.Year, 2023);
    2805           1 :         EXPECT_EQ(sField.Date.Month, 7);
    2806           1 :         EXPECT_EQ(sField.Date.Day, 11);
    2807           1 :         EXPECT_EQ(sField.Date.Hour, 17);
    2808           1 :         EXPECT_EQ(sField.Date.Minute, 27);
    2809           1 :         EXPECT_EQ(sField.Date.Second, 34.0f);
    2810           1 :         EXPECT_EQ(sField.Date.TZFlag, 0);
    2811             :     }
    2812             :     {
    2813             :         // Invalid
    2814           1 :         char szInput[] = "2023-07-11T17:27:3";
    2815             :         OGRField sField;
    2816           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMSSZ(szInput, &sField), false);
    2817             :     }
    2818             :     {
    2819             :         // Invalid
    2820           1 :         char szInput[] = "2023-07-11T17:27:99";
    2821             :         OGRField sField;
    2822           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMSSZ(szInput, &sField), false);
    2823             :     }
    2824           1 : }
    2825             : 
    2826           4 : TEST_F(test_ogr, OGRParseDateTimeYYYYMMDDTHHMMSSsssZ)
    2827             : {
    2828             :     {
    2829           1 :         char szInput[] = "2023-07-11T17:27:34.123Z";
    2830             :         OGRField sField;
    2831           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMSSsssZ(szInput, &sField), true);
    2832           1 :         EXPECT_EQ(sField.Date.Year, 2023);
    2833           1 :         EXPECT_EQ(sField.Date.Month, 7);
    2834           1 :         EXPECT_EQ(sField.Date.Day, 11);
    2835           1 :         EXPECT_EQ(sField.Date.Hour, 17);
    2836           1 :         EXPECT_EQ(sField.Date.Minute, 27);
    2837           1 :         EXPECT_EQ(sField.Date.Second, 34.123f);
    2838           1 :         EXPECT_EQ(sField.Date.TZFlag, 100);
    2839             :     }
    2840             :     {
    2841           1 :         char szInput[] = "2023-07-11T17:27:34.123";
    2842             :         OGRField sField;
    2843           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMSSsssZ(szInput, &sField), true);
    2844           1 :         EXPECT_EQ(sField.Date.Year, 2023);
    2845           1 :         EXPECT_EQ(sField.Date.Month, 7);
    2846           1 :         EXPECT_EQ(sField.Date.Day, 11);
    2847           1 :         EXPECT_EQ(sField.Date.Hour, 17);
    2848           1 :         EXPECT_EQ(sField.Date.Minute, 27);
    2849           1 :         EXPECT_EQ(sField.Date.Second, 34.123f);
    2850           1 :         EXPECT_EQ(sField.Date.TZFlag, 0);
    2851             :     }
    2852             :     {
    2853             :         // Invalid
    2854           1 :         char szInput[] = "2023-07-11T17:27:34.12";
    2855             :         OGRField sField;
    2856           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMSSsssZ(szInput, &sField), false);
    2857             :     }
    2858             :     {
    2859             :         // Invalid
    2860           1 :         char szInput[] = "2023-07-11T17:27:99.123";
    2861             :         OGRField sField;
    2862           1 :         EXPECT_EQ(OGRParseDateTimeYYYYMMDDTHHMMSSsssZ(szInput, &sField), false);
    2863             :     }
    2864           1 : }
    2865             : 
    2866           4 : TEST_F(test_ogr, OGRGetISO8601DateTime)
    2867             : {
    2868             :     OGRField sField;
    2869           1 :     sField.Date.Year = 2023;
    2870           1 :     sField.Date.Month = 7;
    2871           1 :     sField.Date.Day = 11;
    2872           1 :     sField.Date.Hour = 17;
    2873           1 :     sField.Date.Minute = 27;
    2874           1 :     sField.Date.Second = 34.567f;
    2875           1 :     sField.Date.TZFlag = 100;
    2876             :     {
    2877             :         char szResult[OGR_SIZEOF_ISO8601_DATETIME_BUFFER];
    2878             :         OGRISO8601Format sFormat;
    2879           1 :         sFormat.ePrecision = OGRISO8601Precision::AUTO;
    2880           1 :         OGRGetISO8601DateTime(&sField, sFormat, szResult);
    2881           1 :         EXPECT_STREQ(szResult, "2023-07-11T17:27:34.567Z");
    2882             :     }
    2883             :     {
    2884             :         char szResult[OGR_SIZEOF_ISO8601_DATETIME_BUFFER];
    2885             :         OGRISO8601Format sFormat;
    2886           1 :         sFormat.ePrecision = OGRISO8601Precision::MILLISECOND;
    2887           1 :         OGRGetISO8601DateTime(&sField, sFormat, szResult);
    2888           1 :         EXPECT_STREQ(szResult, "2023-07-11T17:27:34.567Z");
    2889             :     }
    2890             :     {
    2891             :         char szResult[OGR_SIZEOF_ISO8601_DATETIME_BUFFER];
    2892             :         OGRISO8601Format sFormat;
    2893           1 :         sFormat.ePrecision = OGRISO8601Precision::SECOND;
    2894           1 :         OGRGetISO8601DateTime(&sField, sFormat, szResult);
    2895           1 :         EXPECT_STREQ(szResult, "2023-07-11T17:27:35Z");
    2896             :     }
    2897             :     {
    2898             :         char szResult[OGR_SIZEOF_ISO8601_DATETIME_BUFFER];
    2899             :         OGRISO8601Format sFormat;
    2900           1 :         sFormat.ePrecision = OGRISO8601Precision::MINUTE;
    2901           1 :         OGRGetISO8601DateTime(&sField, sFormat, szResult);
    2902           1 :         EXPECT_STREQ(szResult, "2023-07-11T17:27Z");
    2903             :     }
    2904           1 :     sField.Date.Second = 34.0f;
    2905             :     {
    2906             :         char szResult[OGR_SIZEOF_ISO8601_DATETIME_BUFFER];
    2907             :         OGRISO8601Format sFormat;
    2908           1 :         sFormat.ePrecision = OGRISO8601Precision::AUTO;
    2909           1 :         OGRGetISO8601DateTime(&sField, sFormat, szResult);
    2910           1 :         EXPECT_STREQ(szResult, "2023-07-11T17:27:34Z");
    2911             :     }
    2912           1 : }
    2913             : 
    2914             : // Test calling importFromWkb() multiple times on the same geometry object
    2915           4 : TEST_F(test_ogr, importFromWkbReuse)
    2916             : {
    2917             :     {
    2918           2 :         OGRPoint oPoint;
    2919             :         {
    2920           1 :             size_t nBytesConsumed = 0;
    2921           1 :             EXPECT_EQ(oPoint.importFromWkb(
    2922             :                           reinterpret_cast<const GByte *>(
    2923             :                               "\x01\x01\x00\x00\x00"                // Point
    2924             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"    // 1.0
    2925             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"),  // 2.0
    2926             :                           21, wkbVariantIso, nBytesConsumed),
    2927             :                       OGRERR_NONE);
    2928           1 :             EXPECT_EQ(nBytesConsumed, 21);
    2929           1 :             EXPECT_EQ(oPoint.getX(), 1.0);
    2930           1 :             EXPECT_EQ(oPoint.getY(), 2.0);
    2931             :         }
    2932             :         {
    2933           1 :             size_t nBytesConsumed = 0;
    2934           1 :             EXPECT_EQ(oPoint.importFromWkb(
    2935             :                           reinterpret_cast<const GByte *>(
    2936             :                               "\x01\x01\x00\x00\x00"              // Point
    2937             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"  // 2.0
    2938             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"  // 1.0
    2939             :                               ),
    2940             :                           21, wkbVariantIso, nBytesConsumed),
    2941             :                       OGRERR_NONE);
    2942           1 :             EXPECT_EQ(nBytesConsumed, 21);
    2943           1 :             EXPECT_EQ(oPoint.getX(), 2.0);
    2944           1 :             EXPECT_EQ(oPoint.getY(), 1.0);
    2945             :         }
    2946             :     }
    2947             : 
    2948             :     {
    2949           1 :         OGRLineString oLS;
    2950             :         {
    2951           1 :             size_t nBytesConsumed = 0;
    2952           1 :             EXPECT_EQ(oLS.importFromWkb(
    2953             :                           reinterpret_cast<const GByte *>(
    2954             :                               "\x01\x02\x00\x00\x00"              // LineString
    2955             :                               "\x01\x00\x00\x00"                  // 1 point
    2956             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"  // 1.0
    2957             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"),  // 2.0
    2958             :                           25, wkbVariantIso, nBytesConsumed),
    2959             :                       OGRERR_NONE);
    2960           1 :             EXPECT_EQ(nBytesConsumed, 25);
    2961           1 :             ASSERT_EQ(oLS.getNumPoints(), 1);
    2962           1 :             EXPECT_EQ(oLS.getX(0), 1.0);
    2963           1 :             EXPECT_EQ(oLS.getY(0), 2.0);
    2964             :         }
    2965             :         {
    2966           1 :             size_t nBytesConsumed = 0;
    2967           1 :             EXPECT_EQ(oLS.importFromWkb(
    2968             :                           reinterpret_cast<const GByte *>(
    2969             :                               "\x01\x02\x00\x00\x00"              // LineString
    2970             :                               "\x02\x00\x00\x00"                  // 2 points
    2971             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"  // 1.0
    2972             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"  // 2.0
    2973             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"  // 2.0
    2974             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"),  // 1.0
    2975             :                           41, wkbVariantIso, nBytesConsumed),
    2976             :                       OGRERR_NONE);
    2977           1 :             EXPECT_EQ(nBytesConsumed, 41);
    2978           1 :             ASSERT_EQ(oLS.getNumPoints(), 2);
    2979           1 :             EXPECT_EQ(oLS.getX(0), 1.0);
    2980           1 :             EXPECT_EQ(oLS.getY(0), 2.0);
    2981           1 :             EXPECT_EQ(oLS.getX(1), 2.0);
    2982           1 :             EXPECT_EQ(oLS.getY(1), 1.0);
    2983             :         }
    2984             :         {
    2985           1 :             size_t nBytesConsumed = 0;
    2986           1 :             EXPECT_EQ(oLS.importFromWkb(
    2987             :                           reinterpret_cast<const GByte *>(
    2988             :                               "\x01\x02\x00\x00\x00"              // LineString
    2989             :                               "\x01\x00\x00\x00"                  // 1 point
    2990             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"  // 2.0
    2991             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"),  // 1.0
    2992             :                           25, wkbVariantIso, nBytesConsumed),
    2993             :                       OGRERR_NONE);
    2994           1 :             EXPECT_EQ(nBytesConsumed, 25);
    2995           1 :             ASSERT_EQ(oLS.getNumPoints(), 1);
    2996           1 :             EXPECT_EQ(oLS.getX(0), 2.0);
    2997           1 :             EXPECT_EQ(oLS.getY(0), 1.0);
    2998             :         }
    2999             :     }
    3000             : 
    3001             :     {
    3002           1 :         OGRPolygon oPoly;
    3003             :         {
    3004           1 :             size_t nBytesConsumed = 0;
    3005           1 :             EXPECT_EQ(oPoly.importFromWkb(
    3006             :                           reinterpret_cast<const GByte *>(
    3007             :                               "\x01\x03\x00\x00\x00"                // Polygon
    3008             :                               "\x01\x00\x00\x00"                    // 1 ring
    3009             :                               "\x01\x00\x00\x00"                    // 1 point
    3010             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"    // 1.0
    3011             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"),  // 2.0
    3012             :                           29, wkbVariantIso, nBytesConsumed),
    3013             :                       OGRERR_NONE);
    3014           1 :             EXPECT_EQ(nBytesConsumed, 29);
    3015           1 :             ASSERT_TRUE(oPoly.getExteriorRing() != nullptr);
    3016           1 :             ASSERT_EQ(oPoly.getNumInteriorRings(), 0);
    3017           1 :             auto poLS = oPoly.getExteriorRing();
    3018           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3019           1 :             EXPECT_EQ(poLS->getX(0), 1.0);
    3020           1 :             EXPECT_EQ(poLS->getY(0), 2.0);
    3021             :         }
    3022             :         {
    3023           1 :             size_t nBytesConsumed = 0;
    3024           1 :             EXPECT_EQ(oPoly.importFromWkb(
    3025             :                           reinterpret_cast<const GByte *>(
    3026             :                               "\x01\x03\x00\x00\x00"                // Polygon
    3027             :                               "\x01\x00\x00\x00"                    // 1 ring
    3028             :                               "\x01\x00\x00\x00"                    // 1 point
    3029             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"    // 2.0
    3030             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"),  // 1.0
    3031             :                           29, wkbVariantIso, nBytesConsumed),
    3032             :                       OGRERR_NONE);
    3033           1 :             EXPECT_EQ(nBytesConsumed, 29);
    3034           1 :             ASSERT_TRUE(oPoly.getExteriorRing() != nullptr);
    3035           1 :             ASSERT_EQ(oPoly.getNumInteriorRings(), 0);
    3036           1 :             auto poLS = oPoly.getExteriorRing();
    3037           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3038           1 :             EXPECT_EQ(poLS->getX(0), 2.0);
    3039           1 :             EXPECT_EQ(poLS->getY(0), 1.0);
    3040             :         }
    3041             :         {
    3042           1 :             size_t nBytesConsumed = 0;
    3043           1 :             EXPECT_EQ(oPoly.importFromWkb(reinterpret_cast<const GByte *>(
    3044             :                                               "\x01\x03\x00\x00\x00"  // Polygon
    3045             :                                               "\x00\x00\x00\x00"),    // 0 ring
    3046             :                                           9, wkbVariantIso, nBytesConsumed),
    3047             :                       OGRERR_NONE);
    3048           1 :             EXPECT_EQ(nBytesConsumed, 9);
    3049           1 :             ASSERT_TRUE(oPoly.getExteriorRing() == nullptr);
    3050           1 :             ASSERT_EQ(oPoly.getNumInteriorRings(), 0);
    3051             :         }
    3052             :         {
    3053           1 :             size_t nBytesConsumed = 0;
    3054           1 :             EXPECT_EQ(oPoly.importFromWkb(
    3055             :                           reinterpret_cast<const GByte *>(
    3056             :                               "\x01\x03\x00\x00\x00"                // Polygon
    3057             :                               "\x01\x00\x00\x00"                    // 1 ring
    3058             :                               "\x01\x00\x00\x00"                    // 1 point
    3059             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"    // 1.0
    3060             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"),  // 2.0
    3061             :                           29, wkbVariantIso, nBytesConsumed),
    3062             :                       OGRERR_NONE);
    3063           1 :             EXPECT_EQ(nBytesConsumed, 29);
    3064           1 :             ASSERT_TRUE(oPoly.getExteriorRing() != nullptr);
    3065           1 :             ASSERT_EQ(oPoly.getNumInteriorRings(), 0);
    3066           1 :             auto poLS = oPoly.getExteriorRing();
    3067           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3068           1 :             EXPECT_EQ(poLS->getX(0), 1.0);
    3069           1 :             EXPECT_EQ(poLS->getY(0), 2.0);
    3070             :         }
    3071             :         {
    3072           1 :             size_t nBytesConsumed = 0;
    3073           1 :             EXPECT_EQ(oPoly.importFromWkb(
    3074             :                           reinterpret_cast<const GByte *>(
    3075             :                               "\x01\x03\x00\x00\x00"                // Polygon
    3076             :                               "\x01\x00\x00\x00"                    // 1 ring
    3077             :                               "\x01\x00\x00\x00"                    // 1 point
    3078             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"    // 2.0
    3079             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"),  // 1.0
    3080             :                           static_cast<size_t>(-1), wkbVariantIso,
    3081             :                           nBytesConsumed),
    3082             :                       OGRERR_NONE);
    3083           1 :             EXPECT_EQ(nBytesConsumed, 29);
    3084           1 :             ASSERT_TRUE(oPoly.getExteriorRing() != nullptr);
    3085           1 :             ASSERT_EQ(oPoly.getNumInteriorRings(), 0);
    3086           1 :             auto poLS = oPoly.getExteriorRing();
    3087           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3088           1 :             EXPECT_EQ(poLS->getX(0), 2.0);
    3089           1 :             EXPECT_EQ(poLS->getY(0), 1.0);
    3090             :         }
    3091             :         {
    3092           1 :             size_t nBytesConsumed = 0;
    3093             :             // Truncated WKB
    3094           1 :             EXPECT_NE(oPoly.importFromWkb(reinterpret_cast<const GByte *>(
    3095             :                                               "\x01\x03\x00\x00\x00"  // Polygon
    3096             :                                               "\x01\x00\x00\x00"      // 1 ring
    3097             :                                               "\x01\x00\x00\x00"),    // 1 point
    3098             :                                           13, wkbVariantIso, nBytesConsumed),
    3099             :                       OGRERR_NONE);
    3100           1 :             ASSERT_TRUE(oPoly.getExteriorRing() == nullptr);
    3101           1 :             ASSERT_EQ(oPoly.getNumInteriorRings(), 0);
    3102             :         }
    3103             :     }
    3104             : 
    3105             :     {
    3106           1 :         OGRMultiLineString oMLS;
    3107             :         {
    3108           1 :             size_t nBytesConsumed = 0;
    3109           1 :             EXPECT_EQ(oMLS.importFromWkb(
    3110             :                           reinterpret_cast<const GByte *>(
    3111             :                               "\x01\x05\x00\x00\x00"  // MultiLineString
    3112             :                               "\x01\x00\x00\x00"      // 1-part
    3113             :                               "\x01\x02\x00\x00\x00"  // LineString
    3114             :                               "\x01\x00\x00\x00"      // 1 point
    3115             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"    // 1.0
    3116             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"),  // 2.0
    3117             :                           34, wkbVariantIso, nBytesConsumed),
    3118             :                       OGRERR_NONE);
    3119           1 :             EXPECT_EQ(nBytesConsumed, 34);
    3120           1 :             ASSERT_EQ(oMLS.getNumGeometries(), 1);
    3121           1 :             auto poLS = oMLS.getGeometryRef(0);
    3122           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3123           1 :             EXPECT_EQ(poLS->getX(0), 1.0);
    3124           1 :             EXPECT_EQ(poLS->getY(0), 2.0);
    3125             :         }
    3126             :         {
    3127           1 :             size_t nBytesConsumed = 0;
    3128           1 :             EXPECT_EQ(oMLS.importFromWkb(
    3129             :                           reinterpret_cast<const GByte *>(
    3130             :                               "\x01\x05\x00\x00\x00"  // MultiLineString
    3131             :                               "\x01\x00\x00\x00"      // 1-part
    3132             :                               "\x01\x02\x00\x00\x00"  // LineString
    3133             :                               "\x01\x00\x00\x00"      // 1 point
    3134             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"    // 2.0
    3135             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"),  // 1.0
    3136             :                           34, wkbVariantIso, nBytesConsumed),
    3137             :                       OGRERR_NONE);
    3138           1 :             EXPECT_EQ(nBytesConsumed, 34);
    3139           1 :             ASSERT_EQ(oMLS.getNumGeometries(), 1);
    3140           1 :             auto poLS = oMLS.getGeometryRef(0);
    3141           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3142           1 :             EXPECT_EQ(poLS->getX(0), 2.0);
    3143           1 :             EXPECT_EQ(poLS->getY(0), 1.0);
    3144             :         }
    3145             :         {
    3146           1 :             size_t nBytesConsumed = 0;
    3147           1 :             EXPECT_EQ(oMLS.importFromWkb(
    3148             :                           reinterpret_cast<const GByte *>(
    3149             :                               "\x01\x05\x00\x00\x00"  // MultiLineString
    3150             :                               "\x01\x00\x00\x00"      // 1-part
    3151             :                               "\x01\x02\x00\x00\x00"  // LineString
    3152             :                               "\x01\x00\x00\x00"      // 1 point
    3153             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"    // 1.0
    3154             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"),  // 2.0
    3155             :                           static_cast<size_t>(-1), wkbVariantIso,
    3156             :                           nBytesConsumed),
    3157             :                       OGRERR_NONE);
    3158           1 :             EXPECT_EQ(nBytesConsumed, 34);
    3159           1 :             ASSERT_EQ(oMLS.getNumGeometries(), 1);
    3160           1 :             auto poLS = oMLS.getGeometryRef(0);
    3161           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3162           1 :             EXPECT_EQ(poLS->getX(0), 1.0);
    3163           1 :             EXPECT_EQ(poLS->getY(0), 2.0);
    3164             :         }
    3165             :         {
    3166           1 :             size_t nBytesConsumed = 0;
    3167             :             // Truncated WKB
    3168           1 :             EXPECT_NE(oMLS.importFromWkb(
    3169             :                           reinterpret_cast<const GByte *>(
    3170             :                               "\x01\x05\x00\x00\x00"  // MultiLineString
    3171             :                               "\x01\x00\x00\x00"      // 1-part
    3172             :                               "\x01\x02\x00\x00\x00"  // LineString
    3173             :                               "\x01\x00\x00\x00"      // 1 point
    3174             :                               ),
    3175             :                           18, wkbVariantIso, nBytesConsumed),
    3176             :                       OGRERR_NONE);
    3177           1 :             ASSERT_EQ(oMLS.getNumGeometries(), 0);
    3178             :         }
    3179             :     }
    3180             : 
    3181             :     {
    3182           1 :         OGRMultiPolygon oMP;
    3183             :         {
    3184           1 :             size_t nBytesConsumed = 0;
    3185           1 :             EXPECT_EQ(oMP.importFromWkb(
    3186             :                           reinterpret_cast<const GByte *>(
    3187             :                               "\x01\x06\x00\x00\x00"  // MultiPolygon
    3188             :                               "\x01\x00\x00\x00"      // 1-part
    3189             :                               "\x01\x03\x00\x00\x00"  // Polygon
    3190             :                               "\x01\x00\x00\x00"      // 1 ring
    3191             :                               "\x01\x00\x00\x00"      // 1 point
    3192             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"    // 1.0
    3193             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"),  // 2.0
    3194             :                           38, wkbVariantIso, nBytesConsumed),
    3195             :                       OGRERR_NONE);
    3196           1 :             EXPECT_EQ(nBytesConsumed, 38);
    3197           1 :             ASSERT_EQ(oMP.getNumGeometries(), 1);
    3198           1 :             auto poPoly = oMP.getGeometryRef(0);
    3199           1 :             ASSERT_TRUE(poPoly->getExteriorRing() != nullptr);
    3200           1 :             ASSERT_EQ(poPoly->getNumInteriorRings(), 0);
    3201           1 :             auto poLS = poPoly->getExteriorRing();
    3202           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3203           1 :             EXPECT_EQ(poLS->getX(0), 1.0);
    3204           1 :             EXPECT_EQ(poLS->getY(0), 2.0);
    3205             :         }
    3206             :         {
    3207           1 :             size_t nBytesConsumed = 0;
    3208           1 :             EXPECT_EQ(oMP.importFromWkb(
    3209             :                           reinterpret_cast<const GByte *>(
    3210             :                               "\x01\x06\x00\x00\x00"  // MultiPolygon
    3211             :                               "\x01\x00\x00\x00"      // 1-part
    3212             :                               "\x01\x03\x00\x00\x00"  // Polygon
    3213             :                               "\x01\x00\x00\x00"      // 1 ring
    3214             :                               "\x01\x00\x00\x00"      // 1 point
    3215             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"    // 2.0
    3216             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"),  // 1.0
    3217             :                           38, wkbVariantIso, nBytesConsumed),
    3218             :                       OGRERR_NONE);
    3219           1 :             EXPECT_EQ(nBytesConsumed, 38);
    3220           1 :             ASSERT_EQ(oMP.getNumGeometries(), 1);
    3221           1 :             auto poPoly = oMP.getGeometryRef(0);
    3222           1 :             ASSERT_TRUE(poPoly->getExteriorRing() != nullptr);
    3223           1 :             ASSERT_EQ(poPoly->getNumInteriorRings(), 0);
    3224           1 :             auto poLS = poPoly->getExteriorRing();
    3225           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3226           1 :             EXPECT_EQ(poLS->getX(0), 2.0);
    3227           1 :             EXPECT_EQ(poLS->getY(0), 1.0);
    3228             :         }
    3229             :         {
    3230           1 :             size_t nBytesConsumed = 0;
    3231           1 :             EXPECT_EQ(oMP.importFromWkb(
    3232             :                           reinterpret_cast<const GByte *>(
    3233             :                               "\x01\x06\x00\x00\x00"  // MultiPolygon
    3234             :                               "\x01\x00\x00\x00"      // 1-part
    3235             :                               "\x01\x03\x00\x00\x00"  // Polygon
    3236             :                               "\x01\x00\x00\x00"      // 1 ring
    3237             :                               "\x01\x00\x00\x00"      // 1 point
    3238             :                               "\x00\x00\x00\x00\x00\x00\xf0\x3f"    // 1.0
    3239             :                               "\x00\x00\x00\x00\x00\x00\x00\x40"),  // 2.0
    3240             :                           static_cast<size_t>(-1), wkbVariantIso,
    3241             :                           nBytesConsumed),
    3242             :                       OGRERR_NONE);
    3243           1 :             EXPECT_EQ(nBytesConsumed, 38);
    3244           1 :             ASSERT_EQ(oMP.getNumGeometries(), 1);
    3245           1 :             auto poPoly = oMP.getGeometryRef(0);
    3246           1 :             ASSERT_TRUE(poPoly->getExteriorRing() != nullptr);
    3247           1 :             ASSERT_EQ(poPoly->getNumInteriorRings(), 0);
    3248           1 :             auto poLS = poPoly->getExteriorRing();
    3249           1 :             ASSERT_EQ(poLS->getNumPoints(), 1);
    3250           1 :             EXPECT_EQ(poLS->getX(0), 1.0);
    3251           1 :             EXPECT_EQ(poLS->getY(0), 2.0);
    3252             :         }
    3253             :         {
    3254           1 :             size_t nBytesConsumed = 0;
    3255             :             // Truncated WKB
    3256           1 :             EXPECT_NE(
    3257             :                 oMP.importFromWkb(reinterpret_cast<const GByte *>(
    3258             :                                       "\x01\x06\x00\x00\x00"  // MultiPolygon
    3259             :                                       "\x01\x00\x00\x00"      // 1-part
    3260             :                                       "\x01\x03\x00\x00\x00"  // Polygon
    3261             :                                       "\x01\x00\x00\x00"      // 1 ring
    3262             :                                       "\x01\x00\x00\x00"      // 1 point
    3263             :                                       ),
    3264             :                                   22, wkbVariantIso, nBytesConsumed),
    3265             :                 OGRERR_NONE);
    3266           1 :             ASSERT_EQ(oMP.getNumGeometries(), 0);
    3267             :         }
    3268             :     }
    3269             : }
    3270             : 
    3271             : // Test sealing functionality on OGRFieldDefn
    3272           4 : TEST_F(test_ogr, OGRFieldDefn_sealing)
    3273             : {
    3274           2 :     OGRFieldDefn oFieldDefn("test", OFTString);
    3275           1 :     oFieldDefn.Seal();
    3276             : 
    3277             :     {
    3278           1 :         CPLErrorReset();
    3279           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3280           1 :         oFieldDefn.SetName("new_name");
    3281           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3282             :     }
    3283             : 
    3284             :     {
    3285           1 :         CPLErrorReset();
    3286           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3287           1 :         oFieldDefn.SetType(OFTInteger);
    3288           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3289             :     }
    3290             : 
    3291             :     {
    3292           1 :         CPLErrorReset();
    3293           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3294           1 :         oFieldDefn.SetSubType(OFSTJSON);
    3295           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3296             :     }
    3297             : 
    3298             :     {
    3299           1 :         CPLErrorReset();
    3300           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3301           1 :         oFieldDefn.SetWidth(1);
    3302           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3303             :     }
    3304             : 
    3305             :     {
    3306           1 :         CPLErrorReset();
    3307           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3308           1 :         oFieldDefn.SetPrecision(1);
    3309           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3310             :     }
    3311             : 
    3312             :     {
    3313           1 :         CPLErrorReset();
    3314           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3315           1 :         oFieldDefn.SetDefault("");
    3316           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3317             :     }
    3318             : 
    3319             :     {
    3320           1 :         CPLErrorReset();
    3321           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3322           1 :         oFieldDefn.SetUnique(true);
    3323           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3324             :     }
    3325             : 
    3326             :     {
    3327           1 :         CPLErrorReset();
    3328           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3329           1 :         oFieldDefn.SetNullable(false);
    3330           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3331             :     }
    3332             : 
    3333             :     {
    3334           1 :         CPLErrorReset();
    3335           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3336           1 :         oFieldDefn.SetComment("");
    3337           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3338             :     }
    3339             : 
    3340             :     {
    3341           1 :         CPLErrorReset();
    3342           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3343           1 :         oFieldDefn.SetAlternativeName("");
    3344           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3345             :     }
    3346             : 
    3347             :     {
    3348           1 :         CPLErrorReset();
    3349           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3350           1 :         oFieldDefn.SetDomainName("");
    3351           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3352             :     }
    3353             : 
    3354             :     {
    3355           1 :         CPLErrorReset();
    3356           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3357           1 :         oFieldDefn.SetTZFlag(0);
    3358           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3359             :     }
    3360             : 
    3361             :     {
    3362           2 :         auto oTemporaryUnsealer(oFieldDefn.GetTemporaryUnsealer());
    3363           1 :         CPLErrorReset();
    3364           1 :         oFieldDefn.SetName("new_name");
    3365           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") == nullptr);
    3366             :     }
    3367             : 
    3368             :     {
    3369           1 :         CPLErrorReset();
    3370           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3371           1 :         oFieldDefn.SetName("new_name");
    3372           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3373             :     }
    3374             : 
    3375             :     {
    3376           1 :         CPLErrorReset();
    3377           1 :         whileUnsealing(&oFieldDefn)->SetName("new_name");
    3378           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3379             :     }
    3380           1 : }
    3381             : 
    3382             : // Test sealing functionality on OGRGeomFieldDefn
    3383           4 : TEST_F(test_ogr, OGRGeomFieldDefn_sealing)
    3384             : {
    3385           2 :     OGRGeomFieldDefn oFieldDefn("test", wkbUnknown);
    3386             : 
    3387           1 :     oFieldDefn.Seal();
    3388             : 
    3389             :     {
    3390           1 :         CPLErrorReset();
    3391           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3392           1 :         oFieldDefn.SetName("new_name");
    3393           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3394             :     }
    3395             : 
    3396             :     {
    3397           1 :         CPLErrorReset();
    3398           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3399           1 :         oFieldDefn.SetType(wkbPoint);
    3400           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3401             :     }
    3402             : 
    3403             :     {
    3404           1 :         CPLErrorReset();
    3405           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3406           1 :         oFieldDefn.SetSpatialRef(nullptr);
    3407           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3408             :     }
    3409             : 
    3410             :     {
    3411           1 :         CPLErrorReset();
    3412           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3413           1 :         oFieldDefn.SetNullable(false);
    3414           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3415             :     }
    3416             : 
    3417             :     {
    3418           2 :         auto oTemporaryUnsealer(oFieldDefn.GetTemporaryUnsealer());
    3419           1 :         CPLErrorReset();
    3420           1 :         oFieldDefn.SetName("new_name");
    3421           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3422             :     }
    3423             : 
    3424             :     {
    3425           1 :         CPLErrorReset();
    3426           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3427           1 :         oFieldDefn.SetName("new_name");
    3428           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3429             :     }
    3430             : 
    3431             :     {
    3432           1 :         CPLErrorReset();
    3433           1 :         whileUnsealing(&oFieldDefn)->SetName("new_name");
    3434           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3435             :     }
    3436           1 : }
    3437             : 
    3438             : // Test sealing functionality on OGRFeatureDefn
    3439           4 : TEST_F(test_ogr, OGRFeatureDefn_sealing)
    3440             : {
    3441           2 :     OGRFeatureDefn oFDefn;
    3442           1 :     CPLErrorReset();
    3443             :     {
    3444           1 :         oFDefn.SetName("new_name");
    3445           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3446             :     }
    3447             :     {
    3448           2 :         OGRFieldDefn oFieldDefn("test", OFTString);
    3449           1 :         oFDefn.AddFieldDefn(&oFieldDefn);
    3450           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3451             :     }
    3452             :     {
    3453           2 :         OGRGeomFieldDefn oFieldDefn("test", wkbUnknown);
    3454           1 :         oFDefn.AddGeomFieldDefn(&oFieldDefn);
    3455           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3456             :     }
    3457             : 
    3458             :     {
    3459           1 :         CPLErrorReset();
    3460           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3461           1 :         oFDefn.Unseal(true);
    3462           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "unsealed") != nullptr);
    3463             : 
    3464           1 :         CPLErrorReset();
    3465           2 :         auto oTemporaryUnsealer1(oFDefn.GetTemporaryUnsealer(false));
    3466           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "unsealed") != nullptr);
    3467           1 :         CPLErrorReset();
    3468           2 :         auto oTemporaryUnsealer2(oFDefn.GetTemporaryUnsealer(false));
    3469           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3470             :     }
    3471             : 
    3472           1 :     oFDefn.Seal(true);
    3473             : 
    3474             :     {
    3475           1 :         CPLErrorReset();
    3476           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3477           1 :         oFDefn.Seal(true);
    3478           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3479             :     }
    3480             : 
    3481             :     {
    3482           1 :         CPLErrorReset();
    3483           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3484           1 :         oFDefn.SetName("new_name");
    3485           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3486             :     }
    3487             : 
    3488             :     {
    3489           1 :         CPLErrorReset();
    3490           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3491           2 :         OGRFieldDefn oFieldDefn("test2", OFTString);
    3492           1 :         oFDefn.AddFieldDefn(&oFieldDefn);
    3493           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3494             :     }
    3495             : 
    3496             :     {
    3497           1 :         CPLErrorReset();
    3498           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3499           1 :         oFDefn.DeleteFieldDefn(0);
    3500           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3501             :     }
    3502             : 
    3503             :     {
    3504           1 :         CPLErrorReset();
    3505           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3506           1 :         int map[] = {0};
    3507           1 :         oFDefn.ReorderFieldDefns(map);
    3508           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3509             :     }
    3510             : 
    3511             :     {
    3512           1 :         CPLErrorReset();
    3513           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3514           2 :         OGRGeomFieldDefn oFieldDefn("test2", wkbUnknown);
    3515           1 :         oFDefn.AddGeomFieldDefn(&oFieldDefn);
    3516           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3517             :     }
    3518             : 
    3519             :     {
    3520           1 :         CPLErrorReset();
    3521           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3522           1 :         oFDefn.DeleteGeomFieldDefn(0);
    3523           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3524             :     }
    3525             : 
    3526             :     {
    3527           2 :         auto oTemporaryUnsealer(oFDefn.GetTemporaryUnsealer(false));
    3528           1 :         CPLErrorReset();
    3529           1 :         oFDefn.SetName("new_name");
    3530           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3531             : 
    3532           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3533           1 :         CPLErrorReset();
    3534           1 :         oFDefn.GetFieldDefn(0)->SetName("new_name");
    3535           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3536             : 
    3537           1 :         CPLErrorReset();
    3538           1 :         oFDefn.GetGeomFieldDefn(0)->SetName("new_name");
    3539           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3540             :     }
    3541             : 
    3542             :     {
    3543             : 
    3544           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3545           1 :         CPLErrorReset();
    3546           1 :         oFDefn.GetFieldDefn(0)->SetName("new_name");
    3547           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3548             : 
    3549           1 :         CPLErrorReset();
    3550           1 :         oFDefn.GetGeomFieldDefn(0)->SetName("new_name");
    3551           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3552             :     }
    3553             : 
    3554             :     {
    3555           2 :         auto oTemporaryUnsealer(oFDefn.GetTemporaryUnsealer(true));
    3556           1 :         CPLErrorReset();
    3557           1 :         oFDefn.SetName("new_name");
    3558           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3559             : 
    3560           1 :         oFDefn.GetFieldDefn(0)->SetName("new_name");
    3561           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3562             : 
    3563           2 :         auto oTemporaryUnsealer2(oFDefn.GetTemporaryUnsealer(true));
    3564             : 
    3565           1 :         oFDefn.GetGeomFieldDefn(0)->SetName("new_name");
    3566           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3567             :     }
    3568             : 
    3569             :     {
    3570             : 
    3571           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    3572           1 :         CPLErrorReset();
    3573           1 :         oFDefn.GetFieldDefn(0)->SetName("new_name");
    3574           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3575             : 
    3576           1 :         CPLErrorReset();
    3577           1 :         oFDefn.GetGeomFieldDefn(0)->SetName("new_name");
    3578           1 :         EXPECT_TRUE(strstr(CPLGetLastErrorMsg(), "sealed") != nullptr);
    3579             :     }
    3580             : 
    3581             :     {
    3582           1 :         CPLErrorReset();
    3583           1 :         whileUnsealing(&oFDefn)->SetName("new_name");
    3584           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_None);
    3585             :     }
    3586           1 : }
    3587             : 
    3588             : // Test wkbExportOptions
    3589           4 : TEST_F(test_ogr, wkbExportOptions_default)
    3590             : {
    3591           1 :     OGRwkbExportOptions *psOptions = OGRwkbExportOptionsCreate();
    3592           1 :     ASSERT_TRUE(psOptions != nullptr);
    3593           2 :     OGRPoint p(1.23456789012345678, 2.23456789012345678, 3);
    3594           2 :     std::vector<GByte> abyWKB(p.WkbSize());
    3595           1 :     OGR_G_ExportToWkbEx(OGRGeometry::ToHandle(&p), &abyWKB[0], psOptions);
    3596           1 :     OGRwkbExportOptionsDestroy(psOptions);
    3597             : 
    3598           2 :     std::vector<GByte> abyRegularWKB(p.WkbSize());
    3599           1 :     OGR_G_ExportToWkb(OGRGeometry::ToHandle(&p), wkbNDR, &abyRegularWKB[0]);
    3600             : 
    3601           1 :     EXPECT_TRUE(abyWKB == abyRegularWKB);
    3602             : }
    3603             : 
    3604             : // Test wkbExportOptions
    3605           4 : TEST_F(test_ogr, wkbExportOptions)
    3606             : {
    3607           1 :     OGRwkbExportOptions *psOptions = OGRwkbExportOptionsCreate();
    3608           1 :     ASSERT_TRUE(psOptions != nullptr);
    3609           1 :     OGRwkbExportOptionsSetByteOrder(psOptions, wkbXDR);
    3610           1 :     OGRwkbExportOptionsSetVariant(psOptions, wkbVariantIso);
    3611             : 
    3612           1 :     auto hPrec = OGRGeomCoordinatePrecisionCreate();
    3613           1 :     OGRGeomCoordinatePrecisionSet(hPrec, 1e-1, 1e-2, 1e-4);
    3614           1 :     OGRwkbExportOptionsSetPrecision(psOptions, hPrec);
    3615           1 :     OGRGeomCoordinatePrecisionDestroy(hPrec);
    3616             : 
    3617             :     OGRPoint p(1.23456789012345678, -1.23456789012345678, 1.23456789012345678,
    3618           1 :                1.23456789012345678);
    3619           1 :     std::vector<GByte> abyWKB(p.WkbSize());
    3620           1 :     OGR_G_ExportToWkbEx(OGRGeometry::ToHandle(&p), &abyWKB[0], psOptions);
    3621           1 :     OGRwkbExportOptionsDestroy(psOptions);
    3622             : 
    3623             :     const std::vector<GByte> expectedWKB{
    3624             :         0x00, 0x00, 0x00, 0x0B, 0xB9, 0x3F, 0xF3, 0x80, 0x00, 0x00,
    3625             :         0x00, 0x00, 0x00, 0xBF, 0xF3, 0x80, 0x00, 0x00, 0x00, 0x00,
    3626             :         0x00, 0x3F, 0xF3, 0xC0, 0x00, 0x00, 0x00, 0x00, 0x00, 0x3F,
    3627           1 :         0xF3, 0xC0, 0xC0, 0x00, 0x00, 0x00, 0x00};
    3628           1 :     EXPECT_TRUE(abyWKB == expectedWKB);
    3629             : 
    3630           1 :     OGRGeometry *poGeom = nullptr;
    3631           1 :     OGRGeometryFactory::createFromWkb(abyWKB.data(), nullptr, &poGeom);
    3632           1 :     ASSERT_NE(poGeom, nullptr);
    3633           1 :     EXPECT_NEAR(poGeom->toPoint()->getX(), 1.2, 1e-1);
    3634           1 :     EXPECT_NEAR(poGeom->toPoint()->getY(), -1.2, 1e-1);
    3635           1 :     EXPECT_NEAR(poGeom->toPoint()->getZ(), 1.23, 1e-2);
    3636           1 :     EXPECT_NEAR(poGeom->toPoint()->getM(), 1.2346, 1e-4);
    3637           1 :     delete poGeom;
    3638             : }
    3639             : 
    3640             : // Test OGRGeometry::roundCoordinatesIEEE754()
    3641           4 : TEST_F(test_ogr, roundCoordinatesIEEE754)
    3642             : {
    3643           2 :     OGRLineString oLS;
    3644           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234, -1.2345678901234, 0.012345);
    3645           1 :     oLS.addPoint(-1.2345678901234, 1.2345678901234, 1.2345678901234, -0.012345);
    3646           1 :     oLS.addPoint(std::numeric_limits<double>::infinity(),
    3647             :                  std::numeric_limits<double>::quiet_NaN());
    3648           1 :     OGRGeomCoordinateBinaryPrecision sBinaryPrecision;
    3649           2 :     OGRGeomCoordinatePrecision sPrecision;
    3650           1 :     sPrecision.dfXYResolution = 1e-10;
    3651           1 :     sPrecision.dfZResolution = 1e-3;
    3652           1 :     sPrecision.dfMResolution = 1e-5;
    3653           1 :     sBinaryPrecision.SetFrom(sPrecision);
    3654           2 :     OGRLineString oLSOri(oLS);
    3655           1 :     oLS.roundCoordinatesIEEE754(sBinaryPrecision);
    3656           1 :     EXPECT_NE(oLS.getX(0), oLSOri.getX(0));
    3657           1 :     EXPECT_NE(oLS.getY(0), oLSOri.getY(0));
    3658           1 :     EXPECT_NE(oLS.getZ(0), oLSOri.getZ(0));
    3659           1 :     EXPECT_NE(oLS.getM(0), oLSOri.getM(0));
    3660           1 :     EXPECT_NEAR(oLS.getX(0), oLSOri.getX(0), sPrecision.dfXYResolution);
    3661           1 :     EXPECT_NEAR(oLS.getY(0), oLSOri.getY(0), sPrecision.dfXYResolution);
    3662           1 :     EXPECT_NEAR(oLS.getZ(0), oLSOri.getZ(0), sPrecision.dfZResolution);
    3663           1 :     EXPECT_NEAR(oLS.getM(0), oLSOri.getM(0), sPrecision.dfMResolution);
    3664           1 :     EXPECT_NEAR(oLS.getX(1), oLSOri.getX(1), sPrecision.dfXYResolution);
    3665           1 :     EXPECT_NEAR(oLS.getY(1), oLSOri.getY(1), sPrecision.dfXYResolution);
    3666           1 :     EXPECT_NEAR(oLS.getZ(1), oLSOri.getZ(1), sPrecision.dfZResolution);
    3667           1 :     EXPECT_NEAR(oLS.getM(1), oLSOri.getM(1), sPrecision.dfMResolution);
    3668           1 :     EXPECT_EQ(oLS.getX(2), std::numeric_limits<double>::infinity());
    3669           1 :     EXPECT_TRUE(std::isnan(oLS.getY(2)));
    3670           1 : }
    3671             : 
    3672             : // Test discarding of bits in WKB export
    3673           4 : TEST_F(test_ogr, wkb_linestring_2d_xy_precision)
    3674             : {
    3675           2 :     OGRLineString oLS;
    3676           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234);
    3677           1 :     oLS.addPoint(-1.2345678901234, 1.2345678901234);
    3678           1 :     OGRwkbExportOptions sOptions;
    3679           2 :     OGRGeomCoordinatePrecision sPrecision;
    3680           1 :     sPrecision.dfXYResolution = 1e-10;
    3681           1 :     sOptions.sPrecision.SetFrom(sPrecision);
    3682           2 :     std::vector<GByte> abyWKB(oLS.WkbSize());
    3683           1 :     oLS.exportToWkb(&abyWKB[0], &sOptions);
    3684           3 :     for (int i = 0; i < oLS.getDimension() * oLS.getNumPoints(); ++i)
    3685             :     {
    3686           2 :         EXPECT_EQ(abyWKB[5 + 4 + 0 + 8 * i], 0);
    3687             :     }
    3688           1 :     OGRGeometry *poGeom = nullptr;
    3689           1 :     OGRGeometryFactory::createFromWkb(abyWKB.data(), nullptr, &poGeom);
    3690           1 :     EXPECT_NE(poGeom->toLineString()->getX(0), oLS.getX(0));
    3691           1 :     EXPECT_NE(poGeom->toLineString()->getY(0), oLS.getY(0));
    3692           1 :     EXPECT_NEAR(poGeom->toLineString()->getX(0), oLS.getX(0),
    3693             :                 sPrecision.dfXYResolution);
    3694           1 :     EXPECT_NEAR(poGeom->toLineString()->getY(0), oLS.getY(0),
    3695             :                 sPrecision.dfXYResolution);
    3696           1 :     EXPECT_NEAR(poGeom->toLineString()->getX(1), oLS.getX(1),
    3697             :                 sPrecision.dfXYResolution);
    3698           1 :     EXPECT_NEAR(poGeom->toLineString()->getY(1), oLS.getY(1),
    3699             :                 sPrecision.dfXYResolution);
    3700           1 :     delete poGeom;
    3701           1 : }
    3702             : 
    3703             : // Test discarding of bits in WKB export
    3704           4 : TEST_F(test_ogr, wkb_linestring_3d_discard_lsb_bits)
    3705             : {
    3706           2 :     OGRLineString oLS;
    3707           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234, -1.2345678901234);
    3708           1 :     oLS.addPoint(-1.2345678901234, 1.2345678901234, 1.2345678901234);
    3709           1 :     OGRwkbExportOptions sOptions;
    3710           2 :     OGRGeomCoordinatePrecision sPrecision;
    3711           1 :     sPrecision.dfXYResolution = 1e-10;
    3712           1 :     sPrecision.dfZResolution = 1e-3;
    3713           1 :     sOptions.sPrecision.SetFrom(sPrecision);
    3714           2 :     std::vector<GByte> abyWKB(oLS.WkbSize());
    3715           1 :     oLS.exportToWkb(&abyWKB[0], &sOptions);
    3716           3 :     for (int i = 0; i < oLS.getDimension() * oLS.getNumPoints(); ++i)
    3717             :     {
    3718           2 :         EXPECT_EQ(abyWKB[5 + 4 + 0 + 8 * i], 0);
    3719             :     }
    3720           1 :     OGRGeometry *poGeom = nullptr;
    3721           1 :     OGRGeometryFactory::createFromWkb(abyWKB.data(), nullptr, &poGeom);
    3722           1 :     EXPECT_NE(poGeom->toLineString()->getX(0), oLS.getX(0));
    3723           1 :     EXPECT_NE(poGeom->toLineString()->getY(0), oLS.getY(0));
    3724           1 :     EXPECT_NE(poGeom->toLineString()->getZ(0), oLS.getZ(0));
    3725           1 :     EXPECT_NEAR(poGeom->toLineString()->getX(0), oLS.getX(0),
    3726             :                 sPrecision.dfXYResolution);
    3727           1 :     EXPECT_NEAR(poGeom->toLineString()->getY(0), oLS.getY(0),
    3728             :                 sPrecision.dfXYResolution);
    3729           1 :     EXPECT_NEAR(poGeom->toLineString()->getZ(0), oLS.getZ(0),
    3730             :                 sPrecision.dfZResolution);
    3731           1 :     EXPECT_NEAR(poGeom->toLineString()->getX(1), oLS.getX(1),
    3732             :                 sPrecision.dfXYResolution);
    3733           1 :     EXPECT_NEAR(poGeom->toLineString()->getY(1), oLS.getY(1),
    3734             :                 sPrecision.dfXYResolution);
    3735           1 :     EXPECT_NEAR(poGeom->toLineString()->getZ(1), oLS.getZ(1),
    3736             :                 sPrecision.dfZResolution);
    3737           1 :     delete poGeom;
    3738           1 : }
    3739             : 
    3740             : // Test discarding of bits in WKB export
    3741           4 : TEST_F(test_ogr, wkb_linestring_xym_discard_lsb_bits)
    3742             : {
    3743           2 :     OGRLineString oLS;
    3744           1 :     oLS.addPointM(1.2345678901234, -1.2345678901234, -1.2345678901234);
    3745           1 :     oLS.addPointM(-1.2345678901234, 1.2345678901234, 1.2345678901234);
    3746           1 :     OGRwkbExportOptions sOptions;
    3747           2 :     OGRGeomCoordinatePrecision sPrecision;
    3748           1 :     sPrecision.dfXYResolution = 1e-10;
    3749           1 :     sPrecision.dfMResolution = 1e-3;
    3750           1 :     sOptions.sPrecision.SetFrom(sPrecision);
    3751           2 :     std::vector<GByte> abyWKB(oLS.WkbSize());
    3752           1 :     oLS.exportToWkb(&abyWKB[0], &sOptions);
    3753           3 :     for (int i = 0; i < oLS.getDimension() * oLS.getNumPoints(); ++i)
    3754             :     {
    3755           2 :         EXPECT_EQ(abyWKB[5 + 4 + 0 + 8 * i], 0);
    3756             :     }
    3757           1 :     OGRGeometry *poGeom = nullptr;
    3758           1 :     OGRGeometryFactory::createFromWkb(abyWKB.data(), nullptr, &poGeom);
    3759           1 :     EXPECT_NE(poGeom->toLineString()->getX(0), oLS.getX(0));
    3760           1 :     EXPECT_NE(poGeom->toLineString()->getY(0), oLS.getY(0));
    3761           1 :     EXPECT_NE(poGeom->toLineString()->getM(0), oLS.getM(0));
    3762           1 :     EXPECT_NEAR(poGeom->toLineString()->getX(0), oLS.getX(0),
    3763             :                 sPrecision.dfXYResolution);
    3764           1 :     EXPECT_NEAR(poGeom->toLineString()->getY(0), oLS.getY(0),
    3765             :                 sPrecision.dfXYResolution);
    3766           1 :     EXPECT_NEAR(poGeom->toLineString()->getM(0), oLS.getM(0),
    3767             :                 sPrecision.dfMResolution);
    3768           1 :     EXPECT_NEAR(poGeom->toLineString()->getX(1), oLS.getX(1),
    3769             :                 sPrecision.dfXYResolution);
    3770           1 :     EXPECT_NEAR(poGeom->toLineString()->getY(1), oLS.getY(1),
    3771             :                 sPrecision.dfXYResolution);
    3772           1 :     EXPECT_NEAR(poGeom->toLineString()->getM(1), oLS.getM(1),
    3773             :                 sPrecision.dfMResolution);
    3774           1 :     delete poGeom;
    3775           1 : }
    3776             : 
    3777             : // Test discarding of bits in WKB export
    3778           4 : TEST_F(test_ogr, wkb_linestring_xyzm_discard_lsb_bits)
    3779             : {
    3780           2 :     OGRLineString oLS;
    3781           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234, -1.2345678901234, 0.012345);
    3782           1 :     oLS.addPoint(-1.2345678901234, 1.2345678901234, 1.2345678901234, 0.012345);
    3783           1 :     OGRwkbExportOptions sOptions;
    3784           2 :     OGRGeomCoordinatePrecision sPrecision;
    3785           1 :     sPrecision.dfXYResolution = 1e-10;
    3786           1 :     sPrecision.dfZResolution = 1e-3;
    3787           1 :     sPrecision.dfMResolution = 1e-5;
    3788           1 :     sOptions.sPrecision.SetFrom(sPrecision);
    3789           2 :     std::vector<GByte> abyWKB(oLS.WkbSize());
    3790           1 :     oLS.exportToWkb(&abyWKB[0], &sOptions);
    3791           3 :     for (int i = 0; i < oLS.getDimension() * oLS.getNumPoints(); ++i)
    3792             :     {
    3793           2 :         EXPECT_EQ(abyWKB[5 + 4 + 0 + 8 * i], 0);
    3794             :     }
    3795           1 :     OGRGeometry *poGeom = nullptr;
    3796           1 :     OGRGeometryFactory::createFromWkb(abyWKB.data(), nullptr, &poGeom);
    3797           1 :     EXPECT_NE(poGeom->toLineString()->getX(0), oLS.getX(0));
    3798           1 :     EXPECT_NE(poGeom->toLineString()->getY(0), oLS.getY(0));
    3799           1 :     EXPECT_NE(poGeom->toLineString()->getZ(0), oLS.getZ(0));
    3800           1 :     EXPECT_NE(poGeom->toLineString()->getM(0), oLS.getM(0));
    3801           1 :     EXPECT_NEAR(poGeom->toLineString()->getX(0), oLS.getX(0),
    3802             :                 sPrecision.dfXYResolution);
    3803           1 :     EXPECT_NEAR(poGeom->toLineString()->getY(0), oLS.getY(0),
    3804             :                 sPrecision.dfXYResolution);
    3805           1 :     EXPECT_NEAR(poGeom->toLineString()->getZ(0), oLS.getZ(0),
    3806             :                 sPrecision.dfZResolution);
    3807           1 :     EXPECT_NEAR(poGeom->toLineString()->getM(0), oLS.getM(0),
    3808             :                 sPrecision.dfMResolution);
    3809           1 :     EXPECT_NEAR(poGeom->toLineString()->getX(1), oLS.getX(1),
    3810             :                 sPrecision.dfXYResolution);
    3811           1 :     EXPECT_NEAR(poGeom->toLineString()->getY(1), oLS.getY(1),
    3812             :                 sPrecision.dfXYResolution);
    3813           1 :     EXPECT_NEAR(poGeom->toLineString()->getZ(1), oLS.getZ(1),
    3814             :                 sPrecision.dfZResolution);
    3815           1 :     EXPECT_NEAR(poGeom->toLineString()->getM(1), oLS.getM(1),
    3816             :                 sPrecision.dfMResolution);
    3817           1 :     delete poGeom;
    3818           1 : }
    3819             : 
    3820             : // Test discarding of bits in WKB export
    3821           4 : TEST_F(test_ogr, wkb_polygon_2d_xy_precision)
    3822             : {
    3823           2 :     OGRLinearRing oLS;
    3824           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234);
    3825           1 :     oLS.addPoint(-1.2345678901234, -1.2345678901234);
    3826           1 :     oLS.addPoint(-2.2345678901234, 1.2345678901234);
    3827           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234);
    3828           2 :     OGRPolygon oPoly;
    3829           1 :     oPoly.addRing(&oLS);
    3830           1 :     OGRwkbExportOptions sOptions;
    3831           2 :     OGRGeomCoordinatePrecision sPrecision;
    3832           1 :     sPrecision.dfXYResolution = 1e-10;
    3833           1 :     sOptions.sPrecision.SetFrom(sPrecision);
    3834           2 :     std::vector<GByte> abyWKB(oPoly.WkbSize());
    3835           1 :     oPoly.exportToWkb(&abyWKB[0], &sOptions);
    3836           5 :     for (int i = 0; i < oLS.getDimension() * oLS.getNumPoints(); ++i)
    3837             :     {
    3838           4 :         EXPECT_EQ(abyWKB[5 + 4 + 4 + 0 + 8 * i], 0);
    3839             :     }
    3840           1 :     OGRGeometry *poGeom = nullptr;
    3841           1 :     OGRGeometryFactory::createFromWkb(abyWKB.data(), nullptr, &poGeom);
    3842           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getX(0), oLS.getX(0));
    3843           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getY(0), oLS.getY(0));
    3844           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getX(0), oLS.getX(0),
    3845             :                 sPrecision.dfXYResolution);
    3846           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getY(0), oLS.getY(0),
    3847             :                 sPrecision.dfXYResolution);
    3848           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getX(1), oLS.getX(1),
    3849             :                 sPrecision.dfXYResolution);
    3850           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getY(1), oLS.getY(1),
    3851             :                 sPrecision.dfXYResolution);
    3852           1 :     delete poGeom;
    3853           1 : }
    3854             : 
    3855             : // Test discarding of bits in WKB export
    3856           4 : TEST_F(test_ogr, wkb_polygon_3d_discard_lsb_bits)
    3857             : {
    3858           2 :     OGRLinearRing oLS;
    3859           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234, 1.2345678901234);
    3860           1 :     oLS.addPoint(-1.2345678901234, -1.2345678901234, -1.2345678901234);
    3861           1 :     oLS.addPoint(-2.2345678901234, 1.2345678901234, -1.2345678901234);
    3862           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234, 1.2345678901234);
    3863           2 :     OGRPolygon oPoly;
    3864           1 :     oPoly.addRing(&oLS);
    3865           2 :     OGRSpatialReference oSRS;
    3866           1 :     oSRS.importFromEPSG(4326);
    3867           1 :     OGRwkbExportOptions sOptions;
    3868           2 :     OGRGeomCoordinatePrecision sPrecision;
    3869           1 :     sPrecision.SetFromMeter(&oSRS, 1e-3, 1e-3, 0);
    3870           1 :     sOptions.sPrecision.SetFrom(sPrecision);
    3871           2 :     std::vector<GByte> abyWKB(oPoly.WkbSize());
    3872           1 :     oPoly.exportToWkb(&abyWKB[0], &sOptions);
    3873           5 :     for (int i = 0; i < oLS.getDimension() * oLS.getNumPoints(); ++i)
    3874             :     {
    3875           4 :         EXPECT_EQ(abyWKB[5 + 4 + 4 + 0 + 8 * i], 0);
    3876             :     }
    3877           1 :     OGRGeometry *poGeom = nullptr;
    3878           1 :     OGRGeometryFactory::createFromWkb(abyWKB.data(), nullptr, &poGeom);
    3879           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getX(0), oLS.getX(0));
    3880           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getY(0), oLS.getY(0));
    3881           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getZ(0), oLS.getZ(0));
    3882           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getX(0), oLS.getX(0),
    3883             :                 8.9e-9);
    3884           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getY(0), oLS.getY(0),
    3885             :                 8.9e-9);
    3886           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getZ(0), oLS.getZ(0),
    3887             :                 1e-3);
    3888           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getX(1), oLS.getX(1),
    3889             :                 8.9e-9);
    3890           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getY(1), oLS.getY(1),
    3891             :                 8.9e-9);
    3892           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getZ(1), oLS.getZ(1),
    3893             :                 1e-3);
    3894           1 :     delete poGeom;
    3895           1 : }
    3896             : 
    3897             : // Test discarding of bits in WKB export
    3898           4 : TEST_F(test_ogr, wkb_polygon_xym_discard_lsb_bits)
    3899             : {
    3900           2 :     OGRLinearRing oLS;
    3901           1 :     oLS.addPointM(1.2345678901234, -1.2345678901234, 1.2345678901234);
    3902           1 :     oLS.addPointM(-1.2345678901234, -1.2345678901234, -1.2345678901234);
    3903           1 :     oLS.addPointM(-2.2345678901234, 1.2345678901234, -1.2345678901234);
    3904           1 :     oLS.addPointM(1.2345678901234, -1.2345678901234, 1.2345678901234);
    3905           2 :     OGRPolygon oPoly;
    3906           1 :     oPoly.addRing(&oLS);
    3907           2 :     OGRSpatialReference oSRS;
    3908           1 :     oSRS.importFromEPSG(4326);
    3909           1 :     OGRwkbExportOptions sOptions;
    3910           2 :     OGRGeomCoordinatePrecision sPrecision;
    3911           1 :     sPrecision.SetFromMeter(&oSRS, 1e-3, 0, 1e-3);
    3912           1 :     sOptions.sPrecision.SetFrom(sPrecision);
    3913           2 :     std::vector<GByte> abyWKB(oPoly.WkbSize());
    3914           1 :     oPoly.exportToWkb(&abyWKB[0], &sOptions);
    3915           5 :     for (int i = 0; i < oLS.getDimension() * oLS.getNumPoints(); ++i)
    3916             :     {
    3917           4 :         EXPECT_EQ(abyWKB[5 + 4 + 4 + 0 + 8 * i], 0);
    3918             :     }
    3919           1 :     OGRGeometry *poGeom = nullptr;
    3920           1 :     OGRGeometryFactory::createFromWkb(abyWKB.data(), nullptr, &poGeom);
    3921           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getX(0), oLS.getX(0));
    3922           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getY(0), oLS.getY(0));
    3923           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getM(0), oLS.getM(0));
    3924           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getX(0), oLS.getX(0),
    3925             :                 8.9e-9);
    3926           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getY(0), oLS.getY(0),
    3927             :                 8.9e-9);
    3928           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getM(0), oLS.getM(0),
    3929             :                 1e-3);
    3930           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getX(1), oLS.getX(1),
    3931             :                 8.9e-9);
    3932           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getY(1), oLS.getY(1),
    3933             :                 8.9e-9);
    3934           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getM(1), oLS.getM(1),
    3935             :                 1e-3);
    3936           1 :     delete poGeom;
    3937           1 : }
    3938             : 
    3939             : // Test discarding of bits in WKB export
    3940           4 : TEST_F(test_ogr, wkb_polygon_xyzm_discard_lsb_bits)
    3941             : {
    3942           2 :     OGRLinearRing oLS;
    3943           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234, 1.2345678901234, 0.012345);
    3944           1 :     oLS.addPoint(-1.2345678901234, -1.2345678901234, -1.2345678901234, 12345);
    3945           1 :     oLS.addPoint(-2.2345678901234, 1.2345678901234, -1.2345678901234, 0.012345);
    3946           1 :     oLS.addPoint(1.2345678901234, -1.2345678901234, 1.2345678901234, 0.012345);
    3947           2 :     OGRPolygon oPoly;
    3948           1 :     oPoly.addRing(&oLS);
    3949           2 :     OGRSpatialReference oSRS;
    3950           1 :     oSRS.importFromEPSG(4326);
    3951           1 :     OGRwkbExportOptions sOptions;
    3952           2 :     OGRGeomCoordinatePrecision sPrecision;
    3953           1 :     sPrecision.SetFromMeter(&oSRS, 1e-3, 1e-3, 1e-4);
    3954           1 :     sOptions.sPrecision.SetFrom(sPrecision);
    3955           2 :     std::vector<GByte> abyWKB(oPoly.WkbSize());
    3956           1 :     oPoly.exportToWkb(&abyWKB[0], &sOptions);
    3957           5 :     for (int i = 0; i < oLS.getDimension() * oLS.getNumPoints(); ++i)
    3958             :     {
    3959           4 :         EXPECT_EQ(abyWKB[5 + 4 + 4 + 0 + 8 * i], 0);
    3960             :     }
    3961           1 :     OGRGeometry *poGeom = nullptr;
    3962           1 :     OGRGeometryFactory::createFromWkb(abyWKB.data(), nullptr, &poGeom);
    3963           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getX(0), oLS.getX(0));
    3964           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getY(0), oLS.getY(0));
    3965           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getZ(0), oLS.getZ(0));
    3966           1 :     EXPECT_NE(poGeom->toPolygon()->getExteriorRing()->getM(0), oLS.getM(0));
    3967           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getX(0), oLS.getX(0),
    3968             :                 8.9e-9);
    3969           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getY(0), oLS.getY(0),
    3970             :                 8.9e-9);
    3971           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getZ(0), oLS.getZ(0),
    3972             :                 1e-3);
    3973           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getM(0), oLS.getM(0),
    3974             :                 1e-4);
    3975           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getX(1), oLS.getX(1),
    3976             :                 8.9e-9);
    3977           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getY(1), oLS.getY(1),
    3978             :                 8.9e-9);
    3979           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getZ(1), oLS.getZ(1),
    3980             :                 1e-3);
    3981           1 :     EXPECT_NEAR(poGeom->toPolygon()->getExteriorRing()->getM(1), oLS.getM(1),
    3982             :                 1e-4);
    3983           1 :     delete poGeom;
    3984           1 : }
    3985             : 
    3986             : // Test OGRFeature::SerializeToBinary() and DeserializeFromBinary();
    3987           4 : TEST_F(test_ogr, OGRFeature_SerializeToBinary)
    3988             : {
    3989             :     {
    3990           2 :         OGRFeatureDefn oFDefn;
    3991           1 :         oFDefn.SetGeomType(wkbNone);
    3992           1 :         oFDefn.Reference();
    3993             : 
    3994             :         {
    3995           2 :             OGRFeature oFeatSrc(&oFDefn);
    3996           1 :             oFeatSrc.SetFID(1);
    3997           2 :             std::vector<GByte> abyBuffer;
    3998             : 
    3999           1 :             EXPECT_TRUE(oFeatSrc.SerializeToBinary(abyBuffer));
    4000           1 :             EXPECT_EQ(abyBuffer.size(), 1);
    4001           1 :             EXPECT_EQ(abyBuffer[0], 1);
    4002             : 
    4003           2 :             OGRFeature oFeatDst(&oFDefn);
    4004           1 :             EXPECT_FALSE(oFeatDst.DeserializeFromBinary(abyBuffer.data(), 0));
    4005           1 :             EXPECT_TRUE(oFeatDst.DeserializeFromBinary(abyBuffer.data(),
    4006             :                                                        abyBuffer.size()));
    4007           1 :             EXPECT_EQ(oFeatDst.GetFID(), 1);
    4008             :         }
    4009             : 
    4010             :         {
    4011           2 :             OGRFeature oFeatSrc(&oFDefn);
    4012           1 :             oFeatSrc.SetFID(static_cast<GIntBig>(-12345678901234));
    4013           2 :             std::vector<GByte> abyBuffer;
    4014             : 
    4015           1 :             EXPECT_TRUE(oFeatSrc.SerializeToBinary(abyBuffer));
    4016             : 
    4017           2 :             OGRFeature oFeatDst(&oFDefn);
    4018             :             // Try truncated buffers
    4019           8 :             for (size_t i = 0; i < abyBuffer.size(); ++i)
    4020             :             {
    4021           7 :                 EXPECT_FALSE(
    4022             :                     oFeatDst.DeserializeFromBinary(abyBuffer.data(), i));
    4023             :             }
    4024           1 :             EXPECT_TRUE(oFeatDst.DeserializeFromBinary(abyBuffer.data(),
    4025             :                                                        abyBuffer.size()));
    4026           1 :             EXPECT_EQ(oFeatDst.GetFID(), static_cast<GIntBig>(-12345678901234));
    4027             :         }
    4028             :     }
    4029             : 
    4030             :     {
    4031           1 :         OGRFeatureDefn oFDefn;
    4032           1 :         oFDefn.Reference();
    4033             :         {
    4034           2 :             OGRFieldDefn oFieldDefn("int", OFTInteger);
    4035           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4036             :         }
    4037             :         {
    4038           2 :             OGRFieldDefn oFieldDefn("int64", OFTInteger64);
    4039           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4040             :         }
    4041             :         {
    4042           2 :             OGRFieldDefn oFieldDefn("real", OFTReal);
    4043           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4044             :         }
    4045             :         {
    4046           2 :             OGRFieldDefn oFieldDefn("str", OFTString);
    4047           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4048             :         }
    4049             :         {
    4050           2 :             OGRFieldDefn oFieldDefn("binary", OFTBinary);
    4051           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4052             :         }
    4053             :         {
    4054           2 :             OGRFieldDefn oFieldDefn("intlist", OFTIntegerList);
    4055           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4056             :         }
    4057             :         {
    4058           2 :             OGRFieldDefn oFieldDefn("int64list", OFTInteger64List);
    4059           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4060             :         }
    4061             :         {
    4062           2 :             OGRFieldDefn oFieldDefn("reallist", OFTRealList);
    4063           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4064             :         }
    4065             :         {
    4066           2 :             OGRFieldDefn oFieldDefn("strlist", OFTStringList);
    4067           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4068             :         }
    4069             :         {
    4070           2 :             OGRFieldDefn oFieldDefn("date", OFTDate);
    4071           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4072             :         }
    4073             :         {
    4074           2 :             OGRFieldDefn oFieldDefn("time", OFTTime);
    4075           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4076             :         }
    4077             :         {
    4078           2 :             OGRFieldDefn oFieldDefn("datetime", OFTDateTime);
    4079           1 :             oFDefn.AddFieldDefn(&oFieldDefn);
    4080             :         }
    4081             : 
    4082             :         {
    4083           2 :             OGRFeature oFeatSrc(&oFDefn);
    4084           2 :             std::vector<GByte> abyBuffer;
    4085             : 
    4086           1 :             EXPECT_TRUE(oFeatSrc.SerializeToBinary(abyBuffer));
    4087           1 :             EXPECT_EQ(abyBuffer.size(), 5);
    4088             : 
    4089           2 :             OGRFeature oFeatDst(&oFDefn);
    4090           6 :             for (size_t i = 0; i < abyBuffer.size(); ++i)
    4091             :             {
    4092           5 :                 EXPECT_FALSE(
    4093             :                     oFeatDst.DeserializeFromBinary(abyBuffer.data(), i));
    4094             :             }
    4095           1 :             EXPECT_TRUE(oFeatDst.DeserializeFromBinary(abyBuffer.data(),
    4096             :                                                        abyBuffer.size()));
    4097           1 :             EXPECT_TRUE(oFeatDst.Equal(&oFeatSrc));
    4098             :         }
    4099             : 
    4100             :         {
    4101           1 :             OGRFeature oFeatSrc(&oFDefn);
    4102           1 :             std::vector<GByte> abyBuffer;
    4103             : 
    4104           1 :             const int iFieldInt = oFDefn.GetFieldIndex("int");
    4105           1 :             ASSERT_TRUE(iFieldInt >= 0);
    4106           1 :             oFeatSrc.SetFieldNull(iFieldInt);
    4107           1 :             EXPECT_TRUE(oFeatSrc.SerializeToBinary(abyBuffer));
    4108           1 :             EXPECT_EQ(abyBuffer.size(), 5);
    4109             : 
    4110           2 :             OGRFeature oFeatDst(&oFDefn);
    4111             : 
    4112             :             // Try truncated buffers
    4113           6 :             for (size_t i = 0; i < abyBuffer.size(); ++i)
    4114             :             {
    4115           5 :                 EXPECT_FALSE(
    4116             :                     oFeatDst.DeserializeFromBinary(abyBuffer.data(), i));
    4117             :             }
    4118             : 
    4119           1 :             EXPECT_TRUE(oFeatDst.DeserializeFromBinary(abyBuffer.data(),
    4120             :                                                        abyBuffer.size()));
    4121           1 :             EXPECT_TRUE(oFeatDst.Equal(&oFeatSrc));
    4122             :         }
    4123             : 
    4124             :         {
    4125           1 :             OGRFeature oFeatSrc(&oFDefn);
    4126           1 :             oFeatSrc.SetFID(1);
    4127           1 :             oFeatSrc.SetField("int", -123);
    4128           1 :             oFeatSrc.SetField("int64", static_cast<GIntBig>(-12345678901234));
    4129           1 :             oFeatSrc.SetField("real", 1.25);
    4130           1 :             oFeatSrc.SetField("str", "foo");
    4131           1 :             const int iFieldBinary = oFDefn.GetFieldIndex("binary");
    4132           1 :             ASSERT_TRUE(iFieldBinary >= 0);
    4133           1 :             oFeatSrc.SetField(iFieldBinary, 3,
    4134             :                               static_cast<const void *>("abc"));
    4135           1 :             oFeatSrc.SetField("intlist", 2,
    4136           2 :                               std::vector<int>{1, -123456}.data());
    4137           1 :             oFeatSrc.SetField("int64list", 2,
    4138           2 :                               std::vector<GIntBig>{1, -12345678901234}.data());
    4139           1 :             oFeatSrc.SetField("reallist", 2,
    4140           2 :                               std::vector<double>{1.5, -2.5}.data());
    4141           2 :             CPLStringList aosList;
    4142           1 :             aosList.AddString("foo");
    4143           1 :             aosList.AddString("barbaz");
    4144           1 :             oFeatSrc.SetField("strlist", aosList.List());
    4145           1 :             oFeatSrc.SetField("date", 2023, 1, 3);
    4146           1 :             oFeatSrc.SetField("time", 0, 0, 0, 12, 34, 56.789f);
    4147           1 :             oFeatSrc.SetField("datetime", 2023, 1, 3, 12, 34, 56.789f);
    4148           2 :             OGRPoint p(1, 2);
    4149           1 :             oFeatSrc.SetGeometry(&p);
    4150           2 :             std::vector<GByte> abyBuffer;
    4151             : 
    4152           1 :             EXPECT_TRUE(oFeatSrc.SerializeToBinary(abyBuffer));
    4153             : 
    4154           2 :             OGRFeature oFeatDst(&oFDefn);
    4155             : 
    4156             :             // Try truncated buffers
    4157         118 :             for (size_t i = 0; i < abyBuffer.size(); ++i)
    4158             :             {
    4159         117 :                 EXPECT_FALSE(
    4160             :                     oFeatDst.DeserializeFromBinary(abyBuffer.data(), i));
    4161             :             }
    4162             : 
    4163             :             // Try corrupted buffers
    4164             :             {
    4165           2 :                 CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    4166         118 :                 for (size_t i = 0; i < abyBuffer.size(); ++i)
    4167             :                 {
    4168             :                     // Might succeed or fail, but shouldn't crash..
    4169         117 :                     const GByte backup = abyBuffer[i];
    4170         117 :                     abyBuffer[i] = static_cast<GByte>(~abyBuffer[i]);
    4171         117 :                     (void)oFeatDst.DeserializeFromBinary(abyBuffer.data(),
    4172             :                                                          abyBuffer.size());
    4173         117 :                     abyBuffer[i] = backup;
    4174             :                 }
    4175             :             }
    4176             : 
    4177           1 :             EXPECT_TRUE(oFeatDst.DeserializeFromBinary(abyBuffer.data(),
    4178             :                                                        abyBuffer.size()));
    4179             :             // oFeatSrc.DumpReadable(stdout);
    4180             :             // oFeatDst.DumpReadable(stdout);
    4181           1 :             EXPECT_TRUE(oFeatDst.Equal(&oFeatSrc));
    4182             :         }
    4183             :     }
    4184             : }
    4185             : 
    4186             : // Test OGRGeometry::IsRectangle()
    4187           4 : TEST_F(test_ogr, OGRGeometry_IsRectangle)
    4188             : {
    4189             :     // Not a polygon
    4190             :     {
    4191           1 :         OGRGeometry *poGeom = nullptr;
    4192           1 :         OGRGeometryFactory::createFromWkt("POINT EMPTY", nullptr, &poGeom);
    4193           1 :         ASSERT_NE(poGeom, nullptr);
    4194           1 :         EXPECT_FALSE(poGeom->IsRectangle());
    4195           1 :         delete poGeom;
    4196             :     }
    4197             :     // Polygon empty
    4198             :     {
    4199           1 :         OGRGeometry *poGeom = nullptr;
    4200           1 :         OGRGeometryFactory::createFromWkt("POLYGON EMPTY", nullptr, &poGeom);
    4201           1 :         ASSERT_NE(poGeom, nullptr);
    4202           1 :         EXPECT_FALSE(poGeom->IsRectangle());
    4203           1 :         delete poGeom;
    4204             :     }
    4205             :     // Polygon with inner ring
    4206             :     {
    4207           1 :         OGRGeometry *poGeom = nullptr;
    4208           1 :         OGRGeometryFactory::createFromWkt(
    4209             :             "POLYGON ((0 0,0 1,1 1,1 0,0 0),(0.2 0.2,0.2 0.8,0.8 0.8,0.8 "
    4210             :             "0.2,0.2 0.2))",
    4211             :             nullptr, &poGeom);
    4212           1 :         ASSERT_NE(poGeom, nullptr);
    4213           1 :         EXPECT_FALSE(poGeom->IsRectangle());
    4214           1 :         delete poGeom;
    4215             :     }
    4216             :     // Polygon with 3 points
    4217             :     {
    4218           1 :         OGRGeometry *poGeom = nullptr;
    4219           1 :         OGRGeometryFactory::createFromWkt("POLYGON ((0 0,0 1,1 1))", nullptr,
    4220             :                                           &poGeom);
    4221           1 :         ASSERT_NE(poGeom, nullptr);
    4222           1 :         EXPECT_FALSE(poGeom->IsRectangle());
    4223           1 :         delete poGeom;
    4224             :     }
    4225             :     // Polygon with 6 points
    4226             :     {
    4227           1 :         OGRGeometry *poGeom = nullptr;
    4228           1 :         OGRGeometryFactory::createFromWkt(
    4229             :             "POLYGON ((0 0,0.1 0,0.2 0,0.3 0,1 1,0 0))", nullptr, &poGeom);
    4230           1 :         ASSERT_NE(poGeom, nullptr);
    4231           1 :         EXPECT_FALSE(poGeom->IsRectangle());
    4232           1 :         delete poGeom;
    4233             :     }
    4234             :     // Polygon with 5 points, but last one not matching first (invalid)
    4235             :     {
    4236           1 :         OGRGeometry *poGeom = nullptr;
    4237           1 :         OGRGeometryFactory::createFromWkt(
    4238             :             "POLYGON ((0 0,0 1,1 1,1 0,-999 -999))", nullptr, &poGeom);
    4239           1 :         ASSERT_NE(poGeom, nullptr);
    4240           1 :         EXPECT_FALSE(poGeom->IsRectangle());
    4241           1 :         delete poGeom;
    4242             :     }
    4243             :     // Polygon with 5 points, but not rectangle
    4244             :     {
    4245           1 :         OGRGeometry *poGeom = nullptr;
    4246           1 :         OGRGeometryFactory::createFromWkt("POLYGON ((0 0,0 1.1,1 1,1 0,0 0))",
    4247             :                                           nullptr, &poGeom);
    4248           1 :         ASSERT_NE(poGeom, nullptr);
    4249           1 :         EXPECT_FALSE(poGeom->IsRectangle());
    4250           1 :         delete poGeom;
    4251             :     }
    4252             :     // Rectangle (type 1)
    4253             :     {
    4254           1 :         OGRGeometry *poGeom = nullptr;
    4255           1 :         OGRGeometryFactory::createFromWkt("POLYGON ((0 0,0 1,1 1,1 0,0 0))",
    4256             :                                           nullptr, &poGeom);
    4257           1 :         ASSERT_NE(poGeom, nullptr);
    4258           1 :         EXPECT_TRUE(poGeom->IsRectangle());
    4259           1 :         delete poGeom;
    4260             :     }
    4261             :     // Rectangle2(type 1)
    4262             :     {
    4263           1 :         OGRGeometry *poGeom = nullptr;
    4264           1 :         OGRGeometryFactory::createFromWkt("POLYGON ((0 0,1 0,1 1,0 1,0 0))",
    4265             :                                           nullptr, &poGeom);
    4266           1 :         ASSERT_NE(poGeom, nullptr);
    4267           1 :         EXPECT_TRUE(poGeom->IsRectangle());
    4268           1 :         delete poGeom;
    4269             :     }
    4270             : }
    4271             : 
    4272             : // Test OGRGeometry::removeEmptyParts()
    4273           4 : TEST_F(test_ogr, OGRGeometry_removeEmptyParts)
    4274             : {
    4275             :     {
    4276           1 :         OGRGeometry *poGeom = nullptr;
    4277           1 :         OGRGeometryFactory::createFromWkt("POINT EMPTY", nullptr, &poGeom);
    4278           1 :         ASSERT_NE(poGeom, nullptr);
    4279           1 :         EXPECT_FALSE(poGeom->hasEmptyParts());
    4280           1 :         poGeom->removeEmptyParts();
    4281           1 :         EXPECT_TRUE(poGeom->IsEmpty());
    4282           1 :         delete poGeom;
    4283             :     }
    4284             :     {
    4285           1 :         OGRGeometry *poGeom = nullptr;
    4286           1 :         OGRGeometryFactory::createFromWkt("POLYGON ((0 0,0 1,1 0,0 0))",
    4287             :                                           nullptr, &poGeom);
    4288           1 :         ASSERT_NE(poGeom, nullptr);
    4289           1 :         EXPECT_FALSE(poGeom->hasEmptyParts());
    4290           1 :         poGeom->removeEmptyParts();
    4291           1 :         EXPECT_NE(poGeom->toPolygon()->getExteriorRing(), nullptr);
    4292           1 :         delete poGeom;
    4293             :     }
    4294             :     {
    4295           1 :         OGRGeometry *poGeom = nullptr;
    4296           1 :         OGRGeometryFactory::createFromWkt("POLYGON ((0 0,0 1,1 0,0 0))",
    4297             :                                           nullptr, &poGeom);
    4298           1 :         ASSERT_NE(poGeom, nullptr);
    4299           1 :         poGeom->toPolygon()->addRingDirectly(new OGRLinearRing());
    4300           1 :         EXPECT_EQ(poGeom->toPolygon()->getNumInteriorRings(), 1);
    4301           1 :         EXPECT_TRUE(poGeom->hasEmptyParts());
    4302           1 :         poGeom->removeEmptyParts();
    4303           1 :         EXPECT_NE(poGeom->toPolygon()->getExteriorRing(), nullptr);
    4304           1 :         EXPECT_EQ(poGeom->toPolygon()->getNumInteriorRings(), 0);
    4305           1 :         EXPECT_FALSE(poGeom->hasEmptyParts());
    4306           1 :         delete poGeom;
    4307             :     }
    4308             :     {
    4309           1 :         OGRGeometry *poGeom = nullptr;
    4310           1 :         OGRGeometryFactory::createFromWkt("COMPOUNDCURVE ((0 0,1 1))", nullptr,
    4311             :                                           &poGeom);
    4312           1 :         ASSERT_NE(poGeom, nullptr);
    4313           1 :         EXPECT_FALSE(poGeom->hasEmptyParts());
    4314           1 :         poGeom->removeEmptyParts();
    4315           1 :         EXPECT_EQ(poGeom->toCompoundCurve()->getNumCurves(), 1);
    4316           1 :         delete poGeom;
    4317             :     }
    4318             :     {
    4319           1 :         OGRGeometry *poGeom = nullptr;
    4320           1 :         OGRGeometryFactory::createFromWkt("COMPOUNDCURVE ((0 0,1 1),(1 1,2 2))",
    4321             :                                           nullptr, &poGeom);
    4322           1 :         ASSERT_NE(poGeom, nullptr);
    4323           1 :         poGeom->toCompoundCurve()->getCurve(1)->empty();
    4324           1 :         EXPECT_EQ(poGeom->toCompoundCurve()->getNumCurves(), 2);
    4325           1 :         EXPECT_TRUE(poGeom->hasEmptyParts());
    4326           1 :         poGeom->removeEmptyParts();
    4327           1 :         EXPECT_FALSE(poGeom->hasEmptyParts());
    4328           1 :         EXPECT_EQ(poGeom->toCompoundCurve()->getNumCurves(), 1);
    4329           1 :         delete poGeom;
    4330             :     }
    4331             :     {
    4332           1 :         OGRGeometry *poGeom = nullptr;
    4333           1 :         OGRGeometryFactory::createFromWkt("GEOMETRYCOLLECTION (POINT(0 0))",
    4334             :                                           nullptr, &poGeom);
    4335           1 :         ASSERT_NE(poGeom, nullptr);
    4336           1 :         EXPECT_FALSE(poGeom->hasEmptyParts());
    4337           1 :         poGeom->removeEmptyParts();
    4338           1 :         EXPECT_EQ(poGeom->toGeometryCollection()->getNumGeometries(), 1);
    4339           1 :         delete poGeom;
    4340             :     }
    4341             :     {
    4342           1 :         OGRGeometry *poGeom = nullptr;
    4343           1 :         OGRGeometryFactory::createFromWkt(
    4344             :             "GEOMETRYCOLLECTION (POINT EMPTY,POINT(0 0),POINT EMPTY)", nullptr,
    4345             :             &poGeom);
    4346           1 :         ASSERT_NE(poGeom, nullptr);
    4347           1 :         EXPECT_EQ(poGeom->toGeometryCollection()->getNumGeometries(), 3);
    4348           1 :         EXPECT_TRUE(poGeom->hasEmptyParts());
    4349           1 :         poGeom->removeEmptyParts();
    4350           1 :         EXPECT_FALSE(poGeom->hasEmptyParts());
    4351           1 :         EXPECT_EQ(poGeom->toGeometryCollection()->getNumGeometries(), 1);
    4352           1 :         delete poGeom;
    4353             :     }
    4354             :     {
    4355           1 :         OGRGeometry *poGeom = nullptr;
    4356           1 :         OGRGeometryFactory::createFromWkt("GEOMETRYCOLLECTION EMPTY", nullptr,
    4357             :                                           &poGeom);
    4358           1 :         ASSERT_NE(poGeom, nullptr);
    4359           1 :         OGRGeometry *poPoly = nullptr;
    4360           1 :         OGRGeometryFactory::createFromWkt("POLYGON ((0 0,0 1,1 0,0 0))",
    4361             :                                           nullptr, &poPoly);
    4362           1 :         EXPECT_NE(poPoly, nullptr);
    4363           1 :         if (poPoly)
    4364             :         {
    4365           1 :             poPoly->toPolygon()->addRingDirectly(new OGRLinearRing());
    4366           1 :             poGeom->toGeometryCollection()->addGeometryDirectly(poPoly);
    4367           1 :             EXPECT_EQ(poGeom->toGeometryCollection()->getNumGeometries(), 1);
    4368           1 :             EXPECT_TRUE(poGeom->hasEmptyParts());
    4369           1 :             poGeom->removeEmptyParts();
    4370           1 :             EXPECT_FALSE(poGeom->hasEmptyParts());
    4371           1 :             EXPECT_EQ(poGeom->toGeometryCollection()->getNumGeometries(), 1);
    4372             :         }
    4373           1 :         delete poGeom;
    4374             :     }
    4375             :     {
    4376           1 :         OGRGeometry *poGeom = nullptr;
    4377           1 :         OGRGeometryFactory::createFromWkt(
    4378             :             "POLYHEDRALSURFACE (((0 0,0 1,1 1,0 0)))", nullptr, &poGeom);
    4379           1 :         ASSERT_NE(poGeom, nullptr);
    4380           1 :         EXPECT_FALSE(poGeom->hasEmptyParts());
    4381           1 :         poGeom->removeEmptyParts();
    4382           1 :         EXPECT_EQ(poGeom->toPolyhedralSurface()->getNumGeometries(), 1);
    4383           1 :         delete poGeom;
    4384             :     }
    4385             :     {
    4386           1 :         OGRGeometry *poGeom = nullptr;
    4387           1 :         OGRGeometryFactory::createFromWkt(
    4388             :             "POLYHEDRALSURFACE (((0 0,0 1,1 1,0 0)))", nullptr, &poGeom);
    4389           1 :         ASSERT_NE(poGeom, nullptr);
    4390           1 :         poGeom->toPolyhedralSurface()->addGeometryDirectly(new OGRPolygon());
    4391           1 :         EXPECT_EQ(poGeom->toPolyhedralSurface()->getNumGeometries(), 2);
    4392           1 :         EXPECT_TRUE(poGeom->hasEmptyParts());
    4393           1 :         poGeom->removeEmptyParts();
    4394           1 :         EXPECT_FALSE(poGeom->hasEmptyParts());
    4395           1 :         EXPECT_EQ(poGeom->toPolyhedralSurface()->getNumGeometries(), 1);
    4396           1 :         delete poGeom;
    4397             :     }
    4398             : }
    4399             : 
    4400             : // Test OGRCurve::reversePoints()
    4401           4 : TEST_F(test_ogr, OGRCurve_reversePoints)
    4402             : {
    4403             :     {
    4404           1 :         OGRGeometry *poGeom = nullptr;
    4405           1 :         OGRGeometryFactory::createFromWkt(
    4406             :             "COMPOUNDCURVE ZM (CIRCULARSTRING ZM (0 0 10 20,1 1 11 21,2 0 12 "
    4407             :             "22),(2 0 12 22,3 0 13 2))",
    4408             :             nullptr, &poGeom);
    4409           1 :         ASSERT_NE(poGeom, nullptr);
    4410           1 :         poGeom->toCurve()->reversePoints();
    4411           1 :         char *pszWKT = nullptr;
    4412           1 :         poGeom->exportToWkt(&pszWKT, wkbVariantIso);
    4413           1 :         EXPECT_TRUE(pszWKT != nullptr);
    4414           1 :         if (pszWKT)
    4415             :         {
    4416           1 :             EXPECT_STREQ(
    4417             :                 pszWKT, "COMPOUNDCURVE ZM ((3 0 13 2,2 0 12 22),CIRCULARSTRING "
    4418             :                         "ZM (2 0 12 22,1 1 11 21,0 0 10 20))");
    4419             :         }
    4420           1 :         CPLFree(pszWKT);
    4421           1 :         delete poGeom;
    4422             :     }
    4423             : }
    4424             : 
    4425             : // Test OGRGeometryFactory::transformWithOptions()
    4426           4 : TEST_F(test_ogr, transformWithOptions)
    4427             : {
    4428             :     // Projected CRS to national geographic CRS (not including poles or antimeridian)
    4429           1 :     auto [poGeom, err] = OGRGeometryFactory::createFromWkt(
    4430           1 :         "LINESTRING(700000 6600000, 700001 6600001)");
    4431           1 :     ASSERT_NE(poGeom, nullptr);
    4432             : 
    4433           1 :     OGRSpatialReference oEPSG_2154;
    4434           1 :     oEPSG_2154.importFromEPSG(2154);  // "RGF93 v1 / Lambert-93"
    4435           1 :     OGRSpatialReference oEPSG_4171;
    4436           1 :     oEPSG_4171.importFromEPSG(4171);  // "RGF93 v1"
    4437           1 :     oEPSG_4171.SetAxisMappingStrategy(OAMS_TRADITIONAL_GIS_ORDER);
    4438             :     auto poCT = std::unique_ptr<OGRCoordinateTransformation>(
    4439           1 :         OGRCreateCoordinateTransformation(&oEPSG_2154, &oEPSG_4171));
    4440           1 :     OGRGeometryFactory::TransformWithOptionsCache oCache;
    4441             :     auto poNewGeom =
    4442             :         std::unique_ptr<OGRGeometry>(OGRGeometryFactory::transformWithOptions(
    4443           1 :             poGeom.get(), poCT.get(), nullptr, oCache));
    4444           1 :     ASSERT_NE(poNewGeom, nullptr);
    4445           1 :     EXPECT_NEAR(poNewGeom->toLineString()->getX(0), 3, 1e-8);
    4446           1 :     EXPECT_NEAR(poNewGeom->toLineString()->getY(0), 46.5, 1e-8);
    4447             : }
    4448             : 
    4449             : #ifdef HAVE_GEOS
    4450             : 
    4451             : // Test OGRGeometryFactory::transformWithOptions()
    4452           4 : TEST_F(test_ogr, transformWithOptions_GEOS)
    4453             : {
    4454             :     // Projected CRS to national geographic CRS including antimeridian
    4455           1 :     auto [poGeom, err] = OGRGeometryFactory::createFromWkt(
    4456           1 :         "LINESTRING(657630.64 4984896.17,815261.43 4990738.26)");
    4457           1 :     ASSERT_NE(poGeom, nullptr);
    4458             : 
    4459           1 :     OGRSpatialReference oEPSG_6329;
    4460           1 :     oEPSG_6329.importFromEPSG(6329);  // "NAD83(2011) / UTM zone 60N"
    4461           1 :     OGRSpatialReference oEPSG_6318;
    4462           1 :     oEPSG_6318.importFromEPSG(6318);  // "NAD83(2011)"
    4463           1 :     oEPSG_6318.SetAxisMappingStrategy(OAMS_TRADITIONAL_GIS_ORDER);
    4464             :     auto poCT = std::unique_ptr<OGRCoordinateTransformation>(
    4465           1 :         OGRCreateCoordinateTransformation(&oEPSG_6329, &oEPSG_6318));
    4466           1 :     OGRGeometryFactory::TransformWithOptionsCache oCache;
    4467             :     auto poNewGeom =
    4468             :         std::unique_ptr<OGRGeometry>(OGRGeometryFactory::transformWithOptions(
    4469           1 :             poGeom.get(), poCT.get(), nullptr, oCache));
    4470           1 :     ASSERT_NE(poNewGeom, nullptr);
    4471           1 :     EXPECT_EQ(poNewGeom->getGeometryType(), wkbMultiLineString);
    4472           1 :     if (poNewGeom->getGeometryType() == wkbMultiLineString)
    4473             :     {
    4474           1 :         const auto poMLS = poNewGeom->toMultiLineString();
    4475           1 :         EXPECT_EQ(poMLS->getNumGeometries(), 2);
    4476           1 :         if (poMLS->getNumGeometries() == 2)
    4477             :         {
    4478           1 :             const auto poLS = poMLS->getGeometryRef(0);
    4479           1 :             EXPECT_EQ(poLS->getNumPoints(), 2);
    4480           1 :             if (poLS->getNumPoints() == 2)
    4481             :             {
    4482           1 :                 EXPECT_NEAR(poLS->getX(0), 179, 1e-6);
    4483           1 :                 EXPECT_NEAR(poLS->getY(0), 45, 1e-6);
    4484           1 :                 EXPECT_NEAR(poLS->getX(1), 180, 1e-6);
    4485           1 :                 EXPECT_NEAR(poLS->getY(1), 45.004384301691303, 1e-6);
    4486             :             }
    4487             :         }
    4488             :     }
    4489             : }
    4490             : #endif
    4491             : 
    4492             : // Test OGRCurvePolygon::addRingDirectly
    4493           4 : TEST_F(test_ogr, OGRCurvePolygon_addRingDirectly)
    4494             : {
    4495           1 :     OGRCurvePolygon cp;
    4496             :     OGRGeometry *ring;
    4497             : 
    4498             :     // closed CircularString
    4499           1 :     OGRGeometryFactory::createFromWkt(
    4500             :         "CIRCULARSTRING (0 0, 1 1, 2 0, 1 -1, 0 0)", nullptr, &ring);
    4501           1 :     ASSERT_TRUE(ring);
    4502           1 :     EXPECT_EQ(cp.addRingDirectly(ring->toCurve()), OGRERR_NONE);
    4503             : 
    4504             :     // open CircularString
    4505           1 :     OGRGeometryFactory::createFromWkt("CIRCULARSTRING (0 0, 1 1, 2 0)", nullptr,
    4506             :                                       &ring);
    4507           1 :     ASSERT_TRUE(ring);
    4508             :     {
    4509             :         CPLConfigOptionSetter oSetter("OGR_GEOMETRY_ACCEPT_UNCLOSED_RING", "NO",
    4510           1 :                                       false);
    4511           1 :         ASSERT_EQ(cp.addRingDirectly(ring->toCurve()),
    4512             :                   OGRERR_UNSUPPORTED_GEOMETRY_TYPE);
    4513             :     }
    4514           1 :     EXPECT_EQ(cp.addRingDirectly(ring->toCurve()), OGRERR_NONE);
    4515             : 
    4516             :     // closed CompoundCurve
    4517           1 :     OGRGeometryFactory::createFromWkt(
    4518             :         "COMPOUNDCURVE( CIRCULARSTRING (0 0, 1 1, 2 0), (2 0, 0 0))", nullptr,
    4519             :         &ring);
    4520           1 :     ASSERT_TRUE(ring);
    4521           1 :     EXPECT_EQ(cp.addRingDirectly(ring->toCurve()), OGRERR_NONE);
    4522             : 
    4523             :     // closed LineString
    4524           1 :     OGRGeometryFactory::createFromWkt("LINESTRING (0 0, 1 0, 1 1, 0 1, 0 0)",
    4525             :                                       nullptr, &ring);
    4526           1 :     ASSERT_TRUE(ring);
    4527           1 :     EXPECT_EQ(cp.addRingDirectly(ring->toCurve()), OGRERR_NONE);
    4528             : 
    4529             :     // LinearRing
    4530           1 :     auto lr = std::make_unique<OGRLinearRing>();
    4531           1 :     lr->addPoint(0, 0);
    4532           1 :     lr->addPoint(1, 0);
    4533           1 :     lr->addPoint(1, 1);
    4534           1 :     lr->addPoint(0, 1);
    4535           1 :     lr->addPoint(0, 0);
    4536           1 :     ASSERT_TRUE(ring);
    4537           1 :     ASSERT_EQ(cp.addRingDirectly(lr.get()), OGRERR_UNSUPPORTED_GEOMETRY_TYPE);
    4538             : }
    4539             : 
    4540             : // Test OGRPolygon::addRingDirectly
    4541           4 : TEST_F(test_ogr, OGRPolygon_addRingDirectly)
    4542             : {
    4543           1 :     OGRPolygon p;
    4544             :     OGRGeometry *ring;
    4545             : 
    4546             :     // closed CircularString
    4547           1 :     OGRGeometryFactory::createFromWkt(
    4548             :         "CIRCULARSTRING (0 0, 1 1, 2 0, 1 -1, 0 0)", nullptr, &ring);
    4549           1 :     ASSERT_TRUE(ring);
    4550           1 :     EXPECT_EQ(p.addRingDirectly(ring->toCurve()),
    4551             :               OGRERR_UNSUPPORTED_GEOMETRY_TYPE);
    4552           1 :     delete ring;
    4553             : 
    4554             :     // closed LineString
    4555           1 :     OGRGeometryFactory::createFromWkt("LINESTRING (0 0, 1 0, 1 1, 0 1, 0 0)",
    4556             :                                       nullptr, &ring);
    4557           1 :     ASSERT_TRUE(ring);
    4558           1 :     EXPECT_EQ(p.addRingDirectly(ring->toCurve()),
    4559             :               OGRERR_UNSUPPORTED_GEOMETRY_TYPE);
    4560           1 :     delete ring;
    4561             : 
    4562             :     // open LineString
    4563           1 :     OGRGeometryFactory::createFromWkt("LINESTRING (0 0, 1 0)", nullptr, &ring);
    4564           1 :     ASSERT_TRUE(ring);
    4565           1 :     EXPECT_EQ(p.addRingDirectly(ring->toCurve()),
    4566             :               OGRERR_UNSUPPORTED_GEOMETRY_TYPE);
    4567           1 :     delete ring;
    4568             : 
    4569             :     // LinearRing
    4570           1 :     auto lr = std::make_unique<OGRLinearRing>();
    4571           1 :     lr->addPoint(0, 0);
    4572           1 :     lr->addPoint(1, 0);
    4573           1 :     lr->addPoint(1, 1);
    4574           1 :     lr->addPoint(0, 1);
    4575           1 :     lr->addPoint(0, 0);
    4576           1 :     ASSERT_EQ(p.addRingDirectly(lr.release()), OGRERR_NONE);
    4577             : }
    4578             : 
    4579           4 : TEST_F(test_ogr, OGRFeature_SetGeometry)
    4580             : {
    4581           1 :     OGRFeatureDefn *poFeatureDefn = new OGRFeatureDefn();
    4582           1 :     poFeatureDefn->Reference();
    4583             : 
    4584           1 :     OGRFeature oFeat(poFeatureDefn);
    4585           1 :     auto [poGeom, err] = OGRGeometryFactory::createFromWkt("POINT (3 7)");
    4586           1 :     ASSERT_EQ(err, OGRERR_NONE);
    4587             : 
    4588           1 :     ASSERT_EQ(oFeat.SetGeometry(std::move(poGeom)), OGRERR_NONE);
    4589           1 :     EXPECT_EQ(oFeat.GetGeometryRef()->toPoint()->getX(), 3);
    4590           1 :     EXPECT_EQ(oFeat.GetGeometryRef()->toPoint()->getY(), 7);
    4591             : 
    4592             :     // set it again to make sure previous feature geometry is freed
    4593           1 :     std::tie(poGeom, err) = OGRGeometryFactory::createFromWkt("POINT (2 8)");
    4594           1 :     ASSERT_EQ(err, OGRERR_NONE);
    4595           1 :     ASSERT_EQ(oFeat.SetGeometry(std::move(poGeom)), OGRERR_NONE);
    4596           1 :     EXPECT_EQ(oFeat.GetGeometryRef()->toPoint()->getX(), 2);
    4597           1 :     EXPECT_EQ(oFeat.GetGeometryRef()->toPoint()->getY(), 8);
    4598             : 
    4599           1 :     poFeatureDefn->Release();
    4600             : }
    4601             : 
    4602           4 : TEST_F(test_ogr, OGRFeature_SetGeomField)
    4603             : {
    4604           1 :     OGRFeatureDefn *poFeatureDefn = new OGRFeatureDefn();
    4605           1 :     poFeatureDefn->Reference();
    4606             : 
    4607           1 :     OGRGeomFieldDefn oGeomField("second", wkbPoint);
    4608           1 :     poFeatureDefn->AddGeomFieldDefn(&oGeomField);
    4609             : 
    4610           1 :     OGRFeature oFeat(poFeatureDefn);
    4611             : 
    4612             :     // failure
    4613             :     {
    4614           1 :         auto [poGeom, err] = OGRGeometryFactory::createFromWkt("POINT (3 7)");
    4615           1 :         ASSERT_EQ(err, OGRERR_NONE);
    4616           1 :         EXPECT_EQ(oFeat.SetGeomField(13, std::move(poGeom)), OGRERR_FAILURE);
    4617             :     }
    4618             : 
    4619             :     // success
    4620             :     {
    4621           1 :         auto [poGeom, err] = OGRGeometryFactory::createFromWkt("POINT (3 7)");
    4622           1 :         ASSERT_EQ(err, OGRERR_NONE);
    4623           1 :         EXPECT_EQ(oFeat.SetGeomField(1, std::move(poGeom)), OGRERR_NONE);
    4624             :     }
    4625             : 
    4626           1 :     poFeatureDefn->Release();
    4627             : }
    4628             : 
    4629           4 : TEST_F(test_ogr, GetArrowStream_DateTime_As_String)
    4630             : {
    4631             :     auto poDS = std::unique_ptr<GDALDataset>(
    4632             :         GetGDALDriverManager()->GetDriverByName("MEM")->Create(
    4633           1 :             "", 0, 0, 0, GDT_Unknown, nullptr));
    4634           1 :     auto poLayer = poDS->CreateLayer("test", nullptr, wkbNone);
    4635           1 :     OGRFieldDefn oFieldDefn("dt", OFTDateTime);
    4636           1 :     EXPECT_EQ(poLayer->CreateField(&oFieldDefn), OGRERR_NONE);
    4637             :     struct ArrowArrayStream stream;
    4638           1 :     CPLStringList aosOptions;
    4639           1 :     aosOptions.SetNameValue("INCLUDE_FID", "NO");
    4640           1 :     aosOptions.SetNameValue("DATETIME_AS_STRING", "YES");
    4641           1 :     ASSERT_TRUE(poLayer->GetArrowStream(&stream, aosOptions.List()));
    4642             :     struct ArrowSchema schema;
    4643           1 :     memset(&schema, 0, sizeof(schema));
    4644           1 :     EXPECT_EQ(stream.get_schema(&stream, &schema), 0);
    4645           1 :     EXPECT_TRUE(schema.n_children == 1 &&
    4646             :                 strcmp(schema.children[0]->format, "u") == 0)
    4647           0 :         << schema.n_children;
    4648           1 :     if (schema.n_children == 1 && strcmp(schema.children[0]->format, "u") == 0)
    4649             :     {
    4650           1 :         EXPECT_TRUE(schema.children[0]->metadata != nullptr);
    4651           1 :         if (schema.children[0]->metadata)
    4652             :         {
    4653             :             auto oMapKeyValue =
    4654           2 :                 OGRParseArrowMetadata(schema.children[0]->metadata);
    4655           1 :             EXPECT_EQ(oMapKeyValue.size(), 1);
    4656           1 :             if (oMapKeyValue.size() == 1)
    4657             :             {
    4658           1 :                 EXPECT_STREQ(oMapKeyValue.begin()->first.c_str(),
    4659             :                              "GDAL:OGR:type");
    4660           1 :                 EXPECT_STREQ(oMapKeyValue.begin()->second.c_str(), "DateTime");
    4661             :             }
    4662             :         }
    4663             :     }
    4664           1 :     schema.release(&schema);
    4665           1 :     stream.release(&stream);
    4666             : }
    4667             : 
    4668             : // Test OGRFeatureDefn::GetFieldSubTypeByName()
    4669           4 : TEST_F(test_ogr, OGRFieldDefnGetFieldSubTypeByName)
    4670             : {
    4671           6 :     for (int i = 0; i < OFSTMaxSubType; i++)
    4672             :     {
    4673             :         const char *pszName =
    4674           5 :             OGRFieldDefn::GetFieldSubTypeName(static_cast<OGRFieldSubType>(i));
    4675           5 :         if (pszName != nullptr)
    4676             :         {
    4677           5 :             EXPECT_EQ(OGRFieldDefn::GetFieldSubTypeByName(pszName), i);
    4678             :         }
    4679             :     }
    4680           1 : }
    4681             : 
    4682             : // Test OGRFeatureDefn::GetFieldTypeByName()
    4683           4 : TEST_F(test_ogr, OGRFieldDefnGetFieldTypeByName)
    4684             : {
    4685          14 :     for (int i = 0; i < OFTMaxType; i++)
    4686             :     {
    4687             :         // deprecated types
    4688          13 :         if (i == OFTWideString || i == OFTWideStringList)
    4689             :         {
    4690           2 :             continue;
    4691             :         }
    4692             :         const char *pszName =
    4693          11 :             OGRFieldDefn::GetFieldTypeName(static_cast<OGRFieldType>(i));
    4694          11 :         if (pszName != nullptr)
    4695             :         {
    4696          11 :             EXPECT_EQ(OGRFieldDefn::GetFieldTypeByName(pszName), i);
    4697             :         }
    4698             :     }
    4699           1 : }
    4700             : 
    4701             : // Test OGRGeometryFactory::GetDefaultArcStepSize()
    4702           4 : TEST_F(test_ogr, GetDefaultArcStepSize)
    4703             : {
    4704           1 :     if (CPLGetConfigOption("OGR_ARC_STEPSIZE", nullptr) == nullptr)
    4705             :     {
    4706           1 :         EXPECT_EQ(OGRGeometryFactory::GetDefaultArcStepSize(), 4.0);
    4707             :     }
    4708             :     {
    4709             :         CPLConfigOptionSetter oSetter("OGR_ARC_STEPSIZE", "0.00001",
    4710           2 :                                       /* bSetOnlyIfUndefined = */ false);
    4711           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    4712           1 :         EXPECT_EQ(OGRGeometryFactory::GetDefaultArcStepSize(), 1e-2);
    4713           1 :         EXPECT_TRUE(
    4714             :             strstr(CPLGetLastErrorMsg(),
    4715             :                    "Too small value for OGR_ARC_STEPSIZE. Clamping it to"));
    4716             :     }
    4717             :     {
    4718             :         CPLConfigOptionSetter oSetter("OGR_ARC_STEPSIZE", "190",
    4719           2 :                                       /* bSetOnlyIfUndefined = */ false);
    4720           2 :         CPLErrorHandlerPusher oErrorHandler(CPLQuietErrorHandler);
    4721           1 :         EXPECT_EQ(OGRGeometryFactory::GetDefaultArcStepSize(), 180);
    4722           1 :         EXPECT_TRUE(
    4723             :             strstr(CPLGetLastErrorMsg(),
    4724             :                    "Too large value for OGR_ARC_STEPSIZE. Clamping it to"));
    4725             :     }
    4726           1 : }
    4727             : 
    4728           4 : TEST_F(test_ogr, OGRPolygon_two_vertex_constructor)
    4729             : {
    4730           2 :     OGRPolygon p(1, 2, 3, 4);
    4731           1 :     char *outWKT = nullptr;
    4732           1 :     p.exportToWkt(&outWKT, wkbVariantIso);
    4733           1 :     EXPECT_STREQ(outWKT, "POLYGON ((1 2,1 4,3 4,3 2,1 2))");
    4734           1 :     CPLFree(outWKT);
    4735           1 : }
    4736             : 
    4737             : // Test OGRGeometryCollection::addComponents()
    4738           4 : TEST_F(test_ogr, OGRGeometryCollection_addComponents)
    4739             : {
    4740           1 :     auto gc1 = std::make_unique<OGRGeometryCollection>();
    4741           1 :     auto gc2 = std::make_unique<OGRGeometryCollection>();
    4742             : 
    4743           3 :     for (const auto &wkt : {"POINT (3 7)", "LINESTRING M (9 3 2, 4 1 9)"})
    4744             :     {
    4745           2 :         auto [poGeom, err] = OGRGeometryFactory::createFromWkt(wkt);
    4746           2 :         ASSERT_EQ(err, OGRERR_NONE);
    4747           2 :         gc1->addGeometry(std::move(poGeom));
    4748             :     }
    4749             : 
    4750           2 :     for (const auto &wkt :
    4751           3 :          {"POINT Z (6 4 8)", "POLYGON ((0 0, 1 1, 0 1, 0 0))"})
    4752             :     {
    4753           2 :         auto [poGeom, err] = OGRGeometryFactory::createFromWkt(wkt);
    4754           2 :         ASSERT_EQ(err, OGRERR_NONE);
    4755           2 :         gc2->addGeometry(std::move(poGeom));
    4756             :     }
    4757             : 
    4758           1 :     ASSERT_EQ(gc1->addGeometryComponents(std::move(gc2)), OGRERR_NONE);
    4759             : 
    4760           1 :     OGRWktOptions wktOptions;
    4761           1 :     wktOptions.variant = wkbVariantIso;
    4762             : 
    4763           1 :     EXPECT_EQ(gc1->getNumGeometries(), 4);
    4764           1 :     EXPECT_EQ(gc1->getGeometryType(), wkbGeometryCollectionZM);
    4765           2 :     EXPECT_EQ(gc1->getGeometryRef(0)->exportToWkt(wktOptions),
    4766             :               "POINT ZM (3 7 0 0)");
    4767           2 :     EXPECT_EQ(gc1->getGeometryRef(1)->exportToWkt(wktOptions),
    4768             :               "LINESTRING ZM (9 3 0 2,4 1 0 9)");
    4769           2 :     EXPECT_EQ(gc1->getGeometryRef(2)->exportToWkt(wktOptions),
    4770             :               "POINT ZM (6 4 8 0)");
    4771           2 :     EXPECT_EQ(gc1->getGeometryRef(3)->exportToWkt(wktOptions),
    4772             :               "POLYGON ZM ((0 0 0 0,1 1 0 0,0 1 0 0,0 0 0 0))");
    4773             : }
    4774             : 
    4775             : struct ProgressData
    4776             : {
    4777             :     std::vector<double> adfRatios;
    4778             :     bool bInterrupt{false};
    4779             : };
    4780             : 
    4781         107 : int CPL_STDCALL MyProgress(double dfComplete, const char *, void *userData)
    4782             : {
    4783         107 :     auto *progressData = static_cast<ProgressData *>(userData);
    4784         107 :     progressData->adfRatios.push_back(dfComplete);
    4785         107 :     return !progressData->bInterrupt;
    4786             : }
    4787             : 
    4788           4 : TEST_F(test_ogr, GEOSProgress)
    4789             : {
    4790           1 :     auto coll = std::make_unique<OGRGeometryCollection>();
    4791             : 
    4792          11 :     for (int i = 0; i < 10; i++)
    4793             :     {
    4794         110 :         for (int j = 0; j < 10; j++)
    4795             :         {
    4796             :             std::unique_ptr<OGRGeometry> g(OGRGeometry::FromHandle(
    4797         100 :                 OGR_G_CreateFromEnvelope(i, j, i + 1, j + 1, nullptr)));
    4798         100 :             coll->addGeometry(std::move(g));
    4799             :         }
    4800             :     }
    4801             : 
    4802             : #if GEOS_VERSION_MAJOR > 3 ||                                                  \
    4803             :     (GEOS_VERSION_MAJOR == 3 && GEOS_VERSION_MINOR >= 15)
    4804             :     // Run to completion
    4805             :     {
    4806           1 :         ProgressData progressData;
    4807             :         std::unique_ptr<OGRGeometry> result(
    4808           1 :             coll->UnaryUnion(MyProgress, &progressData));
    4809             : 
    4810           1 :         ASSERT_TRUE(progressData.adfRatios.size() > 3);
    4811           1 :         EXPECT_TRUE(std::is_sorted(progressData.adfRatios.begin(),
    4812             :                                    progressData.adfRatios.end()));
    4813           1 :         EXPECT_EQ(progressData.adfRatios.back(), 1);
    4814             :     }
    4815             : 
    4816             :     // Interrupt GEOS
    4817             :     {
    4818           2 :         CPLErrorStateBackuper oBackuper(CPLQuietErrorHandler);
    4819             : 
    4820           2 :         ProgressData progressData;
    4821           1 :         progressData.bInterrupt = true;
    4822             :         std::unique_ptr<OGRGeometry> result(
    4823           2 :             coll->UnaryUnion(MyProgress, &progressData));
    4824             : 
    4825           1 :         EXPECT_EQ(CPLGetLastErrorType(), CE_Failure);
    4826           1 :         EXPECT_TRUE(cpl::starts_with(CPLGetLastErrorMsg(), "Interrupted"));
    4827           1 :         EXPECT_EQ(result, nullptr);
    4828             :     }
    4829             : #else
    4830             :     // Run to completion
    4831             :     {
    4832             :         ProgressData progressData;
    4833             :         std::unique_ptr<OGRGeometry> result(
    4834             :             coll->UnaryUnion(&MyProgress, &progressData));
    4835             :         EXPECT_TRUE(progressData.adfRatios.empty());
    4836             :     }
    4837             : #endif
    4838             : }
    4839             : 
    4840             : }  // namespace

Generated by: LCOV version 1.14