LCOV - code coverage report
Current view: top level - frmts/wmts - wmtsdataset.cpp (source / functions) Hit Total Coverage
Test: gdal_filtered.info Lines: 1021 1124 90.8 %
Date: XXXX-XX-XX Functions: 40 42 95.2 %
Branches: 1694 3169 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( oTM.dfPixelSize <= 0.0 )
     723                 :            :             {
     724                 :            :                 CPLError(CE_Failure, CPLE_AppDefined,
     725         [ #  # ]:          0 :                          "Invalid ScaleDenominator");
     726                 :          0 :                 return FALSE;
     727                 :            :             }
     728 [ +  - ][ +  + ]:          1 :             if( oTMS.oSRS.IsGeographic() )
     729                 :          1 :                 oTM.dfPixelSize *= WMTS_WGS84_DEG_PER_METER;
     730         [ +  - ]:          1 :             double dfVal1 = CPLAtof(pszTopLeftCorner);
     731         [ +  - ]:          1 :             double dfVal2 = CPLAtof(strchr(pszTopLeftCorner, ' ')+1);
     732 [ +  + ][ -  + ]:          1 :             if( !bSwap ||
                 [ #  # ]
     733                 :            :                 /* Hack for http://osm.geobretagne.fr/gwc01/service/wmts?request=getcapabilities */
     734                 :          1 :                 ( STARTS_WITH_CI(l_pszIdentifier, "EPSG:4326:") &&
     735                 :            :                   dfVal1 == -180.0 ) )
     736                 :            :             {
     737                 :          1 :                 oTM.dfTLX = dfVal1;
     738                 :          1 :                 oTM.dfTLY = dfVal2;
     739                 :            :             }
     740                 :            :             else
     741                 :            :             {
     742                 :          1 :                 oTM.dfTLX = dfVal2;
     743                 :          1 :                 oTM.dfTLY = dfVal1;
     744                 :            :             }
     745                 :          1 :             oTM.nTileWidth = atoi(pszTileWidth);
     746                 :          1 :             oTM.nTileHeight = atoi(pszTileHeight);
     747 [ +  - ][ +  - ]:          1 :             if( oTM.nTileWidth <= 0 || oTM.nTileWidth > 4096 ||
         [ +  - ][ -  + ]
     748                 :            :                 oTM.nTileHeight <= 0 || oTM.nTileHeight > 4096 )
     749                 :            :             {
     750                 :            :                 CPLError(CE_Failure, CPLE_AppDefined,
     751         [ #  # ]:          0 :                          "Invalid TileWidth/TileHeight element");
     752                 :          0 :                 return FALSE;
     753                 :            :             }
     754                 :          1 :             oTM.nMatrixWidth = atoi(pszMatrixWidth);
     755                 :          1 :             oTM.nMatrixHeight = atoi(pszMatrixHeight);
     756                 :            :             // http://datacarto.geonormandie.fr/mapcache/wmts?SERVICE=WMTS&REQUEST=GetCapabilities
     757                 :            :             // has a TileMatrix 0 with MatrixWidth = MatrixHeight = 0
     758 [ +  - ][ -  + ]:          1 :             if( oTM.nMatrixWidth < 1 || oTM.nMatrixHeight < 1 )
     759                 :          0 :                 continue;
     760         [ +  - ]:          1 :             oTMS.aoTM.push_back(oTM);
     761 [ +  + ][ +  - ]:          1 :             if( (nMaxZoomLevel >= 0 && static_cast<int>(oTMS.aoTM.size())-1
         [ +  + ][ +  + ]
         [ +  + ][ +  + ]
     762                 :            :                                                         == nMaxZoomLevel) ||
     763         [ +  - ]:          1 :                 (!osMaxTileMatrixIdentifier.empty() &&
     764         [ +  - ]:          1 :                  EQUAL(osMaxTileMatrixIdentifier, l_pszIdentifier)) )
     765                 :            :             {
     766                 :          1 :                 bFoundTileMatrix = true;
     767                 :            :                 break;
     768                 :            :             }
     769      [ +  +  - ]:          1 :         }
     770 [ +  + ][ +  + ]:          1 :         if( nMaxZoomLevel >= 0 && !bFoundTileMatrix )
     771                 :            :         {
     772                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
     773                 :            :                      "Cannot find TileMatrix of zoom level %d in TileMatrixSet '%s'",
     774                 :            :                      nMaxZoomLevel,
     775                 :          1 :                      osIdentifier.c_str());
     776                 :          1 :             return FALSE;
     777                 :            :         }
     778 [ +  + ][ +  + ]:          1 :         if( !osMaxTileMatrixIdentifier.empty() && !bFoundTileMatrix )
                 [ +  + ]
     779                 :            :         {
     780                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
     781                 :            :                      "Cannot find TileMatrix '%s' in TileMatrixSet '%s'",
     782                 :            :                      osMaxTileMatrixIdentifier.c_str(),
     783                 :          1 :                      osIdentifier.c_str());
     784                 :          1 :             return FALSE;
     785                 :            :         }
     786         [ +  + ]:          1 :         if( oTMS.aoTM.empty() )
     787                 :            :         {
     788                 :            :             CPLError(CE_Failure, CPLE_AppDefined, "Cannot find TileMatrix in TileMatrixSet '%s'",
     789                 :          1 :                      osIdentifier.c_str());
     790                 :          1 :             return FALSE;
     791                 :            :         }
     792                 :          1 :         return TRUE;
     793                 :            :     }
     794                 :            :     CPLError(CE_Failure, CPLE_AppDefined, "Cannot find TileMatrixSet '%s'",
     795                 :          1 :              osIdentifier.c_str());
     796                 :          1 :     return FALSE;
     797                 :            : }
     798                 :            : 
     799                 :            : /************************************************************************/
     800                 :            : /*                              ReadTMLimits()                          */
     801                 :            : /************************************************************************/
     802                 :            : 
     803                 :          1 : int WMTSDataset::ReadTMLimits(CPLXMLNode* psTMSLimits,
     804                 :            :                               std::map<CPLString, WMTSTileMatrixLimits>& aoMapTileMatrixLimits)
     805                 :            : {
     806         [ +  + ]:          1 :     for(CPLXMLNode* psIter = psTMSLimits->psChild; psIter; psIter = psIter->psNext)
     807                 :            :     {
     808 [ +  - ][ -  + ]:          1 :         if( psIter->eType != CXT_Element || strcmp(psIter->pszValue, "TileMatrixLimits") != 0 )
     809                 :          0 :             continue;
     810                 :          1 :         WMTSTileMatrixLimits oTMLimits;
     811         [ +  - ]:          1 :         const char* pszTileMatrix = CPLGetXMLValue(psIter, "TileMatrix", NULL);
     812         [ +  - ]:          1 :         const char* pszMinTileRow = CPLGetXMLValue(psIter, "MinTileRow", NULL);
     813         [ +  - ]:          1 :         const char* pszMaxTileRow = CPLGetXMLValue(psIter, "MaxTileRow", NULL);
     814         [ +  - ]:          1 :         const char* pszMinTileCol = CPLGetXMLValue(psIter, "MinTileCol", NULL);
     815         [ +  - ]:          1 :         const char* pszMaxTileCol = CPLGetXMLValue(psIter, "MaxTileCol", NULL);
     816 [ +  - ][ +  - ]:          1 :         if( pszTileMatrix == NULL ||
         [ +  - ][ +  - ]
                 [ -  + ]
     817                 :            :             pszMinTileRow == NULL || pszMaxTileRow == NULL ||
     818                 :            :             pszMinTileCol == NULL || pszMaxTileCol == NULL )
     819                 :            :         {
     820                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
     821         [ #  # ]:          0 :                      "Missing required element in TileMatrixLimits element");
     822                 :          0 :             return FALSE;
     823                 :            :         }
     824 [ +  - ][ +  - ]:          1 :         oTMLimits.osIdentifier = pszTileMatrix;
                 [ +  - ]
     825                 :          1 :         oTMLimits.nMinTileRow = atoi(pszMinTileRow);
     826                 :          1 :         oTMLimits.nMaxTileRow = atoi(pszMaxTileRow);
     827                 :          1 :         oTMLimits.nMinTileCol = atoi(pszMinTileCol);
     828                 :          1 :         oTMLimits.nMaxTileCol = atoi(pszMaxTileCol);
     829 [ +  - ][ +  - ]:          1 :         aoMapTileMatrixLimits[pszTileMatrix] = oTMLimits;
         [ +  - ][ +  - ]
                 [ +  - ]
     830                 :          1 :     }
     831                 :          1 :     return TRUE;
     832                 :            : }
     833                 :            : 
     834                 :            : /************************************************************************/
     835                 :            : /*                               Replace()                              */
     836                 :            : /************************************************************************/
     837                 :            : 
     838                 :          1 : CPLString WMTSDataset::Replace(const CPLString& osStr, const char* pszOld,
     839                 :            :                                const char* pszNew)
     840                 :            : {
     841                 :          1 :     size_t nPos = osStr.ifind(pszOld);
     842         [ +  + ]:          1 :     if( nPos == std::string::npos )
     843                 :          1 :         return osStr;
     844 [ +  - ][ +  - ]:          1 :     CPLString osRet(osStr.substr(0, nPos));
     845         [ +  - ]:          1 :     osRet += pszNew;
     846 [ +  - ][ +  - ]:          1 :     osRet += osStr.substr(nPos + strlen(pszOld));
                 [ +  - ]
     847         [ +  - ]:          1 :     return osRet;
     848                 :            : }
     849                 :            : 
     850                 :            : /************************************************************************/
     851                 :            : /*                       GetCapabilitiesResponse()                      */
     852                 :            : /************************************************************************/
     853                 :            : 
     854                 :          1 : CPLXMLNode* WMTSDataset::GetCapabilitiesResponse(const CPLString& osFilename,
     855                 :            :                                                  char** papszHTTPOptions)
     856                 :            : {
     857                 :            :     CPLXMLNode* psXML;
     858                 :            :     VSIStatBufL sStat;
     859         [ +  + ]:          1 :     if( VSIStatL(osFilename, &sStat) == 0 )
     860                 :          1 :         psXML = CPLParseXMLFile(osFilename);
     861                 :            :     else
     862                 :            :     {
     863                 :          1 :         CPLHTTPResult* psResult = CPLHTTPFetch(osFilename, papszHTTPOptions);
     864         [ -  + ]:          1 :         if( psResult == NULL )
     865                 :          0 :             return NULL;
     866         [ +  - ]:          1 :         if( psResult->pabyData == NULL )
     867                 :            :         {
     868                 :          1 :             CPLHTTPDestroyResult(psResult);
     869                 :          1 :             return NULL;
     870                 :            :         }
     871                 :          0 :         psXML = CPLParseXMLString((const char*)psResult->pabyData);
     872                 :          0 :         CPLHTTPDestroyResult(psResult);
     873                 :            :     }
     874                 :          1 :     return psXML;
     875                 :            : }
     876                 :            : 
     877                 :            : /************************************************************************/
     878                 :            : /*                          WMTSAddOtherXML()                           */
     879                 :            : /************************************************************************/
     880                 :            : 
     881                 :          1 : static void WMTSAddOtherXML(CPLXMLNode* psRoot, const char* pszElement,
     882                 :            :                             CPLString& osOtherXML)
     883                 :            : {
     884                 :          1 :     CPLXMLNode* psElement = CPLGetXMLNode(psRoot, pszElement);
     885         [ +  + ]:          1 :     if( psElement )
     886                 :            :     {
     887                 :          1 :         CPLXMLNode* psNext = psElement->psNext;
     888                 :          1 :         psElement->psNext = NULL;
     889                 :          1 :         char* pszTmp = CPLSerializeXMLTree(psElement);
     890                 :          1 :         osOtherXML += pszTmp;
     891                 :          1 :         CPLFree(pszTmp);
     892                 :          1 :         psElement->psNext = psNext;
     893                 :            :     }
     894                 :          1 : }
     895                 :            : 
     896                 :            : /************************************************************************/
     897                 :            : /*                          GetOperationKVPURL()                        */
     898                 :            : /************************************************************************/
     899                 :            : 
     900                 :          1 : CPLString WMTSDataset::GetOperationKVPURL(CPLXMLNode* psXML,
     901                 :            :                                           const char* pszOperation)
     902                 :            : {
     903                 :          1 :     CPLString osRet;
     904         [ +  - ]:          1 :     CPLXMLNode* psOM = CPLGetXMLNode(psXML, "=Capabilities.OperationsMetadata");
     905 [ +  + ][ +  + ]:          1 :     for(CPLXMLNode* psIter = psOM ? psOM->psChild : NULL;
     906                 :            :         psIter != NULL; psIter = psIter->psNext)
     907                 :            :     {
     908 [ +  - ][ +  - ]:          1 :         if( psIter->eType != CXT_Element ||
         [ +  + ][ +  + ]
     909                 :          1 :             strcmp(psIter->pszValue, "Operation") != 0 ||
     910         [ +  - ]:          1 :             !EQUAL(CPLGetXMLValue(psIter, "name", ""), pszOperation) )
     911                 :            :         {
     912                 :          1 :             continue;
     913                 :            :         }
     914         [ +  - ]:          1 :         CPLXMLNode* psHTTP = CPLGetXMLNode(psIter, "DCP.HTTP");
     915 [ +  - ][ +  + ]:          1 :         for(CPLXMLNode* psGet = psHTTP ? psHTTP->psChild : NULL;
     916                 :            :                 psGet != NULL; psGet = psGet->psNext)
     917                 :            :         {
     918 [ +  - ][ -  + ]:          1 :             if( psGet->eType != CXT_Element ||
     919                 :          1 :                 strcmp(psGet->pszValue, "Get") != 0 )
     920                 :            :             {
     921                 :          0 :                 continue;
     922                 :            :             }
     923 [ +  - ][ -  + ]:          1 :             if( !EQUAL(CPLGetXMLValue(psGet, "Constraint.AllowedValues.Value", "KVP"), "KVP") )
     924                 :          0 :                 continue;
     925 [ +  - ][ +  - ]:          1 :             osRet = CPLGetXMLValue(psGet, "href", "");
         [ +  - ][ +  - ]
     926                 :            :         }
     927                 :            :     }
     928                 :          1 :     return osRet;
     929                 :            : }
     930                 :            : 
     931                 :            : /************************************************************************/
     932                 :            : /*                           BuildHTTPRequestOpts()                     */
     933                 :            : /************************************************************************/
     934                 :            : 
     935                 :          1 : char** WMTSDataset::BuildHTTPRequestOpts(CPLString osOtherXML)
     936                 :            : {
     937 [ +  - ][ +  - ]:          1 :     osOtherXML = "<Root>" + osOtherXML + "</Root>";
         [ +  - ][ +  - ]
                 [ +  - ]
     938                 :          1 :     CPLXMLNode* psXML = CPLParseXMLString(osOtherXML);
     939                 :          1 :     char **http_request_opts = NULL;
     940         [ -  + ]:          1 :     if (CPLGetXMLValue(psXML, "Timeout", NULL)) {
     941                 :          0 :         CPLString optstr;
     942 [ #  # ][ #  # ]:          0 :         optstr.Printf("TIMEOUT=%s", CPLGetXMLValue(psXML, "Timeout", NULL));
     943 [ #  # ][ #  # ]:          0 :         http_request_opts = CSLAddString(http_request_opts, optstr.c_str());
     944                 :            :     }
     945         [ -  + ]:          1 :     if (CPLGetXMLValue(psXML, "UserAgent", NULL)) {
     946                 :          0 :         CPLString optstr;
     947 [ #  # ][ #  # ]:          0 :         optstr.Printf("USERAGENT=%s", CPLGetXMLValue(psXML, "UserAgent", NULL));
     948 [ #  # ][ #  # ]:          0 :         http_request_opts = CSLAddString(http_request_opts, optstr.c_str());
     949                 :            :     }
     950         [ -  + ]:          1 :     if (CPLGetXMLValue(psXML, "Referer", NULL)) {
     951                 :          0 :         CPLString optstr;
     952 [ #  # ][ #  # ]:          0 :         optstr.Printf("REFERER=%s", CPLGetXMLValue(psXML, "Referer", NULL));
     953 [ #  # ][ #  # ]:          0 :         http_request_opts = CSLAddString(http_request_opts, optstr.c_str());
     954                 :            :     }
     955         [ +  + ]:          1 :     if (CPLTestBool(CPLGetXMLValue(psXML, "UnsafeSSL", "false"))) {
     956                 :          1 :         http_request_opts = CSLAddString(http_request_opts, "UNSAFESSL=1");
     957                 :            :     }
     958         [ -  + ]:          1 :     if (CPLGetXMLValue(psXML, "UserPwd", NULL)) {
     959                 :          0 :         CPLString optstr;
     960 [ #  # ][ #  # ]:          0 :         optstr.Printf("USERPWD=%s", CPLGetXMLValue(psXML, "UserPwd", NULL));
     961 [ #  # ][ #  # ]:          0 :         http_request_opts = CSLAddString(http_request_opts, optstr.c_str());
     962                 :            :     }
     963                 :          1 :     CPLDestroyXMLNode(psXML);
     964                 :          1 :     return http_request_opts;
     965                 :            : }
     966                 :            : 
     967                 :            : /************************************************************************/
     968                 :            : /*                                Open()                                */
     969                 :            : /************************************************************************/
     970                 :            : 
     971                 :          1 : GDALDataset* WMTSDataset::Open(GDALOpenInfo* poOpenInfo)
     972                 :            : {
     973         [ +  + ]:          1 :     if (!Identify(poOpenInfo))
     974                 :          1 :         return NULL;
     975                 :            : 
     976                 :          1 :     CPLXMLNode* psXML = NULL;
     977                 :          1 :     CPLString osTileFormat;
     978         [ +  - ]:          1 :     CPLString osInfoFormat;
     979                 :            : 
     980                 :            :     CPLString osGetCapabilitiesURL = CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     981 [ +  - ][ +  - ]:          1 :                                                 "URL", "");
     982                 :            :     CPLString osLayer = CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     983 [ +  - ][ +  - ]:          1 :                                     "LAYER", "");
     984                 :            :     CPLString osTMS = CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     985 [ +  - ][ +  - ]:          1 :                                     "TILEMATRIXSET", "");
     986                 :            :     CPLString osMaxTileMatrixIdentifier = CSLFetchNameValueDef(
     987                 :            :                                     poOpenInfo->papszOpenOptions,
     988 [ +  - ][ +  - ]:          1 :                                     "TILEMATRIX", "");
     989                 :            :     int nUserMaxZoomLevel = atoi(CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     990                 :            :                                     "ZOOM_LEVEL",
     991                 :            :                                     CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     992 [ +  - ][ +  - ]:          1 :                                     "ZOOMLEVEL", "-1")));
     993                 :            :     CPLString osStyle = CSLFetchNameValueDef(poOpenInfo->papszOpenOptions,
     994 [ +  - ][ +  - ]:          1 :                                     "STYLE", "");
     995                 :            : 
     996                 :            :     int bExtendBeyondDateLine =
     997                 :            :         CPLFetchBool(poOpenInfo->papszOpenOptions,
     998         [ +  - ]:          1 :                      "EXTENDBEYONDDATELINE", false);
     999                 :            : 
    1000                 :            :     CPLString osOtherXML = "<Cache />"
    1001                 :            :                      "<UnsafeSSL>true</UnsafeSSL>"
    1002                 :            :                      "<ZeroBlockHttpCodes>204,404</ZeroBlockHttpCodes>"
    1003         [ +  - ]:          1 :                      "<ZeroBlockOnServerException>true</ZeroBlockOnServerException>";
    1004                 :            : 
    1005         [ +  + ]:          1 :     if( STARTS_WITH_CI(poOpenInfo->pszFilename, "WMTS:") )
    1006                 :            :     {
    1007                 :            :         char** papszTokens = CSLTokenizeString2( poOpenInfo->pszFilename + 5,
    1008         [ +  - ]:          1 :                                                  ",", CSLT_HONOURSTRINGS );
    1009 [ +  - ][ +  + ]:          1 :         if( papszTokens && papszTokens[0] )
    1010                 :            :         {
    1011 [ +  - ][ +  - ]:          1 :             osGetCapabilitiesURL = papszTokens[0];
                 [ +  - ]
    1012         [ +  + ]:          1 :             for(char** papszIter = papszTokens+1; *papszIter; papszIter++)
    1013                 :            :             {
    1014                 :          1 :                 char* pszKey = NULL;
    1015         [ +  - ]:          1 :                 const char* pszValue = CPLParseNameValue(*papszIter, &pszKey);
    1016 [ +  - ][ +  - ]:          1 :                 if( pszKey && pszValue )
    1017                 :            :                 {
    1018         [ +  + ]:          1 :                     if( EQUAL(pszKey, "layer") )
    1019 [ +  - ][ +  - ]:          1 :                         osLayer = pszValue;
                 [ +  - ]
    1020         [ +  + ]:          1 :                     else if( EQUAL(pszKey, "tilematrixset") )
    1021 [ +  - ][ +  - ]:          1 :                         osTMS = pszValue;
                 [ +  - ]
    1022         [ +  + ]:          1 :                     else if( EQUAL(pszKey, "tilematrix") )
    1023 [ +  - ][ +  - ]:          1 :                         osMaxTileMatrixIdentifier = pszValue;
                 [ +  - ]
    1024 [ +  + ][ -  + ]:          1 :                     else if( EQUAL(pszKey, "zoom_level") ||
    1025                 :          1 :                              EQUAL(pszKey, "zoomlevel") )
    1026                 :          1 :                         nUserMaxZoomLevel = atoi(pszValue);
    1027         [ +  + ]:          1 :                     else if( EQUAL(pszKey, "style") )
    1028 [ +  - ][ +  - ]:          1 :                         osStyle = pszValue;
                 [ +  - ]
    1029         [ +  - ]:          1 :                     else if( EQUAL(pszKey, "extendbeyonddateline") )
    1030         [ +  - ]:          1 :                         bExtendBeyondDateLine = CPLTestBool(pszValue);
    1031                 :            :                     else
    1032                 :            :                         CPLError(CE_Warning, CPLE_AppDefined,
    1033         [ #  # ]:          1 :                                  "Unknown parameter: %s'", pszKey);
    1034                 :            :                 }
    1035         [ +  - ]:          1 :                 CPLFree(pszKey);
    1036                 :            :             }
    1037                 :            :         }
    1038         [ +  - ]:          1 :         CSLDestroy(papszTokens);
    1039                 :            : 
    1040 [ +  - ][ +  - ]:          1 :         char** papszHTTPOptions = BuildHTTPRequestOpts(osOtherXML);
                 [ +  - ]
    1041         [ +  - ]:          1 :         psXML = GetCapabilitiesResponse(osGetCapabilitiesURL, papszHTTPOptions);
    1042         [ +  - ]:          1 :         CSLDestroy(papszHTTPOptions);
    1043                 :            :     }
    1044                 :            : 
    1045                 :          1 :     int bHasAOI = FALSE;
    1046         [ +  - ]:          1 :     OGREnvelope sAOI;
    1047                 :          1 :     int nBands = 4;
    1048         [ +  - ]:          1 :     CPLString osProjection;
    1049                 :            : 
    1050 [ +  + ][ +  - ]:          1 :     if( (psXML != NULL && CPLGetXMLNode(psXML, "=GDAL_WMTS") != NULL ) ||
         [ +  + ][ +  + ]
           [ +  +  +  + ]
                 [ +  + ]
    1051                 :          1 :         STARTS_WITH_CI(poOpenInfo->pszFilename, "<GDAL_WMTS") ||
    1052                 :            :         (poOpenInfo->nHeaderBytes > 0 &&
    1053                 :          1 :          strstr((const char*)poOpenInfo->pabyHeader, "<GDAL_WMTS")) )
    1054                 :            :     {
    1055                 :            :         CPLXMLNode* psGDALWMTS;
    1056 [ +  + ][ +  - ]:          1 :         if( psXML != NULL && CPLGetXMLNode(psXML, "=GDAL_WMTS") != NULL )
         [ +  - ][ +  + ]
    1057         [ +  - ]:          1 :             psGDALWMTS = CPLCloneXMLTree(psXML);
    1058         [ +  + ]:          1 :         else if( STARTS_WITH_CI(poOpenInfo->pszFilename, "<GDAL_WMTS") )
    1059         [ +  - ]:          1 :             psGDALWMTS = CPLParseXMLString(poOpenInfo->pszFilename);
    1060                 :            :         else
    1061         [ +  - ]:          1 :             psGDALWMTS = CPLParseXMLFile(poOpenInfo->pszFilename);
    1062         [ +  + ]:          1 :         if( psGDALWMTS == NULL )
    1063                 :          1 :             return NULL;
    1064         [ +  - ]:          1 :         CPLXMLNode* psRoot = CPLGetXMLNode(psGDALWMTS, "=GDAL_WMTS");
    1065         [ +  + ]:          1 :         if( psRoot == NULL )
    1066                 :            :         {
    1067         [ +  - ]:          1 :             CPLError(CE_Failure, CPLE_AppDefined, "Cannot find root <GDAL_WMTS>");
    1068         [ +  - ]:          1 :             CPLDestroyXMLNode(psGDALWMTS);
    1069                 :          1 :             return NULL;
    1070                 :            :         }
    1071 [ +  - ][ +  - ]:          1 :         osGetCapabilitiesURL = CPLGetXMLValue(psRoot, "GetCapabilitiesUrl", "");
         [ +  - ][ +  - ]
    1072 [ +  - ][ +  + ]:          1 :         if( osGetCapabilitiesURL.empty() )
    1073                 :            :         {
    1074         [ +  - ]:          1 :             CPLError(CE_Failure, CPLE_AppDefined, "Missing <GetCapabilitiesUrl>");
    1075         [ +  - ]:          1 :             CPLDestroyXMLNode(psGDALWMTS);
    1076                 :          1 :             return NULL;
    1077                 :            :         }
    1078                 :            : 
    1079 [ +  - ][ +  - ]:          1 :         osLayer = CPLGetXMLValue(psRoot, "Layer", osLayer);
         [ +  - ][ +  - ]
                 [ +  - ]
    1080 [ +  - ][ +  - ]:          1 :         osTMS = CPLGetXMLValue(psRoot, "TileMatrixSet", osTMS);
         [ +  - ][ +  - ]
                 [ +  - ]
    1081                 :            :         osMaxTileMatrixIdentifier = CPLGetXMLValue(psRoot, "TileMatrix",
    1082 [ +  - ][ +  - ]:          1 :                                                    osMaxTileMatrixIdentifier);
         [ +  - ][ +  - ]
                 [ +  - ]
    1083                 :            :         nUserMaxZoomLevel = atoi(CPLGetXMLValue(psRoot, "ZoomLevel",
    1084 [ +  - ][ +  - ]:          1 :                                        CPLSPrintf("%d", nUserMaxZoomLevel)));
    1085 [ +  - ][ +  - ]:          1 :         osStyle = CPLGetXMLValue(psRoot, "Style", osStyle);
         [ +  - ][ +  - ]
                 [ +  - ]
    1086 [ +  - ][ +  - ]:          1 :         osTileFormat = CPLGetXMLValue(psRoot, "Format", osTileFormat);
         [ +  - ][ +  - ]
                 [ +  - ]
    1087 [ +  - ][ +  - ]:          1 :         osInfoFormat = CPLGetXMLValue(psRoot, "InfoFormat", osInfoFormat);
         [ +  - ][ +  - ]
                 [ +  - ]
    1088 [ +  - ][ +  - ]:          1 :         osProjection = CPLGetXMLValue(psRoot, "Projection", osProjection);
         [ +  - ][ +  - ]
                 [ +  - ]
    1089                 :            :         bExtendBeyondDateLine = CPLTestBool(CPLGetXMLValue(psRoot, "ExtendBeyondDateLine",
    1090 [ -  + ][ +  - ]:          1 :                                             (bExtendBeyondDateLine) ? "true": "false"));
                 [ +  - ]
    1091                 :            : 
    1092 [ +  - ][ +  - ]:          1 :         osOtherXML = "";
                 [ +  - ]
    1093         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "Cache", osOtherXML);
    1094         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "MaxConnections", osOtherXML);
    1095         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "Timeout", osOtherXML);
    1096         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "OfflineMode", osOtherXML);
    1097         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "MaxConnections", osOtherXML);
    1098         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "UserAgent", osOtherXML);
    1099         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "UserPwd", osOtherXML);
    1100         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "UnsafeSSL", osOtherXML);
    1101         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "Referer", osOtherXML);
    1102         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "ZeroBlockHttpCodes", osOtherXML);
    1103         [ +  - ]:          1 :         WMTSAddOtherXML(psRoot, "ZeroBlockOnServerException", osOtherXML);
    1104                 :            : 
    1105         [ +  - ]:          1 :         nBands = atoi(CPLGetXMLValue(psRoot, "BandsCount", "4"));
    1106                 :            : 
    1107         [ +  - ]:          1 :         const char* pszULX = CPLGetXMLValue(psRoot, "DataWindow.UpperLeftX", NULL);
    1108         [ +  - ]:          1 :         const char* pszULY = CPLGetXMLValue(psRoot, "DataWindow.UpperLeftY", NULL);
    1109         [ +  - ]:          1 :         const char* pszLRX = CPLGetXMLValue(psRoot, "DataWindow.LowerRightX", NULL);
    1110         [ +  - ]:          1 :         const char* pszLRY = CPLGetXMLValue(psRoot, "DataWindow.LowerRightY", NULL);
    1111 [ +  + ][ +  - ]:          1 :         if( pszULX && pszULY && pszLRX && pszLRY )
         [ +  - ][ +  - ]
    1112                 :            :         {
    1113         [ +  - ]:          1 :             sAOI.MinX = CPLAtof(pszULX);
    1114         [ +  - ]:          1 :             sAOI.MaxY = CPLAtof(pszULY);
    1115         [ +  - ]:          1 :             sAOI.MaxX = CPLAtof(pszLRX);
    1116         [ +  - ]:          1 :             sAOI.MinY = CPLAtof(pszLRY);
    1117                 :          1 :             bHasAOI = TRUE;
    1118                 :            :         }
    1119                 :            : 
    1120         [ +  - ]:          1 :         CPLDestroyXMLNode(psGDALWMTS);
    1121                 :            : 
    1122         [ +  - ]:          1 :         CPLDestroyXMLNode(psXML);
    1123 [ +  - ][ +  - ]:          1 :         char** papszHTTPOptions = BuildHTTPRequestOpts(osOtherXML);
                 [ +  - ]
    1124         [ +  - ]:          1 :         psXML = GetCapabilitiesResponse(osGetCapabilitiesURL, papszHTTPOptions);
    1125         [ +  - ]:          1 :         CSLDestroy(papszHTTPOptions);
    1126                 :            :     }
    1127         [ +  + ]:          1 :     else if( !STARTS_WITH_CI(poOpenInfo->pszFilename, "WMTS:") )
    1128                 :            :     {
    1129 [ +  - ][ +  - ]:          1 :         osGetCapabilitiesURL = poOpenInfo->pszFilename;
                 [ +  - ]
    1130         [ +  - ]:          1 :         psXML = CPLParseXMLFile(poOpenInfo->pszFilename);
    1131                 :            :     }
    1132         [ +  + ]:          1 :     if( psXML == NULL )
    1133                 :          1 :         return NULL;
    1134         [ +  - ]:          1 :     CPLStripXMLNamespace(psXML, NULL, TRUE);
    1135                 :            : 
    1136         [ +  - ]:          1 :     CPLXMLNode* psContents = CPLGetXMLNode(psXML, "=Capabilities.Contents");
    1137         [ +  + ]:          1 :     if( psContents == NULL )
    1138                 :            :     {
    1139         [ +  - ]:          1 :         CPLError(CE_Failure, CPLE_AppDefined, "Missing Capabilities.Contents element");
    1140         [ +  - ]:          1 :         CPLDestroyXMLNode(psXML);
    1141                 :          1 :         return NULL;
    1142                 :            :     }
    1143                 :            : 
    1144 [ +  - ][ +  - ]:          1 :     if( STARTS_WITH(osGetCapabilitiesURL, "/vsimem/") )
    1145                 :            :     {
    1146 [ +  - ][ +  + ]:          1 :         if( CPLGetXMLValue(psXML, "=Capabilities.ServiceMetadataURL.href", NULL) )
    1147 [ +  - ][ +  - ]:          1 :             osGetCapabilitiesURL = CPLGetXMLValue(psXML, "=Capabilities.ServiceMetadataURL.href", NULL);
         [ +  - ][ +  - ]
    1148                 :            :         else
    1149                 :            :         {
    1150 [ +  - ][ +  - ]:          1 :             osGetCapabilitiesURL = GetOperationKVPURL(psXML, "GetCapabilities");
                 [ +  - ]
    1151 [ +  - ][ +  + ]:          1 :             if( !osGetCapabilitiesURL.empty() )
    1152                 :            :             {
    1153 [ +  - ][ +  - ]:          1 :                 osGetCapabilitiesURL = CPLURLAddKVP(osGetCapabilitiesURL, "service", "WMTS");
         [ +  - ][ +  - ]
    1154 [ +  - ][ +  - ]:          1 :                 osGetCapabilitiesURL = CPLURLAddKVP(osGetCapabilitiesURL, "request", "GetCapabilities");
         [ +  - ][ +  - ]
    1155                 :            :             }
    1156                 :            :         }
    1157                 :            :     }
    1158         [ +  - ]:          1 :     CPLString osCapabilitiesFilename(osGetCapabilitiesURL);
    1159 [ +  - ][ +  - ]:          1 :     if( !STARTS_WITH_CI(osCapabilitiesFilename, "WMTS:") )
    1160 [ +  - ][ +  - ]:          1 :         osCapabilitiesFilename = "WMTS:" + osGetCapabilitiesURL;
         [ +  - ][ +  - ]
                 [ +  - ]
    1161                 :            : 
    1162                 :          1 :     int nLayerCount = 0;
    1163         [ +  - ]:          1 :     CPLStringList aosSubDatasets;
    1164 [ +  - ][ +  - ]:          1 :     CPLString osSelectLayer(osLayer), osSelectTMS(osTMS), osSelectStyle(osStyle);
                 [ +  - ]
    1165 [ +  - ][ +  - ]:          1 :     CPLString osSelectLayerTitle, osSelectLayerAbstract;
    1166 [ +  - ][ +  - ]:          1 :     CPLString osSelectTileFormat(osTileFormat), osSelectInfoFormat(osInfoFormat);
    1167                 :          1 :     int nCountTileFormat = 0;
    1168                 :          1 :     int nCountInfoFormat = 0;
    1169         [ +  - ]:          1 :     CPLString osURLTileTemplate;
    1170         [ +  - ]:          1 :     CPLString osURLFeatureInfoTemplate;
    1171         [ +  - ]:          1 :     std::set<CPLString> aoSetLayers;
    1172         [ +  - ]:          1 :     std::map<CPLString, OGREnvelope> aoMapBoundingBox;
    1173         [ +  - ]:          1 :     std::map<CPLString, WMTSTileMatrixLimits> aoMapTileMatrixLimits;
    1174         [ +  - ]:          1 :     std::map<CPLString, CPLString> aoMapDimensions;
    1175                 :            : 
    1176         [ +  + ]:          1 :     for(CPLXMLNode* psIter = psContents->psChild; psIter != NULL; psIter = psIter->psNext )
    1177                 :            :     {
    1178 [ +  - ][ +  + ]:          1 :         if( psIter->eType != CXT_Element || strcmp(psIter->pszValue, "Layer") != 0 )
    1179                 :          1 :             continue;
    1180         [ +  - ]:          1 :         const char* pszIdentifier = CPLGetXMLValue(psIter, "Identifier", "");
    1181 [ +  - ][ +  - ]:          1 :         if( aoSetLayers.find(pszIdentifier) != aoSetLayers.end() )
         [ +  - ][ +  - ]
         [ +  - ][ -  + ]
    1182                 :            :         {
    1183                 :            :             CPLError(CE_Warning, CPLE_AppDefined,
    1184                 :            :                      "Several layers with identifier '%s'. Only first one kept",
    1185         [ #  # ]:          0 :                      pszIdentifier);
    1186                 :            :         }
    1187 [ +  - ][ +  - ]:          1 :         aoSetLayers.insert(pszIdentifier);
                 [ +  - ]
    1188 [ +  - ][ +  + ]:          1 :         if( !osLayer.empty() && strcmp(osLayer, pszIdentifier) != 0 )
         [ +  - ][ +  + ]
                 [ +  + ]
    1189                 :          1 :             continue;
    1190         [ +  - ]:          1 :         const char* pszTitle = CPLGetXMLValue(psIter, "Title", NULL);
    1191 [ +  - ][ +  + ]:          1 :         if( osSelectLayer.empty() )
    1192                 :            :         {
    1193 [ +  - ][ +  - ]:          1 :             osSelectLayer = pszIdentifier;
                 [ +  - ]
    1194                 :            :         }
    1195 [ +  - ][ +  - ]:          1 :         if( strcmp(osSelectLayer, pszIdentifier) == 0 )
    1196                 :            :         {
    1197         [ +  + ]:          1 :             if( pszTitle != NULL )
    1198 [ +  - ][ +  - ]:          1 :                 osSelectLayerTitle = pszTitle;
                 [ +  - ]
    1199         [ +  - ]:          1 :             const char* pszAbstract = CPLGetXMLValue(psIter, "Abstract", NULL);
    1200         [ +  + ]:          1 :             if( pszAbstract != NULL )
    1201 [ +  - ][ +  - ]:          1 :                 osSelectLayerAbstract = pszAbstract;
                 [ +  - ]
    1202                 :            :         }
    1203                 :            : 
    1204         [ +  - ]:          1 :         std::vector<CPLString> aosTMS;
    1205         [ +  - ]:          1 :         std::vector<CPLString> aosStylesIdentifier;
    1206         [ +  - ]:          1 :         std::vector<CPLString> aosStylesTitle;
    1207                 :            : 
    1208                 :          1 :         CPLXMLNode* psSubIter = psIter->psChild;
    1209         [ +  + ]:          1 :         for(; psSubIter != NULL; psSubIter = psSubIter->psNext )
    1210                 :            :         {
    1211         [ +  + ]:          1 :             if( psSubIter->eType != CXT_Element )
    1212                 :          1 :                 continue;
    1213 [ +  - ][ +  - ]:          1 :             if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1214                 :          1 :                 strcmp(psSubIter->pszValue, "Format") == 0 )
    1215                 :            :             {
    1216         [ +  - ]:          1 :                 const char* pszValue = CPLGetXMLValue(psSubIter, "", "");
    1217 [ +  - ][ +  + ]:          1 :                 if( !osTileFormat.empty() && strcmp(osTileFormat, pszValue) != 0 )
         [ +  - ][ +  + ]
                 [ +  + ]
    1218                 :          1 :                     continue;
    1219                 :          1 :                 nCountTileFormat ++;
    1220 [ +  - ][ +  + ]:          1 :                 if( osSelectTileFormat.empty() ||
         [ +  - ][ +  - ]
    1221                 :          1 :                     EQUAL(pszValue, "image/png") )
    1222                 :            :                 {
    1223 [ +  - ][ +  - ]:          1 :                     osSelectTileFormat = pszValue;
                 [ +  - ]
    1224                 :            :                 }
    1225                 :            :             }
    1226 [ +  - ][ +  - ]:          1 :             else if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1227                 :          1 :                      strcmp(psSubIter->pszValue, "InfoFormat") == 0 )
    1228                 :            :             {
    1229         [ +  - ]:          1 :                 const char* pszValue = CPLGetXMLValue(psSubIter, "", "");
    1230 [ +  - ][ -  + ]:          1 :                 if( !osInfoFormat.empty() && strcmp(osInfoFormat, pszValue) != 0 )
         [ #  # ][ #  # ]
                 [ -  + ]
    1231                 :          0 :                     continue;
    1232                 :          1 :                 nCountInfoFormat ++;
    1233 [ +  - ][ -  + ]:          1 :                 if( osSelectInfoFormat.empty() ||
         [ #  # ][ #  # ]
         [ #  # ][ +  - ]
    1234                 :          0 :                     (EQUAL(pszValue, "application/vnd.ogc.gml") &&
    1235         [ #  # ]:          0 :                      !EQUAL(osSelectInfoFormat, "application/vnd.ogc.gml/3.1.1")) ||
    1236                 :          0 :                     EQUAL(pszValue, "application/vnd.ogc.gml/3.1.1") )
    1237                 :            :                 {
    1238 [ +  - ][ +  - ]:          1 :                     osSelectInfoFormat = pszValue;
                 [ +  - ]
    1239                 :            :                 }
    1240                 :            :             }
    1241 [ +  - ][ +  - ]:          1 :             else if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1242                 :          1 :                      strcmp(psSubIter->pszValue, "Dimension") == 0 )
    1243                 :            :             {
    1244                 :            :                 /* Cf http://wmts.geo.admin.ch/1.0.0/WMTSCapabilities.xml */
    1245         [ +  - ]:          1 :                 const char* pszDimensionIdentifier = CPLGetXMLValue(psSubIter, "Identifier", NULL);
    1246         [ +  - ]:          1 :                 const char* pszDefault = CPLGetXMLValue(psSubIter, "Default", "");
    1247         [ +  - ]:          1 :                 if( pszDimensionIdentifier != NULL )
    1248 [ +  - ][ +  - ]:          1 :                     aoMapDimensions[pszDimensionIdentifier] = pszDefault;
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    1249                 :            :             }
    1250         [ +  + ]:          1 :             else if( strcmp(psSubIter->pszValue, "TileMatrixSetLink") == 0 )
    1251                 :            :             {
    1252                 :            :                 const char* pszTMS = CPLGetXMLValue(
    1253         [ +  - ]:          1 :                                             psSubIter, "TileMatrixSet", "");
    1254 [ +  - ][ +  + ]:          1 :                 if( !osTMS.empty() && strcmp(osTMS, pszTMS) != 0 )
         [ +  - ][ +  + ]
                 [ +  + ]
    1255                 :          1 :                     continue;
    1256 [ +  - ][ +  - ]:          1 :                 if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1257         [ +  - ]:          1 :                     osSelectTMS.empty() )
    1258                 :            :                 {
    1259 [ +  - ][ +  - ]:          1 :                     osSelectTMS = pszTMS;
                 [ +  - ]
    1260                 :            :                 }
    1261 [ +  - ][ +  - ]:          1 :                 if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1262         [ +  - ]:          1 :                     strcmp(osSelectTMS, pszTMS) == 0 )
    1263                 :            :                 {
    1264                 :            :                     CPLXMLNode* psTMSLimits = CPLGetXMLNode(
    1265         [ +  - ]:          1 :                                         psSubIter, "TileMatrixSetLimits");
    1266         [ +  + ]:          1 :                     if( psTMSLimits )
    1267         [ +  - ]:          1 :                         ReadTMLimits(psTMSLimits, aoMapTileMatrixLimits);
    1268                 :            :                 }
    1269 [ +  - ][ +  - ]:          1 :                 aosTMS.push_back(pszTMS);
                 [ +  - ]
    1270                 :            :             }
    1271         [ +  + ]:          1 :             else if( strcmp(psSubIter->pszValue, "Style") == 0 )
    1272                 :            :             {
    1273                 :            :                 int bIsDefault = CPLTestBool(CPLGetXMLValue(
    1274 [ +  - ][ +  - ]:          1 :                                             psSubIter, "isDefault", "false"));
    1275                 :            :                 const char* l_pszIdentifier = CPLGetXMLValue(
    1276         [ +  - ]:          1 :                                             psSubIter, "Identifier", "");
    1277 [ +  - ][ +  + ]:          1 :                 if( !osStyle.empty() && strcmp(osStyle, l_pszIdentifier) != 0 )
         [ +  - ][ +  + ]
                 [ +  + ]
    1278                 :          1 :                     continue;
    1279                 :            :                 const char* pszStyleTitle = CPLGetXMLValue(
    1280         [ +  - ]:          1 :                                         psSubIter, "Title", l_pszIdentifier);
    1281         [ +  + ]:          1 :                 if( bIsDefault )
    1282                 :            :                 {
    1283                 :            :                     aosStylesIdentifier.insert(aosStylesIdentifier.begin(),
    1284 [ +  - ][ +  - ]:          1 :                                                 CPLString(l_pszIdentifier));
         [ +  - ][ +  - ]
    1285                 :            :                     aosStylesTitle.insert(aosStylesTitle.begin(),
    1286 [ +  - ][ +  - ]:          1 :                                             CPLString(pszStyleTitle));
         [ +  - ][ +  - ]
    1287 [ +  - ][ -  + ]:          1 :                     if( strcmp(osSelectLayer, l_pszIdentifier) == 0 &&
         [ #  # ][ -  + ]
    1288         [ #  # ]:          0 :                         osSelectStyle.empty() )
    1289                 :            :                     {
    1290 [ #  # ][ #  # ]:          0 :                         osSelectStyle = l_pszIdentifier;
                 [ #  # ]
    1291                 :            :                     }
    1292                 :            :                 }
    1293                 :            :                 else
    1294                 :            :                 {
    1295 [ +  - ][ +  - ]:          1 :                     aosStylesIdentifier.push_back(l_pszIdentifier);
                 [ +  - ]
    1296 [ +  - ][ +  - ]:          1 :                     aosStylesTitle.push_back(pszStyleTitle);
                 [ +  - ]
    1297                 :            :                 }
    1298                 :            :             }
    1299 [ +  - ][ +  - ]:          1 :             else if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
                 [ +  + ]
    1300                 :          1 :                      (strcmp(psSubIter->pszValue, "BoundingBox") == 0 ||
    1301                 :          1 :                       strcmp(psSubIter->pszValue, "WGS84BoundingBox") == 0) )
    1302                 :            :             {
    1303 [ +  - ][ +  - ]:          1 :                 CPLString osCRS = CPLGetXMLValue(psSubIter, "crs", "");
    1304 [ +  - ][ +  + ]:          1 :                 if( osCRS.empty() )
    1305                 :            :                 {
    1306         [ +  - ]:          1 :                     if( strcmp(psSubIter->pszValue, "WGS84BoundingBox") == 0 )
    1307                 :            :                     {
    1308 [ +  - ][ +  - ]:          1 :                         osCRS = "EPSG:4326";
                 [ +  - ]
    1309                 :            :                     }
    1310                 :            :                     else
    1311                 :            :                     {
    1312                 :          0 :                         int nCountTileMatrixSet = 0;
    1313         [ #  # ]:          0 :                         CPLString osSingleTileMatrixSet;
    1314         [ #  # ]:          0 :                         for(CPLXMLNode* psIter3 = psContents->psChild; psIter3 != NULL; psIter3 = psIter3->psNext )
    1315                 :            :                         {
    1316 [ #  # ][ #  # ]:          0 :                             if( psIter3->eType != CXT_Element || strcmp(psIter3->pszValue, "TileMatrixSet") != 0 )
    1317                 :          0 :                                 continue;
    1318                 :          0 :                             nCountTileMatrixSet ++;
    1319         [ #  # ]:          0 :                             if( nCountTileMatrixSet == 1 )
    1320 [ #  # ][ #  # ]:          0 :                                 osSingleTileMatrixSet = CPLGetXMLValue(psIter3, "Identifier", "");
         [ #  # ][ #  # ]
    1321                 :            :                         }
    1322         [ #  # ]:          0 :                         if( nCountTileMatrixSet == 1 )
    1323                 :            :                         {
    1324                 :            :                             // For 13-082_WMTS_Simple_Profile/schemas/wmts/1.0/profiles/WMTSSimple/examples/wmtsGetCapabilities_response_OSM.xml
    1325         [ #  # ]:          0 :                             WMTSTileMatrixSet oTMS;
    1326         [ #  # ]:          0 :                             if( ReadTMS(psContents, osSingleTileMatrixSet,
    1327 [ #  # ][ #  # ]:          0 :                                         CPLString(), -1, oTMS) )
                 [ #  # ]
    1328                 :            :                             {
    1329         [ #  # ]:          0 :                                 osCRS = oTMS.osSRS;
    1330         [ #  # ]:          0 :                             }
    1331         [ #  # ]:          0 :                         }
    1332                 :            :                     }
    1333                 :            :                 }
    1334 [ +  - ][ +  - ]:          1 :                 CPLString osLowerCorner = CPLGetXMLValue(psSubIter, "LowerCorner", "");
    1335 [ +  - ][ +  - ]:          1 :                 CPLString osUpperCorner = CPLGetXMLValue(psSubIter, "UpperCorner", "");
    1336         [ +  - ]:          1 :                 OGRSpatialReference oSRS;
    1337 [ +  - ][ +  - ]:          1 :                 if( !osCRS.empty() && !osLowerCorner.empty() && !osUpperCorner.empty() &&
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    1338 [ +  - ][ +  - ]:          1 :                     oSRS.SetFromUserInput(FixCRSName(osCRS)) == OGRERR_NONE )
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ #  # ]
    1339                 :            :                 {
    1340         [ +  - ]:          1 :                     int bSwap = oSRS.EPSGTreatsAsLatLong() ||
    1341 [ +  + ][ +  - ]:          1 :                                 oSRS.EPSGTreatsAsNorthingEasting();
                 [ -  + ]
    1342 [ +  - ][ +  - ]:          1 :                     char** papszLC = CSLTokenizeString(osLowerCorner);
    1343 [ +  - ][ +  - ]:          1 :                     char** papszUC = CSLTokenizeString(osUpperCorner);
    1344 [ +  - ][ +  - ]:          1 :                     if( CSLCount(papszLC) == 2 && CSLCount(papszUC) == 2 )
         [ +  - ][ +  - ]
                 [ +  - ]
    1345                 :            :                     {
    1346         [ +  - ]:          1 :                         OGREnvelope sEnvelope;
    1347 [ +  + ][ +  - ]:          1 :                         sEnvelope.MinX = CPLAtof(papszLC[(bSwap)? 1 : 0]);
    1348 [ +  + ][ +  - ]:          1 :                         sEnvelope.MinY = CPLAtof(papszLC[(bSwap)? 0 : 1]);
    1349 [ +  + ][ +  - ]:          1 :                         sEnvelope.MaxX = CPLAtof(papszUC[(bSwap)? 1 : 0]);
    1350 [ +  + ][ +  - ]:          1 :                         sEnvelope.MaxY = CPLAtof(papszUC[(bSwap)? 0 : 1]);
    1351         [ +  - ]:          1 :                         aoMapBoundingBox[osCRS] = sEnvelope;
    1352                 :            :                     }
    1353         [ +  - ]:          1 :                     CSLDestroy(papszLC);
    1354         [ +  - ]:          1 :                     CSLDestroy(papszUC);
    1355 [ +  - ][ +  - ]:          1 :                 }
         [ +  - ][ +  - ]
    1356                 :            :             }
    1357 [ +  - ][ +  - ]:          1 :             else if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
    1358                 :          1 :                      strcmp(psSubIter->pszValue, "ResourceURL") == 0 )
    1359                 :            :             {
    1360 [ +  - ][ +  + ]:          1 :                 if( EQUAL(CPLGetXMLValue(psSubIter, "resourceType", ""), "tile") )
    1361                 :            :                 {
    1362         [ +  - ]:          1 :                     const char* pszFormat = CPLGetXMLValue(psSubIter, "format", "");
    1363 [ +  - ][ -  + ]:          1 :                     if( !osTileFormat.empty() && strcmp(osTileFormat, pszFormat) != 0 )
         [ #  # ][ #  # ]
                 [ -  + ]
    1364                 :          0 :                         continue;
    1365 [ +  - ][ +  - ]:          1 :                     if( osURLTileTemplate.empty() )
    1366 [ +  - ][ +  - ]:          1 :                         osURLTileTemplate = CPLGetXMLValue(psSubIter, "template", "");
         [ +  - ][ +  - ]
    1367                 :            :                 }
    1368 [ +  - ][ +  - ]:          1 :                 else if( EQUAL(CPLGetXMLValue(psSubIter, "resourceType", ""), "FeatureInfo") )
    1369                 :            :                 {
    1370         [ +  - ]:          1 :                     const char* pszFormat = CPLGetXMLValue(psSubIter, "format", "");
    1371 [ +  - ][ -  + ]:          1 :                     if( !osInfoFormat.empty() && strcmp(osInfoFormat, pszFormat) != 0 )
         [ #  # ][ #  # ]
                 [ -  + ]
    1372                 :          0 :                         continue;
    1373 [ +  - ][ +  - ]:          1 :                     if( osURLFeatureInfoTemplate.empty() )
    1374 [ +  - ][ +  - ]:          1 :                         osURLFeatureInfoTemplate = CPLGetXMLValue(psSubIter, "template", "");
         [ +  - ][ +  - ]
    1375                 :            :                 }
    1376                 :            :             }
    1377                 :            :         }
    1378 [ +  - ][ +  - ]:          1 :         if( strcmp(osSelectLayer, pszIdentifier) == 0 &&
         [ +  + ][ +  + ]
                 [ +  + ]
    1379 [ +  - ][ +  - ]:          1 :             osSelectStyle.empty() && !aosStylesIdentifier.empty() )
    1380                 :            :         {
    1381 [ +  - ][ +  - ]:          1 :             osSelectStyle = aosStylesIdentifier[0];
    1382                 :            :         }
    1383 [ +  - ][ +  + ]:          1 :         for(size_t i=0;i<aosTMS.size();i++)
    1384                 :            :         {
    1385 [ +  - ][ +  + ]:          1 :             for(size_t j=0;j<aosStylesIdentifier.size();j++)
    1386                 :            :             {
    1387         [ +  - ]:          1 :                 int nIdx = 1 + aosSubDatasets.size() / 2;
    1388         [ +  - ]:          1 :                 CPLString osName(osCapabilitiesFilename);
    1389         [ +  - ]:          1 :                 osName += ",layer=";
    1390 [ +  - ][ +  - ]:          1 :                 osName += QuoteIfNecessary(pszIdentifier);
                 [ +  - ]
    1391 [ +  - ][ +  + ]:          1 :                 if( aosTMS.size() > 1 )
    1392                 :            :                 {
    1393         [ +  - ]:          1 :                     osName += ",tilematrixset=";
    1394 [ +  - ][ +  - ]:          1 :                     osName += QuoteIfNecessary(aosTMS[i]);
         [ +  - ][ +  - ]
                 [ +  - ]
    1395                 :            :                 }
    1396 [ +  - ][ +  + ]:          1 :                 if( aosStylesIdentifier.size() > 1 )
    1397                 :            :                 {
    1398         [ +  - ]:          1 :                     osName += ",style=";
    1399 [ +  - ][ +  - ]:          1 :                     osName += QuoteIfNecessary(aosStylesIdentifier[j]);
         [ +  - ][ +  - ]
                 [ +  - ]
    1400                 :            :                 }
    1401                 :            :                 aosSubDatasets.AddNameValue(
    1402 [ +  - ][ +  - ]:          1 :                     CPLSPrintf("SUBDATASET_%d_NAME", nIdx), osName);
                 [ +  - ]
    1403                 :            : 
    1404         [ +  - ]:          1 :                 CPLString osDesc("Layer ");
    1405 [ +  + ][ +  - ]:          1 :                 osDesc += pszTitle ? pszTitle : pszIdentifier;
    1406 [ +  - ][ +  + ]:          1 :                 if( aosTMS.size() > 1 )
    1407                 :            :                 {
    1408         [ +  - ]:          1 :                     osDesc += ", tile matrix set ";
    1409 [ +  - ][ +  - ]:          1 :                     osDesc += aosTMS[i];
    1410                 :            :                 }
    1411 [ +  - ][ +  + ]:          1 :                 if( aosStylesIdentifier.size() > 1 )
    1412                 :            :                 {
    1413         [ +  - ]:          1 :                     osDesc += ", style ";
    1414 [ +  - ][ +  - ]:          1 :                     osDesc += QuoteIfNecessary(aosStylesTitle[j]);
         [ +  - ][ +  - ]
                 [ +  - ]
    1415                 :            :                 }
    1416                 :            :                 aosSubDatasets.AddNameValue(
    1417 [ +  - ][ +  - ]:          1 :                     CPLSPrintf("SUBDATASET_%d_DESC", nIdx), osDesc);
                 [ +  - ]
    1418 [ +  - ][ +  - ]:          1 :             }
    1419                 :            :         }
    1420 [ +  - ][ +  + ]:          1 :         if( !aosTMS.empty() && !aosStylesIdentifier.empty() )
         [ +  - ][ +  + ]
                 [ +  + ]
    1421                 :          1 :             nLayerCount ++;
    1422                 :            :         else
    1423         [ +  - ]:          1 :             CPLError(CE_Failure, CPLE_AppDefined, "Missing TileMatrixSetLink and/or Style");
    1424 [ +  - ][ +  - ]:          1 :     }
                 [ +  - ]
    1425                 :            : 
    1426         [ +  + ]:          1 :     if( nLayerCount == 0 )
    1427                 :            :     {
    1428         [ +  - ]:          1 :         CPLDestroyXMLNode(psXML);
    1429                 :          1 :         return NULL;
    1430                 :            :     }
    1431                 :            : 
    1432 [ +  - ][ +  - ]:          1 :     WMTSDataset* poDS = new WMTSDataset();
    1433                 :            : 
    1434 [ +  - ][ +  + ]:          1 :     if( aosSubDatasets.size() > 2 )
    1435 [ +  - ][ +  - ]:          1 :         poDS->SetMetadata(aosSubDatasets.List(), "SUBDATASETS");
    1436                 :            : 
    1437         [ +  - ]:          1 :     if( nLayerCount == 1 )
    1438                 :            :     {
    1439 [ +  - ][ +  + ]:          1 :         if( !osSelectLayerTitle.empty() )
    1440 [ +  - ][ +  - ]:          1 :             poDS->SetMetadataItem("TITLE", osSelectLayerTitle);
    1441 [ +  - ][ +  + ]:          1 :         if( !osSelectLayerAbstract.empty() )
    1442 [ +  - ][ +  - ]:          1 :             poDS->SetMetadataItem("ABSTRACT", osSelectLayerAbstract);
    1443                 :            : 
    1444 [ +  - ][ +  - ]:          1 :         poDS->papszHTTPOptions = BuildHTTPRequestOpts(osOtherXML);
                 [ +  - ]
    1445         [ +  - ]:          1 :         poDS->osLayer = osSelectLayer;
    1446         [ +  - ]:          1 :         poDS->osTMS = osSelectTMS;
    1447                 :            : 
    1448         [ +  - ]:          1 :         WMTSTileMatrixSet oTMS;
    1449         [ +  + ]:          1 :         if( !ReadTMS(psContents, osSelectTMS, osMaxTileMatrixIdentifier,
    1450         [ +  - ]:          1 :                      nUserMaxZoomLevel, oTMS) )
    1451                 :            :         {
    1452         [ +  - ]:          1 :             CPLDestroyXMLNode(psXML);
    1453 [ +  - ][ +  - ]:          1 :             delete poDS;
    1454                 :          1 :             return NULL;
    1455                 :            :         }
    1456                 :            : 
    1457                 :            :         const char* pszExtentMethod = CSLFetchNameValueDef(
    1458         [ +  - ]:          1 :             poOpenInfo->papszOpenOptions, "EXTENT_METHOD", "AUTO");
    1459                 :          1 :         ExtentMethod eExtentMethod = AUTO;
    1460         [ -  + ]:          1 :         if( EQUAL(pszExtentMethod, "LAYER_BBOX") )
    1461                 :          0 :             eExtentMethod = LAYER_BBOX;
    1462         [ -  + ]:          1 :         else if( EQUAL(pszExtentMethod, "TILE_MATRIX_SET") )
    1463                 :          0 :             eExtentMethod = TILE_MATRIX_SET;
    1464         [ -  + ]:          1 :         else if( EQUAL(pszExtentMethod, "MOST_PRECISE_TILE_MATRIX") )
    1465                 :          0 :             eExtentMethod = MOST_PRECISE_TILE_MATRIX;
    1466                 :            : 
    1467                 :            :         // Use in priority layer bounding box expressed in the SRS of the TMS
    1468 [ +  + ][ -  + ]:          1 :         if( (!bHasAOI || bExtendBeyondDateLine) &&
         [ -  + ][ #  # ]
         [ +  + ][ +  + ]
    1469                 :            :             (eExtentMethod == AUTO || eExtentMethod == LAYER_BBOX) &&
    1470 [ +  - ][ +  - ]:          1 :             aoMapBoundingBox.find(oTMS.osSRS) != aoMapBoundingBox.end() )
                 [ +  - ]
    1471                 :            :         {
    1472         [ +  - ]:          1 :             if( !bHasAOI )
    1473                 :            :             {
    1474         [ +  - ]:          1 :                 sAOI = aoMapBoundingBox[oTMS.osSRS];
    1475                 :          1 :                 bHasAOI = TRUE;
    1476                 :            :             }
    1477                 :            : 
    1478                 :          1 :             int bRecomputeAOI = FALSE;
    1479         [ +  + ]:          1 :             if( bExtendBeyondDateLine )
    1480                 :            :             {
    1481                 :          1 :                 bExtendBeyondDateLine = FALSE;
    1482                 :            : 
    1483         [ +  - ]:          1 :                 OGRSpatialReference oWGS84;
    1484         [ +  - ]:          1 :                     oWGS84.SetFromUserInput(SRS_WKT_WGS84);
    1485                 :            :                 OGRCoordinateTransformation* poCT =
    1486         [ +  - ]:          1 :                     OGRCreateCoordinateTransformation(&oTMS.oSRS, &oWGS84);
    1487         [ +  - ]:          1 :                 if( poCT != NULL )
    1488                 :            :                 {
    1489                 :          1 :                     double dfX1 = sAOI.MinX;
    1490                 :          1 :                     double dfY1 = sAOI.MinY;
    1491                 :          1 :                     double dfX2 = sAOI.MaxX;
    1492                 :          1 :                     double dfY2 = sAOI.MaxY;
    1493 [ +  - ][ +  - ]:          1 :                     if( poCT->Transform(1, &dfX1, &dfY1) &&
         [ +  - ][ +  - ]
    1494         [ +  - ]:          1 :                         poCT->Transform(1, &dfX2, &dfY2) )
    1495                 :            :                     {
    1496 [ +  - ][ +  - ]:          1 :                         if( fabs(dfX1 + 180) < 1e-8 &&
    1497                 :          1 :                             fabs(dfX2 - 180) < 1e-8 )
    1498                 :            :                         {
    1499                 :          1 :                             bExtendBeyondDateLine = TRUE;
    1500                 :          1 :                             bRecomputeAOI = TRUE;
    1501                 :            :                         }
    1502         [ #  # ]:          0 :                         else if( dfX2 < dfX1 )
    1503                 :            :                         {
    1504                 :          0 :                             bExtendBeyondDateLine = TRUE;
    1505                 :            :                         }
    1506                 :            :                         else
    1507                 :            :                         {
    1508                 :            :                             CPLError(CE_Warning, CPLE_AppDefined,
    1509                 :            :                                     "ExtendBeyondDateLine disabled, since longitudes of %s "
    1510                 :            :                                     "BoundingBox do not span from -180 to 180 but from %.16g to %.16g, "
    1511                 :            :                                     "or longitude of upper right corner is not lesser than the one of lower left corner",
    1512 [ #  # ][ #  # ]:          1 :                                     oTMS.osSRS.c_str(), dfX1, dfX2);
    1513                 :            :                         }
    1514                 :            :                     }
    1515 [ +  - ][ +  - ]:          1 :                     delete poCT;
    1516         [ +  - ]:          1 :                 }
    1517                 :            :             }
    1518 [ +  + ][ +  - ]:          1 :             if( bExtendBeyondDateLine && bRecomputeAOI )
    1519                 :            :             {
    1520                 :          1 :                 bExtendBeyondDateLine = FALSE;
    1521                 :            : 
    1522         [ +  - ]:          1 :                 std::map<CPLString, OGREnvelope>::iterator oIter = aoMapBoundingBox.begin();
    1523 [ +  - ][ +  - ]:          1 :                 for(; oIter != aoMapBoundingBox.end(); ++oIter )
         [ +  - ][ +  - ]
    1524                 :            :                 {
    1525         [ +  - ]:          1 :                     OGRSpatialReference oSRS;
    1526 [ +  - ][ +  - ]:          1 :                     if( oSRS.SetFromUserInput(FixCRSName(oIter->first)) == OGRERR_NONE )
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1527                 :            :                     {
    1528         [ +  - ]:          1 :                         OGRSpatialReference oWGS84;
    1529         [ +  - ]:          1 :                         oWGS84.SetFromUserInput(SRS_WKT_WGS84);
    1530                 :            :                         OGRCoordinateTransformation* poCT =
    1531         [ +  - ]:          1 :                             OGRCreateCoordinateTransformation(&oSRS, &oWGS84);
    1532         [ +  - ]:          1 :                         double dfX1 = oIter->second.MinX;
    1533         [ +  - ]:          1 :                         double dfY1 = oIter->second.MinY;
    1534         [ +  - ]:          1 :                         double dfX2 = oIter->second.MaxX;
    1535         [ +  - ]:          1 :                         double dfY2 = oIter->second.MaxY;
    1536 [ +  - ][ +  - ]:          1 :                         if( poCT != NULL &&
         [ +  - ][ +  + ]
                 [ +  + ]
    1537         [ +  - ]:          1 :                             poCT->Transform(1, &dfX1, &dfY1) &&
    1538         [ +  - ]:          1 :                             poCT->Transform(1, &dfX2, &dfY2) &&
    1539                 :            :                             dfX2 < dfX1 )
    1540                 :            :                         {
    1541 [ +  - ][ +  - ]:          1 :                             delete poCT;
    1542                 :          1 :                             dfX2 += 360;
    1543                 :          1 :                             char* pszProj4 = NULL;
    1544         [ +  - ]:          1 :                             oTMS.oSRS.exportToProj4(&pszProj4);
    1545 [ +  - ][ +  - ]:          1 :                             oSRS.SetFromUserInput(CPLSPrintf("%s +over +wktext", pszProj4));
    1546         [ +  - ]:          1 :                             CPLFree(pszProj4);
    1547         [ +  - ]:          1 :                             poCT = OGRCreateCoordinateTransformation(&oWGS84, &oSRS);
    1548 [ +  - ][ +  - ]:          1 :                             if( poCT &&
         [ +  - ][ +  - ]
    1549         [ +  - ]:          1 :                                 poCT->Transform(1, &dfX1, &dfY1) &&
    1550         [ +  - ]:          1 :                                 poCT->Transform(1, &dfX2, &dfY2) )
    1551                 :            :                             {
    1552                 :          1 :                                 bExtendBeyondDateLine = TRUE;
    1553         [ +  - ]:          1 :                                 sAOI.MinX = std::min(dfX1, dfX2);
    1554         [ +  - ]:          1 :                                 sAOI.MinY = std::min(dfY1, dfY2);
    1555         [ +  - ]:          1 :                                 sAOI.MaxX = std::max(dfX1, dfX2);
    1556         [ +  - ]:          1 :                                 sAOI.MaxY = std::max(dfY1, dfY2);
    1557                 :            :                                 CPLDebug("WMTS",
    1558                 :            :                                          "ExtendBeyondDateLine using %s bounding box",
    1559 [ +  - ][ +  - ]:          1 :                                          oIter->first.c_str());
                 [ +  - ]
    1560                 :            :                             }
    1561 [ +  - ][ +  - ]:          1 :                             delete poCT;
    1562                 :            :                             break;
    1563                 :            :                         }
    1564 [ +  - ][ +  - ]:          1 :                         delete poCT;
         [ +  - ][ +  + ]
    1565                 :            :                     }
    1566 [ +  - ][ +  + ]:          1 :                 }
    1567                 :            :             }
    1568                 :            :         }
    1569                 :            :         else
    1570                 :            :         {
    1571         [ -  + ]:          1 :             if( bExtendBeyondDateLine )
    1572                 :            :             {
    1573                 :            :                 CPLError(CE_Warning, CPLE_AppDefined,
    1574                 :            :                          "ExtendBeyondDateLine disabled, since BoundingBox of %s is missing",
    1575 [ #  # ][ #  # ]:          0 :                          oTMS.osSRS.c_str());
    1576                 :          0 :                 bExtendBeyondDateLine = FALSE;
    1577                 :            :             }
    1578                 :            :         }
    1579                 :            : 
    1580                 :            :         // Otherwise default to reproject a layer bounding box expressed in
    1581                 :            :         // another SRS
    1582 [ +  + ][ +  - ]:          1 :         if( !bHasAOI && !aoMapBoundingBox.empty() &&
         [ +  + ][ -  + ]
         [ #  # ][ +  + ]
    1583                 :            :             (eExtentMethod == AUTO || eExtentMethod == LAYER_BBOX) )
    1584                 :            :         {
    1585         [ +  - ]:          1 :             std::map<CPLString, OGREnvelope>::iterator oIter = aoMapBoundingBox.begin();
    1586 [ #  # ][ +  - ]:          1 :             for(; oIter != aoMapBoundingBox.end(); ++oIter )
         [ +  - ][ +  - ]
    1587                 :            :             {
    1588         [ +  - ]:          1 :                 OGRSpatialReference oSRS;
    1589 [ +  - ][ +  - ]:          1 :                 if( oSRS.SetFromUserInput(FixCRSName(oIter->first)) == OGRERR_NONE )
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1590                 :            :                 {
    1591                 :            :                     // Check if this doesn't match the most precise tile matrix
    1592                 :            :                     // by densifying its contour
    1593         [ +  - ]:          1 :                     const WMTSTileMatrix& oTM = oTMS.aoTM.back();
    1594                 :            : 
    1595                 :          1 :                     bool bMatchFound = false;
    1596         [ +  - ]:          1 :                     const char *pszProjectionTMS = oTMS.oSRS.GetAttrValue("PROJECTION");
    1597         [ +  - ]:          1 :                     const char *pszProjectionBBOX = oSRS.GetAttrValue("PROJECTION");
    1598                 :            :                     const bool bIsTMerc = (pszProjectionTMS != NULL &&
    1599                 :          1 :                         EQUAL(pszProjectionTMS, SRS_PT_TRANSVERSE_MERCATOR)) ||
    1600                 :            :                         (pszProjectionBBOX != NULL &&
    1601 [ +  + ][ +  + ]:          1 :                         EQUAL(pszProjectionBBOX, SRS_PT_TRANSVERSE_MERCATOR));
         [ -  + ][ #  # ]
    1602                 :            :                     // If one of the 2 SRS is a TMerc, try with classical tmerc
    1603                 :            :                     // or etmerc.
    1604 [ +  + ][ +  + ]:          1 :                     for( int j = 0; j < (bIsTMerc ? 2 : 1); j++ )
    1605                 :            :                     {
    1606                 :            :                         CPLString osOldVal =
    1607 [ +  - ][ +  - ]:          1 :                             CPLGetThreadLocalConfigOption("OSR_USE_ETMERC", "");
    1608         [ +  + ]:          1 :                         if( bIsTMerc )
    1609                 :            :                         {
    1610                 :            :                             CPLSetThreadLocalConfigOption("OSR_USE_ETMERC",
    1611 [ +  - ][ +  - ]:          1 :                                                       (j==0) ? "NO" : "YES");
    1612                 :            :                         }
    1613                 :            :                         OGRCoordinateTransformation* poRevCT =
    1614         [ +  - ]:          1 :                             OGRCreateCoordinateTransformation(&oTMS.oSRS, &oSRS);
    1615         [ +  + ]:          1 :                         if( bIsTMerc )
    1616                 :            :                         {
    1617                 :            :                             CPLSetThreadLocalConfigOption("OSR_USE_ETMERC",
    1618 [ +  - ][ +  - ]:          1 :                                 osOldVal.empty() ? NULL : osOldVal.c_str());
         [ #  # ][ +  - ]
    1619                 :            :                         }
    1620         [ +  - ]:          1 :                         if( poRevCT != NULL )
    1621                 :            :                         {
    1622                 :          1 :                             const double dfX0 = oTM.dfTLX;
    1623                 :          1 :                             const double dfY1 = oTM.dfTLY;
    1624                 :            :                             const double dfX1 = oTM.dfTLX +
    1625                 :          1 :                             oTM.nMatrixWidth  * oTM.dfPixelSize * oTM.nTileWidth;
    1626                 :            :                             const double dfY0 = oTM.dfTLY -
    1627                 :          1 :                             oTM.nMatrixHeight * oTM.dfPixelSize * oTM.nTileHeight;
    1628                 :          1 :                             double dfXMin = std::numeric_limits<double>::infinity();
    1629                 :          1 :                             double dfYMin = std::numeric_limits<double>::infinity();
    1630                 :          1 :                             double dfXMax = -std::numeric_limits<double>::infinity();
    1631                 :          1 :                             double dfYMax = -std::numeric_limits<double>::infinity();
    1632                 :            : 
    1633                 :          1 :                             const int NSTEPS = 20;
    1634         [ +  + ]:          1 :                             for(int i=0;i<=NSTEPS;i++)
    1635                 :            :                             {
    1636                 :          1 :                                 double dfX = dfX0 + (dfX1 - dfX0) * i / NSTEPS;
    1637                 :          1 :                                 double dfY = dfY0;
    1638 [ +  - ][ +  - ]:          1 :                                 if( poRevCT->Transform(1, &dfX, &dfY) )
    1639                 :            :                                 {
    1640         [ +  - ]:          1 :                                     dfXMin = std::min(dfXMin, dfX);
    1641         [ +  - ]:          1 :                                     dfYMin = std::min(dfYMin, dfY);
    1642         [ +  - ]:          1 :                                     dfXMax = std::max(dfXMax, dfX);
    1643         [ +  - ]:          1 :                                     dfYMax = std::max(dfYMax, dfY);
    1644                 :            :                                 }
    1645                 :            : 
    1646                 :          1 :                                 dfX = dfX0 + (dfX1 - dfX0) * i / NSTEPS;
    1647                 :          1 :                                 dfY = dfY1;
    1648 [ +  - ][ +  - ]:          1 :                                 if( poRevCT->Transform(1, &dfX, &dfY) )
    1649                 :            :                                 {
    1650         [ +  - ]:          1 :                                     dfXMin = std::min(dfXMin, dfX);
    1651         [ +  - ]:          1 :                                     dfYMin = std::min(dfYMin, dfY);
    1652         [ +  - ]:          1 :                                     dfXMax = std::max(dfXMax, dfX);
    1653         [ +  - ]:          1 :                                     dfYMax = std::max(dfYMax, dfY);
    1654                 :            :                                 }
    1655                 :            : 
    1656                 :          1 :                                 dfX = dfX0;
    1657                 :          1 :                                 dfY = dfY0 + (dfY1 - dfY0) * i / NSTEPS;
    1658 [ +  - ][ +  - ]:          1 :                                 if( poRevCT->Transform(1, &dfX, &dfY) )
    1659                 :            :                                 {
    1660         [ +  - ]:          1 :                                     dfXMin = std::min(dfXMin, dfX);
    1661         [ +  - ]:          1 :                                     dfYMin = std::min(dfYMin, dfY);
    1662         [ +  - ]:          1 :                                     dfXMax = std::max(dfXMax, dfX);
    1663         [ +  - ]:          1 :                                     dfYMax = std::max(dfYMax, dfY);
    1664                 :            :                                 }
    1665                 :            : 
    1666                 :          1 :                                 dfX = dfX1;
    1667                 :          1 :                                 dfY = dfY0 + (dfY1 - dfY0) * i / NSTEPS;
    1668 [ +  - ][ +  - ]:          1 :                                 if( poRevCT->Transform(1, &dfX, &dfY) )
    1669                 :            :                                 {
    1670         [ +  - ]:          1 :                                     dfXMin = std::min(dfXMin, dfX);
    1671         [ +  - ]:          1 :                                     dfYMin = std::min(dfYMin, dfY);
    1672         [ +  - ]:          1 :                                     dfXMax = std::max(dfXMax, dfX);
    1673         [ +  - ]:          1 :                                     dfYMax = std::max(dfYMax, dfY);
    1674                 :            :                                 }
    1675                 :            :                             }
    1676                 :            : 
    1677 [ +  - ][ +  - ]:          1 :                             delete poRevCT;
    1678                 :            : #ifdef DEBUG_VERBOSE
    1679                 :            :                             CPLDebug("WMTS", "Reprojected densified bbox of most "
    1680                 :            :                                     "precise tile matrix in %s: %.8g %8g %8g %8g",
    1681                 :            :                                     oIter->first.c_str(),
    1682                 :            :                                     dfXMin, dfYMin, dfXMax, dfYMax);
    1683                 :            : #endif
    1684 [ +  - ][ +  + ]:          1 :                             if( fabs(oIter->second.MinX - dfXMin) < 1e-5 *
         [ +  - ][ +  - ]
         [ +  - ][ +  + ]
    1685 [ +  - ][ +  - ]:          1 :                                 std::max(fabs(oIter->second.MinX),fabs(dfXMin)) &&
    1686         [ +  - ]:          1 :                                 fabs(oIter->second.MinY - dfYMin) < 1e-5 *
    1687 [ +  - ][ +  - ]:          1 :                                 std::max(fabs(oIter->second.MinY),fabs(dfYMin)) &&
    1688         [ +  - ]:          1 :                                 fabs(oIter->second.MaxX - dfXMax) < 1e-5 *
    1689 [ +  - ][ +  - ]:          1 :                                 std::max(fabs(oIter->second.MaxX),fabs(dfXMax)) &&
    1690         [ +  - ]:          1 :                                 fabs(oIter->second.MaxY - dfYMax) < 1e-5 *
    1691 [ +  - ][ +  - ]:          1 :                                 std::max(fabs(oIter->second.MaxY),fabs(dfYMax)) )
    1692                 :            :                             {
    1693                 :          1 :                                 bMatchFound = true;
    1694                 :            : #ifdef DEBUG_VERBOSE
    1695                 :            :                                 CPLDebug("WMTS", "Matches layer bounding box, so "
    1696                 :            :                                         "that one is not significant");
    1697                 :            : #endif
    1698                 :            :                                 break;
    1699                 :            :                             }
    1700                 :            :                         }
    1701 [ +  - ][ +  + ]:          1 :                     }
    1702                 :            : 
    1703         [ +  + ]:          1 :                     if( bMatchFound )
    1704                 :            :                     {
    1705         [ -  + ]:          1 :                         if( eExtentMethod == LAYER_BBOX )
    1706                 :          1 :                             eExtentMethod = MOST_PRECISE_TILE_MATRIX;
    1707                 :            :                         break;
    1708                 :            :                     }
    1709                 :            : 
    1710                 :            :                     OGRCoordinateTransformation* poCT =
    1711         [ +  - ]:          1 :                         OGRCreateCoordinateTransformation(&oSRS, &oTMS.oSRS);
    1712         [ +  - ]:          1 :                     if( poCT != NULL )
    1713                 :            :                     {
    1714         [ +  - ]:          1 :                         double dfX1 = oIter->second.MinX;
    1715         [ +  - ]:          1 :                         double dfY1 = oIter->second.MinY;
    1716         [ +  - ]:          1 :                         double dfX2 = oIter->second.MaxX;
    1717         [ +  - ]:          1 :                         double dfY2 = oIter->second.MinY;
    1718         [ +  - ]:          1 :                         double dfX3 = oIter->second.MaxX;
    1719         [ +  - ]:          1 :                         double dfY3 = oIter->second.MaxY;
    1720         [ +  - ]:          1 :                         double dfX4 = oIter->second.MinX;
    1721         [ +  - ]:          1 :                         double dfY4 = oIter->second.MaxY;
    1722 [ +  - ][ +  - ]:          1 :                         if( poCT->Transform(1, &dfX1, &dfY1) &&
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    1723         [ +  - ]:          1 :                             poCT->Transform(1, &dfX2, &dfY2) &&
    1724         [ +  - ]:          1 :                             poCT->Transform(1, &dfX3, &dfY3) &&
    1725         [ +  - ]:          1 :                             poCT->Transform(1, &dfX4, &dfY4) )
    1726                 :            :                         {
    1727         [ +  - ]:          1 :                             sAOI.MinX = std::min(std::min(dfX1, dfX2),
    1728 [ +  - ][ +  - ]:          1 :                                                  std::min(dfX3, dfX4));
    1729         [ +  - ]:          1 :                             sAOI.MinY = std::min(std::min(dfY1, dfY2),
    1730 [ +  - ][ +  - ]:          1 :                                                  std::min(dfY3, dfY4));
    1731         [ +  - ]:          1 :                             sAOI.MaxX = std::max(std::max(dfX1, dfX2),
    1732 [ +  - ][ +  - ]:          1 :                                                  std::max(dfX3, dfX4));
    1733         [ +  - ]:          1 :                             sAOI.MaxY = std::max(std::max(dfY1, dfY2),
    1734 [ +  - ][ +  - ]:          1 :                                                  std::max(dfY3, dfY4));
    1735                 :          1 :                             bHasAOI = TRUE;
    1736                 :            :                         }
    1737 [ +  - ][ +  - ]:          1 :                         delete poCT;
    1738                 :            :                     }
    1739                 :            :                     break;
    1740                 :            :                 }
    1741 [ +  - ][ -  + ]:          1 :             }
    1742                 :            :         }
    1743                 :            : 
    1744                 :            :         // Otherwise default to BoundingBox of the TMS
    1745 [ +  + ][ +  + ]:          1 :         if( !bHasAOI && oTMS.bBoundingBoxValid &&
         [ -  + ][ #  # ]
    1746                 :            :             (eExtentMethod == AUTO || eExtentMethod == TILE_MATRIX_SET) )
    1747                 :            :         {
    1748         [ +  - ]:          1 :             CPLDebug("WMTS", "Using TMS bounding box");
    1749                 :          1 :             sAOI = oTMS.sBoundingBox;
    1750                 :          1 :             bHasAOI = TRUE;
    1751                 :            :         }
    1752                 :            : 
    1753                 :            :         // Otherwise default to implied BoundingBox of the most precise TM
    1754 [ +  + ][ -  + ]:          1 :         if( !bHasAOI &&
                 [ #  # ]
    1755                 :            :             (eExtentMethod == AUTO || eExtentMethod == MOST_PRECISE_TILE_MATRIX) )
    1756                 :            :         {
    1757         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM.back();
    1758 [ +  - ][ +  - ]:          1 :             CPLDebug("WMTS", "Using TM level %s bounding box", oTM.osIdentifier.c_str() );
    1759                 :            : 
    1760                 :          1 :             sAOI.MinX = oTM.dfTLX;
    1761                 :          1 :             sAOI.MaxY = oTM.dfTLY;
    1762                 :          1 :             sAOI.MaxX = oTM.dfTLX + oTM.nMatrixWidth  * oTM.dfPixelSize * oTM.nTileWidth;
    1763                 :          1 :             sAOI.MinY = oTM.dfTLY - oTM.nMatrixHeight * oTM.dfPixelSize * oTM.nTileHeight;
    1764                 :          1 :             bHasAOI = TRUE;
    1765                 :            :         }
    1766                 :            : 
    1767         [ -  + ]:          1 :         if( !bHasAOI )
    1768                 :            :         {
    1769                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
    1770         [ #  # ]:          0 :                      "Could not determine raster extent");
    1771         [ #  # ]:          0 :             CPLDestroyXMLNode(psXML);
    1772 [ #  # ][ #  # ]:          0 :             delete poDS;
    1773                 :          0 :             return NULL;
    1774                 :            :         }
    1775                 :            : 
    1776                 :            :         {
    1777                 :            :             // Clip with implied BoundingBox of the most precise TM
    1778                 :            :             // Useful for http://tileserver.maptiler.com/wmts
    1779         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM.back();
    1780                 :            : 
    1781                 :            :             // For https://data.linz.govt.nz/services;key=XXXXXXXX/wmts/1.0.0/set/69/WMTSCapabilities.xml
    1782                 :            :             // only clip in Y since there's a warp over dateline.
    1783                 :            :             // Update: it sems that the content of the server has changed since
    1784                 :            :             // initial coding. So do X clipping in default mode.
    1785         [ +  + ]:          1 :             if( !bExtendBeyondDateLine )
    1786                 :            :             {
    1787         [ +  - ]:          1 :                 sAOI.MinX = std::max(sAOI.MinX, oTM.dfTLX);
    1788                 :            :                 sAOI.MaxX = std::min(sAOI.MaxX,
    1789                 :            :                     oTM.dfTLX +
    1790         [ +  - ]:          1 :                     oTM.nMatrixWidth  * oTM.dfPixelSize * oTM.nTileWidth);
    1791                 :            :             }
    1792         [ +  - ]:          1 :             sAOI.MaxY = std::min(sAOI.MaxY, oTM.dfTLY);
    1793                 :            :             sAOI.MinY =
    1794                 :            :                 std::max(sAOI.MinY,
    1795                 :            :                          oTM.dfTLY -
    1796         [ +  - ]:          1 :                          oTM.nMatrixHeight * oTM.dfPixelSize * oTM.nTileHeight);
    1797                 :            :         }
    1798                 :            : 
    1799                 :            :         // Clip with limits of most precise TM when available
    1800                 :            :         {
    1801         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM.back();
    1802 [ +  - ][ +  - ]:          1 :             if( aoMapTileMatrixLimits.find(oTM.osIdentifier) != aoMapTileMatrixLimits.end() )
         [ +  - ][ +  + ]
    1803                 :            :             {
    1804         [ +  - ]:          1 :                 const WMTSTileMatrixLimits& oTMLimits = aoMapTileMatrixLimits[oTM.osIdentifier];
    1805                 :          1 :                 double dfTileWidthUnits = oTM.dfPixelSize * oTM.nTileWidth;
    1806                 :          1 :                 double dfTileHeightUnits = oTM.dfPixelSize * oTM.nTileHeight;
    1807         [ +  - ]:          1 :                 sAOI.MinX = std::max(sAOI.MinX, oTM.dfTLX + oTMLimits.nMinTileCol * dfTileWidthUnits);
    1808         [ +  - ]:          1 :                 sAOI.MaxY = std::min(sAOI.MaxY, oTM.dfTLY - oTMLimits.nMinTileRow * dfTileHeightUnits);
    1809         [ +  - ]:          1 :                 sAOI.MaxX = std::min(sAOI.MaxX, oTM.dfTLX + (oTMLimits.nMaxTileCol + 1) * dfTileWidthUnits);
    1810         [ +  - ]:          1 :                 sAOI.MinY = std::max(sAOI.MinY, oTM.dfTLY - (oTMLimits.nMaxTileRow + 1) * dfTileHeightUnits);
    1811                 :            :             }
    1812                 :            :         }
    1813                 :            : 
    1814                 :            :         // Establish raster dimension and extent
    1815         [ +  - ]:          1 :         int nMaxZoomLevel = (int)oTMS.aoTM.size()-1;
    1816         [ +  - ]:          1 :         while(nMaxZoomLevel >= 0)
    1817                 :            :         {
    1818         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM[nMaxZoomLevel];
    1819                 :          1 :             double dfRasterXSize = (sAOI.MaxX - sAOI.MinX) / oTM.dfPixelSize;
    1820                 :          1 :             double dfRasterYSize = (sAOI.MaxY - sAOI.MinY) / oTM.dfPixelSize;
    1821 [ +  + ][ +  - ]:          1 :             if( dfRasterXSize < INT_MAX && dfRasterYSize < INT_MAX )
    1822                 :            :             {
    1823 [ +  - ][ +  + ]:          1 :                 if( nMaxZoomLevel != (int)oTMS.aoTM.size()-1 )
    1824                 :            :                 {
    1825                 :            :                     CPLDebug("WMTS", "Using zoom level %s instead of %s to avoid int overflow",
    1826         [ +  - ]:          1 :                              oTMS.aoTM[nMaxZoomLevel].osIdentifier.c_str(),
    1827 [ +  - ][ +  - ]:          1 :                              oTMS.aoTM.back().osIdentifier.c_str());
         [ +  - ][ +  - ]
    1828                 :            :                 }
    1829                 :            : 
    1830                 :            :                 // Align AOI on pixel boundaries with respect to TopLeftCorner of
    1831                 :            :                 // this tile matrix
    1832                 :          1 :                 poDS->adfGT[0] = oTM.dfTLX + floor((sAOI.MinX - oTM.dfTLX) / oTM.dfPixelSize+1e-10) * oTM.dfPixelSize;
    1833                 :          1 :                 poDS->adfGT[1] = oTM.dfPixelSize;
    1834                 :          1 :                 poDS->adfGT[2] = 0.0;
    1835                 :          1 :                 poDS->adfGT[3] = oTM.dfTLY + ceil((sAOI.MaxY - oTM.dfTLY) / oTM.dfPixelSize-1e-10) * oTM.dfPixelSize;
    1836                 :          1 :                 poDS->adfGT[4] = 0.0;
    1837                 :          1 :                 poDS->adfGT[5] = -oTM.dfPixelSize;
    1838                 :          1 :                 poDS->nRasterXSize = int(0.5 + (sAOI.MaxX - poDS->adfGT[0]) / oTM.dfPixelSize);
    1839                 :          1 :                 poDS->nRasterYSize = int(0.5 + (poDS->adfGT[3] - sAOI.MinY) / oTM.dfPixelSize);
    1840                 :          1 :                 break;
    1841                 :            :             }
    1842                 :          1 :             nMaxZoomLevel --;
    1843                 :            :         }
    1844         [ -  + ]:          1 :         if( nMaxZoomLevel < 0 )
    1845                 :            :         {
    1846                 :            :             CPLError(CE_Failure, CPLE_AppDefined,
    1847         [ #  # ]:          0 :                      "No zoom level in tile matrix set found");
    1848         [ #  # ]:          0 :             CPLDestroyXMLNode(psXML);
    1849 [ #  # ][ #  # ]:          0 :             delete poDS;
    1850                 :          0 :             return NULL;
    1851                 :            :         }
    1852                 :            :         CPLDebug("WMTS", "Using tilematrix=%s (zoom level %d)",
    1853 [ +  - ][ +  - ]:          1 :                  oTMS.aoTM[nMaxZoomLevel].osIdentifier.c_str(), nMaxZoomLevel);
                 [ +  - ]
    1854 [ +  - ][ +  - ]:          1 :         oTMS.aoTM.resize(1 + nMaxZoomLevel);
                 [ +  - ]
    1855         [ +  - ]:          1 :         poDS->oTMS = oTMS;
    1856                 :            : 
    1857 [ +  - ][ -  + ]:          1 :         if( !osProjection.empty() )
    1858                 :            :         {
    1859         [ #  # ]:          0 :             OGRSpatialReference oSRS;
    1860 [ #  # ][ #  # ]:          0 :             if( oSRS.SetFromUserInput(osProjection) == OGRERR_NONE )
                 [ #  # ]
    1861                 :            :             {
    1862                 :          0 :                 char* pszWKT = NULL;
    1863         [ #  # ]:          0 :                 oSRS.exportToWkt(&pszWKT);
    1864 [ #  # ][ #  # ]:          0 :                 poDS->osProjection = pszWKT;
                 [ #  # ]
    1865         [ #  # ]:          0 :                 CPLFree(pszWKT);
    1866         [ #  # ]:          0 :             }
    1867                 :            :         }
    1868 [ +  - ][ +  - ]:          1 :         if( poDS->osProjection.empty() )
    1869                 :            :         {
    1870                 :            :             // Strip AXIS
    1871         [ +  - ]:          1 :             OGR_SRSNode *poGEOGCS = oTMS.oSRS.GetAttrNode( "GEOGCS" );
    1872         [ +  - ]:          1 :             if( poGEOGCS != NULL )
    1873         [ +  - ]:          1 :                 poGEOGCS->StripNodes( "AXIS" );
    1874                 :            : 
    1875         [ +  - ]:          1 :             OGR_SRSNode *poPROJCS = oTMS.oSRS.GetAttrNode( "PROJCS" );
    1876 [ +  + ][ +  - ]:          1 :             if (poPROJCS != NULL && oTMS.oSRS.EPSGTreatsAsNorthingEasting())
         [ -  + ][ -  + ]
    1877         [ #  # ]:          0 :                 poPROJCS->StripNodes( "AXIS" );
    1878                 :            : 
    1879                 :          1 :             char* pszWKT = NULL;
    1880         [ +  - ]:          1 :             oTMS.oSRS.exportToWkt(&pszWKT);
    1881 [ +  - ][ +  - ]:          1 :             poDS->osProjection = pszWKT;
                 [ +  - ]
    1882         [ +  - ]:          1 :             CPLFree(pszWKT);
    1883                 :            :         }
    1884                 :            : 
    1885 [ +  - ][ +  + ]:          1 :         if( osURLTileTemplate.empty() )
    1886                 :            :         {
    1887 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = GetOperationKVPURL(psXML, "GetTile");
                 [ +  - ]
    1888 [ +  - ][ +  + ]:          1 :             if( osURLTileTemplate.empty() )
    1889                 :            :             {
    1890                 :            :                 CPLError(CE_Failure, CPLE_AppDefined,
    1891         [ +  - ]:          1 :                          "No RESTful nor KVP GetTile operation found");
    1892         [ +  - ]:          1 :                 CPLDestroyXMLNode(psXML);
    1893 [ +  - ][ +  - ]:          1 :                 delete poDS;
    1894                 :          1 :                 return NULL;
    1895                 :            :             }
    1896 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "service", "WMTS");
         [ +  - ][ +  - ]
    1897 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "request", "GetTile");
         [ +  - ][ +  - ]
    1898 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "version", "1.0.0");
         [ +  - ][ +  - ]
    1899 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "layer", osSelectLayer);
         [ +  - ][ +  - ]
                 [ +  - ]
    1900 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "style", osSelectStyle);
         [ +  - ][ +  - ]
                 [ +  - ]
    1901 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "format", osSelectTileFormat);
         [ +  - ][ +  - ]
                 [ +  - ]
    1902 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate, "TileMatrixSet", osSelectTMS);
         [ +  - ][ +  - ]
                 [ +  - ]
    1903         [ +  - ]:          1 :             osURLTileTemplate += "&TileMatrix={TileMatrix}";
    1904         [ +  - ]:          1 :             osURLTileTemplate += "&TileRow=${y}";
    1905         [ +  - ]:          1 :             osURLTileTemplate += "&TileCol=${x}";
    1906                 :            : 
    1907         [ +  - ]:          1 :             std::map<CPLString,CPLString>::iterator oIter = aoMapDimensions.begin();
    1908 [ +  - ][ +  - ]:          1 :             for(; oIter != aoMapDimensions.end(); ++oIter )
         [ +  - ][ +  + ]
    1909                 :            :             {
    1910                 :            :                 osURLTileTemplate = CPLURLAddKVP(osURLTileTemplate,
    1911         [ +  - ]:          1 :                                                  oIter->first,
    1912 [ +  - ][ +  - ]:          1 :                                                  oIter->second);
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1913                 :            :             }
    1914                 :            :             //CPLDebug("WMTS", "osURLTileTemplate = %s", osURLTileTemplate.c_str());
    1915                 :            :         }
    1916                 :            :         else
    1917                 :            :         {
    1918 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = Replace(osURLTileTemplate, "{Style}", osSelectStyle);
         [ +  - ][ +  - ]
    1919 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = Replace(osURLTileTemplate, "{TileMatrixSet}", osSelectTMS);
         [ +  - ][ +  - ]
    1920 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = Replace(osURLTileTemplate, "{TileCol}", "${x}");
                 [ +  - ]
    1921 [ +  - ][ +  - ]:          1 :             osURLTileTemplate = Replace(osURLTileTemplate, "{TileRow}", "${y}");
                 [ +  - ]
    1922                 :            : 
    1923         [ +  - ]:          1 :             std::map<CPLString,CPLString>::iterator oIter = aoMapDimensions.begin();
    1924 [ +  - ][ +  - ]:          1 :             for(; oIter != aoMapDimensions.end(); ++oIter )
         [ +  - ][ +  + ]
    1925                 :            :             {
    1926                 :            :                 osURLTileTemplate = Replace(osURLTileTemplate,
    1927         [ +  - ]:          1 :                                             CPLSPrintf("{%s}", oIter->first.c_str()),
    1928 [ +  - ][ +  - ]:          1 :                                             oIter->second);
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1929                 :            :             }
    1930                 :            :         }
    1931                 :            : 
    1932 [ +  - ][ +  + ]:          1 :         if( osURLFeatureInfoTemplate.empty() && !osSelectInfoFormat.empty() )
         [ +  - ][ +  + ]
                 [ +  + ]
    1933                 :            :         {
    1934 [ +  - ][ +  - ]:          1 :             osURLFeatureInfoTemplate = GetOperationKVPURL(psXML, "GetFeatureInfo");
                 [ +  - ]
    1935 [ +  - ][ +  - ]:          1 :             if( !osURLFeatureInfoTemplate.empty() )
    1936                 :            :             {
    1937 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "service", "WMTS");
         [ +  - ][ +  - ]
    1938 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "request", "GetFeatureInfo");
         [ +  - ][ +  - ]
    1939 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "version", "1.0.0");
         [ +  - ][ +  - ]
    1940 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "layer", osSelectLayer);
         [ +  - ][ +  - ]
                 [ +  - ]
    1941 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "style", osSelectStyle);
         [ +  - ][ +  - ]
                 [ +  - ]
    1942                 :            :                 //osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "format", osSelectTileFormat);
    1943 [ +  - ][ +  - ]:          1 :                 osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate, "InfoFormat", osSelectInfoFormat);
         [ +  - ][ +  - ]
                 [ +  - ]
    1944         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&TileMatrixSet={TileMatrixSet}";
    1945         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&TileMatrix={TileMatrix}";
    1946         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&TileRow={TileRow}";
    1947         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&TileCol={TileCol}";
    1948         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&J={J}";
    1949         [ +  - ]:          1 :                 osURLFeatureInfoTemplate += "&I={I}";
    1950                 :            : 
    1951         [ +  - ]:          1 :                 std::map<CPLString,CPLString>::iterator oIter = aoMapDimensions.begin();
    1952 [ +  - ][ +  - ]:          1 :                 for(; oIter != aoMapDimensions.end(); ++oIter )
         [ +  - ][ +  + ]
    1953                 :            :                 {
    1954                 :            :                     osURLFeatureInfoTemplate = CPLURLAddKVP(osURLFeatureInfoTemplate,
    1955         [ +  - ]:          1 :                                                             oIter->first,
    1956 [ +  - ][ +  - ]:          1 :                                                             oIter->second);
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1957                 :            :                 }
    1958                 :            :                 //CPLDebug("WMTS", "osURLFeatureInfoTemplate = %s", osURLFeatureInfoTemplate.c_str());
    1959                 :            :             }
    1960                 :            :         }
    1961                 :            :         else
    1962                 :            :         {
    1963 [ +  - ][ +  - ]:          1 :              osURLFeatureInfoTemplate = Replace(osURLFeatureInfoTemplate, "{Style}", osSelectStyle);
         [ +  - ][ +  - ]
    1964                 :            : 
    1965         [ +  - ]:          1 :             std::map<CPLString,CPLString>::iterator oIter = aoMapDimensions.begin();
    1966 [ +  - ][ +  - ]:          1 :             for(; oIter != aoMapDimensions.end(); ++oIter )
         [ +  - ][ +  + ]
    1967                 :            :             {
    1968                 :            :                 osURLFeatureInfoTemplate = Replace(osURLFeatureInfoTemplate,
    1969         [ +  - ]:          1 :                                             CPLSPrintf("{%s}", oIter->first.c_str()),
    1970 [ +  - ][ +  - ]:          1 :                                             oIter->second);
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    1971                 :            :             }
    1972                 :            :         }
    1973         [ +  - ]:          1 :         poDS->osURLFeatureInfoTemplate = osURLFeatureInfoTemplate;
    1974                 :            : 
    1975                 :            :         // Build all TMS datasets, wrapped in VRT datasets
    1976         [ +  + ]:          1 :         for(int i=nMaxZoomLevel;i>=0;i--)
    1977                 :            :         {
    1978         [ +  - ]:          1 :             const WMTSTileMatrix& oTM = oTMS.aoTM[i];
    1979                 :          1 :             int nRasterXSize = int(0.5 + poDS->nRasterXSize / oTM.dfPixelSize * poDS->adfGT[1]);
    1980                 :          1 :             int nRasterYSize = int(0.5 + poDS->nRasterYSize / oTM.dfPixelSize * poDS->adfGT[1]);
    1981 [ +  - ][ +  + ]:          1 :             if( !poDS->apoDatasets.empty() &&
         [ +  - ][ +  + ]
                 [ +  + ]
    1982                 :            :                 (nRasterXSize < 128 || nRasterYSize < 128) )
    1983                 :            :             {
    1984                 :          1 :                 break;
    1985                 :            :             }
    1986 [ +  - ][ +  - ]:          1 :             CPLString osURL(Replace(osURLTileTemplate, "{TileMatrix}", oTM.osIdentifier));
    1987                 :            : 
    1988                 :          1 :             double dfTileWidthUnits = oTM.dfPixelSize * oTM.nTileWidth;
    1989                 :          1 :             double dfTileHeightUnits = oTM.dfPixelSize * oTM.nTileHeight;
    1990                 :            : 
    1991                 :            :             // Compute the shift in terms of tiles between AOI and TM origin
    1992                 :          1 :             int nTileX = (int)(floor(poDS->adfGT[0] - oTM.dfTLX + 1e-10) / dfTileWidthUnits);
    1993                 :          1 :             int nTileY = (int)(floor(oTM.dfTLY - poDS->adfGT[3] + 1e-10) / dfTileHeightUnits);
    1994                 :            : 
    1995                 :            :             // Compute extent of this zoom level slightly larger than the AOI and
    1996                 :            :             // aligned on tile boundaries at this TM
    1997                 :          1 :             double dfULX = oTM.dfTLX + nTileX * dfTileWidthUnits;
    1998                 :          1 :             double dfULY = oTM.dfTLY - nTileY * dfTileHeightUnits;
    1999                 :          1 :             double dfLRX = poDS->adfGT[0] + poDS->nRasterXSize * poDS->adfGT[1];
    2000                 :          1 :             double dfLRY = poDS->adfGT[3] + poDS->nRasterYSize * poDS->adfGT[5];
    2001                 :          1 :             dfLRX = dfULX + ceil((dfLRX - dfULX) / dfTileWidthUnits - 1e-10) * dfTileWidthUnits;
    2002                 :          1 :             dfLRY = dfULY + floor((dfLRY - dfULY) / dfTileHeightUnits + 1e-10) * dfTileHeightUnits;
    2003                 :            : 
    2004                 :          1 :             int nSizeX = int(0.5+(dfLRX - dfULX) / oTM.dfPixelSize);
    2005                 :          1 :             int nSizeY = int(0.5+(dfULY - dfLRY) / oTM.dfPixelSize);
    2006                 :            : 
    2007                 :          1 :             double dfDateLineX = oTM.dfTLX + oTM.nMatrixWidth * dfTileWidthUnits;
    2008                 :          1 :             int nSizeX1 = int(0.5+(dfDateLineX - dfULX) / oTM.dfPixelSize);
    2009                 :          1 :             int nSizeX2 = int(0.5+(dfLRX - dfDateLineX) / oTM.dfPixelSize);
    2010 [ +  + ][ -  + ]:          1 :             if( bExtendBeyondDateLine && dfDateLineX > dfLRX )
    2011                 :            :             {
    2012         [ #  # ]:          0 :                 CPLDebug("WMTS", "ExtendBeyondDateLine ignored in that case");
    2013                 :          0 :                 bExtendBeyondDateLine = FALSE;
    2014                 :            :             }
    2015                 :            : 
    2016                 :            : #define WMS_TMS_TEMPLATE \
    2017                 :            :     "<GDAL_WMS>" \
    2018                 :            :     "<Service name=\"TMS\">" \
    2019                 :            :     "    <ServerUrl>%s</ServerUrl>" \
    2020                 :            :     "</Service>" \
    2021                 :            :     "<DataWindow>" \
    2022                 :            :     "    <UpperLeftX>%.16g</UpperLeftX>" \
    2023                 :            :     "    <UpperLeftY>%.16g</UpperLeftY>" \
    2024                 :            :     "    <LowerRightX>%.16g</LowerRightX>" \
    2025                 :            :     "    <LowerRightY>%.16g</LowerRightY>" \
    2026                 :            :     "    <TileLevel>0</TileLevel>" \
    2027                 :            :     "    <TileX>%d</TileX>" \
    2028                 :            :     "    <TileY>%d</TileY>" \
    2029                 :            :     "    <SizeX>%d</SizeX>" \
    2030                 :            :     "    <SizeY>%d</SizeY>" \
    2031                 :            :     "    <YOrigin>top</YOrigin>" \
    2032                 :            :     "</DataWindow>" \
    2033                 :            :     "<BlockSizeX>%d</BlockSizeX>" \
    2034                 :            :     "<BlockSizeY>%d</BlockSizeY>" \
    2035                 :            :     "<BandsCount>%d</BandsCount>" \
    2036                 :            :     "%s" \
    2037                 :            : "</GDAL_WMS>"
    2038                 :            : 
    2039                 :            :             CPLString osStr(CPLSPrintf( WMS_TMS_TEMPLATE,
    2040                 :            :                 WMTSEscapeXML(osURL).c_str(),
    2041                 :            :                 dfULX, dfULY, (bExtendBeyondDateLine) ? dfDateLineX : dfLRX, dfLRY,
    2042                 :            :                 nTileX, nTileY, (bExtendBeyondDateLine) ? nSizeX1 : nSizeX, nSizeY,
    2043                 :            :                 oTM.nTileWidth, oTM.nTileHeight, nBands,
    2044 [ +  - ][ +  + ]:          1 :                 osOtherXML.c_str()));
         [ +  + ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    2045                 :            :             GDALDataset* poWMSDS = (GDALDataset*)GDALOpenEx(
    2046 [ +  - ][ +  - ]:          1 :                 osStr, GDAL_OF_RASTER | GDAL_OF_SHARED | GDAL_OF_VERBOSE_ERROR, NULL, NULL, NULL);
    2047         [ -  + ]:          1 :             if( poWMSDS == NULL )
    2048                 :            :             {
    2049         [ #  # ]:          0 :                 CPLDestroyXMLNode(psXML);
    2050 [ #  # ][ #  # ]:          0 :                 delete poDS;
    2051                 :          0 :                 return NULL;
    2052                 :            :             }
    2053                 :            : 
    2054         [ +  - ]:          1 :             VRTDatasetH hVRTDS = VRTCreate( nRasterXSize, nRasterYSize );
    2055         [ +  + ]:          1 :             for(int iBand=1;iBand<=nBands;iBand++)
    2056                 :            :             {
    2057         [ +  - ]:          1 :                 VRTAddBand( hVRTDS, GDT_Byte, NULL );
    2058                 :            :             }
    2059                 :            : 
    2060                 :            :             int nSrcXOff, nSrcYOff, nDstXOff, nDstYOff;
    2061                 :            : 
    2062                 :          1 :             nSrcXOff = 0;
    2063                 :          1 :             nDstXOff = (int)(0.5 + (dfULX - poDS->adfGT[0]) / oTM.dfPixelSize);
    2064                 :            : 
    2065                 :          1 :             nSrcYOff = 0;
    2066                 :          1 :             nDstYOff = (int)(0.5 + (poDS->adfGT[3] - dfULY) / oTM.dfPixelSize);
    2067                 :            : 
    2068         [ +  + ]:          1 :             if( bExtendBeyondDateLine )
    2069                 :            :             {
    2070                 :            :                 int nSrcXOff2, nDstXOff2;
    2071                 :            : 
    2072                 :          1 :                 nSrcXOff2 = 0;
    2073                 :          1 :                 nDstXOff2 = (int)(0.5 + (dfDateLineX - poDS->adfGT[0]) / oTM.dfPixelSize);
    2074                 :            : 
    2075                 :            :                 osStr = CPLSPrintf( WMS_TMS_TEMPLATE,
    2076                 :            :                     WMTSEscapeXML(osURL).c_str(),
    2077                 :            :                     -dfDateLineX, dfULY, dfLRX - 2 * dfDateLineX, dfLRY,
    2078                 :            :                     0, nTileY, nSizeX2, nSizeY,
    2079                 :            :                     oTM.nTileWidth, oTM.nTileHeight, nBands,
    2080 [ +  - ][ +  - ]:          1 :                     osOtherXML.c_str());
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    2081                 :            : 
    2082                 :            :                 GDALDataset* poWMSDS2 = (GDALDataset*)GDALOpenEx(
    2083 [ +  - ][ +  - ]:          1 :                     osStr, GDAL_OF_RASTER | GDAL_OF_SHARED, NULL, NULL, NULL);
    2084         [ -  + ]:          1 :                 CPLAssert(poWMSDS2);
    2085                 :            : 
    2086         [ +  + ]:          1 :                 for(int iBand=1;iBand<=nBands;iBand++)
    2087                 :            :                 {
    2088                 :            :                     VRTSourcedRasterBandH hVRTBand =
    2089         [ +  - ]:          1 :                         (VRTSourcedRasterBandH) GDALGetRasterBand(hVRTDS, iBand);
    2090                 :            :                     VRTAddSimpleSource( hVRTBand, GDALGetRasterBand(poWMSDS, iBand),
    2091                 :            :                                         nSrcXOff, nSrcYOff, nSizeX1, nSizeY,
    2092                 :            :                                         nDstXOff, nDstYOff, nSizeX1, nSizeY,
    2093 [ +  - ][ +  - ]:          1 :                                         "NEAR", VRT_NODATA_UNSET);
    2094                 :            :                     VRTAddSimpleSource( hVRTBand, GDALGetRasterBand(poWMSDS2, iBand),
    2095                 :            :                                         nSrcXOff2, nSrcYOff, nSizeX2, nSizeY,
    2096                 :            :                                         nDstXOff2, nDstYOff, nSizeX2, nSizeY,
    2097 [ +  - ][ +  - ]:          1 :                                         "NEAR", VRT_NODATA_UNSET);
    2098                 :            :                 }
    2099                 :            : 
    2100         [ +  - ]:          1 :                 poWMSDS2->Dereference();
    2101                 :            :             }
    2102                 :            :             else
    2103                 :            :             {
    2104         [ +  + ]:          1 :                 for(int iBand=1;iBand<=nBands;iBand++)
    2105                 :            :                 {
    2106                 :            :                     VRTSourcedRasterBandH hVRTBand =
    2107         [ +  - ]:          1 :                         (VRTSourcedRasterBandH) GDALGetRasterBand(hVRTDS, iBand);
    2108                 :            :                     VRTAddSimpleSource( hVRTBand, GDALGetRasterBand(poWMSDS, iBand),
    2109                 :            :                                         nSrcXOff, nSrcYOff, nSizeX, nSizeY,
    2110                 :            :                                         nDstXOff, nDstYOff, nSizeX, nSizeY,
    2111 [ +  - ][ +  - ]:          1 :                                         "NEAR", VRT_NODATA_UNSET);
    2112                 :            :                 }
    2113                 :            :             }
    2114                 :            : 
    2115         [ +  - ]:          1 :             poWMSDS->Dereference();
    2116                 :            : 
    2117 [ +  - ][ -  + ]:          1 :             poDS->apoDatasets.push_back((GDALDataset*)hVRTDS);
    2118 [ +  - ][ +  - ]:          1 :         }
                 [ -  + ]
    2119                 :            : 
    2120 [ +  - ][ -  + ]:          1 :         if( poDS->apoDatasets.empty() )
    2121                 :            :         {
    2122         [ #  # ]:          0 :             CPLError(CE_Failure, CPLE_AppDefined, "No zoom level found");
    2123         [ #  # ]:          0 :             CPLDestroyXMLNode(psXML);
    2124 [ #  # ][ #  # ]:          0 :             delete poDS;
    2125                 :          0 :             return NULL;
    2126                 :            :         }
    2127                 :            : 
    2128         [ +  - ]:          1 :         poDS->SetMetadataItem("INTERLEAVE", "PIXEL", "IMAGE_STRUCTURE");
    2129         [ +  + ]:          1 :         for(int i=0;i<nBands;i++)
    2130 [ +  - ][ +  - ]:          1 :             poDS->SetBand(i+1, new WMTSBand(poDS,i+1));
                 [ +  - ]
    2131                 :            : 
    2132 [ +  - ][ +  - ]:          1 :         poDS->osXML = "<GDAL_WMTS>\n";
                 [ +  - ]
    2133                 :            :         poDS->osXML += "  <GetCapabilitiesUrl>" +
    2134                 :            :                      WMTSEscapeXML(osGetCapabilitiesURL) +
    2135 [ +  - ][ +  - ]:          1 :                      "</GetCapabilitiesUrl>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2136 [ +  - ][ +  + ]:          1 :         if( !osSelectLayer.empty() )
    2137 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <Layer>" + WMTSEscapeXML(osSelectLayer) + "</Layer>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2138 [ +  - ][ +  + ]:          1 :         if( !osSelectStyle.empty() )
    2139 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <Style>" + WMTSEscapeXML(osSelectStyle) + "</Style>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2140 [ +  - ][ +  + ]:          1 :         if( !osSelectTMS.empty() )
    2141 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <TileMatrixSet>" + WMTSEscapeXML(osSelectTMS) + "</TileMatrixSet>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2142 [ +  - ][ +  + ]:          1 :         if( !osMaxTileMatrixIdentifier.empty() )
    2143 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <TileMatrix>" + WMTSEscapeXML(osMaxTileMatrixIdentifier) + "</TileMatrix>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2144         [ +  + ]:          1 :         if( nUserMaxZoomLevel >= 0 )
    2145 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <ZoomLevel>" + CPLString().Printf("%d", nUserMaxZoomLevel) + "</ZoomLevel>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2146 [ +  + ][ +  - ]:          1 :         if( nCountTileFormat > 1 && !osSelectTileFormat.empty() )
         [ +  - ][ +  + ]
    2147 [ +  - ][ +  - ]:          1 :             poDS->osXML += "  <Format>" + WMTSEscapeXML(osSelectTileFormat) + "</Format>\n";
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
    2148 [ -  + ][ #  # ]:          1 :         if( nCountInfoFormat > 1 && !osSelectInfoFormat.empty() )
         [ #  # ][ -  + ]
    2149 [ #  # ][ #  # ]:          0 :             poDS->osXML += "  <InfoFormat>" + WMTSEscapeXML(osSelectInfoFormat) + "</InfoFormat>\n";
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
    2150         [ +  - ]:          1 :         poDS->osXML += "  <DataWindow>\n";
    2151                 :            :         poDS->osXML += CPLSPrintf("    <UpperLeftX>%.16g</UpperLeftX>\n",
    2152 [ +  - ][ +  - ]:          1 :                              poDS->adfGT[0]);
    2153                 :            :         poDS->osXML += CPLSPrintf("    <UpperLeftY>%.16g</UpperLeftY>\n",
    2154 [ +  - ][ +  - ]:          1 :                              poDS->adfGT[3]);
    2155                 :            :         poDS->osXML += CPLSPrintf("    <LowerRightX>%.16g</LowerRightX>\n",
    2156 [ +  - ][ +  - ]:          1 :                              poDS->adfGT[0] +  poDS->adfGT[1] *  poDS->nRasterXSize);
    2157                 :            :         poDS->osXML += CPLSPrintf("    <LowerRightY>%.16g</LowerRightY>\n",
    2158 [ +  - ][ +  - ]:          1 :                              poDS->adfGT[3] +  poDS->adfGT[5] *  poDS->nRasterYSize);
    2159         [ +  - ]:          1 :         poDS->osXML += "  </DataWindow>\n";
    2160         [ +  + ]:          1 :         if( bExtendBeyondDateLine )
    2161         [ +  - ]:          1 :             poDS->osXML += "  <ExtendBeyondDateLine>true</ExtendBeyondDateLine>\n";
    2162 [ +  - ][ +  - ]:          1 :         poDS->osXML += CPLSPrintf("  <BandsCount>%d</BandsCount>\n", nBands);
    2163         [ +  - ]:          1 :         poDS->osXML += "  <Cache />\n";
    2164         [ +  - ]:          1 :         poDS->osXML += "  <UnsafeSSL>true</UnsafeSSL>\n";
    2165         [ +  - ]:          1 :         poDS->osXML += "  <ZeroBlockHttpCodes>204,404</ZeroBlockHttpCodes>\n";
    2166         [ +  - ]:          1 :         poDS->osXML += "  <ZeroBlockOnServerException>true</ZeroBlockOnServerException>\n";
    2167 [ +  - ][ +  - ]:          1 :         poDS->osXML += "</GDAL_WMTS>\n";
                 [ +  + ]
    2168                 :            :     }
    2169                 :            : 
    2170         [ +  - ]:          1 :     CPLDestroyXMLNode(psXML);
    2171                 :            : 
    2172         [ +  - ]:          1 :     poDS->SetPamFlags(0);
    2173 [ +  - ][ +  - ]:          1 :     return poDS;
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
         [ +  - ][ +  - ]
                 [ +  - ]
    2174                 :            : }
    2175                 :            : /************************************************************************/
    2176                 :            : /*                             CreateCopy()                             */
    2177                 :            : /************************************************************************/
    2178                 :            : 
    2179                 :          1 : GDALDataset *WMTSDataset::CreateCopy( const char * pszFilename,
    2180                 :            :                                          GDALDataset *poSrcDS,
    2181                 :            :                                          CPL_UNUSED int bStrict,
    2182                 :            :                                          CPL_UNUSED char ** papszOptions,
    2183                 :            :                                          CPL_UNUSED GDALProgressFunc pfnProgress,
    2184                 :            :                                          CPL_UNUSED void * pProgressData )
    2185                 :            : {
    2186   [ +  -  +  + ]:          1 :     if( poSrcDS->GetDriver() == NULL ||
                 [ +  + ]
    2187                 :          1 :         poSrcDS->GetDriver() != GDALGetDriverByName("WMTS") )
    2188                 :            :     {
    2189                 :            :         CPLError(CE_Failure, CPLE_NotSupported,
    2190                 :          1 :                  "Source dataset must be a WMTS dataset");
    2191                 :          1 :         return NULL;
    2192                 :            :     }
    2193                 :            : 
    2194                 :          1 :     const char* pszXML = poSrcDS->GetMetadataItem("XML", "WMTS");
    2195         [ -  + ]:          1 :     if (pszXML == NULL)
    2196                 :            :     {
    2197                 :            :         CPLError(CE_Failure, CPLE_AppDefined,
    2198                 :          0 :                  "Cannot get XML definition of source WMTS dataset");
    2199                 :          0 :         return NULL;
    2200                 :            :     }
    2201                 :            : 
    2202                 :          1 :     VSILFILE* fp = VSIFOpenL(pszFilename, "wb");
    2203         [ -  + ]:          1 :     if (fp == NULL)
    2204                 :          0 :         return NULL;
    2205                 :            : 
    2206                 :          1 :     VSIFWriteL(pszXML, 1, strlen(pszXML), fp);
    2207                 :          1 :     VSIFCloseL(fp);
    2208                 :            : 
    2209                 :          1 :     GDALOpenInfo oOpenInfo(pszFilename, GA_ReadOnly);
    2210         [ +  - ]:          1 :     return Open(&oOpenInfo);
    2211                 :            : }
    2212                 :            : 
    2213                 :            : /************************************************************************/
    2214                 :            : /*                       GDALRegister_WMTS()                            */
    2215                 :            : /************************************************************************/
    2216                 :            : 
    2217                 :          1 : void GDALRegister_WMTS()
    2218                 :            : 
    2219                 :            : {
    2220         [ -  + ]:          1 :     if( !GDAL_CHECK_VERSION( "WMTS driver" ) )
    2221                 :          0 :         return;
    2222                 :            : 
    2223         [ +  + ]:          1 :     if( GDALGetDriverByName( "WMTS" ) != NULL )
    2224                 :          1 :         return;
    2225                 :            : 
    2226         [ +  - ]:          1 :     GDALDriver *poDriver = new GDALDriver();
    2227                 :            : 
    2228                 :          1 :     poDriver->SetDescription( "WMTS" );
    2229                 :          1 :     poDriver->SetMetadataItem( GDAL_DCAP_RASTER, "YES" );
    2230                 :          1 :     poDriver->SetMetadataItem( GDAL_DMD_LONGNAME, "OGC Web Map Tile Service" );
    2231                 :          1 :     poDriver->SetMetadataItem( GDAL_DMD_HELPTOPIC, "frmt_wmts.html" );
    2232                 :            : 
    2233                 :          1 :     poDriver->SetMetadataItem( GDAL_DMD_CONNECTION_PREFIX, "WMTS:" );
    2234                 :            : 
    2235                 :          1 :     poDriver->SetMetadataItem( GDAL_DCAP_VIRTUALIO, "YES" );
    2236                 :            : 
    2237                 :            :     poDriver->SetMetadataItem( GDAL_DMD_OPENOPTIONLIST,
    2238                 :            : "<OpenOptionList>"
    2239                 :            : "  <Option name='URL' type='string' description='URL that points to GetCapabilities response' required='YES'/>"
    2240                 :            : "  <Option name='LAYER' type='string' description='Layer identifier'/>"
    2241                 :            : "  <Option name='TILEMATRIXSET' alias='TMS' type='string' description='Tile matrix set identifier'/>"
    2242                 :            : "  <Option name='TILEMATRIX' type='string' description='Tile matrix identifier of maximum zoom level. Exclusive with ZOOM_LEVEL.'/>"
    2243                 :            : "  <Option name='ZOOM_LEVEL' alias='ZOOMLEVEL' type='int' description='Maximum zoom level. Exclusive with TILEMATRIX.'/>"
    2244                 :            : "  <Option name='STYLE' type='string' description='Style identifier'/>"
    2245                 :            : "  <Option name='EXTENDBEYONDDATELINE' type='boolean' description='Whether to enable extend-beyond-dateline behaviour' default='NO'/>"
    2246                 :            : "  <Option name='EXTENT_METHOD' type='string-select' description='How the raster extent is computed' default='AUTO'>"
    2247                 :            : "       <Value>AUTO</Value>"
    2248                 :            : "       <Value>LAYER_BBOX</Value>"
    2249                 :            : "       <Value>TILE_MATRIX_SET</Value>"
    2250                 :            : "       <Value>MOST_PRECISE_TILE_MATRIX</Value>"
    2251                 :            : "  </Option>"
    2252                 :          1 : "</OpenOptionList>");
    2253                 :            : 
    2254                 :          1 :     poDriver->pfnOpen = WMTSDataset::Open;
    2255                 :          1 :     poDriver->pfnIdentify = WMTSDataset::Identify;
    2256                 :          1 :     poDriver->pfnCreateCopy = WMTSDataset::CreateCopy;
    2257                 :            : 
    2258                 :          1 :     GetGDALDriverManager()->RegisterDriver( poDriver );
    2259                 :            : }

Generated by: LCOV version 1.9