LCOV - code coverage report
Current view: top level - frmts/wmts - wmtsdataset.cpp (source / functions) Hit Total Coverage
Test: gdal_filtered.info Lines: 1020 1121 91.0 %
Date: XXXX-XX-XX Functions: 40 42 95.2 %
Branches: 1693 3165 53.5 %

           Branch data     Line data    Source code
       1                 :            : /******************************************************************************
       2                 :            :  *
       3                 :            :  * Project:  GDAL WMTS driver
       4                 :            :  * Purpose:  Implement GDAL WMTS support
       5                 :            :  * Author:   Even Rouault, <even dot rouault at spatialys dot com>
       6                 :            :  * Funded by Land Information New Zealand (LINZ)
       7                 :            :  *
       8                 :            :  **********************************************************************
       9                 :            :  * Copyright (c) 2015, Even Rouault <even dot rouault at spatialys dot com>
      10                 :            :  *
      11                 :            :  * Permission is hereby granted, free of charge, to any person obtaining a
      12                 :            :  * copy of this software and associated documentation files (the "Software"),
      13                 :            :  * to deal in the Software without restriction, including without limitation
      14                 :            :  * the rights to use, copy, modify, merge, publish, distribute, sublicense,
      15                 :            :  * and/or sell copies of the Software, and to permit persons to whom the
      16                 :            :  * Software is furnished to do so, subject to the following conditions:
      17                 :            :  *
      18                 :            :  * The above copyright notice and this permission notice shall be included
      19                 :            :  * in all copies or substantial portions of the Software.
      20                 :            :  *
      21                 :            :  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
      22                 :            :  * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
      23                 :            :  * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
      24                 :            :  * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
      25                 :            :  * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
      26                 :            :  * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
      27                 :            :  * DEALINGS IN THE SOFTWARE.
      28                 :            :  ****************************************************************************/
      29                 :            : 
      30                 :            : #include "cpl_http.h"
      31                 :            : #include "cpl_minixml.h"
      32                 :            : #include "gdal_frmts.h"
      33                 :            : #include "gdal_pam.h"
      34                 :            : #include "ogr_spatialref.h"
      35                 :            : #include "../vrt/gdal_vrt.h"
      36                 :            : 
      37                 :            : #include <algorithm>
      38                 :            : #include <map>
      39                 :            : #include <set>
      40                 :            : #include <vector>
      41                 :            : #include <limits>
      42                 :            : 
      43                 :            : extern "C" void GDALRegister_WMTS();
      44                 :            : 
      45                 :            : // g++ -g -Wall -fPIC frmts/wmts/wmtsdataset.cpp -shared -o gdal_WMTS.so -Iport -Igcore -Iogr -Iogr/ogrsf_frmts -L. -lgdal
      46                 :            : 
      47                 :            : /* Set in stone by WMTS spec. In pixel/meter */
      48                 :            : #define WMTS_PITCH                      0.00028
      49                 :            : 
      50                 :            : #define WMTS_WGS84_DEG_PER_METER    (180 / M_PI / SRS_WGS84_SEMIMAJOR)
      51                 :            : 
      52                 :            : CPL_CVSID("$Id$");
      53                 :            : 
      54                 :            : typedef enum
      55                 :            : {
      56                 :            :     AUTO,
      57                 :            :     LAYER_BBOX,
      58                 :            :     TILE_MATRIX_SET,
      59                 :            :     MOST_PRECISE_TILE_MATRIX
      60                 :            : } ExtentMethod;
      61                 :            : 
      62                 :            : /************************************************************************/
      63                 :            : /* ==================================================================== */
      64                 :            : /*                            WMTSTileMatrix                            */
      65                 :            : /* ==================================================================== */
      66                 :            : /************************************************************************/
      67                 :            : 
      68                 :          1 : class WMTSTileMatrix
      69                 :            : {
      70                 :            :     public:
      71                 :            :         CPLString osIdentifier;
      72                 :            :         double    dfScaleDenominator;
      73                 :            :         double    dfPixelSize;
      74                 :            :         double    dfTLX;
      75                 :            :         double    dfTLY;
      76                 :            :         int       nTileWidth;
      77                 :            :         int       nTileHeight;
      78                 :            :         int       nMatrixWidth;
      79                 :            :         int       nMatrixHeight;
      80                 :            : };
      81                 :            : 
      82                 :            : /************************************************************************/
      83                 :            : /* ==================================================================== */
      84                 :            : /*                          WMTSTileMatrixLimits                        */
      85                 :            : /* ==================================================================== */
      86                 :            : /************************************************************************/
      87                 :            : 
      88                 :          1 : class WMTSTileMatrixLimits
      89                 :            : {
      90                 :            :     public:
      91                 :            :         CPLString osIdentifier;
      92                 :            :         int nMinTileRow;
      93                 :            :         int nMaxTileRow;
      94                 :            :         int nMinTileCol;
      95                 :            :         int nMaxTileCol;
      96                 :            : };
      97                 :            : 
      98                 :            : /************************************************************************/
      99                 :            : /* ==================================================================== */
     100                 :            : /*                          WMTSTileMatrixSet                           */
     101                 :            : /* ==================================================================== */
     102                 :            : /************************************************************************/
     103                 :            : 
     104 [ +  - ][ +  - ]:          1 : class WMTSTileMatrixSet
     105                 :            : {
     106                 :            :     public:
     107                 :            :         OGRSpatialReference         oSRS;
     108                 :            :         CPLString                   osSRS;
     109                 :            :         bool                        bBoundingBoxValid;
     110                 :            :         OGREnvelope                 sBoundingBox; /* expressed in TMS SRS */
     111                 :            :         std::vector<WMTSTileMatrix> aoTM;
     112                 :            : 
     113                 :          1 :         WMTSTileMatrixSet() :
     114                 :            :             oSRS( OGRSpatialReference() ),
     115 [ +  - ][ +  - ]:          1 :             bBoundingBoxValid(false)
                 [ +  - ]
     116                 :            :         {
     117                 :          1 :         }
     118                 :            : };
     119                 :            : 
     120                 :            : /************************************************************************/
     121                 :            : /* ==================================================================== */
     122                 :            : /*                              WMTSDataset                             */
     123                 :            : /* ==================================================================== */
     124                 :            : /************************************************************************/
     125                 :            : 
     126                 :            : class WMTSDataset : public GDALPamDataset
     127                 :            : {
     128                 :            :   friend class WMTSBand;
     129                 :            : 
     130                 :            :     CPLString                 osLayer;
     131                 :            :     CPLString                 osTMS;
     132                 :            :     CPLString                 osXML;
     133                 :            :     CPLString                 osURLFeatureInfoTemplate;
     134                 :            :     WMTSTileMatrixSet         oTMS;
     135                 :            : 
     136                 :            :     char                    **papszHTTPOptions;
     137                 :            : 
     138                 :            :     std::vector<GDALDataset*> apoDatasets;
     139                 :            :     CPLString                 osProjection;
     140                 :            :     double                    adfGT[6];
     141                 :            : 
     142                 :            :     CPLString                 osLastGetFeatureInfoURL;
     143                 :            :     CPLString                 osMetadataItemGetFeatureInfo;
     144                 :            : 
     145                 :            :     static char**       BuildHTTPRequestOpts(CPLString osOtherXML);
     146                 :            :     static CPLXMLNode*  GetCapabilitiesResponse(const CPLString& osFilename,
     147                 :            :                                                 char** papszHTTPOptions);
     148                 :            :     static CPLString    FixCRSName(const char* pszCRS);
     149                 :            :     static CPLString    Replace(const CPLString& osStr, const char* pszOld, const char* pszNew);
     150                 :            :     static CPLString    GetOperationKVPURL(CPLXMLNode* psXML,
     151                 :            :                                            const char* pszOperation);
     152                 :            :     static int          ReadTMS(CPLXMLNode* psContents,
     153                 :            :                                 const CPLString& osIdentifier,
     154                 :            :                                 const CPLString& osMaxTileMatrixIdentifier,
     155                 :            :                                 int nMaxZoomLevel,
     156                 :            :                                 WMTSTileMatrixSet& oTMS);
     157                 :            :     static int          ReadTMLimits(CPLXMLNode* psTMSLimits,
     158                 :            :                                      std::map<CPLString, WMTSTileMatrixLimits>& aoMapTileMatrixLimits);
     159                 :            : 
     160                 :            :   public:
     161                 :            :                  WMTSDataset();
     162                 :            :     virtual     ~WMTSDataset();
     163                 :            : 
     164                 :            :     virtual CPLErr GetGeoTransform(double* padfGT) override;
     165                 :            :     virtual const char* GetProjectionRef() override;
     166                 :            :     virtual const char* GetMetadataItem(const char* pszName,
     167                 :            :                                         const char* pszDomain) override;
     168                 :            : 
     169                 :            :     static GDALDataset *Open( GDALOpenInfo * );
     170                 :            :     static int          Identify( GDALOpenInfo * );
     171                 :            :     static GDALDataset *CreateCopy( const char * pszFilename,
     172                 :            :                                          GDALDataset *poSrcDS,
     173                 :            :                                          CPL_UNUSED int bStrict,
     174                 :            :                                          CPL_UNUSED char ** papszOptions,
     175                 :            :                                          CPL_UNUSED GDALProgressFunc pfnProgress,
     176                 :            :                                          CPL_UNUSED void * pProgressData );
     177                 :            : 
     178                 :            :   protected:
     179                 :            :     virtual int         CloseDependentDatasets() override;
     180                 :            : 
     181                 :            :     virtual CPLErr  IRasterIO( GDALRWFlag eRWFlag,
     182                 :            :                                int nXOff, int nYOff, int nXSize, int nYSize,
     183                 :            :                                void * pData, int nBufXSize, int nBufYSize,
     184                 :            :                                GDALDataType eBufType,
     185                 :            :                                int nBandCount, int *panBandMap,
     186                 :            :                                GSpacing nPixelSpace, GSpacing nLineSpace,
     187                 :            :                                GSpacing nBandSpace,
     188                 :            :                                GDALRasterIOExtraArg* psExtraArg) override;
     189                 :            : };
     190                 :            : 
     191                 :            : /************************************************************************/
     192                 :            : /* ==================================================================== */
     193                 :            : /*                               WMTSBand                               */
     194                 :            : /* ==================================================================== */
     195                 :            : /************************************************************************/
     196                 :            : 
     197         [ -  + ]:          1 : class WMTSBand : public GDALPamRasterBand
     198                 :            : {
     199                 :            :   public:
     200                 :            :                   WMTSBand(WMTSDataset* poDS, int nBand);
     201                 :            : 
     202                 :            :     virtual GDALRasterBand* GetOverview(int nLevel) override;
     203                 :            :     virtual int GetOverviewCount() override;
     204                 :            :     virtual GDALColorInterp GetColorInterpretation() override;
     205                 :            :     virtual const char* GetMetadataItem(const char* pszName,
     206                 :            :                                         const char* pszDomain) override;
     207                 :            : 
     208                 :            :   protected:
     209                 :            :     virtual CPLErr IReadBlock( int nBlockXOff, int nBlockYOff, void * pImage) override;
     210                 :            :     virtual CPLErr IRasterIO( GDALRWFlag, int, int, int, int,
     211                 :            :                               void *, int, int, GDALDataType,
     212                 :            :                               GSpacing, GSpacing,
     213                 :            :                               GDALRasterIOExtraArg* psExtraArg ) override;
     214                 :            : };
     215                 :            : 
     216                 :            : /************************************************************************/
     217                 :            : /*                            WMTSBand()                                */
     218                 :            : /************************************************************************/
     219                 :            : 
     220                 :          1 : WMTSBand::WMTSBand( WMTSDataset* poDSIn, int nBandIn )
     221                 :            : {
     222                 :          1 :     poDS = poDSIn;
     223                 :          1 :     nBand = nBandIn;
     224                 :          1 :     eDataType = GDT_Byte;
     225         [ +  - ]:          1 :     poDSIn->apoDatasets[0]->GetRasterBand(1)->
     226 [ +  - ][ +  - ]:          1 :         GetBlockSize(&nBlockXSize, &nBlockYSize);
     227                 :          1 : }
     228                 :            : 
     229                 :            : /************************************************************************/
     230                 :            : /*                            IReadBlock()                              */
     231                 :            : /************************************************************************/
     232                 :            : 
     233                 :          0 : CPLErr WMTSBand::IReadBlock( int nBlockXOff, int nBlockYOff, void * pImage)
     234                 :            : {
     235                 :          0 :     WMTSDataset* poGDS = (WMTSDataset*) poDS;
     236                 :          0 :     return poGDS->apoDatasets[0]->GetRasterBand(nBand)->ReadBlock(nBlockXOff, nBlockYOff, pImage);
     237                 :            : }
     238                 :            : 
     239                 :            : /************************************************************************/
     240                 :            : /*                             IRasterIO()                              */
     241                 :            : /************************************************************************/
     242                 :            : 
     243                 :          1 : CPLErr WMTSBand::IRasterIO( GDALRWFlag eRWFlag,
     244                 :            :                             int nXOff, int nYOff, int nXSize, int nYSize,
     245                 :            :                             void * pData, int nBufXSize, int nBufYSize,
     246                 :            :                             GDALDataType eBufType,
     247                 :            :                             GSpacing nPixelSpace, GSpacing nLineSpace,
     248                 :            :                             GDALRasterIOExtraArg* psExtraArg )
     249                 :            : {
     250                 :          1 :     WMTSDataset* poGDS = (WMTSDataset*) poDS;
     251                 :            : 
     252         [ +  + ]:          1 :     if( (nBufXSize < nXSize || nBufYSize < nYSize)
           [ -  +  +  - ]
         [ +  - ][ +  + ]
     253                 :          1 :         && poGDS->apoDatasets.size() > 1 && eRWFlag == GF_Read )
     254                 :            :     {
     255                 :            :         int bTried;
     256                 :            :         CPLErr eErr = TryOverviewRasterIO( eRWFlag,
     257                 :            :                                     nXOff, nYOff, nXSize, nYSize,
     258                 :            :                                     pData, nBufXSize, nBufYSize,
     259                 :            :                                     eBufType,
     260                 :            :                                     nPixelSpace, nLineSpace,
     261                 :            :                                     psExtraArg,
     262                 :          1 :                                     &bTried );
     263         [ +  - ]:          1 :         if( bTried )
     264                 :          1 :             return eErr;
     265                 :            :     }
     266                 :            : 
     267                 :          1 :     return poGDS->apoDatasets[0]->GetRasterBand(nBand)->RasterIO(
     268                 :            :                                          eRWFlag, nXOff, nYOff, nXSize, nYSize,
     269                 :            :                                          pData, nBufXSize, nBufYSize, eBufType,
     270                 :          1 :                                          nPixelSpace, nLineSpace, psExtraArg );
     271                 :            : }
     272                 :            : 
     273                 :            : /************************************************************************/
     274                 :            : /*                         GetOverviewCount()                           */
     275                 :            : /************************************************************************/
     276                 :            : 
     277                 :          1 : int WMTSBand::GetOverviewCount()
     278                 :            : {
     279                 :          1 :     WMTSDataset* poGDS = (WMTSDataset*) poDS;
     280                 :            : 
     281         [ +  + ]:          1 :     if( poGDS->apoDatasets.size() > 1 )
     282                 :          1 :         return (int)poGDS->apoDatasets.size() - 1;
     283                 :            :     else
     284                 :          1 :         return 0;
     285                 :            : }
     286                 :            : 
     287                 :            : /************************************************************************/
     288                 :            : /*                              GetOverview()                           */
     289                 :            : /************************************************************************/
     290                 :            : 
     291                 :          1 : GDALRasterBand* WMTSBand::GetOverview(int nLevel)
     292                 :            : {
     293                 :          1 :     WMTSDataset* poGDS = (WMTSDataset*) poDS;
     294                 :            : 
     295 [ +  - ][ +  + ]:          1 :     if (nLevel < 0 || nLevel >= GetOverviewCount())
                 [ +  + ]
     296                 :          1 :         return NULL;
     297                 :            : 
     298                 :          1 :     GDALDataset* poOvrDS = poGDS->apoDatasets[nLevel+1];
     299         [ +  - ]:          1 :     if (poOvrDS)
     300                 :          1 :         return poOvrDS->GetRasterBand(nBand);
     301                 :            :     else
     302                 :          1 :         return NULL;
     303                 :            : }
     304                 :            : 
     305                 :            : /************************************************************************/
     306                 :            : /*                   GetColorInterpretation()                           */
     307                 :            : /************************************************************************/
     308                 :            : 
     309                 :          1 : GDALColorInterp WMTSBand::GetColorInterpretation()
     310                 :            : {
     311                 :          1 :     WMTSDataset* poGDS = (WMTSDataset*) poDS;
     312         [ -  + ]:          1 :     if (poGDS->nBands == 1)
     313                 :            :     {
     314                 :          0 :         return GCI_GrayIndex;
     315                 :            :     }
     316 [ +  - ][ +  - ]:          1 :     else if (poGDS->nBands == 3 || poGDS->nBands == 4)
     317                 :            :     {
     318         [ +  + ]:          1 :         if (nBand == 1)
     319                 :          1 :             return GCI_RedBand;
     320         [ +  + ]:          1 :         else if (nBand == 2)
     321                 :          1 :             return GCI_GreenBand;
     322         [ +  + ]:          1 :         else if (nBand == 3)
     323                 :          1 :             return GCI_BlueBand;
     324         [ +  - ]:          1 :         else if (nBand == 4)
     325                 :          1 :             return GCI_AlphaBand;
     326                 :            :     }
     327                 :            : 
     328                 :          1 :     return GCI_Undefined;
     329                 :            : }
     330                 :            : 
     331                 :            : /************************************************************************/
     332                 :            : /*                         GetMetadataItem()                            */
     333                 :            : /************************************************************************/
     334                 :            : 
     335                 :          1 : const char *WMTSBand::GetMetadataItem( const char * pszName,
     336                 :            :                                        const char * pszDomain )
     337                 :            : {
     338                 :          1 :     WMTSDataset* poGDS = (WMTSDataset*) poDS;
     339                 :            : 
     340                 :            : /* ==================================================================== */
     341                 :            : /*      LocationInfo handling.                                          */
     342                 :            : /* ==================================================================== */
     343 [ +  - ][ +  + ]:          1 :     if( pszDomain != NULL && EQUAL(pszDomain,"LocationInfo") &&
         [ +  - ][ +  -  
             +  -  +  + ]
                 [ +  + ]
     344                 :          1 :         pszName != NULL && STARTS_WITH_CI(pszName, "Pixel_") &&
     345                 :          1 :         !poGDS->oTMS.aoTM.empty() &&
     346                 :          1 :         !poGDS->osURLFeatureInfoTemplate.empty() )
     347                 :            :     {
     348                 :            :         int iPixel, iLine;
     349                 :            : 
     350                 :            : /* -------------------------------------------------------------------- */
     351                 :            : /*      What pixel are we aiming at?                                    */
     352                 :            : /* -------------------------------------------------------------------- */
     353         [ -  + ]:          1 :         if( sscanf( pszName+6, "%d_%d", &iPixel, &iLine ) != 2 )
     354                 :          0 :             return NULL;
     355                 :            : 
     356                 :          1 :         const WMTSTileMatrix& oTM = poGDS->oTMS.aoTM.back();
     357                 :            : 
     358                 :          1 :         iPixel += (int)floor(0.5 + (poGDS->adfGT[0] - oTM.dfTLX) / oTM.dfPixelSize);
     359                 :          1 :         iLine += (int)floor(0.5 + (oTM.dfTLY - poGDS->adfGT[3]) / oTM.dfPixelSize);
     360                 :            : 
     361                 :          1 :         CPLString osURL(poGDS->osURLFeatureInfoTemplate);
     362 [ +  - ][ +  - ]:          1 :         osURL = WMTSDataset::Replace(osURL, "{TileMatrixSet}", poGDS->osTMS);
         [ +  - ][ +  - ]
     363 [ +  - ][ +  - ]:          1 :         osURL = WMTSDataset::Replace(osURL, "{TileMatrix}", oTM.osIdentifier);
         [ +  - ][ +  - ]
     364                 :            :         osURL = WMTSDataset::Replace(osURL, "{TileCol}",
     365 [ +  - ][ +  - ]:          1 :                                      CPLSPrintf("%d", iPixel / oTM.nTileWidth));
         [ +  - ][ +  - ]
     366                 :            :         osURL = WMTSDataset::Replace(osURL, "{TileRow}",
     367 [ +  - ][ +  - ]:          1 :                                      CPLSPrintf("%d", iLine / oTM.nTileHeight));
         [ +  - ][ +  - ]
     368                 :            :         osURL = WMTSDataset::Replace(osURL, "{I}",
     369 [ +  - ][ +  - ]:          1 :                                      CPLSPrintf("%d", iPixel % oTM.nTileWidth));
         [ +  - ][ +  - ]
     370                 :            :         osURL = WMTSDataset::Replace(osURL, "{J}",
     371 [ +  - ][ +  - ]:          1 :                                      CPLSPrintf("%d", iLine % oTM.nTileHeight));
         [ +  - ][ +  - ]
     372                 :            : 
     373 [ +  - ][ +  + ]:          1 :         if( poGDS->osLastGetFeatureInfoURL.compare(osURL) != 0 )
     374                 :            :         {
     375         [ +  - ]:          1 :             poGDS->osLastGetFeatureInfoURL = osURL;
     376 [ +  - ][ +  - ]:          1 :             poGDS->osMetadataItemGetFeatureInfo = "";
                 [ +  - ]
     377                 :          1 :             char* pszRes = NULL;
     378 [ +  - ][ +  - ]:          1 :             CPLHTTPResult* psResult = CPLHTTPFetch( osURL, poGDS->papszHTTPOptions);
     379 [ +  - ][ +  + ]:          1 :             if( psResult && psResult->nStatus == 0 && psResult->pabyData )
                 [ +  - ]
     380         [ +  - ]:          1 :                 pszRes = CPLStrdup((const char*) psResult->pabyData);
     381         [ +  - ]:          1 :             CPLHTTPDestroyResult(psResult);
     382                 :            : 
     383         [ +  + ]:          1 :             if (pszRes)
     384                 :            :             {
     385 [ +  - ][ +  - ]:          1 :                 poGDS->osMetadataItemGetFeatureInfo = "<LocationInfo>";
                 [ +  - ]
     386         [ +  - ]:          1 :                 CPLPushErrorHandler(CPLQuietErrorHandler);
     387         [ +  - ]:          1 :                 CPLXMLNode* psXML = CPLParseXMLString(pszRes);
     388         [ +  - ]:          1 :                 CPLPopErrorHandler();
     389 [ +  - ][ +  + ]:          1 :                 if (psXML != NULL && psXML->eType == CXT_Element)
     390                 :            :                 {
     391         [ +  - ]:          1 :                     if (strcmp(psXML->pszValue, "?xml") == 0)
     392                 :            :                     {
     393         [ +  - ]:          1 :                         if (psXML->psNext)
     394                 :            :                         {
     395         [ +  - ]:          1 :                             char* pszXML = CPLSerializeXMLTree(psXML->psNext);
     396         [ +  - ]:          1 :                             poGDS->osMetadataItemGetFeatureInfo += pszXML;
     397         [ +  - ]:          1 :                             CPLFree(pszXML);
     398                 :            :                         }
     399                 :            :                     }
     400                 :            :                     else
     401                 :            :                     {
     402         [ #  # ]:          0 :                         poGDS->osMetadataItemGetFeatureInfo += pszRes;
     403                 :          1 :                     }
     404                 :            :                 }
     405                 :            :                 else
     406                 :            :                 {
     407         [ +  - ]:          1 :                     char* pszEscapedXML = CPLEscapeString(pszRes, -1, CPLES_XML_BUT_QUOTES);
     408         [ +  - ]:          1 :                     poGDS->osMetadataItemGetFeatureInfo += pszEscapedXML;
     409         [ +  - ]:          1 :                     CPLFree(pszEscapedXML);
     410                 :            :                 }
     411         [ +  - ]:          1 :                 if (psXML != NULL)
     412         [ +  - ]:          1 :                     CPLDestroyXMLNode(psXML);
     413                 :            : 
     414         [ +  - ]:          1 :                 poGDS->osMetadataItemGetFeatureInfo += "</LocationInfo>";
     415         [ +  - ]:          1 :                 CPLFree(pszRes);
     416                 :            :             }
     417                 :            :         }
     418         [ +  - ]:          1 :         return poGDS->osMetadataItemGetFeatureInfo.c_str();
     419                 :            :     }
     420                 :            : 
     421                 :          1 :     return GDALPamRasterBand::GetMetadataItem(pszName, pszDomain);
     422                 :            : }
     423                 :            : 
     424                 :            : /************************************************************************/
     425                 :            : /*                          WMTSDataset()                               */
     426                 :            : /************************************************************************/
     427                 :            : 
     428                 :          1 : WMTSDataset::WMTSDataset() :
     429 [ +  - ][ +  - ]:          1 :     papszHTTPOptions(NULL)
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
     430                 :            : {
     431                 :          1 :     adfGT[0] = 0;
     432                 :          1 :     adfGT[1] = 1;
     433                 :          1 :     adfGT[2] = 0;
     434                 :          1 :     adfGT[3] = 0;
     435                 :          1 :     adfGT[4] = 0;
     436                 :          1 :     adfGT[5] = 1;
     437                 :          1 : }
     438                 :            : 
     439                 :            : /************************************************************************/
     440                 :            : /*                        ~WMTSDataset()                                */
     441                 :            : /************************************************************************/
     442                 :            : 
     443 [ +  - ][ +  - ]:          1 : WMTSDataset::~WMTSDataset()
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
     444                 :            : {
     445         [ +  - ]:          1 :     CloseDependentDatasets();
     446         [ +  - ]:          1 :     CSLDestroy(papszHTTPOptions);
     447         [ -  + ]:          1 : }
     448                 :            : 
     449                 :            : /************************************************************************/
     450                 :            : /*                      CloseDependentDatasets()                        */
     451                 :            : /************************************************************************/
     452                 :            : 
     453                 :          1 : int WMTSDataset::CloseDependentDatasets()
     454                 :            : {
     455                 :          1 :     int bRet = GDALPamDataset::CloseDependentDatasets();
     456         [ +  + ]:          1 :     if( !apoDatasets.empty() )
     457                 :            :     {
     458         [ +  + ]:          1 :         for(size_t i=0;i<apoDatasets.size();i++)
     459         [ +  - ]:          1 :             delete apoDatasets[i];
     460                 :          1 :         apoDatasets.resize(0);
     461                 :          1 :         bRet = TRUE;
     462                 :            :     }
     463                 :          1 :     return bRet;
     464                 :            : }
     465                 :            : 
     466                 :            : /************************************************************************/
     467                 :            : /*                             IRasterIO()                              */
     468                 :            : /************************************************************************/
     469                 :            : 
     470                 :          1 : CPLErr  WMTSDataset::IRasterIO( GDALRWFlag eRWFlag,
     471                 :            :                                int nXOff, int nYOff, int nXSize, int nYSize,
     472                 :            :                                void * pData, int nBufXSize, int nBufYSize,
     473                 :            :                                GDALDataType eBufType,
     474                 :            :                                int nBandCount, int *panBandMap,
     475                 :            :                                GSpacing nPixelSpace, GSpacing nLineSpace,
     476                 :            :                                GSpacing nBandSpace,
     477                 :            :                                GDALRasterIOExtraArg* psExtraArg)
     478                 :            : {
     479         [ +  + ]:          1 :     if( (nBufXSize < nXSize || nBufYSize < nYSize)
           [ -  +  +  - ]
         [ +  - ][ +  + ]
     480                 :          1 :         && apoDatasets.size() > 1 && eRWFlag == GF_Read )
     481                 :            :     {
     482                 :            :         int bTried;
     483                 :            :         CPLErr eErr = TryOverviewRasterIO( eRWFlag,
     484                 :            :                                     nXOff, nYOff, nXSize, nYSize,
     485                 :            :                                     pData, nBufXSize, nBufYSize,
     486                 :            :                                     eBufType,
     487                 :            :                                     nBandCount, panBandMap,
     488                 :            :                                     nPixelSpace, nLineSpace,
     489                 :            :                                     nBandSpace,
     490                 :            :                                     psExtraArg,
     491                 :          1 :                                     &bTried );
     492         [ +  - ]:          1 :         if( bTried )
     493                 :          1 :             return eErr;
     494                 :            :     }
     495                 :            : 
     496                 :          1 :     return apoDatasets[0]->RasterIO( eRWFlag, nXOff, nYOff, nXSize, nYSize,
     497                 :            :                                   pData, nBufXSize, nBufYSize,
     498                 :            :                                   eBufType, nBandCount, panBandMap,
     499                 :            :                                   nPixelSpace, nLineSpace, nBandSpace,
     500                 :          1 :                                   psExtraArg );
     501                 :            : }
     502                 :            : 
     503                 :            : /************************************************************************/
     504                 :            : /*                          GetGeoTransform()                           */
     505                 :            : /************************************************************************/
     506                 :            : 
     507                 :          1 : CPLErr WMTSDataset::GetGeoTransform(double* padfGT)
     508                 :            : {
     509                 :          1 :     memcpy(padfGT, adfGT, 6 * sizeof(double));
     510                 :          1 :     return CE_None;
     511                 :            : }
     512                 :            : 
     513                 :            : /************************************************************************/
     514                 :            : /*                         GetProjectionRef()                           */
     515                 :            : /************************************************************************/
     516                 :            : 
     517                 :          1 : const char* WMTSDataset::GetProjectionRef()
     518                 :            : {
     519                 :          1 :     return osProjection.c_str();
     520                 :            : }
     521                 :            : 
     522                 :            : /************************************************************************/
     523                 :            : /*                          WMTSEscapeXML()                             */
     524                 :            : /************************************************************************/
     525                 :            : 
     526                 :          1 : static CPLString WMTSEscapeXML(const char* pszUnescapedXML)
     527                 :            : {
     528                 :          1 :     CPLString osRet;
     529         [ +  - ]:          1 :     char* pszTmp = CPLEscapeString(pszUnescapedXML, -1, CPLES_XML);
     530 [ +  - ][ +  - ]:          1 :     osRet = pszTmp;
                 [ +  - ]
     531         [ +  - ]:          1 :     CPLFree(pszTmp);
     532                 :          1 :     return osRet;
     533                 :            : }
     534                 :            : 
     535                 :            : /************************************************************************/
     536                 :            : /*                         GetMetadataItem()                            */
     537                 :            : /************************************************************************/
     538                 :            : 
     539                 :          1 : const char* WMTSDataset::GetMetadataItem(const char* pszName,
     540                 :            :                                          const char* pszDomain)
     541                 :            : {
     542 [ +  - ][ +  - ]:          1 :     if( pszName != NULL && EQUAL(pszName, "XML") &&
         [ +  - ][ +  - ]
     543                 :          1 :         pszDomain != NULL && EQUAL(pszDomain, "WMTS") )
     544                 :            :     {
     545                 :          1 :         return osXML.c_str();
     546                 :            :     }
     547                 :            : 
     548                 :          1 :     return GDALPamDataset::GetMetadataItem(pszName, pszDomain);
     549                 :            : }
     550                 :            : 
     551                 :            : /************************************************************************/
     552                 :            : /*                             Identify()                               */
     553                 :            : /************************************************************************/
     554                 :            : 
     555                 :          1 : int WMTSDataset::Identify(GDALOpenInfo* poOpenInfo)
     556                 :            : {
     557         [ +  + ]:          1 :     if( STARTS_WITH_CI(poOpenInfo->pszFilename, "WMTS:") )
     558                 :          1 :         return TRUE;
     559                 :            : 
     560         [ +  + ]:          1 :     if( STARTS_WITH_CI(poOpenInfo->pszFilename, "<GDAL_WMTS") )
     561                 :          1 :         return TRUE;
     562                 :            : 
     563         [ +  + ]:          1 :     if( poOpenInfo->nHeaderBytes == 0 )
     564                 :          1 :         return FALSE;
     565                 :            : 
     566         [ +  + ]:          1 :     if( strstr((const char*)poOpenInfo->pabyHeader, "<GDAL_WMTS") )
     567                 :          1 :         return TRUE;
     568                 :            : 
     569                 :            :     return (strstr((const char*)poOpenInfo->pabyHeader,
     570                 :          1 :                   "<Capabilities") != NULL ||
     571                 :            :             strstr((const char*)poOpenInfo->pabyHeader,
     572                 :          1 :                   "<wmts:Capabilities") != NULL) &&
     573                 :            :             strstr((const char*)poOpenInfo->pabyHeader,
     574   [ +  +  -  + ]:          1 :                     "http://www.opengis.net/wmts/1.0") != NULL;
                 [ +  - ]
     575                 :            : }
     576                 :            : 
     577                 :            : /************************************************************************/
     578                 :            : /*                          QuoteIfNecessary()                          */
     579                 :            : /************************************************************************/
     580                 :            : 
     581                 :          1 : static CPLString QuoteIfNecessary(const char* pszVal)
     582                 :            : {
     583 [ +  + ][ +  - ]:          1 :     if( strchr(pszVal, ' ') || strchr(pszVal, ',') || strchr(pszVal, '=') )
                 [ +  + ]
     584                 :            :     {
     585                 :          1 :         CPLString osVal;
     586         [ +  - ]:          1 :         osVal += "\"";
     587         [ +  - ]:          1 :         osVal += pszVal;
     588         [ +  - ]:          1 :         osVal += "\"";
     589         [ +  - ]:          1 :         return osVal;
     590                 :            :     }
     591                 :            :     else
     592                 :          1 :         return pszVal;
     593                 :            : }
     594                 :            : 
     595                 :            : /************************************************************************/
     596                 :            : /*                             FixCRSName()                             */
     597                 :            : /************************************************************************/
     598                 :            : 
     599                 :          1 : CPLString WMTSDataset::FixCRSName(const char* pszCRS)
     600                 :            : {
     601 [ +  - ][ +  - ]:          1 :     while( *pszCRS == ' ' || *pszCRS == '\r' || *pszCRS == '\n' )
         [ -  + ][ -  + ]
     602                 :          0 :         pszCRS ++;
     603                 :            : 
     604                 :            :     /* http://maps.wien.gv.at/wmts/1.0.0/WMTSCapabilities.xml uses urn:ogc:def:crs:EPSG:6.18:3:3857 */
     605                 :            :     /* instead of urn:ogc:def:crs:EPSG:6.18.3:3857. Coming from an incorrect example of URN in WMTS spec */
     606                 :            :     /* https://portal.opengeospatial.org/files/?artifact_id=50398 */
     607         [ +  + ]:          1 :     if( STARTS_WITH_CI(pszCRS, "urn:ogc:def:crs:EPSG:6.18:3:") )    {
     608                 :            :         return CPLSPrintf("urn:ogc:def:crs:EPSG::%s",
     609                 :          1 :                           pszCRS + strlen("urn:ogc:def:crs:EPSG:6.18:3:"));
     610                 :            :     }
     611                 :            : 
     612         [ -  + ]:          1 :     if( EQUAL(pszCRS, "urn:ogc:def:crs:EPSG::102100") )
     613                 :          0 :         return "EPSG:3857";
     614                 :            : 
     615                 :          1 :     CPLString osRet(pszCRS);
     616 [ +  - ][ +  - ]:          1 :     while( osRet.size() &&
         [ +  - ][ +  - ]
         [ -  + ][ -  + ]
     617 [ +  - ][ +  - ]:          1 :            (osRet.back() == ' ' || osRet.back() == '\r' || osRet.back() == '\n') )
                 [ +  - ]
     618                 :            :     {
     619 [ #  # ][ #  # ]:          0 :         osRet.resize(osRet.size() - 1);
     620                 :            :     }
     621         [ +  - ]:          1 :     return osRet;
     622                 :            : }
     623                 :            : 
     624                 :            : /************************************************************************/
     625                 :            : /*                              ReadTMS()                               */
     626                 :            : /************************************************************************/
     627                 :            : 
     628                 :          1 : int WMTSDataset::ReadTMS(CPLXMLNode* psContents,
     629                 :            :                          const CPLString& osIdentifier,
     630                 :            :                          const CPLString& osMaxTileMatrixIdentifier,
     631                 :            :                          int nMaxZoomLevel,
     632                 :            :                          WMTSTileMatrixSet& oTMS)
     633                 :            : {
     634         [ +  + ]:          1 :     for(CPLXMLNode* psIter = psContents->psChild; psIter != NULL; psIter = psIter->psNext )
     635                 :            :     {
     636 [ +  - ][ +  + ]:          1 :         if( psIter->eType != CXT_Element || strcmp(psIter->pszValue, "TileMatrixSet") != 0 )
     637                 :          1 :             continue;
     638                 :          1 :         const char* pszIdentifier = CPLGetXMLValue(psIter, "Identifier", "");
     639         [ -  + ]:          1 :         if( !EQUAL(osIdentifier, pszIdentifier) )
     640                 :          0 :             continue;
     641                 :          1 :         const char* pszSupportedCRS = CPLGetXMLValue(psIter, "SupportedCRS", NULL);
     642         [ +  + ]:          1 :         if( pszSupportedCRS == NULL )
     643                 :            :         {
     644                 :          1 :             CPLError(CE_Failure, CPLE_AppDefined, "Missing SupportedCRS");
     645                 :          1 :             return FALSE;
     646                 :            :         }
     647         [ +  - ]:          1 :         oTMS.osSRS = pszSupportedCRS;
     648 [ +  - ][ +  - ]:          1 :         if( oTMS.oSRS.SetFromUserInput(FixCRSName(pszSupportedCRS)) != OGRERR_NONE )
                 [ -  + ]
     649                 :            :         {
     650                 :            :             CPLError(CE_Failure, CPLE_AppDefined, "Cannot parse CRS '%s'",
     651                 :          0 :                      pszSupportedCRS);
     652                 :          0 :             return FALSE;
     653                 :            :         }
     654 [ +  + ][ -  + ]:          1 :         int bSwap = oTMS.oSRS.EPSGTreatsAsLatLong() || oTMS.oSRS.EPSGTreatsAsNorthingEasting();
     655                 :          1 :         CPLXMLNode* psBB = CPLGetXMLNode(psIter, "BoundingBox");
     656                 :          1 :         oTMS.bBoundingBoxValid = false;
     657         [ +  + ]:          1 :         if( psBB != NULL )
     658                 :            :         {
     659                 :          1 :             CPLString osCRS = CPLGetXMLValue(psBB, "crs", "");
     660 [ +  - ][ +  - ]:          1 :             if( EQUAL(osCRS, "") || EQUAL(osCRS, pszSupportedCRS) )
         [ +  - ][ +  - ]
                 [ +  - ]
     661                 :            :             {
     662 [ +  - ][ +  - ]:          1 :                 CPLString osLowerCorner = CPLGetXMLValue(psBB, "LowerCorner", "");
     663 [ +  - ][ +  - ]:          1 :                 CPLString osUpperCorner = CPLGetXMLValue(psBB, "UpperCorner", "");
     664 [ +  - ][ +  - ]:          1 :                 if( !osLowerCorner.empty() && !osUpperCorner.empty() )
         [ +  - ][ +  - ]
                 [ +  - ]
     665                 :            :                 {
     666 [ +  - ][ +  - ]:          1 :                     char** papszLC = CSLTokenizeString(osLowerCorner);
     667 [ +  - ][ +  - ]:          1 :                     char** papszUC = CSLTokenizeString(osUpperCorner);
     668 [ +  - ][ +  - ]:          1 :                     if( CSLCount(papszLC) == 2 && CSLCount(papszUC) == 2 )
         [ +  - ][ +  - ]
                 [ +  - ]
     669                 :            :                     {
     670 [ +  + ][ +  - ]:          1 :                         oTMS.sBoundingBox.MinX = CPLAtof(papszLC[(bSwap)? 1 : 0]);
     671 [ +  + ][ +  - ]:          1 :                         oTMS.sBoundingBox.MinY = CPLAtof(papszLC[(bSwap)? 0 : 1]);
     672 [ +  + ][ +  - ]:          1 :                         oTMS.sBoundingBox.MaxX = CPLAtof(papszUC[(bSwap)? 1 : 0]);
     673 [ +  + ][ +  - ]:          1 :                         oTMS.sBoundingBox.MaxY = CPLAtof(papszUC[(bSwap)? 0 : 1]);
     674                 :          1 :                         oTMS.bBoundingBoxValid = true;
     675                 :            :                     }
     676         [ +  - ]:          1 :                     CSLDestroy(papszLC);
     677         [ +  - ]:          1 :                     CSLDestroy(papszUC);
     678 [ +  - ][ +  - ]:          1 :                 }
     679                 :          1 :             }
     680                 :            :         }
     681                 :            :         else
     682                 :            :         {
     683                 :          1 :             const char* pszWellKnownScaleSet = CPLGetXMLValue(psIter, "WellKnownScaleSet", "");
     684 [ +  - ][ +  - ]:          1 :             if( EQUAL(pszIdentifier, "GoogleCRS84Quad") ||
         [ +  - ][ -  + ]
     685                 :          1 :                 EQUAL(pszWellKnownScaleSet, "urn:ogc:def:wkss:OGC:1.0:GoogleCRS84Quad") ||
     686                 :          1 :                 EQUAL(pszIdentifier, "GlobalCRS84Scale") ||
     687                 :          1 :                 EQUAL(pszWellKnownScaleSet, "urn:ogc:def:wkss:OGC:1.0:GlobalCRS84Scale") )
     688                 :            :             {
     689                 :          0 :                 oTMS.sBoundingBox.MinX = -180;
     690                 :          0 :                 oTMS.sBoundingBox.MinY = -90;
     691                 :          0 :                 oTMS.sBoundingBox.MaxX = 180;
     692                 :          0 :                 oTMS.sBoundingBox.MaxY = 90;
     693                 :          0 :                 oTMS.bBoundingBoxValid = true;
     694                 :            :             }
     695                 :            :         }
     696                 :            : 
     697                 :          1 :         bool bFoundTileMatrix = false;
     698         [ +  + ]:          1 :         for(CPLXMLNode* psSubIter = psIter->psChild; psSubIter != NULL; psSubIter = psSubIter->psNext )
     699                 :            :         {
     700 [ +  - ][ +  + ]:          1 :             if( psSubIter->eType != CXT_Element || strcmp(psSubIter->pszValue, "TileMatrix") != 0 )
     701                 :          1 :                 continue;
     702                 :          1 :             const char* l_pszIdentifier = CPLGetXMLValue(psSubIter, "Identifier", NULL);
     703                 :          1 :             const char* pszScaleDenominator = CPLGetXMLValue(psSubIter, "ScaleDenominator", NULL);
     704                 :          1 :             const char* pszTopLeftCorner = CPLGetXMLValue(psSubIter, "TopLeftCorner", NULL);
     705                 :          1 :             const char* pszTileWidth = CPLGetXMLValue(psSubIter, "TileWidth", NULL);
     706                 :          1 :             const char* pszTileHeight = CPLGetXMLValue(psSubIter, "TileHeight", NULL);
     707                 :          1 :             const char* pszMatrixWidth = CPLGetXMLValue(psSubIter, "MatrixWidth", NULL);
     708                 :          1 :             const char* pszMatrixHeight = CPLGetXMLValue(psSubIter, "MatrixHeight", NULL);
     709 [ +  + ][ +  - ]:          1 :             if( l_pszIdentifier == NULL || pszScaleDenominator == NULL ||
           [ +  -  +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ -  + ]
     710                 :          1 :                 pszTopLeftCorner == NULL || strchr(pszTopLeftCorner, ' ') == NULL ||
     711                 :            :                 pszTileWidth == NULL || pszTileHeight == NULL ||
     712                 :            :                 pszMatrixWidth == NULL || pszMatrixHeight == NULL )
     713                 :            :             {
     714                 :            :                 CPLError(CE_Failure, CPLE_AppDefined,
     715                 :          1 :                          "Missing required element in TileMatrix element");
     716                 :          1 :                 return FALSE;
     717                 :            :             }
     718                 :          1 :             WMTSTileMatrix oTM;
     719 [ +  - ][ +  - ]:          1 :             oTM.osIdentifier = l_pszIdentifier;
                 [ +  - ]
     720         [ +  - ]:          1 :             oTM.dfScaleDenominator = CPLAtof(pszScaleDenominator);
     721                 :          1 :             oTM.dfPixelSize = oTM.dfScaleDenominator * WMTS_PITCH;
     722 [ +  - ][ +  + ]:          1 :             if( oTMS.oSRS.IsGeographic() )
     723                 :          1 :                 oTM.dfPixelSize *= WMTS_WGS84_DEG_PER_METER;
     724         [ +  - ]:          1 :             double dfVal1 = CPLAtof(pszTopLeftCorner);
     725         [ +  - ]:          1 :             double dfVal2 = CPLAtof(strchr(pszTopLeftCorner, ' ')+1);
     726 [ +  + ][ -  + ]:          1 :             if( !bSwap ||
                 [ #  # ]
     727                 :            :                 /* Hack for http://osm.geobretagne.fr/gwc01/service/wmts?request=getcapabilities */
     728                 :          1 :                 ( STARTS_WITH_CI(l_pszIdentifier, "EPSG:4326:") &&
     729                 :            :                   dfVal1 == -180.0 ) )
     730                 :            :             {
     731                 :          1 :                 oTM.dfTLX = dfVal1;
     732                 :          1 :                 oTM.dfTLY = dfVal2;
     733                 :            :             }
     734                 :            :             else
     735                 :            :             {
     736                 :          1 :                 oTM.dfTLX = dfVal2;
     737                 :          1 :                 oTM.dfTLY = dfVal1;
     738                 :            :             }
     739                 :          1 :             oTM.nTileWidth = atoi(pszTileWidth);
     740                 :          1 :             oTM.nTileHeight = atoi(pszTileHeight);
     741 [ +  - ][ +  - ]:          1 :             if( oTM.nTileWidth <= 0 || oTM.nTileWidth > 4096 ||
         [ +  - ][ -  + ]
     742                 :            :                 oTM.nTileHeight <= 0 || oTM.nTileHeight > 4096 )
     743                 :            :             {
     744                 :            :                 CPLError(CE_Failure, CPLE_AppDefined,
     745         [ #  # ]:          0 :                          "Invalid TileWidth/TileHeight element");
     746                 :          0 :                 return FALSE;
     747                 :            :             }
     748                 :          1 :             oTM.nMatrixWidth = atoi(pszMatrixWidth);
     749                 :          1 :             oTM.nMatrixHeight = atoi(pszMatrixHeight);
     750                 :            :             // http://datacarto.geonormandie.fr/mapcache/wmts?SERVICE=WMTS&REQUEST=GetCapabilities
     751                 :            :             // has a TileMatrix 0 with MatrixWidth = MatrixHeight = 0
     752 [ +  - ][ -  + ]:          1 :             if( oTM.nMatrixWidth < 1 || oTM.nMatrixHeight < 1 )
     753                 :          0 :                 continue;
     754         [ +  - ]:          1 :             oTMS.aoTM.push_back(oTM);
     755 [ +  + ][ +  - ]:          1 :             if( (nMaxZoomLevel >= 0 && static_cast<int>(oTMS.aoTM.size())-1
         [ +  + ][ +  + ]
         [ +  + ][ +  + ]
     756                 :            :                                                         == nMaxZoomLevel) ||
     757         [ +  - ]:          1 :                 (!osMaxTileMatrixIdentifier.empty() &&
     758         [ +  - ]:          1 :                  EQUAL(osMaxTileMatrixIdentifier, l_pszIdentifier)) )
     759                 :            :             {
     760                 :          1 :                 bFoundTileMatrix = true;
     761                 :            :                 break;
     762                 :            :             }
     763      [ +  +  - ]:          1 :         }
     764 [ +  + ][ +  + ]:          1 :         if( nMaxZoomLevel >= 0 && !bFoundTileMatrix )
     765                 :            :         {
     766                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
     767                 :            :                      "Cannot find TileMatrix of zoom level %d in TileMatrixSet '%s'",
     768                 :            :                      nMaxZoomLevel,
     769                 :          1 :                      osIdentifier.c_str());
     770                 :          1 :             return FALSE;
     771                 :            :         }
     772 [ +  + ][ +  + ]:          1 :         if( !osMaxTileMatrixIdentifier.empty() && !bFoundTileMatrix )
                 [ +  + ]
     773                 :            :         {
     774                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
     775                 :            :                      "Cannot find TileMatrix '%s' in TileMatrixSet '%s'",
     776                 :            :                      osMaxTileMatrixIdentifier.c_str(),
     777                 :          1 :                      osIdentifier.c_str());
     778                 :          1 :             return FALSE;
     779                 :            :         }
     780         [ +  + ]:          1 :         if( oTMS.aoTM.empty() )
     781                 :            :         {
     782                 :            :             CPLError(CE_Failure, CPLE_AppDefined, "Cannot find TileMatrix in TileMatrixSet '%s'",
     783                 :          1 :                      osIdentifier.c_str());
     784                 :          1 :             return FALSE;
     785                 :            :         }
     786                 :          1 :         return TRUE;
     787                 :            :     }
     788                 :            :     CPLError(CE_Failure, CPLE_AppDefined, "Cannot find TileMatrixSet '%s'",
     789                 :          1 :              osIdentifier.c_str());
     790                 :          1 :     return FALSE;
     791                 :            : }
     792                 :            : 
     793                 :            : /************************************************************************/
     794                 :            : /*                              ReadTMLimits()                          */
     795                 :            : /************************************************************************/
     796                 :            : 
     797                 :          1 : int WMTSDataset::ReadTMLimits(CPLXMLNode* psTMSLimits,
     798                 :            :                               std::map<CPLString, WMTSTileMatrixLimits>& aoMapTileMatrixLimits)
     799                 :            : {
     800         [ +  + ]:          1 :     for(CPLXMLNode* psIter = psTMSLimits->psChild; psIter; psIter = psIter->psNext)
     801                 :            :     {
     802 [ +  - ][ -  + ]:          1 :         if( psIter->eType != CXT_Element || strcmp(psIter->pszValue, "TileMatrixLimits") != 0 )
     803                 :          0 :             continue;
     804                 :          1 :         WMTSTileMatrixLimits oTMLimits;
     805         [ +  - ]:          1 :         const char* pszTileMatrix = CPLGetXMLValue(psIter, "TileMatrix", NULL);
     806         [ +  - ]:          1 :         const char* pszMinTileRow = CPLGetXMLValue(psIter, "MinTileRow", NULL);
     807         [ +  - ]:          1 :         const char* pszMaxTileRow = CPLGetXMLValue(psIter, "MaxTileRow", NULL);
     808         [ +  - ]:          1 :         const char* pszMinTileCol = CPLGetXMLValue(psIter, "MinTileCol", NULL);
     809         [ +  - ]:          1 :         const char* pszMaxTileCol = CPLGetXMLValue(psIter, "MaxTileCol", NULL);
     810 [ +  - ][ +  - ]:          1 :         if( pszTileMatrix == NULL ||
         [ +  - ][ +  - ]
                 [ -  + ]
     811                 :            :             pszMinTileRow == NULL || pszMaxTileRow == NULL ||
     812                 :            :             pszMinTileCol == NULL || pszMaxTileCol == NULL )
     813                 :            :         {
     814                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
     815         [ #  # ]:          0 :                      "Missing required element in TileMatrixLimits element");
     816                 :          0 :             return FALSE;
     817                 :            :         }
     818 [ +  - ][ +  - ]:          1 :         oTMLimits.osIdentifier = pszTileMatrix;
                 [ +  - ]
     819                 :          1 :         oTMLimits.nMinTileRow = atoi(pszMinTileRow);
     820                 :          1 :         oTMLimits.nMaxTileRow = atoi(pszMaxTileRow);
     821                 :          1 :         oTMLimits.nMinTileCol = atoi(pszMinTileCol);
     822                 :          1 :         oTMLimits.nMaxTileCol = atoi(pszMaxTileCol);
     823 [ +  - ][ +  - ]:          1 :         aoMapTileMatrixLimits[pszTileMatrix] = oTMLimits;
         [ +  - ][ +  - ]
                 [ +  - ]
     824                 :          1 :     }
     825                 :          1 :     return TRUE;
     826                 :            : }
     827                 :            : 
     828                 :            : /************************************************************************/
     829                 :            : /*                               Replace()                              */
     830                 :            : /************************************************************************/
     831                 :            : 
     832                 :          1 : CPLString WMTSDataset::Replace(const CPLString& osStr, const char* pszOld,
     833                 :            :                                const char* pszNew)
     834                 :            : {
     835                 :          1 :     size_t nPos = osStr.ifind(pszOld);
     836         [ +  + ]:          1 :     if( nPos == std::string::npos )
     837                 :          1 :         return osStr;
     838 [ +  - ][ +  - ]:          1 :     CPLString osRet(osStr.substr(0, nPos));
     839         [ +  - ]:          1 :     osRet += pszNew;
     840 [ +  - ][ +  - ]:          1 :     osRet += osStr.substr(nPos + strlen(pszOld));
                 [ +  - ]
     841         [ +  - ]:          1 :     return osRet;
     842                 :            : }
     843                 :            : 
     844                 :            : /************************************************************************/
     845                 :            : /*                       GetCapabilitiesResponse()                      */
     846                 :            : /************************************************************************/
     847                 :            : 
     848                 :          1 : CPLXMLNode* WMTSDataset::GetCapabilitiesResponse(const CPLString& osFilename,
     849                 :            :                                                  char** papszHTTPOptions)
     850                 :            : {
     851                 :            :     CPLXMLNode* psXML;
     852                 :            :     VSIStatBufL sStat;
     853         [ +  + ]:          1 :     if( VSIStatL(osFilename, &sStat) == 0 )
     854                 :          1 :         psXML = CPLParseXMLFile(osFilename);
     855                 :            :     else
     856                 :            :     {
     857                 :          1 :         CPLHTTPResult* psResult = CPLHTTPFetch(osFilename, papszHTTPOptions);
     858         [ -  + ]:          1 :         if( psResult == NULL )
     859                 :          0 :             return NULL;
     860         [ +  - ]:          1 :         if( psResult->pabyData == NULL )
     861                 :            :         {
     862                 :          1 :             CPLHTTPDestroyResult(psResult);
     863                 :          1 :             return NULL;
     864                 :            :         }
     865                 :          0 :         psXML = CPLParseXMLString((const char*)psResult->pabyData);
     866                 :          0 :         CPLHTTPDestroyResult(psResult);
     867                 :            :     }
     868                 :          1 :     return psXML;
     869                 :            : }
     870                 :            : 
     871                 :            : /************************************************************************/
     872                 :            : /*                          WMTSAddOtherXML()                           */
     873                 :            : /************************************************************************/
     874                 :            : 
     875                 :          1 : static void WMTSAddOtherXML(CPLXMLNode* psRoot, const char* pszElement,
     876                 :            :                             CPLString& osOtherXML)
     877                 :            : {
     878                 :          1 :     CPLXMLNode* psElement = CPLGetXMLNode(psRoot, pszElement);
     879         [ +  + ]:          1 :     if( psElement )
     880                 :            :     {
     881                 :          1 :         CPLXMLNode* psNext = psElement->psNext;
     882                 :          1 :         psElement->psNext = NULL;
     883                 :          1 :         char* pszTmp = CPLSerializeXMLTree(psElement);
     884                 :          1 :         osOtherXML += pszTmp;
     885                 :          1 :         CPLFree(pszTmp);
     886                 :          1 :         psElement->psNext = psNext;
     887                 :            :     }
     888                 :          1 : }
     889                 :            : 
     890                 :            : /************************************************************************/
     891                 :            : /*                          GetOperationKVPURL()                        */
     892                 :            : /************************************************************************/
     893                 :            : 
     894                 :          1 : CPLString WMTSDataset::GetOperationKVPURL(CPLXMLNode* psXML,
     895                 :            :                                           const char* pszOperation)
     896                 :            : {
     897                 :          1 :     CPLString osRet;
     898         [ +  - ]:          1 :     CPLXMLNode* psOM = CPLGetXMLNode(psXML, "=Capabilities.OperationsMetadata");
     899 [ +  + ][ +  + ]:          1 :     for(CPLXMLNode* psIter = psOM ? psOM->psChild : NULL;
     900                 :            :         psIter != NULL; psIter = psIter->psNext)
     901                 :            :     {
     902 [ +  - ][ +  - ]:          1 :         if( psIter->eType != CXT_Element ||
         [ +  + ][ +  + ]
     903                 :          1 :             strcmp(psIter->pszValue, "Operation") != 0 ||
     904         [ +  - ]:          1 :             !EQUAL(CPLGetXMLValue(psIter, "name", ""), pszOperation) )
     905                 :            :         {
     906                 :          1 :             continue;
     907                 :            :         }
     908         [ +  - ]:          1 :         CPLXMLNode* psHTTP = CPLGetXMLNode(psIter, "DCP.HTTP");
     909 [ +  - ][ +  + ]:          1 :         for(CPLXMLNode* psGet = psHTTP ? psHTTP->psChild : NULL;
     910                 :            :                 psGet != NULL; psGet = psGet->psNext)
     911                 :            :         {
     912 [ +  - ][ -  + ]:          1 :             if( psGet->eType != CXT_Element ||
     913                 :          1 :                 strcmp(psGet->pszValue, "Get") != 0 )
     914                 :            :             {
     915                 :          0 :                 continue;
     916                 :            :             }
     917 [ +  - ][ -  + ]:          1 :             if( !EQUAL(CPLGetXMLValue(psGet, "Constraint.AllowedValues.Value", "KVP"), "KVP") )
     918                 :          0 :                 continue;
     919 [ +  - ][ +  - ]:          1 :             osRet = CPLGetXMLValue(psGet, "href", "");
         [ +  - ][ +  - ]
     920                 :            :         }
     921                 :            :     }
     922                 :          1 :     return osRet;
     923                 :            : }
     924                 :            : 
     925                 :            : /************************************************************************/
     926                 :            : /*                           BuildHTTPRequestOpts()                     */
     927                 :            : /************************************************************************/
     928                 :            : 
     929                 :          1 : char** WMTSDataset::BuildHTTPRequestOpts(CPLString osOtherXML)
     930                 :            : {
     931 [ +  - ][ +  - ]:          1 :     osOtherXML = "<Root>" + osOtherXML + "</Root>";
         [ +  - ][ +  - ]
                 [ +  - ]
     932                 :          1 :     CPLXMLNode* psXML = CPLParseXMLString(osOtherXML);
     933                 :          1 :     char **http_request_opts = NULL;
     934         [ -  + ]:          1 :     if (CPLGetXMLValue(psXML, "Timeout", NULL)) {
     935                 :          0 :         CPLString optstr;
     936 [ #  # ][ #  # ]:          0 :         optstr.Printf("TIMEOUT=%s", CPLGetXMLValue(psXML, "Timeout", NULL));
     937 [ #  # ][ #  # ]:          0 :         http_request_opts = CSLAddString(http_request_opts, optstr.c_str());
     938                 :            :     }
     939         [ -  + ]:          1 :     if (CPLGetXMLValue(psXML, "UserAgent", NULL)) {
     940                 :          0 :         CPLString optstr;
     941 [ #  # ][ #  # ]:          0 :         optstr.Printf("USERAGENT=%s", CPLGetXMLValue(psXML, "UserAgent", NULL));
     942 [ #  # ][ #  # ]:          0 :         http_request_opts = CSLAddString(http_request_opts, optstr.c_str());
     943                 :            :     }
     944         [ -  + ]:          1 :     if (CPLGetXMLValue(psXML, "Referer", NULL)) {
     945                 :          0 :         CPLString optstr;
     946 [ #  # ][ #  # ]:          0 :         optstr.Printf("REFERER=%s", CPLGetXMLValue(psXML, "Referer", NULL));
     947 [ #  # ][ #  # ]:          0 :         http_request_opts = CSLAddString(http_request_opts, optstr.c_str());
     948                 :            :     }
     949         [ +  + ]:          1 :     if (CPLTestBool(CPLGetXMLValue(psXML, "UnsafeSSL", "false"))) {
     950                 :          1 :         http_request_opts = CSLAddString(http_request_opts, "UNSAFESSL=1");
     951                 :            :     }
     952         [ -  + ]:          1 :     if (CPLGetXMLValue(psXML, "UserPwd", NULL)) {
     953                 :          0 :         CPLString optstr;
     954 [ #  # ][ #  # ]:          0 :         optstr.Printf("USERPWD=%s", CPLGetXMLValue(psXML, "UserPwd", NULL));
     955 [ #  # ][ #  # ]:          0 :         http_request_opts = CSLAddString(http_request_opts, optstr.c_str());
     956                 :            :     }
     957                 :          1 :     CPLDestroyXMLNode(psXML);
     958                 :          1 :     return http_request_opts;
     959                 :            : }
     960                 :            : 
     961                 :            : /************************************************************************/
     962                 :            : /*                                Open()                                */
     963                 :            : /************************************************************************/
     964                 :            : 
     965                 :          1 : GDALDataset* WMTSDataset::Open(GDALOpenInfo* poOpenInfo)
     966                 :            : {
     967         [ +  + ]:          1 :     if (!Identify(poOpenInfo))
     968                 :          1 :         return NULL;
     969                 :            : 
     970                 :          1 :     CPLXMLNode* psXML = NULL;
     971                 :          1 :     CPLString osTileFormat;
     972         [ +  - ]:          1 :     CPLString osInfoFormat;
     973                 :            : 
     974                 :            :     CPLString osGetCapabilitiesURL = CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     975 [ +  - ][ +  - ]:          1 :                                                 "URL", "");
     976                 :            :     CPLString osLayer = CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     977 [ +  - ][ +  - ]:          1 :                                     "LAYER", "");
     978                 :            :     CPLString osTMS = CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     979 [ +  - ][ +  - ]:          1 :                                     "TILEMATRIXSET", "");
     980                 :            :     CPLString osMaxTileMatrixIdentifier = CSLFetchNameValueDef(
     981                 :            :                                     poOpenInfo->papszOpenOptions,
     982 [ +  - ][ +  - ]:          1 :                                     "TILEMATRIX", "");
     983                 :            :     int nUserMaxZoomLevel = atoi(CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     984                 :            :                                     "ZOOM_LEVEL",
     985                 :            :                                     CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     986 [ +  - ][ +  - ]:          1 :                                     "ZOOMLEVEL", "-1")));
     987                 :            :     CPLString osStyle = CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     988 [ +  - ][ +  - ]:          1 :                                     "STYLE", "");
     989                 :            : 
     990                 :            :     int bExtendBeyondDateLine =
     991                 :            :         CPLFetchBool(poOpenInfo->papszOpenOptions,
     992         [ +  - ]:          1 :                      "EXTENDBEYONDDATELINE", false);
     993                 :            : 
     994                 :            :     CPLString osOtherXML = "<Cache />"
     995                 :            :                      "<UnsafeSSL>true</UnsafeSSL>"
     996                 :            :                      "<ZeroBlockHttpCodes>204,404</ZeroBlockHttpCodes>"
     997         [ +  - ]:          1 :                      "<ZeroBlockOnServerException>true</ZeroBlockOnServerException>";
     998                 :            : 
     999         [ +  + ]:          1 :     if( STARTS_WITH_CI(poOpenInfo->pszFilename, "WMTS:") )
    1000                 :            :     {
    1001                 :            :         char** papszTokens = CSLTokenizeString2( poOpenInfo->pszFilename + 5,
    1002         [ +  - ]:          1 :                                                  ",", CSLT_HONOURSTRINGS );
    1003 [ +  - ][ +  + ]:          1 :         if( papszTokens && papszTokens[0] )
    1004                 :            :         {
    1005 [ +  - ][ +  - ]:          1 :             osGetCapabilitiesURL = papszTokens[0];
                 [ +  - ]
    1006         [ +  + ]:          1 :             for(char** papszIter = papszTokens+1; *papszIter; papszIter++)
    1007                 :            :             {
    1008                 :          1 :                 char* pszKey = NULL;
    1009         [ +  - ]:          1 :                 const char* pszValue = CPLParseNameValue(*papszIter, &pszKey);
    1010 [ +  - ][ +  - ]:          1 :                 if( pszKey && pszValue )
    1011                 :            :                 {
    1012         [ +  + ]:          1 :                     if( EQUAL(pszKey, "layer") )
    1013 [ +  - ][ +  - ]:          1 :                         osLayer = pszValue;
                 [ +  - ]
    1014         [ +  + ]:          1 :                     else if( EQUAL(pszKey, "tilematrixset") )
    1015 [ +  - ][ +  - ]:          1 :                         osTMS = pszValue;
                 [ +  - ]
    1016         [ +  + ]:          1 :                     else if( EQUAL(pszKey, "tilematrix") )
    1017 [ +  - ][ +  - ]:          1 :                         osMaxTileMatrixIdentifier = pszValue;
                 [ +  - ]
    1018 [ +  + ][ -  + ]:          1 :                     else if( EQUAL(pszKey, "zoom_level") ||
    1019                 :          1 :                              EQUAL(pszKey, "zoomlevel") )
    1020                 :          1 :                         nUserMaxZoomLevel = atoi(pszValue);
    1021         [ +  + ]:          1 :                     else if( EQUAL(pszKey, "style") )
    1022 [ +  - ][ +  - ]:          1 :                         osStyle = pszValue;
                 [ +  - ]
    1023         [ +  - ]:          1 :                     else if( EQUAL(pszKey, "extendbeyonddateline") )
    1024         [ +  - ]:          1 :                         bExtendBeyondDateLine = CPLTestBool(pszValue);
    1025                 :            :                     else
    1026                 :            :                         CPLError(CE_Warning, CPLE_AppDefined,
    1027         [ #  # ]:          1 :                                  "Unknown parameter: %s'", pszKey);
    1028                 :            :                 }
    1029         [ +  - ]:          1 :                 CPLFree(pszKey);
    1030                 :            :             }
    1031                 :            :         }
    1032         [ +  - ]:          1 :         CSLDestroy(papszTokens);
    1033                 :            : 
    1034 [ +  - ][ +  - ]:          1 :         char** papszHTTPOptions = BuildHTTPRequestOpts(osOtherXML);
                 [ +  - ]
    1035         [ +  - ]:          1 :         psXML = GetCapabilitiesResponse(osGetCapabilitiesURL, papszHTTPOptions);
    1036         [ +  - ]:          1 :         CSLDestroy(papszHTTPOptions);
    1037                 :            :     }
    1038                 :            : 
    1039                 :          1 :     int bHasAOI = FALSE;
    1040         [ +  - ]:          1 :     OGREnvelope sAOI;
    1041                 :          1 :     int nBands = 4;
    1042         [ +  - ]:          1 :     CPLString osProjection;
    1043                 :            : 
    1044 [ +  + ][ +  - ]:          1 :     if( (psXML != NULL && CPLGetXMLNode(psXML, "=GDAL_WMTS") != NULL ) ||
         [ +  + ][ +  + ]
           [ +  +  +  + ]
                 [ +  + ]
    1045                 :          1 :         STARTS_WITH_CI(poOpenInfo->pszFilename, "<GDAL_WMTS") ||
    1046                 :            :         (poOpenInfo->nHeaderBytes > 0 &&
    1047                 :          1 :          strstr((const char*)poOpenInfo->pabyHeader, "<GDAL_WMTS")) )
    1048                 :            :     {
    1049                 :            :         CPLXMLNode* psGDALWMTS;
    1050 [ +  + ][ +  - ]:          1 :         if( psXML != NULL && CPLGetXMLNode(psXML, "=GDAL_WMTS") != NULL )
         [ +  - ][ +  + ]
    1051         [ +  - ]:          1 :             psGDALWMTS = CPLCloneXMLTree(psXML);
    1052         [ +  + ]:          1 :         else if( STARTS_WITH_CI(poOpenInfo->pszFilename, "<GDAL_WMTS") )
    1053         [ +  - ]:          1 :             psGDALWMTS = CPLParseXMLString(poOpenInfo->pszFilename);
    1054                 :            :         else
    1055         [ +  - ]:          1 :             psGDALWMTS = CPLParseXMLFile(poOpenInfo->pszFilename);
    1056         [ +  + ]:          1 :         if( psGDALWMTS == NULL )
    1057                 :          1 :             return NULL;
    1058         [ +  - ]:          1 :         CPLXMLNode* psRoot = CPLGetXMLNode(psGDALWMTS, "=GDAL_WMTS");
    1059         [ +  + ]:          1 :         if( psRoot == NULL )
    1060                 :            :         {
    1061         [ +  - ]:          1 :             CPLError(CE_Failure, CPLE_AppDefined, "Cannot find root <GDAL_WMTS>");
    1062         [ +  - ]:          1 :             CPLDestroyXMLNode(psGDALWMTS);
    1063                 :          1 :             return NULL;
    1064                 :            :         }
    1065 [ +  - ][ +  - ]:          1 :         osGetCapabilitiesURL = CPLGetXMLValue(psRoot, "GetCapabilitiesUrl", "");
         [ +  - ][ +  - ]
    1066 [ +  - ][ +  + ]:          1 :         if( osGetCapabilitiesURL.empty() )
    1067                 :            :         {
    1068         [ +  - ]:          1 :             CPLError(CE_Failure, CPLE_AppDefined, "Missing <GetCapabilitiesUrl>");
    1069         [ +  - ]:          1 :             CPLDestroyXMLNode(psGDALWMTS);
    1070                 :          1 :             return NULL;
    1071                 :            :         }
    1072                 :            : 
    1073 [ +  - ][ +  - ]:          1 :         osLayer = CPLGetXMLValue(psRoot, "Layer", osLayer);
         [ +  - ][ +  - ]
                 [ +  - ]
    1074 [ +  - ][ +  - ]:          1 :         osTMS = CPLGetXMLValue(psRoot, "TileMatrixSet", osTMS);
         [ +  - ][ +  - ]
                 [ +  - ]
    1075                 :            :         osMaxTileMatrixIdentifier = CPLGetXMLValue(psRoot, "TileMatrix",
    1076 [ +  - ][ +  - ]:          1 :                                                    osMaxTileMatrixIdentifier);
         [ +  - ][ +  - ]
                 [ +  - ]
    1077                 :            :         nUserMaxZoomLevel = atoi(CPLGetXMLValue(psRoot, "ZoomLevel",
    1078 [ +  - ][ +  - ]:          1 :                                        CPLSPrintf("%d", nUserMaxZoomLevel)));
    1079 [ +  - ][ +  - ]:          1 :         osStyle = CPLGetXMLValue(psRoot, "Style", osStyle);
         [ +  - ][ +  - ]
                 [ +  - ]
    1080 [ +  - ][ +  - ]:          1 :         osTileFormat = CPLGetXMLValue(psRoot, "Format", osTileFormat);
         [ +  - ][ +  - ]
                 [ +  - ]
    1081 [ +  - ][ +  - ]:          1 :         osInfoFormat = CPLGetXMLValue(psRoot, "InfoFormat", osInfoFormat);
         [ +  - ][ +  - ]
                 [ +  - ]
    1082 [ +  - ][ +  - ]:          1 :         osProjection = CPLGetXMLValue(psRoot, "Projection", osProjection);
         [ +  - ][ +  - ]
                 [ +  - ]
    1083                 :            :         bExtendBeyondDateLine = CPLTestBool(CPLGetXMLValue(psRoot, "ExtendBeyondDateLine",
    1084 [ -  + ][ +  - ]:          1 :                                             (bExtendBeyondDateLine) ? "true": "false"));
                 [ +  - ]
    1085                 :            : 
    1086 [ +  - ][ +  - ]:          1 :         osOtherXML = "";
                 [ +  - ]
    1087         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "Cache", osOtherXML);
    1088         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "MaxConnections", osOtherXML);
    1089         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "Timeout", osOtherXML);
    1090         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "OfflineMode", osOtherXML);
    1091         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "MaxConnections", osOtherXML);
    1092         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "UserAgent", osOtherXML);
    1093         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "UserPwd", osOtherXML);
    1094         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "UnsafeSSL", osOtherXML);
    1095         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "Referer", osOtherXML);
    1096         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "ZeroBlockHttpCodes", osOtherXML);
    1097         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "ZeroBlockOnServerException", osOtherXML);
    1098                 :            : 
    1099         [ +  - ]:          1 :         nBands = atoi(CPLGetXMLValue(psRoot, "BandsCount", "4"));
    1100                 :            : 
    1101         [ +  - ]:          1 :         const char* pszULX = CPLGetXMLValue(psRoot, "DataWindow.UpperLeftX", NULL);
    1102         [ +  - ]:          1 :         const char* pszULY = CPLGetXMLValue(psRoot, "DataWindow.UpperLeftY", NULL);
    1103         [ +  - ]:          1 :         const char* pszLRX = CPLGetXMLValue(psRoot, "DataWindow.LowerRightX", NULL);
    1104         [ +  - ]:          1 :         const char* pszLRY = CPLGetXMLValue(psRoot, "DataWindow.LowerRightY", NULL);
    1105 [ +  + ][ +  - ]:          1 :         if( pszULX && pszULY && pszLRX && pszLRY )
         [ +  - ][ +  - ]
    1106                 :            :         {
    1107         [ +  - ]:          1 :             sAOI.MinX = CPLAtof(pszULX);
    1108         [ +  - ]:          1 :             sAOI.MaxY = CPLAtof(pszULY);
    1109         [ +  - ]:          1 :             sAOI.MaxX = CPLAtof(pszLRX);
    1110         [ +  - ]:          1 :             sAOI.MinY = CPLAtof(pszLRY);
    1111                 :          1 :             bHasAOI = TRUE;
    1112                 :            :         }
    1113                 :            : 
    1114         [ +  - ]:          1 :         CPLDestroyXMLNode(psGDALWMTS);
    1115                 :            : 
    1116         [ +  - ]:          1 :         CPLDestroyXMLNode(psXML);
    1117 [ +  - ][ +  - ]:          1 :         char** papszHTTPOptions = BuildHTTPRequestOpts(osOtherXML);
                 [ +  - ]
    1118         [ +  - ]:          1 :         psXML = GetCapabilitiesResponse(osGetCapabilitiesURL, papszHTTPOptions);
    1119         [ +  - ]:          1 :         CSLDestroy(papszHTTPOptions);
    1120                 :            :     }
    1121         [ +  + ]:          1 :     else if( !STARTS_WITH_CI(poOpenInfo->pszFilename, "WMTS:") )
    1122                 :            :     {
    1123 [ +  - ][ +  - ]:          1 :         osGetCapabilitiesURL = poOpenInfo->pszFilename;
                 [ +  - ]
    1124         [ +  - ]:          1 :         psXML = CPLParseXMLFile(poOpenInfo->pszFilename);
    1125                 :            :     }
    1126         [ +  + ]:          1 :     if( psXML == NULL )
    1127                 :          1 :         return NULL;
    1128         [ +  - ]:          1 :     CPLStripXMLNamespace(psXML, NULL, TRUE);
    1129                 :            : 
    1130         [ +  - ]:          1 :     CPLXMLNode* psContents = CPLGetXMLNode(psXML, "=Capabilities.Contents");
    1131         [ +  + ]:          1 :     if( psContents == NULL )
    1132                 :            :     {
    1133         [ +  - ]:          1 :         CPLError(CE_Failure, CPLE_AppDefined, "Missing Capabilities.Contents element");
    1134         [ +  - ]:          1 :         CPLDestroyXMLNode(psXML);
    1135                 :          1 :         return NULL;
    1136                 :            :     }
    1137                 :            : 
    1138 [ +  - ][ +  - ]:          1 :     if( STARTS_WITH(osGetCapabilitiesURL, "/vsimem/") )
    1139                 :            :     {
    1140 [ +  - ][ +  + ]:          1 :         if( CPLGetXMLValue(psXML, "=Capabilities.ServiceMetadataURL.href", NULL) )
    1141 [ +  - ][ +  - ]:          1 :             osGetCapabilitiesURL = CPLGetXMLValue(psXML, "=Capabilities.ServiceMetadataURL.href", NULL);
         [ +  - ][ +  - ]
    1142                 :            :         else
    1143                 :            :         {
    1144 [ +  - ][ +  - ]:          1 :             osGetCapabilitiesURL = GetOperationKVPURL(psXML, "GetCapabilities");
                 [ +  - ]
    1145 [ +  - ][ +  + ]:          1 :             if( !osGetCapabilitiesURL.empty() )
    1146                 :            :             {
    1147 [ +  - ][ +  - ]:          1 :                 osGetCapabilitiesURL = CPLURLAddKVP(osGetCapabilitiesURL, "service", "WMTS");
         [ +  - ][ +  - ]
    1148 [ +  - ][ +  - ]:          1 :                 osGetCapabilitiesURL = CPLURLAddKVP(osGetCapabilitiesURL, "request", "GetCapabilities");
         [ +  - ][ +  - ]
    1149                 :            :             }
    1150                 :            :         }
    1151                 :            :     }
    1152         [ +  - ]:          1 :     CPLString osCapabilitiesFilename(osGetCapabilitiesURL);
    1153 [ +  - ][ +  - ]:          1 :     if( !STARTS_WITH_CI(osCapabilitiesFilename, "WMTS:") )
    1154 [ +  - ][ +  - ]:          1 :         osCapabilitiesFilename = "WMTS:" + osGetCapabilitiesURL;
         [ +  - ][ +  - ]
                 [ +  - ]
    1155                 :            : 
    1156                 :          1 :     int nLayerCount = 0;
    1157         [ +  - ]:          1 :     CPLStringList aosSubDatasets;
    1158 [ +  - ][ +  - ]:          1 :     CPLString osSelectLayer(osLayer), osSelectTMS(osTMS), osSelectStyle(osStyle);
                 [ +  - ]
    1159 [ +  - ][ +  - ]:          1 :     CPLString osSelectLayerTitle, osSelectLayerAbstract;
    1160 [ +  - ][ +  - ]:          1 :     CPLString osSelectTileFormat(osTileFormat), osSelectInfoFormat(osInfoFormat);
    1161                 :          1 :     int nCountTileFormat = 0;
    1162                 :          1 :     int nCountInfoFormat = 0;
    1163         [ +  - ]:          1 :     CPLString osURLTileTemplate;
    1164         [ +  - ]:          1 :     CPLString osURLFeatureInfoTemplate;
    1165         [ +  - ]:          1 :     std::set<CPLString> aoSetLayers;
    1166         [ +  - ]:          1 :     std::map<CPLString, OGREnvelope> aoMapBoundingBox;
    1167         [ +  - ]:          1 :     std::map<CPLString, WMTSTileMatrixLimits> aoMapTileMatrixLimits;
    1168         [ +  - ]:          1 :     std::map<CPLString, CPLString> aoMapDimensions;
    1169                 :            : 
    1170         [ +  + ]:          1 :     for(CPLXMLNode* psIter = psContents->psChild; psIter != NULL; psIter = psIter->psNext )
    1171                 :            :     {
    1172 [ +  - ][ +  + ]:          1 :         if( psIter->eType != CXT_Element || strcmp(psIter->pszValue, "Layer") != 0 )
    1173                 :          1 :             continue;
    1174         [ +  - ]:          1 :         const char* pszIdentifier = CPLGetXMLValue(psIter, "Identifier", "");
    1175 [ +  - ][ +  - ]:          1 :         if( aoSetLayers.find(pszIdentifier) != aoSetLayers.end() )
         [ +  - ][ +  - ]
         [ +  - ][ -  + ]
    1176                 :            :         {
    1177                 :            :             CPLError(CE_Warning, CPLE_AppDefined,
    1178                 :            :                      "Several layers with identifier '%s'. Only first one kept",
    1179         [ #  # ]:          0 :                      pszIdentifier);
    1180                 :            :         }
    1181 [ +  - ][ +  - ]:          1 :         aoSetLayers.insert(pszIdentifier);
                 [ +  - ]
    1182 [ +  - ][ +  + ]:          1 :         if( !osLayer.empty() && strcmp(osLayer, pszIdentifier) != 0 )
         [ +  - ][ +  + ]
                 [ +  + ]
    1183                 :          1 :             continue;
    1184         [ +  - ]:          1 :         const char* pszTitle = CPLGetXMLValue(psIter, "Title", NULL);
    1185 [ +  - ][ +  + ]:          1 :         if( osSelectLayer.empty() )
    1186                 :            :         {
    1187 [ +  - ][ +  - ]:          1 :             osSelectLayer = pszIdentifier;
                 [ +  - ]
    1188                 :            :         }
    1189 [ +  - ][ +  - ]:          1 :         if( strcmp(osSelectLayer, pszIdentifier) == 0 )
    1190                 :            :         {
    1191         [ +  + ]:          1 :             if( pszTitle != NULL )
    1192 [ +  - ][ +  - ]:          1 :                 osSelectLayerTitle = pszTitle;
                 [ +  - ]
    1193         [ +  - ]:          1 :             const char* pszAbstract = CPLGetXMLValue(psIter, "Abstract", NULL);
    1194         [ +  + ]:          1 :             if( pszAbstract != NULL )
    1195 [ +  - ][ +  - ]:          1 :                 osSelectLayerAbstract = pszAbstract;
                 [ +  - ]
    1196                 :            :         }
    1197                 :            : 
    1198         [ +  - ]:          1 :         std::vector<CPLString> aosTMS;
    1199         [ +  - ]:          1 :         std::vector<CPLString> aosStylesIdentifier;
    1200         [ +  - ]:          1 :         std::vector<CPLString> aosStylesTitle;
    1201                 :            : 
    1202                 :          1 :         CPLXMLNode* psSubIter = psIter->psChild;
    1203         [ +  + ]:          1 :         for(; psSubIter != NULL; psSubIter = psSubIter->psNext )
    1204                 :            :         {
    1205         [ +  + ]:          1 :             if( psSubIter->eType != CXT_Element )
    1206                 :          1 :                 continue;
    1207 [ +  - ][ +  - ]:          1 :             if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1208                 :          1 :                 strcmp(psSubIter->pszValue, "Format") == 0 )
    1209                 :            :             {
    1210         [ +  - ]:          1 :                 const char* pszValue = CPLGetXMLValue(psSubIter, "", "");
    1211 [ +  - ][ +  + ]:          1 :                 if( !osTileFormat.empty() && strcmp(osTileFormat, pszValue) != 0 )
         [ +  - ][ +  + ]
                 [ +  + ]
    1212                 :          1 :                     continue;
    1213                 :          1 :                 nCountTileFormat ++;
    1214 [ +  - ][ +  + ]:          1 :                 if( osSelectTileFormat.empty() ||
         [ +  - ][ +  - ]
    1215                 :          1 :                     EQUAL(pszValue, "image/png") )
    1216                 :            :                 {
    1217 [ +  - ][ +  - ]:          1 :                     osSelectTileFormat = pszValue;
                 [ +  - ]
    1218                 :            :                 }
    1219                 :            :             }
    1220 [ +  - ][ +  - ]:          1 :             else if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1221                 :          1 :                      strcmp(psSubIter->pszValue, "InfoFormat") == 0 )
    1222                 :            :             {
    1223         [ +  - ]:          1 :                 const char* pszValue = CPLGetXMLValue(psSubIter, "", "");
    1224 [ +  - ][ -  + ]:          1 :                 if( !osInfoFormat.empty() && strcmp(osInfoFormat, pszValue) != 0 )
         [ #  # ][ #  # ]
                 [ -  + ]
    1225                 :          0 :                     continue;
    1226                 :          1 :                 nCountInfoFormat ++;
    1227 [ +  - ][ -  + ]:          1 :                 if( osSelectInfoFormat.empty() ||
         [ #  # ][ #  # ]
         [ #  # ][ +  - ]
    1228                 :          0 :                     (EQUAL(pszValue, "application/vnd.ogc.gml") &&
    1229         [ #  # ]:          0 :                      !EQUAL(osSelectInfoFormat, "application/vnd.ogc.gml/3.1.1")) ||
    1230                 :          0 :                     EQUAL(pszValue, "application/vnd.ogc.gml/3.1.1") )
    1231                 :            :                 {
    1232 [ +  - ][ +  - ]:          1 :                     osSelectInfoFormat = pszValue;
                 [ +  - ]
    1233                 :            :                 }
    1234                 :            :             }
    1235 [ +  - ][ +  - ]:          1 :             else if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1236                 :          1 :                      strcmp(psSubIter->pszValue, "Dimension") == 0 )
    1237                 :            :             {
    1238                 :            :                 /* Cf http://wmts.geo.admin.ch/1.0.0/WMTSCapabilities.xml */
    1239         [ +  - ]:          1 :                 const char* pszDimensionIdentifier = CPLGetXMLValue(psSubIter, "Identifier", NULL);
    1240         [ +  - ]:          1 :                 const char* pszDefault = CPLGetXMLValue(psSubIter, "Default", "");
    1241         [ +  - ]:          1 :                 if( pszDimensionIdentifier != NULL )
    1242 [ +  - ][ +  - ]:          1 :                     aoMapDimensions[pszDimensionIdentifier] = pszDefault;
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    1243                 :            :             }
    1244         [ +  + ]:          1 :             else if( strcmp(psSubIter->pszValue, "TileMatrixSetLink") == 0 )
    1245                 :            :             {
    1246                 :            :                 const char* pszTMS = CPLGetXMLValue(
    1247         [ +  - ]:          1 :                                             psSubIter, "TileMatrixSet", "");
    1248 [ +  - ][ +  + ]:          1 :                 if( !osTMS.empty() && strcmp(osTMS, pszTMS) != 0 )
         [ +  - ][ +  + ]
                 [ +  + ]
    1249                 :          1 :                     continue;
    1250 [ +  - ][ +  - ]:          1 :                 if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1251         [ +  - ]:          1 :                     osSelectTMS.empty() )
    1252                 :            :                 {
    1253 [ +  - ][ +  - ]:          1 :                     osSelectTMS = pszTMS;
                 [ +  - ]
    1254                 :            :                 }
    1255 [ +  - ][ +  - ]:          1 :                 if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1256         [ +  - ]:          1 :                     strcmp(osSelectTMS, pszTMS) == 0 )
    1257                 :            :                 {
    1258                 :            :                     CPLXMLNode* psTMSLimits = CPLGetXMLNode(
    1259         [ +  - ]:          1 :                                         psSubIter, "TileMatrixSetLimits");
    1260         [ +  + ]:          1 :                     if( psTMSLimits )
    1261         [ +  - ]:          1 :                         ReadTMLimits(psTMSLimits, aoMapTileMatrixLimits);
    1262                 :            :                 }
    1263 [ +  - ][ +  - ]:          1 :                 aosTMS.push_back(pszTMS);
                 [ +  - ]
    1264                 :            :             }
    1265         [ +  + ]:          1 :             else if( strcmp(psSubIter->pszValue, "Style") == 0 )
    1266                 :            :             {
    1267                 :            :                 int bIsDefault = CPLTestBool(CPLGetXMLValue(
    1268 [ +  - ][ +  - ]:          1 :                                             psSubIter, "isDefault", "false"));
    1269                 :            :                 const char* l_pszIdentifier = CPLGetXMLValue(
    1270         [ +  - ]:          1 :                                             psSubIter, "Identifier", "");
    1271 [ +  - ][ +  + ]:          1 :                 if( !osStyle.empty() && strcmp(osStyle, l_pszIdentifier) != 0 )
         [ +  - ][ +  + ]
                 [ +  + ]
    1272                 :          1 :                     continue;
    1273                 :            :                 const char* pszStyleTitle = CPLGetXMLValue(
    1274         [ +  - ]:          1 :                                         psSubIter, "Title", l_pszIdentifier);
    1275         [ +  + ]:          1 :                 if( bIsDefault )
    1276                 :            :                 {
    1277                 :            :                     aosStylesIdentifier.insert(aosStylesIdentifier.begin(),
    1278 [ +  - ][ +  - ]:          1 :                                                 CPLString(l_pszIdentifier));
         [ +  - ][ +  - ]
    1279                 :            :                     aosStylesTitle.insert(aosStylesTitle.begin(),
    1280 [ +  - ][ +  - ]:          1 :                                             CPLString(pszStyleTitle));
         [ +  - ][ +  - ]
    1281 [ +  - ][ -  + ]:          1 :                     if( strcmp(osSelectLayer, l_pszIdentifier) == 0 &&
         [ #  # ][ -  + ]
    1282         [ #  # ]:          0 :                         osSelectStyle.empty() )
    1283                 :            :                     {
    1284 [ #  # ][ #  # ]:          0 :                         osSelectStyle = l_pszIdentifier;
                 [ #  # ]
    1285                 :            :                     }
    1286                 :            :                 }
    1287                 :            :                 else
    1288                 :            :                 {
    1289 [ +  - ][ +  - ]:          1 :                     aosStylesIdentifier.push_back(l_pszIdentifier);
                 [ +  - ]
    1290 [ +  - ][ +  - ]:          1 :                     aosStylesTitle.push_back(pszStyleTitle);
                 [ +  - ]
    1291                 :            :                 }
    1292                 :            :             }
    1293 [ +  - ][ +  - ]:          1 :             else if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
                 [ +  + ]
    1294                 :          1 :                      (strcmp(psSubIter->pszValue, "BoundingBox") == 0 ||
    1295                 :          1 :                       strcmp(psSubIter->pszValue, "WGS84BoundingBox") == 0) )
    1296                 :            :             {
    1297 [ +  - ][ +  - ]:          1 :                 CPLString osCRS = CPLGetXMLValue(psSubIter, "crs", "");
    1298 [ +  - ][ +  + ]:          1 :                 if( osCRS.empty() )
    1299                 :            :                 {
    1300         [ +  - ]:          1 :                     if( strcmp(psSubIter->pszValue, "WGS84BoundingBox") == 0 )
    1301                 :            :                     {
    1302 [ +  - ][ +  - ]:          1 :                         osCRS = "EPSG:4326";
                 [ +  - ]
    1303                 :            :                     }
    1304                 :            :                     else
    1305                 :            :                     {
    1306                 :          0 :                         int nCountTileMatrixSet = 0;
    1307         [ #  # ]:          0 :                         CPLString osSingleTileMatrixSet;
    1308         [ #  # ]:          0 :                         for(CPLXMLNode* psIter3 = psContents->psChild; psIter3 != NULL; psIter3 = psIter3->psNext )
    1309                 :            :                         {
    1310 [ #  # ][ #  # ]:          0 :                             if( psIter3->eType != CXT_Element || strcmp(psIter3->pszValue, "TileMatrixSet") != 0 )
    1311                 :          0 :                                 continue;
    1312                 :          0 :                             nCountTileMatrixSet ++;
    1313         [ #  # ]:          0 :                             if( nCountTileMatrixSet == 1 )
    1314 [ #  # ][ #  # ]:          0 :                                 osSingleTileMatrixSet = CPLGetXMLValue(psIter3, "Identifier", "");
         [ #  # ][ #  # ]
    1315                 :            :                         }
    1316         [ #  # ]:          0 :                         if( nCountTileMatrixSet == 1 )
    1317                 :            :                         {
    1318                 :            :                             // For 13-082_WMTS_Simple_Profile/schemas/wmts/1.0/profiles/WMTSSimple/examples/wmtsGetCapabilities_response_OSM.xml
    1319         [ #  # ]:          0 :                             WMTSTileMatrixSet oTMS;
    1320         [ #  # ]:          0 :                             if( ReadTMS(psContents, osSingleTileMatrixSet,
    1321 [ #  # ][ #  # ]:          0 :                                         CPLString(), -1, oTMS) )
                 [ #  # ]
    1322                 :            :                             {
    1323         [ #  # ]:          0 :                                 osCRS = oTMS.osSRS;
    1324         [ #  # ]:          0 :                             }
    1325         [ #  # ]:          0 :                         }
    1326                 :            :                     }
    1327                 :            :                 }
    1328 [ +  - ][ +  - ]:          1 :                 CPLString osLowerCorner = CPLGetXMLValue(psSubIter, "LowerCorner", "");
    1329 [ +  - ][ +  - ]:          1 :                 CPLString osUpperCorner = CPLGetXMLValue(psSubIter, "UpperCorner", "");
    1330         [ +  - ]:          1 :                 OGRSpatialReference oSRS;
    1331 [ +  - ][ +  - ]:          1 :                 if( !osCRS.empty() && !osLowerCorner.empty() && !osUpperCorner.empty() &&
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    1332 [ +  - ][ +  - ]:          1 :                     oSRS.SetFromUserInput(FixCRSName(osCRS)) == OGRERR_NONE )
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ #  # ]
    1333                 :            :                 {
    1334         [ +  - ]:          1 :                     int bSwap = oSRS.EPSGTreatsAsLatLong() ||
    1335 [ +  + ][ +  - ]:          1 :                                 oSRS.EPSGTreatsAsNorthingEasting();
                 [ -  + ]
    1336 [ +  - ][ +  - ]:          1 :                     char** papszLC = CSLTokenizeString(osLowerCorner);
    1337 [ +  - ][ +  - ]:          1 :                     char** papszUC = CSLTokenizeString(osUpperCorner);
    1338 [ +  - ][ +  - ]:          1 :                     if( CSLCount(papszLC) == 2 && CSLCount(papszUC) == 2 )
         [ +  - ][ +  - ]
                 [ +  - ]
    1339                 :            :                     {
    1340         [ +  - ]:          1 :                         OGREnvelope sEnvelope;
    1341 [ +  + ][ +  - ]:          1 :                         sEnvelope.MinX = CPLAtof(papszLC[(bSwap)? 1 : 0]);
    1342 [ +  + ][ +  - ]:          1 :                         sEnvelope.MinY = CPLAtof(papszLC[(bSwap)? 0 : 1]);
    1343 [ +  + ][ +  - ]:          1 :                         sEnvelope.MaxX = CPLAtof(papszUC[(bSwap)? 1 : 0]);
    1344 [ +  + ][ +  - ]:          1 :                         sEnvelope.MaxY = CPLAtof(papszUC[(bSwap)? 0 : 1]);
    1345         [ +  - ]:          1 :                         aoMapBoundingBox[osCRS] = sEnvelope;
    1346                 :            :                     }
    1347         [ +  - ]:          1 :                     CSLDestroy(papszLC);
    1348         [ +  - ]:          1 :                     CSLDestroy(papszUC);
    1349 [ +  - ][ +  - ]:          1 :                 }
         [ +  - ][ +  - ]
    1350                 :            :             }
    1351 [ +  - ][ +  - ]:          1 :             else if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1352                 :          1 :                      strcmp(psSubIter->pszValue, "ResourceURL") == 0 )
    1353                 :            :             {
    1354 [ +  - ][ +  + ]:          1 :                 if( EQUAL(CPLGetXMLValue(psSubIter, "resourceType", ""), "tile") )
    1355                 :            :                 {
    1356         [ +  - ]:          1 :                     const char* pszFormat = CPLGetXMLValue(psSubIter, "format", "");
    1357 [ +  - ][ -  + ]:          1 :                     if( !osTileFormat.empty() && strcmp(osTileFormat, pszFormat) != 0 )
         [ #  # ][ #  # ]
                 [ -  + ]
    1358                 :          0 :                         continue;
    1359 [ +  - ][ +  - ]:          1 :                     if( osURLTileTemplate.empty() )
    1360 [ +  - ][ +  - ]:          1 :                         osURLTileTemplate = CPLGetXMLValue(psSubIter, "template", "");
         [ +  - ][ +  - ]
    1361                 :            :                 }
    1362 [ +  - ][ +  - ]:          1 :                 else if( EQUAL(CPLGetXMLValue(psSubIter, "resourceType", ""), "FeatureInfo") )
    1363                 :            :                 {
    1364         [ +  - ]:          1 :                     const char* pszFormat = CPLGetXMLValue(psSubIter, "format", "");
    1365 [ +  - ][ -  + ]:          1 :                     if( !osInfoFormat.empty() && strcmp(osInfoFormat, pszFormat) != 0 )
         [ #  # ][ #  # ]
                 [ -  + ]
    1366                 :          0 :                         continue;
    1367 [ +  - ][ +  - ]:          1 :                     if( osURLFeatureInfoTemplate.empty() )
    1368 [ +  - ][ +  - ]:          1 :                         osURLFeatureInfoTemplate = CPLGetXMLValue(psSubIter, "template", "");
         [ +  - ][ +  - ]
    1369                 :            :                 }
    1370                 :            :             }
    1371                 :            :         }
    1372 [ +  - ][ +  - ]:          1 :         if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
                 [ +  + ]
    1373 [ +  - ][ +  - ]:          1 :             osSelectStyle.empty() && !aosStylesIdentifier.empty() )
    1374                 :            :         {
    1375 [ +  - ][ +  - ]:          1 :             osSelectStyle = aosStylesIdentifier[0];
    1376                 :            :         }
    1377 [ +  - ][ +  + ]:          1 :         for(size_t i=0;i<aosTMS.size();i++)
    1378                 :            :         {
    1379 [ +  - ][ +  + ]:          1 :             for(size_t j=0;j<aosStylesIdentifier.size();j++)
    1380                 :            :             {
    1381         [ +  - ]:          1 :                 int nIdx = 1 + aosSubDatasets.size() / 2;
    1382         [ +  - ]:          1 :                 CPLString osName(osCapabilitiesFilename);
    1383         [ +  - ]:          1 :                 osName += ",layer=";
    1384 [ +  - ][ +  - ]:          1 :                 osName += QuoteIfNecessary(pszIdentifier);
                 [ +  - ]
    1385 [ +  - ][ +  + ]:          1 :                 if( aosTMS.size() > 1 )
    1386                 :            :                 {
    1387         [ +  - ]:          1 :                     osName += ",tilematrixset=";
    1388 [ +  - ][ +  - ]:          1 :                     osName += QuoteIfNecessary(aosTMS[i]);
         [ +  - ][ +  - ]
                 [ +  - ]
    1389                 :            :                 }
    1390 [ +  - ][ +  + ]:          1 :                 if( aosStylesIdentifier.size() > 1 )
    1391                 :            :                 {
    1392         [ +  - ]:          1 :                     osName += ",style=";
    1393 [ +  - ][ +  - ]:          1 :                     osName += QuoteIfNecessary(aosStylesIdentifier[j]);
         [ +  - ][ +  - ]
                 [ +  - ]
    1394                 :            :                 }
    1395                 :            :                 aosSubDatasets.AddNameValue(
    1396 [ +  - ][ +  - ]:          1 :                     CPLSPrintf("SUBDATASET_%d_NAME", nIdx), osName);
                 [ +  - ]
    1397                 :            : 
    1398         [ +  - ]:          1 :                 CPLString osDesc("Layer ");
    1399 [ +  + ][ +  - ]:          1 :                 osDesc += pszTitle ? pszTitle : pszIdentifier;
    1400 [ +  - ][ +  + ]:          1 :                 if( aosTMS.size() > 1 )
    1401                 :            :                 {
    1402         [ +  - ]:          1 :                     osDesc += ", tile matrix set ";
    1403 [ +  - ][ +  - ]:          1 :                     osDesc += aosTMS[i];
    1404                 :            :                 }
    1405 [ +  - ][ +  + ]:          1 :                 if( aosStylesIdentifier.size() > 1 )
    1406                 :            :                 {
    1407         [ +  - ]:          1 :                     osDesc += ", style ";
    1408 [ +  - ][ +  - ]:          1 :                     osDesc += QuoteIfNecessary(aosStylesTitle[j]);
         [ +  - ][ +  - ]
                 [ +  - ]
    1409                 :            :                 }
    1410                 :            :                 aosSubDatasets.AddNameValue(
    1411 [ +  - ][ +  - ]:          1 :                     CPLSPrintf("SUBDATASET_%d_DESC", nIdx), osDesc);
                 [ +  - ]
    1412 [ +  - ][ +  - ]:          1 :             }
    1413                 :            :         }
    1414 [ +  - ][ +  + ]:          1 :         if( !aosTMS.empty() && !aosStylesIdentifier.empty() )
         [ +  - ][ +  + ]
                 [ +  + ]
    1415                 :          1 :             nLayerCount ++;
    1416                 :            :         else
    1417         [ +  - ]:          1 :             CPLError(CE_Failure, CPLE_AppDefined, "Missing TileMatrixSetLink and/or Style");
    1418 [ +  - ][ +  - ]:          1 :     }
                 [ +  - ]
    1419                 :            : 
    1420         [ +  + ]:          1 :     if( nLayerCount == 0 )
    1421                 :            :     {
    1422         [ +  - ]:          1 :         CPLDestroyXMLNode(psXML);
    1423                 :          1 :         return NULL;
    1424                 :            :     }
    1425                 :            : 
    1426 [ +  - ][ +  - ]:          1 :     WMTSDataset* poDS = new WMTSDataset();
    1427                 :            : 
    1428 [ +  - ][ +  + ]:          1 :     if( aosSubDatasets.size() > 2 )
    1429 [ +  - ][ +  - ]:          1 :         poDS->SetMetadata(aosSubDatasets.List(), "SUBDATASETS");
    1430                 :            : 
    1431         [ +  - ]:          1 :     if( nLayerCount == 1 )
    1432                 :            :     {
    1433 [ +  - ][ +  + ]:          1 :         if( !osSelectLayerTitle.empty() )
    1434 [ +  - ][ +  - ]:          1 :             poDS->SetMetadataItem("TITLE", osSelectLayerTitle);
    1435 [ +  - ][ +  + ]:          1 :         if( !osSelectLayerAbstract.empty() )
    1436 [ +  - ][ +  - ]:          1 :             poDS->SetMetadataItem("ABSTRACT", osSelectLayerAbstract);
    1437                 :            : 
    1438 [ +  - ][ +  - ]:          1 :         poDS->papszHTTPOptions = BuildHTTPRequestOpts(osOtherXML);
                 [ +  - ]
    1439         [ +  - ]:          1 :         poDS->osLayer = osSelectLayer;
    1440         [ +  - ]:          1 :         poDS->osTMS = osSelectTMS;
    1441                 :            : 
    1442         [ +  - ]:          1 :         WMTSTileMatrixSet oTMS;
    1443         [ +  + ]:          1 :         if( !ReadTMS(psContents, osSelectTMS, osMaxTileMatrixIdentifier,
    1444         [ +  - ]:          1 :                      nUserMaxZoomLevel, oTMS) )
    1445                 :            :         {
    1446         [ +  - ]:          1 :             CPLDestroyXMLNode(psXML);
    1447 [ +  - ][ +  - ]:          1 :             delete poDS;
    1448                 :          1 :             return NULL;
    1449                 :            :         }
    1450                 :            : 
    1451                 :            :         const char* pszExtentMethod = CSLFetchNameValueDef(
    1452         [ +  - ]:          1 :             poOpenInfo->papszOpenOptions, "EXTENT_METHOD", "AUTO");
    1453                 :          1 :         ExtentMethod eExtentMethod = AUTO;
    1454         [ -  + ]:          1 :         if( EQUAL(pszExtentMethod, "LAYER_BBOX") )
    1455                 :          0 :             eExtentMethod = LAYER_BBOX;
    1456         [ -  + ]:          1 :         else if( EQUAL(pszExtentMethod, "TILE_MATRIX_SET") )
    1457                 :          0 :             eExtentMethod = TILE_MATRIX_SET;
    1458         [ -  + ]:          1 :         else if( EQUAL(pszExtentMethod, "MOST_PRECISE_TILE_MATRIX") )
    1459                 :          0 :             eExtentMethod = MOST_PRECISE_TILE_MATRIX;
    1460                 :            : 
    1461                 :            :         // Use in priority layer bounding box expressed in the SRS of the TMS
    1462 [ +  + ][ -  + ]:          1 :         if( (!bHasAOI || bExtendBeyondDateLine) &&
         [ -  + ][ #  # ]
         [ +  + ][ +  + ]
    1463                 :            :             (eExtentMethod == AUTO || eExtentMethod == LAYER_BBOX) &&
    1464 [ +  - ][ +  - ]:          1 :             aoMapBoundingBox.find(oTMS.osSRS) != aoMapBoundingBox.end() )
                 [ +  - ]
    1465                 :            :         {
    1466         [ +  - ]:          1 :             if( !bHasAOI )
    1467                 :            :             {
    1468         [ +  - ]:          1 :                 sAOI = aoMapBoundingBox[oTMS.osSRS];
    1469                 :          1 :                 bHasAOI = TRUE;
    1470                 :            :             }
    1471                 :            : 
    1472                 :          1 :             int bRecomputeAOI = FALSE;
    1473         [ +  + ]:          1 :             if( bExtendBeyondDateLine )
    1474                 :            :             {
    1475                 :          1 :                 bExtendBeyondDateLine = FALSE;
    1476                 :            : 
    1477         [ +  - ]:          1 :                 OGRSpatialReference oWGS84;
    1478         [ +  - ]:          1 :                     oWGS84.SetFromUserInput(SRS_WKT_WGS84);
    1479                 :            :                 OGRCoordinateTransformation* poCT =
    1480         [ +  - ]:          1 :                     OGRCreateCoordinateTransformation(&oTMS.oSRS, &oWGS84);
    1481         [ +  - ]:          1 :                 if( poCT != NULL )
    1482                 :            :                 {
    1483                 :          1 :                     double dfX1 = sAOI.MinX;
    1484                 :          1 :                     double dfY1 = sAOI.MinY;
    1485                 :          1 :                     double dfX2 = sAOI.MaxX;
    1486                 :          1 :                     double dfY2 = sAOI.MaxY;
    1487 [ +  - ][ +  - ]:          1 :                     if( poCT->Transform(1, &dfX1, &dfY1) &&
         [ +  - ][ +  - ]
    1488         [ +  - ]:          1 :                         poCT->Transform(1, &dfX2, &dfY2) )
    1489                 :            :                     {
    1490 [ +  - ][ +  - ]:          1 :                         if( fabs(dfX1 + 180) < 1e-8 &&
    1491                 :          1 :                             fabs(dfX2 - 180) < 1e-8 )
    1492                 :            :                         {
    1493                 :          1 :                             bExtendBeyondDateLine = TRUE;
    1494                 :          1 :                             bRecomputeAOI = TRUE;
    1495                 :            :                         }
    1496         [ #  # ]:          0 :                         else if( dfX2 < dfX1 )
    1497                 :            :                         {
    1498                 :          0 :                             bExtendBeyondDateLine = TRUE;
    1499                 :            :                         }
    1500                 :            :                         else
    1501                 :            :                         {
    1502                 :            :                             CPLError(CE_Warning, CPLE_AppDefined,
    1503                 :            :                                     "ExtendBeyondDateLine disabled, since longitudes of %s "
    1504                 :            :                                     "BoundingBox do not span from -180 to 180 but from %.16g to %.16g, "
    1505                 :            :                                     "or longitude of upper right corner is not lesser than the one of lower left corner",
    1506 [ #  # ][ #  # ]:          1 :                                     oTMS.osSRS.c_str(), dfX1, dfX2);
    1507                 :            :                         }
    1508                 :            :                     }
    1509 [ +  - ][ +  - ]:          1 :                     delete poCT;
    1510         [ +  - ]:          1 :                 }
    1511                 :            :             }
    1512 [ +  + ][ +  - ]:          1 :             if( bExtendBeyondDateLine && bRecomputeAOI )
    1513                 :            :             {
    1514                 :          1 :                 bExtendBeyondDateLine = FALSE;
    1515                 :            : 
    1516         [ +  - ]:          1 :                 std::map<CPLString, OGREnvelope>::iterator oIter = aoMapBoundingBox.begin();
    1517 [ +  - ][ +  - ]:          1 :                 for(; oIter != aoMapBoundingBox.end(); ++oIter )
         [ +  - ][ +  - ]
    1518                 :            :                 {
    1519         [ +  - ]:          1 :                     OGRSpatialReference oSRS;
    1520 [ +  - ][ +  - ]:          1 :                     if( oSRS.SetFromUserInput(FixCRSName(oIter->first)) == OGRERR_NONE )
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1521                 :            :                     {
    1522         [ +  - ]:          1 :                         OGRSpatialReference oWGS84;
    1523         [ +  - ]:          1 :                         oWGS84.SetFromUserInput(SRS_WKT_WGS84);
    1524                 :            :                         OGRCoordinateTransformation* poCT =
    1525         [ +  - ]:          1 :                             OGRCreateCoordinateTransformation(&oSRS, &oWGS84);
    1526         [ +  - ]:          1 :                         double dfX1 = oIter->second.MinX;
    1527         [ +  - ]:          1 :                         double dfY1 = oIter->second.MinY;
    1528         [ +  - ]:          1 :                         double dfX2 = oIter->second.MaxX;
    1529         [ +  - ]:          1 :                         double dfY2 = oIter->second.MaxY;
    1530 [ +  - ][ +  - ]:          1 :                         if( poCT != NULL &&
         [ +  - ][ +  + ]
                 [ +  + ]
    1531         [ +  - ]:          1 :                             poCT->Transform(1, &dfX1, &dfY1) &&
    1532         [ +  - ]:          1 :                             poCT->Transform(1, &dfX2, &dfY2) &&
    1533                 :            :                             dfX2 < dfX1 )
    1534                 :            :                         {
    1535 [ +  - ][ +  - ]:          1 :                             delete poCT;
    1536                 :          1 :                             dfX2 += 360;
    1537                 :          1 :                             char* pszProj4 = NULL;
    1538         [ +  - ]:          1 :                             oTMS.oSRS.exportToProj4(&pszProj4);
    1539 [ +  - ][ +  - ]:          1 :                             oSRS.SetFromUserInput(CPLSPrintf("%s +over +wktext", pszProj4));
    1540         [ +  - ]:          1 :                             CPLFree(pszProj4);
    1541         [ +  - ]:          1 :                             poCT = OGRCreateCoordinateTransformation(&oWGS84, &oSRS);
    1542 [ +  - ][ +  - ]:          1 :                             if( poCT &&
         [ +  - ][ +  - ]
    1543         [ +  - ]:          1 :                                 poCT->Transform(1, &dfX1, &dfY1) &&
    1544         [ +  - ]:          1 :                                 poCT->Transform(1, &dfX2, &dfY2) )
    1545                 :            :                             {
    1546                 :          1 :                                 bExtendBeyondDateLine = TRUE;
    1547         [ +  - ]:          1 :                                 sAOI.MinX = std::min(dfX1, dfX2);
    1548         [ +  - ]:          1 :                                 sAOI.MinY = std::min(dfY1, dfY2);
    1549         [ +  - ]:          1 :                                 sAOI.MaxX = std::max(dfX1, dfX2);
    1550         [ +  - ]:          1 :                                 sAOI.MaxY = std::max(dfY1, dfY2);
    1551                 :            :                                 CPLDebug("WMTS",
    1552                 :            :                                          "ExtendBeyondDateLine using %s bounding box",
    1553 [ +  - ][ +  - ]:          1 :                                          oIter->first.c_str());
                 [ +  - ]
    1554                 :            :                             }
    1555 [ +  - ][ +  - ]:          1 :                             delete poCT;
    1556                 :            :                             break;
    1557                 :            :                         }
    1558 [ +  - ][ +  - ]:          1 :                         delete poCT;
         [ +  - ][ +  + ]
    1559                 :            :                     }
    1560 [ +  - ][ +  + ]:          1 :                 }
    1561                 :            :             }
    1562                 :            :         }
    1563                 :            :         else
    1564                 :            :         {
    1565         [ -  + ]:          1 :             if( bExtendBeyondDateLine )
    1566                 :            :             {
    1567                 :            :                 CPLError(CE_Warning, CPLE_AppDefined,
    1568                 :            :                          "ExtendBeyondDateLine disabled, since BoundingBox of %s is missing",
    1569 [ #  # ][ #  # ]:          0 :                          oTMS.osSRS.c_str());
    1570                 :          0 :                 bExtendBeyondDateLine = FALSE;
    1571                 :            :             }
    1572                 :            :         }
    1573                 :            : 
    1574                 :            :         // Otherwise default to reproject a layer bounding box expressed in
    1575                 :            :         // another SRS
    1576 [ +  + ][ +  - ]:          1 :         if( !bHasAOI && !aoMapBoundingBox.empty() &&
         [ +  + ][ -  + ]
         [ #  # ][ +  + ]
    1577                 :            :             (eExtentMethod == AUTO || eExtentMethod == LAYER_BBOX) )
    1578                 :            :         {
    1579         [ +  - ]:          1 :             std::map<CPLString, OGREnvelope>::iterator oIter = aoMapBoundingBox.begin();
    1580 [ #  # ][ +  - ]:          1 :             for(; oIter != aoMapBoundingBox.end(); ++oIter )
         [ +  - ][ +  - ]
    1581                 :            :             {
    1582         [ +  - ]:          1 :                 OGRSpatialReference oSRS;
    1583 [ +  - ][ +  - ]:          1 :                 if( oSRS.SetFromUserInput(FixCRSName(oIter->first)) == OGRERR_NONE )
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1584                 :            :                 {
    1585                 :            :                     // Check if this doesn't match the most precise tile matrix
    1586                 :            :                     // by densifying its contour
    1587         [ +  - ]:          1 :                     const WMTSTileMatrix& oTM = oTMS.aoTM.back();
    1588                 :            : 
    1589                 :          1 :                     bool bMatchFound = false;
    1590         [ +  - ]:          1 :                     const char *pszProjectionTMS = oTMS.oSRS.GetAttrValue("PROJECTION");
    1591         [ +  - ]:          1 :                     const char *pszProjectionBBOX = oSRS.GetAttrValue("PROJECTION");
    1592                 :            :                     const bool bIsTMerc = (pszProjectionTMS != NULL &&
    1593                 :          1 :                         EQUAL(pszProjectionTMS, SRS_PT_TRANSVERSE_MERCATOR)) ||
    1594                 :            :                         (pszProjectionBBOX != NULL &&
    1595 [ +  + ][ +  + ]:          1 :                         EQUAL(pszProjectionBBOX, SRS_PT_TRANSVERSE_MERCATOR));
         [ -  + ][ #  # ]
    1596                 :            :                     // If one of the 2 SRS is a TMerc, try with classical tmerc
    1597                 :            :                     // or etmerc.
    1598 [ +  + ][ +  + ]:          1 :                     for( int j = 0; j < (bIsTMerc ? 2 : 1); j++ )
    1599                 :            :                     {
    1600                 :            :                         CPLString osOldVal =
    1601 [ +  - ][ +  - ]:          1 :                             CPLGetThreadLocalConfigOption("OSR_USE_ETMERC", "");
    1602         [ +  + ]:          1 :                         if( bIsTMerc )
    1603                 :            :                         {
    1604                 :            :                             CPLSetThreadLocalConfigOption("OSR_USE_ETMERC",
    1605 [ +  - ][ +  - ]:          1 :                                                       (j==0) ? "NO" : "YES");
    1606                 :            :                         }
    1607                 :            :                         OGRCoordinateTransformation* poRevCT =
    1608         [ +  - ]:          1 :                             OGRCreateCoordinateTransformation(&oTMS.oSRS, &oSRS);
    1609         [ +  + ]:          1 :                         if( bIsTMerc )
    1610                 :            :                         {
    1611                 :            :                             CPLSetThreadLocalConfigOption("OSR_USE_ETMERC",
    1612 [ +  - ][ +  - ]:          1 :                                 osOldVal.empty() ? NULL : osOldVal.c_str());
         [ #  # ][ +  - ]
    1613                 :            :                         }
    1614         [ +  - ]:          1 :                         if( poRevCT != NULL )
    1615                 :            :                         {
    1616                 :          1 :                             const double dfX0 = oTM.dfTLX;
    1617                 :          1 :                             const double dfY1 = oTM.dfTLY;
    1618                 :            :                             const double dfX1 = oTM.dfTLX +
    1619                 :          1 :                             oTM.nMatrixWidth  * oTM.dfPixelSize * oTM.nTileWidth;
    1620                 :            :                             const double dfY0 = oTM.dfTLY -
    1621                 :          1 :                             oTM.nMatrixHeight * oTM.dfPixelSize * oTM.nTileHeight;
    1622                 :          1 :                             double dfXMin = std::numeric_limits<double>::infinity();
    1623                 :          1 :                             double dfYMin = std::numeric_limits<double>::infinity();
    1624                 :          1 :                             double dfXMax = -std::numeric_limits<double>::infinity();
    1625                 :          1 :                             double dfYMax = -std::numeric_limits<double>::infinity();
    1626                 :            : 
    1627                 :          1 :                             const int NSTEPS = 20;
    1628         [ +  + ]:          1 :                             for(int i=0;i<=NSTEPS;i++)
    1629                 :            :                             {
    1630                 :          1 :                                 double dfX = dfX0 + (dfX1 - dfX0) * i / NSTEPS;
    1631                 :          1 :                                 double dfY = dfY0;
    1632 [ +  - ][ +  - ]:          1 :                                 if( poRevCT->Transform(1, &dfX, &dfY) )
    1633                 :            :                                 {
    1634         [ +  - ]:          1 :                                     dfXMin = std::min(dfXMin, dfX);
    1635         [ +  - ]:          1 :                                     dfYMin = std::min(dfYMin, dfY);
    1636         [ +  - ]:          1 :                                     dfXMax = std::max(dfXMax, dfX);
    1637         [ +  - ]:          1 :                                     dfYMax = std::max(dfYMax, dfY);
    1638                 :            :                                 }
    1639                 :            : 
    1640                 :          1 :                                 dfX = dfX0 + (dfX1 - dfX0) * i / NSTEPS;
    1641                 :          1 :                                 dfY = dfY1;
    1642 [ +  - ][ +  - ]:          1 :                                 if( poRevCT->Transform(1, &dfX, &dfY) )
    1643                 :            :                                 {
    1644         [ +  - ]:          1 :                                     dfXMin = std::min(dfXMin, dfX);
    1645         [ +  - ]:          1 :                                     dfYMin = std::min(dfYMin, dfY);
    1646         [ +  - ]:          1 :                                     dfXMax = std::max(dfXMax, dfX);
    1647         [ +  - ]:          1 :                                     dfYMax = std::max(dfYMax, dfY);
    1648                 :            :                                 }
    1649                 :            : 
    1650                 :          1 :                                 dfX = dfX0;
    1651                 :          1 :                                 dfY = dfY0 + (dfY1 - dfY0) * i / NSTEPS;
    1652 [ +  - ][ +  - ]:          1 :                                 if( poRevCT->Transform(1, &dfX, &dfY) )
    1653                 :            :                                 {
    1654         [ +  - ]:          1 :                                     dfXMin = std::min(dfXMin, dfX);
    1655         [ +  - ]:          1 :                                     dfYMin = std::min(dfYMin, dfY);
    1656         [ +  - ]:          1 :                                     dfXMax = std::max(dfXMax, dfX);
    1657         [ +  - ]:          1 :                                     dfYMax = std::max(dfYMax, dfY);
    1658                 :            :                                 }
    1659                 :            : 
    1660                 :          1 :                                 dfX = dfX1;
    1661                 :          1 :                                 dfY = dfY0 + (dfY1 - dfY0) * i / NSTEPS;
    1662 [ +  - ][ +  - ]:          1 :                                 if( poRevCT->Transform(1, &dfX, &dfY) )
    1663                 :            :                                 {
    1664         [ +  - ]:          1 :                                     dfXMin = std::min(dfXMin, dfX);
    1665         [ +  - ]:          1 :                                     dfYMin = std::min(dfYMin, dfY);
    1666         [ +  - ]:          1 :                                     dfXMax = std::max(dfXMax, dfX);
    1667         [ +  - ]:          1 :                                     dfYMax = std::max(dfYMax, dfY);
    1668                 :            :                                 }
    1669                 :            :                             }
    1670                 :            : 
    1671 [ +  - ][ +  - ]:          1 :                             delete poRevCT;
    1672                 :            : #ifdef DEBUG_VERBOSE
    1673                 :            :                             CPLDebug("WMTS", "Reprojected densified bbox of most "
    1674                 :            :                                     "precise tile matrix in %s: %.8g %8g %8g %8g",
    1675                 :            :                                     oIter->first.c_str(),
    1676                 :            :                                     dfXMin, dfYMin, dfXMax, dfYMax);
    1677                 :            : #endif
    1678 [ +  - ][ +  + ]:          1 :                             if( fabs(oIter->second.MinX - dfXMin) < 1e-5 *
         [ +  - ][ +  - ]
         [ +  - ][ +  + ]
    1679 [ +  - ][ +  - ]:          1 :                                 std::max(fabs(oIter->second.MinX),fabs(dfXMin)) &&
    1680         [ +  - ]:          1 :                                 fabs(oIter->second.MinY - dfYMin) < 1e-5 *
    1681 [ +  - ][ +  - ]:          1 :                                 std::max(fabs(oIter->second.MinY),fabs(dfYMin)) &&
    1682         [ +  - ]:          1 :                                 fabs(oIter->second.MaxX - dfXMax) < 1e-5 *
    1683 [ +  - ][ +  - ]:          1 :                                 std::max(fabs(oIter->second.MaxX),fabs(dfXMax)) &&
    1684         [ +  - ]:          1 :                                 fabs(oIter->second.MaxY - dfYMax) < 1e-5 *
    1685 [ +  - ][ +  - ]:          1 :                                 std::max(fabs(oIter->second.MaxY),fabs(dfYMax)) )
    1686                 :            :                             {
    1687                 :          1 :                                 bMatchFound = true;
    1688                 :            : #ifdef DEBUG_VERBOSE
    1689                 :            :                                 CPLDebug("WMTS", "Matches layer bounding box, so "
    1690                 :            :                                         "that one is not significant");
    1691                 :            : #endif
    1692                 :            :                                 break;
    1693                 :            :                             }
    1694                 :            :                         }
    1695 [ +  - ][ +  + ]:          1 :                     }
    1696                 :            : 
    1697         [ +  + ]:          1 :                     if( bMatchFound )
    1698                 :            :                     {
    1699         [ -  + ]:          1 :                         if( eExtentMethod == LAYER_BBOX )
    1700                 :          1 :                             eExtentMethod = MOST_PRECISE_TILE_MATRIX;
    1701                 :            :                         break;
    1702                 :            :                     }
    1703                 :            : 
    1704                 :            :                     OGRCoordinateTransformation* poCT =
    1705         [ +  - ]:          1 :                         OGRCreateCoordinateTransformation(&oSRS, &oTMS.oSRS);
    1706         [ +  - ]:          1 :                     if( poCT != NULL )
    1707                 :            :                     {
    1708         [ +  - ]:          1 :                         double dfX1 = oIter->second.MinX;
    1709         [ +  - ]:          1 :                         double dfY1 = oIter->second.MinY;
    1710         [ +  - ]:          1 :                         double dfX2 = oIter->second.MaxX;
    1711         [ +  - ]:          1 :                         double dfY2 = oIter->second.MinY;
    1712         [ +  - ]:          1 :                         double dfX3 = oIter->second.MaxX;
    1713         [ +  - ]:          1 :                         double dfY3 = oIter->second.MaxY;
    1714         [ +  - ]:          1 :                         double dfX4 = oIter->second.MinX;
    1715         [ +  - ]:          1 :                         double dfY4 = oIter->second.MaxY;
    1716 [ +  - ][ +  - ]:          1 :                         if( poCT->Transform(1, &dfX1, &dfY1) &&
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    1717         [ +  - ]:          1 :                             poCT->Transform(1, &dfX2, &dfY2) &&
    1718         [ +  - ]:          1 :                             poCT->Transform(1, &dfX3, &dfY3) &&
    1719         [ +  - ]:          1 :                             poCT->Transform(1, &dfX4, &dfY4) )
    1720                 :            :                         {
    1721         [ +  - ]:          1 :                             sAOI.MinX = std::min(std::min(dfX1, dfX2),
    1722 [ +  - ][ +  - ]:          1 :                                                  std::min(dfX3, dfX4));
    1723         [ +  - ]:          1 :                             sAOI.MinY = std::min(std::min(dfY1, dfY2),
    1724 [ +  - ][ +  - ]:          1 :                                                  std::min(dfY3, dfY4));
    1725         [ +  - ]:          1 :                             sAOI.MaxX = std::max(std::max(dfX1, dfX2),
    1726 [ +  - ][ +  - ]:          1 :                                                  std::max(dfX3, dfX4));
    1727         [ +  - ]:          1 :                             sAOI.MaxY = std::max(std::max(dfY1, dfY2),
    1728 [ +  - ][ +  - ]:          1 :                                                  std::max(dfY3, dfY4));
    1729                 :          1 :                             bHasAOI = TRUE;
    1730                 :            :                         }
    1731 [ +  - ][ +  - ]:          1 :                         delete poCT;
    1732                 :            :                     }
    1733                 :            :                     break;
    1734                 :            :                 }
    1735 [ +  - ][ -  + ]:          1 :             }
    1736                 :            :         }
    1737                 :            : 
    1738                 :            :         // Otherwise default to BoundingBox of the TMS
    1739 [ +  + ][ +  + ]:          1 :         if( !bHasAOI && oTMS.bBoundingBoxValid &&
         [ -  + ][ #  # ]
    1740                 :            :             (eExtentMethod == AUTO || eExtentMethod == TILE_MATRIX_SET) )
    1741                 :            :         {
    1742         [ +  - ]:          1 :             CPLDebug("WMTS", "Using TMS bounding box");
    1743                 :          1 :             sAOI = oTMS.sBoundingBox;
    1744                 :          1 :             bHasAOI = TRUE;
    1745                 :            :         }
    1746                 :            : 
    1747                 :            :         // Otherwise default to implied BoundingBox of the most precise TM
    1748 [ +  + ][ -  + ]:          1 :         if( !bHasAOI &&
                 [ #  # ]
    1749                 :            :             (eExtentMethod == AUTO || eExtentMethod == MOST_PRECISE_TILE_MATRIX) )
    1750                 :            :         {
    1751         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM.back();
    1752 [ +  - ][ +  - ]:          1 :             CPLDebug("WMTS", "Using TM level %s bounding box", oTM.osIdentifier.c_str() );
    1753                 :            : 
    1754                 :          1 :             sAOI.MinX = oTM.dfTLX;
    1755                 :          1 :             sAOI.MaxY = oTM.dfTLY;
    1756                 :          1 :             sAOI.MaxX = oTM.dfTLX + oTM.nMatrixWidth  * oTM.dfPixelSize * oTM.nTileWidth;
    1757                 :          1 :             sAOI.MinY = oTM.dfTLY - oTM.nMatrixHeight * oTM.dfPixelSize * oTM.nTileHeight;
    1758                 :          1 :             bHasAOI = TRUE;
    1759                 :            :         }
    1760                 :            : 
    1761         [ -  + ]:          1 :         if( !bHasAOI )
    1762                 :            :         {
    1763                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
    1764         [ #  # ]:          0 :                      "Could not determine raster extent");
    1765         [ #  # ]:          0 :             CPLDestroyXMLNode(psXML);
    1766 [ #  # ][ #  # ]:          0 :             delete poDS;
    1767                 :          0 :             return NULL;
    1768                 :            :         }
    1769                 :            : 
    1770                 :            :         {
    1771                 :            :             // Clip with implied BoundingBox of the most precise TM
    1772                 :            :             // Useful for http://tileserver.maptiler.com/wmts
    1773         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM.back();
    1774                 :            : 
    1775                 :            :             // For https://data.linz.govt.nz/services;key=XXXXXXXX/wmts/1.0.0/set/69/WMTSCapabilities.xml
    1776                 :            :             // only clip in Y since there's a warp over dateline.
    1777                 :            :             // Update: it sems that the content of the server has changed since
    1778                 :            :             // initial coding. So do X clipping in default mode.
    1779         [ +  + ]:          1 :             if( !bExtendBeyondDateLine )
    1780                 :            :             {
    1781         [ +  - ]:          1 :                 sAOI.MinX = std::max(sAOI.MinX, oTM.dfTLX);
    1782                 :            :                 sAOI.MaxX = std::min(sAOI.MaxX,
    1783                 :            :                     oTM.dfTLX +
    1784         [ +  - ]:          1 :                     oTM.nMatrixWidth  * oTM.dfPixelSize * oTM.nTileWidth);
    1785                 :            :             }
    1786         [ +  - ]:          1 :             sAOI.MaxY = std::min(sAOI.MaxY, oTM.dfTLY);
    1787                 :            :             sAOI.MinY =
    1788                 :            :                 std::max(sAOI.MinY,
    1789                 :            :                          oTM.dfTLY -
    1790         [ +  - ]:          1 :                          oTM.nMatrixHeight * oTM.dfPixelSize * oTM.nTileHeight);
    1791                 :            :         }
    1792                 :            : 
    1793                 :            :         // Clip with limits of most precise TM when available
    1794                 :            :         {
    1795         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM.back();
    1796 [ +  - ][ +  - ]:          1 :             if( aoMapTileMatrixLimits.find(oTM.osIdentifier) != aoMapTileMatrixLimits.end() )
         [ +  - ][ +  + ]
    1797                 :            :             {
    1798         [ +  - ]:          1 :                 const WMTSTileMatrixLimits& oTMLimits = aoMapTileMatrixLimits[oTM.osIdentifier];
    1799                 :          1 :                 double dfTileWidthUnits = oTM.dfPixelSize * oTM.nTileWidth;
    1800                 :          1 :                 double dfTileHeightUnits = oTM.dfPixelSize * oTM.nTileHeight;
    1801         [ +  - ]:          1 :                 sAOI.MinX = std::max(sAOI.MinX, oTM.dfTLX + oTMLimits.nMinTileCol * dfTileWidthUnits);
    1802         [ +  - ]:          1 :                 sAOI.MaxY = std::min(sAOI.MaxY, oTM.dfTLY - oTMLimits.nMinTileRow * dfTileHeightUnits);
    1803         [ +  - ]:          1 :                 sAOI.MaxX = std::min(sAOI.MaxX, oTM.dfTLX + (oTMLimits.nMaxTileCol + 1) * dfTileWidthUnits);
    1804         [ +  - ]:          1 :                 sAOI.MinY = std::max(sAOI.MinY, oTM.dfTLY - (oTMLimits.nMaxTileRow + 1) * dfTileHeightUnits);
    1805                 :            :             }
    1806                 :            :         }
    1807                 :            : 
    1808                 :            :         // Establish raster dimension and extent
    1809         [ +  - ]:          1 :         int nMaxZoomLevel = (int)oTMS.aoTM.size()-1;
    1810         [ +  - ]:          1 :         while(nMaxZoomLevel >= 0)
    1811                 :            :         {
    1812         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM[nMaxZoomLevel];
    1813                 :          1 :             double dfRasterXSize = (sAOI.MaxX - sAOI.MinX) / oTM.dfPixelSize;
    1814                 :          1 :             double dfRasterYSize = (sAOI.MaxY - sAOI.MinY) / oTM.dfPixelSize;
    1815 [ +  + ][ +  - ]:          1 :             if( dfRasterXSize < INT_MAX && dfRasterYSize < INT_MAX )
    1816                 :            :             {
    1817 [ +  - ][ +  + ]:          1 :                 if( nMaxZoomLevel != (int)oTMS.aoTM.size()-1 )
    1818                 :            :                 {
    1819                 :            :                     CPLDebug("WMTS", "Using zoom level %s instead of %s to avoid int overflow",
    1820         [ +  - ]:          1 :                              oTMS.aoTM[nMaxZoomLevel].osIdentifier.c_str(),
    1821 [ +  - ][ +  - ]:          1 :                              oTMS.aoTM.back().osIdentifier.c_str());
         [ +  - ][ +  - ]
    1822                 :            :                 }
    1823                 :            : 
    1824                 :            :                 // Align AOI on pixel boundaries with respect to TopLeftCorner of
    1825                 :            :                 // this tile matrix
    1826                 :          1 :                 poDS->adfGT[0] = oTM.dfTLX + floor((sAOI.MinX - oTM.dfTLX) / oTM.dfPixelSize+1e-10) * oTM.dfPixelSize;
    1827                 :          1 :                 poDS->adfGT[1] = oTM.dfPixelSize;
    1828                 :          1 :                 poDS->adfGT[2] = 0.0;
    1829                 :          1 :                 poDS->adfGT[3] = oTM.dfTLY + ceil((sAOI.MaxY - oTM.dfTLY) / oTM.dfPixelSize-1e-10) * oTM.dfPixelSize;
    1830                 :          1 :                 poDS->adfGT[4] = 0.0;
    1831                 :          1 :                 poDS->adfGT[5] = -oTM.dfPixelSize;
    1832                 :          1 :                 poDS->nRasterXSize = int(0.5 + (sAOI.MaxX - poDS->adfGT[0]) / oTM.dfPixelSize);
    1833                 :          1 :                 poDS->nRasterYSize = int(0.5 + (poDS->adfGT[3] - sAOI.MinY) / oTM.dfPixelSize);
    1834                 :          1 :                 break;
    1835                 :            :             }
    1836                 :          1 :             nMaxZoomLevel --;
    1837                 :            :         }
    1838         [ -  + ]:          1 :         if( nMaxZoomLevel < 0 )
    1839                 :            :         {
    1840                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
    1841         [ #  # ]:          0 :                      "No zoom level in tile matrix set found");
    1842         [ #  # ]:          0 :             CPLDestroyXMLNode(psXML);
    1843 [ #  # ][ #  # ]:          0 :             delete poDS;
    1844                 :          0 :             return NULL;
    1845                 :            :         }
    1846                 :            :         CPLDebug("WMTS", "Using tilematrix=%s (zoom level %d)",
    1847 [ +  - ][ +  - ]:          1 :                  oTMS.aoTM[nMaxZoomLevel].osIdentifier.c_str(), nMaxZoomLevel);
                 [ +  - ]
    1848 [ +  - ][ +  - ]:          1 :         oTMS.aoTM.resize(1 + nMaxZoomLevel);
                 [ +  - ]
    1849         [ +  - ]:          1 :         poDS->oTMS = oTMS;
    1850                 :            : 
    1851 [ +  - ][ -  + ]:          1 :         if( !osProjection.empty() )
    1852                 :            :         {
    1853         [ #  # ]:          0 :             OGRSpatialReference oSRS;
    1854 [ #  # ][ #  # ]:          0 :             if( oSRS.SetFromUserInput(osProjection) == OGRERR_NONE )
                 [ #  # ]
    1855                 :            :             {
    1856                 :          0 :                 char* pszWKT = NULL;
    1857         [ #  # ]:          0 :                 oSRS.exportToWkt(&pszWKT);
    1858 [ #  # ][ #  # ]:          0 :                 poDS->osProjection = pszWKT;
                 [ #  # ]
    1859         [ #  # ]:          0 :                 CPLFree(pszWKT);
    1860         [ #  # ]:          0 :             }
    1861                 :            :         }
    1862 [ +  - ][ +  - ]:          1 :         if( poDS->osProjection.empty() )
    1863                 :            :         {
    1864                 :            :             // Strip AXIS
    1865         [ +  - ]:          1 :             OGR_SRSNode *poGEOGCS = oTMS.oSRS.GetAttrNode( "GEOGCS" );
    1866         [ +  - ]:          1 :             if( poGEOGCS != NULL )
    1867         [ +  - ]:          1 :                 poGEOGCS->StripNodes( "AXIS" );
    1868                 :            : 
    1869         [ +  - ]:          1 :             OGR_SRSNode *poPROJCS = oTMS.oSRS.GetAttrNode( "PROJCS" );
    1870 [ +  + ][ +  - ]:          1 :             if (poPROJCS != NULL && oTMS.oSRS.EPSGTreatsAsNorthingEasting())
         [ -  + ][ -  + ]
    1871         [ #  # ]:          0 :                 poPROJCS->StripNodes( "AXIS" );
    1872                 :            : 
    1873                 :          1 :             char* pszWKT = NULL;
    1874         [ +  - ]:          1 :             oTMS.oSRS.exportToWkt(&pszWKT);
    1875 [ +  - ][ +  - ]:          1 :             poDS->osProjection = pszWKT;
                 [ +  - ]
    1876         [ +  - ]:          1 :             CPLFree(pszWKT);
    1877                 :            :         }
    1878                 :            : 
    1879 [ +  - ][ +  + ]:          1 :         if( osURLTileTemplate.empty() )
    1880                 :            :         {
    1881 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = GetOperationKVPURL(psXML, "GetTile");
                 [ +  - ]
    1882 [ +  - ][ +  + ]:          1 :             if( osURLTileTemplate.empty() )
    1883                 :            :             {
    1884                 :            :                 CPLError(CE_Failure, CPLE_AppDefined,
    1885         [ +  - ]:          1 :                          "No RESTful nor KVP GetTile operation found");
    1886         [ +  - ]:          1 :                 CPLDestroyXMLNode(psXML);
    1887 [ +  - ][ +  - ]:          1 :                 delete poDS;
    1888                 :          1 :                 return NULL;
    1889                 :            :             }
    1890 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "service", "WMTS");
         [ +  - ][ +  - ]
    1891 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "request", "GetTile");
         [ +  - ][ +  - ]
    1892 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "version", "1.0.0");
         [ +  - ][ +  - ]
    1893 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "layer", osSelectLayer);
         [ +  - ][ +  - ]
                 [ +  - ]
    1894 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "style", osSelectStyle);
         [ +  - ][ +  - ]
                 [ +  - ]
    1895 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "format", osSelectTileFormat);
         [ +  - ][ +  - ]
                 [ +  - ]
    1896 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "TileMatrixSet", osSelectTMS);
         [ +  - ][ +  - ]
                 [ +  - ]
    1897         [ +  - ]:          1 :             osURLTileTemplate += "&TileMatrix={TileMatrix}";
    1898         [ +  - ]:          1 :             osURLTileTemplate += "&TileRow=${y}";
    1899         [ +  - ]:          1 :             osURLTileTemplate += "&TileCol=${x}";
    1900                 :            : 
    1901         [ +  - ]:          1 :             std::map<CPLString,CPLString>::iterator oIter = aoMapDimensions.begin();
    1902 [ +  - ][ +  - ]:          1 :             for(; oIter != aoMapDimensions.end(); ++oIter )
         [ +  - ][ +  + ]
    1903                 :            :             {
    1904                 :            :                 osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate,
    1905         [ +  - ]:          1 :                                                  oIter->first,
    1906 [ +  - ][ +  - ]:          1 :                                                  oIter->second);
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1907                 :            :             }
    1908                 :            :             //CPLDebug("WMTS", "osURLTileTemplate = %s", osURLTileTemplate.c_str());
    1909                 :            :         }
    1910                 :            :         else
    1911                 :            :         {
    1912 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = Replace(osURLTileTemplate, "{Style}", osSelectStyle);
         [ +  - ][ +  - ]
    1913 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = Replace(osURLTileTemplate, "{TileMatrixSet}", osSelectTMS);
         [ +  - ][ +  - ]
    1914 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = Replace(osURLTileTemplate, "{TileCol}", "${x}");
                 [ +  - ]
    1915 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = Replace(osURLTileTemplate, "{TileRow}", "${y}");
                 [ +  - ]
    1916                 :            : 
    1917         [ +  - ]:          1 :             std::map<CPLString,CPLString>::iterator oIter = aoMapDimensions.begin();
    1918 [ +  - ][ +  - ]:          1 :             for(; oIter != aoMapDimensions.end(); ++oIter )
         [ +  - ][ +  + ]
    1919                 :            :             {
    1920                 :            :                 osURLTileTemplate = Replace(osURLTileTemplate,
    1921         [ +  - ]:          1 :                                             CPLSPrintf("{%s}", oIter->first.c_str()),
    1922 [ +  - ][ +  - ]:          1 :                                             oIter->second);
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1923                 :            :             }
    1924                 :            :         }
    1925                 :            : 
    1926 [ +  - ][ +  + ]:          1 :         if( osURLFeatureInfoTemplate.empty() && !osSelectInfoFormat.empty() )
         [ +  - ][ +  + ]
                 [ +  + ]
    1927                 :            :         {
    1928 [ +  - ][ +  - ]:          1 :             osURLFeatureInfoTemplate = GetOperationKVPURL(psXML, "GetFeatureInfo");
                 [ +  - ]
    1929 [ +  - ][ +  - ]:          1 :             if( !osURLFeatureInfoTemplate.empty() )
    1930                 :            :             {
    1931 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "service", "WMTS");
         [ +  - ][ +  - ]
    1932 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "request", "GetFeatureInfo");
         [ +  - ][ +  - ]
    1933 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "version", "1.0.0");
         [ +  - ][ +  - ]
    1934 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "layer", osSelectLayer);
         [ +  - ][ +  - ]
                 [ +  - ]
    1935 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "style", osSelectStyle);
         [ +  - ][ +  - ]
                 [ +  - ]
    1936                 :            :                 //osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "format", osSelectTileFormat);
    1937 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "InfoFormat", osSelectInfoFormat);
         [ +  - ][ +  - ]
                 [ +  - ]
    1938         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&TileMatrixSet={TileMatrixSet}";
    1939         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&TileMatrix={TileMatrix}";
    1940         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&TileRow={TileRow}";
    1941         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&TileCol={TileCol}";
    1942         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&J={J}";
    1943         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&I={I}";
    1944                 :            : 
    1945         [ +  - ]:          1 :                 std::map<CPLString,CPLString>::iterator oIter = aoMapDimensions.begin();
    1946 [ +  - ][ +  - ]:          1 :                 for(; oIter != aoMapDimensions.end(); ++oIter )
         [ +  - ][ +  + ]
    1947                 :            :                 {
    1948                 :            :                     osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate,
    1949         [ +  - ]:          1 :                                                             oIter->first,
    1950 [ +  - ][ +  - ]:          1 :                                                             oIter->second);
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1951                 :            :                 }
    1952                 :            :                 //CPLDebug("WMTS", "osURLFeatureInfoTemplate = %s", osURLFeatureInfoTemplate.c_str());
    1953                 :            :             }
    1954                 :            :         }
    1955                 :            :         else
    1956                 :            :         {
    1957 [ +  - ][ +  - ]:          1 :              osURLFeatureInfoTemplate = Replace(osURLFeatureInfoTemplate, "{Style}", osSelectStyle);
         [ +  - ][ +  - ]
    1958                 :            : 
    1959         [ +  - ]:          1 :             std::map<CPLString,CPLString>::iterator oIter = aoMapDimensions.begin();
    1960 [ +  - ][ +  - ]:          1 :             for(; oIter != aoMapDimensions.end(); ++oIter )
         [ +  - ][ +  + ]
    1961                 :            :             {
    1962                 :            :                 osURLFeatureInfoTemplate = Replace(osURLFeatureInfoTemplate,
    1963         [ +  - ]:          1 :                                             CPLSPrintf("{%s}", oIter->first.c_str()),
    1964 [ +  - ][ +  - ]:          1 :                                             oIter->second);
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1965                 :            :             }
    1966                 :            :         }
    1967         [ +  - ]:          1 :         poDS->osURLFeatureInfoTemplate = osURLFeatureInfoTemplate;
    1968                 :            : 
    1969                 :            :         // Build all TMS datasets, wrapped in VRT datasets
    1970         [ +  + ]:          1 :         for(int i=nMaxZoomLevel;i>=0;i--)
    1971                 :            :         {
    1972         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM[i];
    1973                 :          1 :             int nRasterXSize = int(0.5 + poDS->nRasterXSize / oTM.dfPixelSize * poDS->adfGT[1]);
    1974                 :          1 :             int nRasterYSize = int(0.5 + poDS->nRasterYSize / oTM.dfPixelSize * poDS->adfGT[1]);
    1975 [ +  - ][ +  + ]:          1 :             if( !poDS->apoDatasets.empty() &&
         [ +  - ][ +  + ]
                 [ +  + ]
    1976                 :            :                 (nRasterXSize < 128 || nRasterYSize < 128) )
    1977                 :            :             {
    1978                 :          1 :                 break;
    1979                 :            :             }
    1980 [ +  - ][ +  - ]:          1 :             CPLString osURL(Replace(osURLTileTemplate, "{TileMatrix}", oTM.osIdentifier));
    1981                 :            : 
    1982                 :          1 :             double dfTileWidthUnits = oTM.dfPixelSize * oTM.nTileWidth;
    1983                 :          1 :             double dfTileHeightUnits = oTM.dfPixelSize * oTM.nTileHeight;
    1984                 :            : 
    1985                 :            :             // Compute the shift in terms of tiles between AOI and TM origin
    1986                 :          1 :             int nTileX = (int)(floor(poDS->adfGT[0] - oTM.dfTLX + 1e-10) / dfTileWidthUnits);
    1987                 :          1 :             int nTileY = (int)(floor(oTM.dfTLY - poDS->adfGT[3] + 1e-10) / dfTileHeightUnits);
    1988                 :            : 
    1989                 :            :             // Compute extent of this zoom level slightly larger than the AOI and
    1990                 :            :             // aligned on tile boundaries at this TM
    1991                 :          1 :             double dfULX = oTM.dfTLX + nTileX * dfTileWidthUnits;
    1992                 :          1 :             double dfULY = oTM.dfTLY - nTileY * dfTileHeightUnits;
    1993                 :          1 :             double dfLRX = poDS->adfGT[0] + poDS->nRasterXSize * poDS->adfGT[1];
    1994                 :          1 :             double dfLRY = poDS->adfGT[3] + poDS->nRasterYSize * poDS->adfGT[5];
    1995                 :          1 :             dfLRX = dfULX + ceil((dfLRX - dfULX) / dfTileWidthUnits - 1e-10) * dfTileWidthUnits;
    1996                 :          1 :             dfLRY = dfULY + floor((dfLRY - dfULY) / dfTileHeightUnits + 1e-10) * dfTileHeightUnits;
    1997                 :            : 
    1998                 :          1 :             int nSizeX = int(0.5+(dfLRX - dfULX) / oTM.dfPixelSize);
    1999                 :          1 :             int nSizeY = int(0.5+(dfULY - dfLRY) / oTM.dfPixelSize);
    2000                 :            : 
    2001                 :          1 :             double dfDateLineX = oTM.dfTLX + oTM.nMatrixWidth * dfTileWidthUnits;
    2002                 :          1 :             int nSizeX1 = int(0.5+(dfDateLineX - dfULX) / oTM.dfPixelSize);
    2003                 :          1 :             int nSizeX2 = int(0.5+(dfLRX - dfDateLineX) / oTM.dfPixelSize);
    2004 [ +  + ][ -  + ]:          1 :             if( bExtendBeyondDateLine && dfDateLineX > dfLRX )
    2005                 :            :             {
    2006         [ #  # ]:          0 :                 CPLDebug("WMTS", "ExtendBeyondDateLine ignored in that case");
    2007                 :          0 :                 bExtendBeyondDateLine = FALSE;
    2008                 :            :             }
    2009                 :            : 
    2010                 :            : #define WMS_TMS_TEMPLATE \
    2011                 :            :     "<GDAL_WMS>" \
    2012                 :            :     "<Service name=\"TMS\">" \
    2013                 :            :     "    <ServerUrl>%s</ServerUrl>" \
    2014                 :            :     "</Service>" \
    2015                 :            :     "<DataWindow>" \
    2016                 :            :     "    <UpperLeftX>%.16g</UpperLeftX>" \
    2017                 :            :     "    <UpperLeftY>%.16g</UpperLeftY>" \
    2018                 :            :     "    <LowerRightX>%.16g</LowerRightX>" \
    2019                 :            :     "    <LowerRightY>%.16g</LowerRightY>" \
    2020                 :            :     "    <TileLevel>0</TileLevel>" \
    2021                 :            :     "    <TileX>%d</TileX>" \
    2022                 :            :     "    <TileY>%d</TileY>" \
    2023                 :            :     "    <SizeX>%d</SizeX>" \
    2024                 :            :     "    <SizeY>%d</SizeY>" \
    2025                 :            :     "    <YOrigin>top</YOrigin>" \
    2026                 :            :     "</DataWindow>" \
    2027                 :            :     "<BlockSizeX>%d</BlockSizeX>" \
    2028                 :            :     "<BlockSizeY>%d</BlockSizeY>" \
    2029                 :            :     "<BandsCount>%d</BandsCount>" \
    2030                 :            :     "%s" \
    2031                 :            : "</GDAL_WMS>"
    2032                 :            : 
    2033                 :            :             CPLString osStr(CPLSPrintf( WMS_TMS_TEMPLATE,
    2034                 :            :                 WMTSEscapeXML(osURL).c_str(),
    2035                 :            :                 dfULX, dfULY, (bExtendBeyondDateLine) ? dfDateLineX : dfLRX, dfLRY,
    2036                 :            :                 nTileX, nTileY, (bExtendBeyondDateLine) ? nSizeX1 : nSizeX, nSizeY,
    2037                 :            :                 oTM.nTileWidth, oTM.nTileHeight, nBands,
    2038 [ +  - ][ +  + ]:          1 :                 osOtherXML.c_str()));
         [ +  + ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    2039                 :            :             GDALDataset* poWMSDS = (GDALDataset*)GDALOpenEx(
    2040 [ +  - ][ +  - ]:          1 :                 osStr, GDAL_OF_RASTER | GDAL_OF_SHARED | GDAL_OF_VERBOSE_ERROR, NULL, NULL, NULL);
    2041         [ -  + ]:          1 :             if( poWMSDS == NULL )
    2042                 :            :             {
    2043         [ #  # ]:          0 :                 CPLDestroyXMLNode(psXML);
    2044 [ #  # ][ #  # ]:          0 :                 delete poDS;
    2045                 :          0 :                 return NULL;
    2046                 :            :             }
    2047                 :            : 
    2048         [ +  - ]:          1 :             VRTDatasetH hVRTDS = VRTCreate( nRasterXSize, nRasterYSize );
    2049         [ +  + ]:          1 :             for(int iBand=1;iBand<=nBands;iBand++)
    2050                 :            :             {
    2051         [ +  - ]:          1 :                 VRTAddBand( hVRTDS, GDT_Byte, NULL );
    2052                 :            :             }
    2053                 :            : 
    2054                 :            :             int nSrcXOff, nSrcYOff, nDstXOff, nDstYOff;
    2055                 :            : 
    2056                 :          1 :             nSrcXOff = 0;
    2057                 :          1 :             nDstXOff = (int)(0.5 + (dfULX - poDS->adfGT[0]) / oTM.dfPixelSize);
    2058                 :            : 
    2059                 :          1 :             nSrcYOff = 0;
    2060                 :          1 :             nDstYOff = (int)(0.5 + (poDS->adfGT[3] - dfULY) / oTM.dfPixelSize);
    2061                 :            : 
    2062         [ +  + ]:          1 :             if( bExtendBeyondDateLine )
    2063                 :            :             {
    2064                 :            :                 int nSrcXOff2, nDstXOff2;
    2065                 :            : 
    2066                 :          1 :                 nSrcXOff2 = 0;
    2067                 :          1 :                 nDstXOff2 = (int)(0.5 + (dfDateLineX - poDS->adfGT[0]) / oTM.dfPixelSize);
    2068                 :            : 
    2069                 :            :                 osStr = CPLSPrintf( WMS_TMS_TEMPLATE,
    2070                 :            :                     WMTSEscapeXML(osURL).c_str(),
    2071                 :            :                     -dfDateLineX, dfULY, dfLRX - 2 * dfDateLineX, dfLRY,
    2072                 :            :                     0, nTileY, nSizeX2, nSizeY,
    2073                 :            :                     oTM.nTileWidth, oTM.nTileHeight, nBands,
    2074 [ +  - ][ +  - ]:          1 :                     osOtherXML.c_str());
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    2075                 :            : 
    2076                 :            :                 GDALDataset* poWMSDS2 = (GDALDataset*)GDALOpenEx(
    2077 [ +  - ][ +  - ]:          1 :                     osStr, GDAL_OF_RASTER | GDAL_OF_SHARED, NULL, NULL, NULL);
    2078         [ -  + ]:          1 :                 CPLAssert(poWMSDS2);
    2079                 :            : 
    2080         [ +  + ]:          1 :                 for(int iBand=1;iBand<=nBands;iBand++)
    2081                 :            :                 {
    2082                 :            :                     VRTSourcedRasterBandH hVRTBand =
    2083         [ +  - ]:          1 :                         (VRTSourcedRasterBandH) GDALGetRasterBand(hVRTDS, iBand);
    2084                 :            :                     VRTAddSimpleSource( hVRTBand, GDALGetRasterBand(poWMSDS, iBand),
    2085                 :            :                                         nSrcXOff, nSrcYOff, nSizeX1, nSizeY,
    2086                 :            :                                         nDstXOff, nDstYOff, nSizeX1, nSizeY,
    2087 [ +  - ][ +  - ]:          1 :                                         "NEAR", VRT_NODATA_UNSET);
    2088                 :            :                     VRTAddSimpleSource( hVRTBand, GDALGetRasterBand(poWMSDS2, iBand),
    2089                 :            :                                         nSrcXOff2, nSrcYOff, nSizeX2, nSizeY,
    2090                 :            :                                         nDstXOff2, nDstYOff, nSizeX2, nSizeY,
    2091 [ +  - ][ +  - ]:          1 :                                         "NEAR", VRT_NODATA_UNSET);
    2092                 :            :                 }
    2093                 :            : 
    2094         [ +  - ]:          1 :                 poWMSDS2->Dereference();
    2095                 :            :             }
    2096                 :            :             else
    2097                 :            :             {
    2098         [ +  + ]:          1 :                 for(int iBand=1;iBand<=nBands;iBand++)
    2099                 :            :                 {
    2100                 :            :                     VRTSourcedRasterBandH hVRTBand =
    2101         [ +  - ]:          1 :                         (VRTSourcedRasterBandH) GDALGetRasterBand(hVRTDS, iBand);
    2102                 :            :                     VRTAddSimpleSource( hVRTBand, GDALGetRasterBand(poWMSDS, iBand),
    2103                 :            :                                         nSrcXOff, nSrcYOff, nSizeX, nSizeY,
    2104                 :            :                                         nDstXOff, nDstYOff, nSizeX, nSizeY,
    2105 [ +  - ][ +  - ]:          1 :                                         "NEAR", VRT_NODATA_UNSET);
    2106                 :            :                 }
    2107                 :            :             }
    2108                 :            : 
    2109         [ +  - ]:          1 :             poWMSDS->Dereference();
    2110                 :            : 
    2111 [ +  - ][ -  + ]:          1 :             poDS->apoDatasets.push_back((GDALDataset*)hVRTDS);
    2112 [ +  - ][ +  - ]:          1 :         }
                 [ -  + ]
    2113                 :            : 
    2114 [ +  - ][ -  + ]:          1 :         if( poDS->apoDatasets.empty() )
    2115                 :            :         {
    2116         [ #  # ]:          0 :             CPLError(CE_Failure, CPLE_AppDefined, "No zoom level found");
    2117         [ #  # ]:          0 :             CPLDestroyXMLNode(psXML);
    2118 [ #  # ][ #  # ]:          0 :             delete poDS;
    2119                 :          0 :             return NULL;
    2120                 :            :         }
    2121                 :            : 
    2122         [ +  - ]:          1 :         poDS->SetMetadataItem("INTERLEAVE", "PIXEL", "IMAGE_STRUCTURE");
    2123         [ +  + ]:          1 :         for(int i=0;i<nBands;i++)
    2124 [ +  - ][ +  - ]:          1 :             poDS->SetBand(i+1, new WMTSBand(poDS,i+1));
                 [ +  - ]
    2125                 :            : 
    2126 [ +  - ][ +  - ]:          1 :         poDS->osXML = "<GDAL_WMTS>\n";
                 [ +  - ]
    2127                 :            :         poDS->osXML += "  <GetCapabilitiesUrl>" +
    2128                 :            :                      WMTSEscapeXML(osGetCapabilitiesURL) +
    2129 [ +  - ][ +  - ]:          1 :                      "</GetCapabilitiesUrl>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2130 [ +  - ][ +  + ]:          1 :         if( !osSelectLayer.empty() )
    2131 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <Layer>" + WMTSEscapeXML(osSelectLayer) + "</Layer>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2132 [ +  - ][ +  + ]:          1 :         if( !osSelectStyle.empty() )
    2133 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <Style>" + WMTSEscapeXML(osSelectStyle) + "</Style>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2134 [ +  - ][ +  + ]:          1 :         if( !osSelectTMS.empty() )
    2135 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <TileMatrixSet>" + WMTSEscapeXML(osSelectTMS) + "</TileMatrixSet>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2136 [ +  - ][ +  + ]:          1 :         if( !osMaxTileMatrixIdentifier.empty() )
    2137 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <TileMatrix>" + WMTSEscapeXML(osMaxTileMatrixIdentifier) + "</TileMatrix>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2138         [ +  + ]:          1 :         if( nUserMaxZoomLevel >= 0 )
    2139 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <ZoomLevel>" + CPLString().Printf("%d", nUserMaxZoomLevel) + "</ZoomLevel>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2140 [ +  + ][ +  - ]:          1 :         if( nCountTileFormat > 1 && !osSelectTileFormat.empty() )
         [ +  - ][ +  + ]
    2141 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <Format>" + WMTSEscapeXML(osSelectTileFormat) + "</Format>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2142 [ -  + ][ #  # ]:          1 :         if( nCountInfoFormat > 1 && !osSelectInfoFormat.empty() )
         [ #  # ][ -  + ]
    2143 [ #  # ][ #  # ]:          0 :             poDS->osXML += "  <InfoFormat>" + WMTSEscapeXML(osSelectInfoFormat) + "</InfoFormat>\n";
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
    2144         [ +  - ]:          1 :         poDS->osXML += "  <DataWindow>\n";
    2145                 :            :         poDS->osXML += CPLSPrintf("    <UpperLeftX>%.16g</UpperLeftX>\n",
    2146 [ +  - ][ +  - ]:          1 :                              poDS->adfGT[0]);
    2147                 :            :         poDS->osXML += CPLSPrintf("    <UpperLeftY>%.16g</UpperLeftY>\n",
    2148 [ +  - ][ +  - ]:          1 :                              poDS->adfGT[3]);
    2149                 :            :         poDS->osXML += CPLSPrintf("    <LowerRightX>%.16g</LowerRightX>\n",
    2150 [ +  - ][ +  - ]:          1 :                              poDS->adfGT[0] +  poDS->adfGT[1] *  poDS->nRasterXSize);
    2151                 :            :         poDS->osXML += CPLSPrintf("    <LowerRightY>%.16g</LowerRightY>\n",
    2152 [ +  - ][ +  - ]:          1 :                              poDS->adfGT[3] +  poDS->adfGT[5] *  poDS->nRasterYSize);
    2153         [ +  - ]:          1 :         poDS->osXML += "  </DataWindow>\n";
    2154         [ +  + ]:          1 :         if( bExtendBeyondDateLine )
    2155         [ +  - ]:          1 :             poDS->osXML += "  <ExtendBeyondDateLine>true</ExtendBeyondDateLine>\n";
    2156 [ +  - ][ +  - ]:          1 :         poDS->osXML += CPLSPrintf("  <BandsCount>%d</BandsCount>\n", nBands);
    2157         [ +  - ]:          1 :         poDS->osXML += "  <Cache />\n";
    2158         [ +  - ]:          1 :         poDS->osXML += "  <UnsafeSSL>true</UnsafeSSL>\n";
    2159         [ +  - ]:          1 :         poDS->osXML += "  <ZeroBlockHttpCodes>204,404</ZeroBlockHttpCodes>\n";
    2160         [ +  - ]:          1 :         poDS->osXML += "  <ZeroBlockOnServerException>true</ZeroBlockOnServerException>\n";
    2161 [ +  - ][ +  - ]:          1 :         poDS->osXML += "</GDAL_WMTS>\n";
                 [ +  + ]
    2162                 :            :     }
    2163                 :            : 
    2164         [ +  - ]:          1 :     CPLDestroyXMLNode(psXML);
    2165                 :            : 
    2166         [ +  - ]:          1 :     poDS->SetPamFlags(0);
    2167 [ +  - ][ +  - ]:          1 :     return poDS;
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    2168                 :            : }
    2169                 :            : /************************************************************************/
    2170                 :            : /*                             CreateCopy()                             */
    2171                 :            : /************************************************************************/
    2172                 :            : 
    2173                 :          1 : GDALDataset *WMTSDataset::CreateCopy( const char * pszFilename,
    2174                 :            :                                          GDALDataset *poSrcDS,
    2175                 :            :                                          CPL_UNUSED int bStrict,
    2176                 :            :                                          CPL_UNUSED char ** papszOptions,
    2177                 :            :                                          CPL_UNUSED GDALProgressFunc pfnProgress,
    2178                 :            :                                          CPL_UNUSED void * pProgressData )
    2179                 :            : {
    2180   [ +  -  +  + ]:          1 :     if( poSrcDS->GetDriver() == NULL ||
                 [ +  + ]
    2181                 :          1 :         poSrcDS->GetDriver() != GDALGetDriverByName("WMTS") )
    2182                 :            :     {
    2183                 :            :         CPLError(CE_Failure, CPLE_NotSupported,
    2184                 :          1 :                  "Source dataset must be a WMTS dataset");
    2185                 :          1 :         return NULL;
    2186                 :            :     }
    2187                 :            : 
    2188                 :          1 :     const char* pszXML = poSrcDS->GetMetadataItem("XML", "WMTS");
    2189         [ -  + ]:          1 :     if (pszXML == NULL)
    2190                 :            :     {
    2191                 :            :         CPLError(CE_Failure, CPLE_AppDefined,
    2192                 :          0 :                  "Cannot get XML definition of source WMTS dataset");
    2193                 :          0 :         return NULL;
    2194                 :            :     }
    2195                 :            : 
    2196                 :          1 :     VSILFILE* fp = VSIFOpenL(pszFilename, "wb");
    2197         [ -  + ]:          1 :     if (fp == NULL)
    2198                 :          0 :         return NULL;
    2199                 :            : 
    2200                 :          1 :     VSIFWriteL(pszXML, 1, strlen(pszXML), fp);
    2201                 :          1 :     VSIFCloseL(fp);
    2202                 :            : 
    2203                 :          1 :     GDALOpenInfo oOpenInfo(pszFilename, GA_ReadOnly);
    2204         [ +  - ]:          1 :     return Open(&oOpenInfo);
    2205                 :            : }
    2206                 :            : 
    2207                 :            : /************************************************************************/
    2208                 :            : /*                       GDALRegister_WMTS()                            */
    2209                 :            : /************************************************************************/
    2210                 :            : 
    2211                 :          1 : void GDALRegister_WMTS()
    2212                 :            : 
    2213                 :            : {
    2214         [ -  + ]:          1 :     if( !GDAL_CHECK_VERSION( "WMTS driver" ) )
    2215                 :          0 :         return;
    2216                 :            : 
    2217         [ +  + ]:          1 :     if( GDALGetDriverByName( "WMTS" ) != NULL )
    2218                 :          1 :         return;
    2219                 :            : 
    2220         [ +  - ]:          1 :     GDALDriver *poDriver = new GDALDriver();
    2221                 :            : 
    2222                 :          1 :     poDriver->SetDescription( "WMTS" );
    2223                 :          1 :     poDriver->SetMetadataItem( GDAL_DCAP_RASTER, "YES" );
    2224                 :          1 :     poDriver->SetMetadataItem( GDAL_DMD_LONGNAME, "OGC Web Mab Tile Service" );
    2225                 :          1 :     poDriver->SetMetadataItem( GDAL_DMD_HELPTOPIC, "frmt_wmts.html" );
    2226                 :            : 
    2227                 :          1 :     poDriver->SetMetadataItem( GDAL_DMD_CONNECTION_PREFIX, "WMTS:" );
    2228                 :            : 
    2229                 :          1 :     poDriver->SetMetadataItem( GDAL_DCAP_VIRTUALIO, "YES" );
    2230                 :            : 
    2231                 :            :     poDriver->SetMetadataItem( GDAL_DMD_OPENOPTIONLIST,
    2232                 :            : "<OpenOptionList>"
    2233                 :            : "  <Option name='URL' type='string' description='URL that points to GetCapabilities response' required='YES'/>"
    2234                 :            : "  <Option name='LAYER' type='string' description='Layer identifier'/>"
    2235                 :            : "  <Option name='TILEMATRIXSET' alias='TMS' type='string' description='Tile matrix set identifier'/>"
    2236                 :            : "  <Option name='TILEMATRIX' type='string' description='Tile matrix identifier of maximum zoom level. Exclusive with ZOOM_LEVEL.'/>"
    2237                 :            : "  <Option name='ZOOM_LEVEL' alias='ZOOMLEVEL' type='int' description='Maximum zoom level. Exclusive with TILEMATRIX.'/>"
    2238                 :            : "  <Option name='STYLE' type='string' description='Style identifier'/>"
    2239                 :            : "  <Option name='EXTENDBEYONDDATELINE' type='boolean' description='Whether to enable extend-beyond-dateline behaviour' default='NO'/>"
    2240                 :            : "  <Option name='EXTENT_METHOD' type='string-select' description='How the raster extent is computed' default='AUTO'>"
    2241                 :            : "       <Value>AUTO</Value>"
    2242                 :            : "       <Value>LAYER_BBOX</Value>"
    2243                 :            : "       <Value>TILE_MATRIX_SET</Value>"
    2244                 :            : "       <Value>MOST_PRECISE_TILE_MATRIX</Value>"
    2245                 :            : "  </Option>"
    2246                 :          1 : "</OpenOptionList>");
    2247                 :            : 
    2248                 :          1 :     poDriver->pfnOpen = WMTSDataset::Open;
    2249                 :          1 :     poDriver->pfnIdentify = WMTSDataset::Identify;
    2250                 :          1 :     poDriver->pfnCreateCopy = WMTSDataset::CreateCopy;
    2251                 :            : 
    2252                 :          1 :     GetGDALDriverManager()->RegisterDriver( poDriver );
    2253                 :            : }

Generated by: LCOV version 1.9