25 constexpr double z0xMin = -20037508.3427892;
26 constexpr double z0yMax = 20037508.3427892;
32 const QgsPointXY &z0TopLeftPoint,
double z0Dimension,
int z0MatrixWidth,
int z0MatrixHeight )
35 double z0xMin = z0TopLeftPoint.
x();
36 double z0yMax = z0TopLeftPoint.
y();
37 double z0xMax = z0xMin + z0MatrixWidth * z0Dimension;
38 double z0yMin = z0yMax - z0MatrixHeight * z0Dimension;
41 constexpr double TILE_SIZE = 256.0;
42 constexpr double PIXELS_TO_M = 2.8 / 10000.0;
45 const double scaleDenom0 = ( z0Dimension / TILE_SIZE ) * ( unitToMeters / PIXELS_TO_M );
47 int numTiles =
static_cast<int>( pow( 2,
zoomLevel ) );
52 tm.mMatrixWidth = z0MatrixWidth * numTiles;
53 tm.mMatrixHeight = z0MatrixHeight * numTiles;
54 tm.mTileXSpan = ( z0xMax - z0xMin ) / tm.mMatrixWidth;
55 tm.mTileYSpan = ( z0yMax - z0yMin ) / tm.mMatrixHeight;
56 tm.mExtent =
QgsRectangle( z0xMin, z0yMin, z0xMax, z0yMax );
57 tm.mScaleDenom = scaleDenom0 / pow( 2,
zoomLevel );
64 int numTiles =
static_cast<int>( pow( 2,
zoomLevel ) );
66 int aNumTiles =
static_cast<int>( pow( 2, aZoomLevel ) );
70 tm.mCrs = tileMatrix.
crs();
72 tm.mMatrixWidth = aMatrixWidth * numTiles / aNumTiles;
73 tm.mMatrixHeight = aMatrixHeight * numTiles / aNumTiles;
74 tm.mTileXSpan = aExtent.
width() / tm.mMatrixWidth;
75 tm.mTileYSpan = aExtent.
height() / tm.mMatrixHeight;
77 tm.mScaleDenom = tileMatrix.
scale() * pow( 2, aZoomLevel ) / pow( 2,
zoomLevel );
83 double xMin = mExtent.
xMinimum() + mTileXSpan *
id.column();
84 double xMax = xMin + mTileXSpan;
85 double yMax = mExtent.
yMaximum() - mTileYSpan *
id.row();
86 double yMin = yMax - mTileYSpan;
92 double x = mExtent.
xMinimum() + mTileXSpan *
id.column() + mTileXSpan / 2;
93 double y = mExtent.
yMaximum() - mTileYSpan *
id.row() - mTileYSpan / 2;
103 if ( x0 >= x1 || y0 >= y1 )
106 double tileX1 = ( x0 - mExtent.
xMinimum() ) / mTileXSpan;
107 double tileX2 = ( x1 - mExtent.
xMinimum() ) / mTileXSpan;
108 double tileY1 = ( mExtent.
yMaximum() - y1 ) / mTileYSpan;
109 double tileY2 = ( mExtent.
yMaximum() - y0 ) / mTileYSpan;
111 QgsDebugMsgLevel( QStringLiteral(
"Tile range of edges [%1,%2] - [%3,%4]" ).arg( tileX1 ).arg( tileY1 ).arg( tileX2 ).arg( tileY2 ), 2 );
114 int startColumn = std::clamp(
static_cast<int>( floor( tileX1 ) ), 0, mMatrixWidth - 1 );
115 int endColumn = std::clamp(
static_cast<int>( floor( tileX2 ) ), 0, mMatrixWidth - 1 );
116 int startRow = std::clamp(
static_cast<int>( floor( tileY1 ) ), 0, mMatrixHeight - 1 );
117 int endRow = std::clamp(
static_cast<int>( floor( tileY2 ) ), 0, mMatrixHeight - 1 );
118 return QgsTileRange( startColumn, endColumn, startRow, endRow );
123 double dx = mapPoint.
x() - mExtent.
xMinimum();
124 double dy = mExtent.
yMaximum() - mapPoint.
y();
125 return QPointF( dx / mTileXSpan, dy / mTileYSpan );
181 if ( res == -1 || it->zoomLevel() < res )
192 if ( res == -1 || it->zoomLevel() > res )
193 res = it->zoomLevel();
230 double scaleUnder = 0;
231 double scaleOver = 0;
248 if ( it->scale() > scale && ( zoomUnder == -1 || zoomUnder < it->zoomLevel() ) )
250 zoomUnder = it->zoomLevel();
251 scaleUnder = it->scale();
253 if ( it->scale() < scale && ( zoomOver == -1 || zoomOver > it->zoomLevel() ) )
255 zoomOver = it->zoomLevel();
256 scaleOver = it->scale();
269 if ( scaleOver < scale && scale < scaleUnder )
271 return ( scaleUnder - scale ) / ( scaleUnder - scaleOver ) * ( zoomOver - zoomUnder ) + zoomUnder;
273 scaleUnder = scaleOver;
274 zoomUnder = zoomOver;
275 scaleOver = scaleOver / 2;
280 return ( scaleUnder - scale ) / ( scaleUnder - scaleOver ) * ( zoomOver - zoomUnder ) + zoomUnder;
289 tileZoom =
static_cast<int>( round(
scaleToZoom( scale ) ) );
292 tileZoom =
static_cast<int>( floor(
scaleToZoom( scale ) ) );
305 context.
painter()->device()->logicalDpiX() );
314 return actualMapScale;
321 constexpr double METERS_PER_DEGREE = M_PI / 180.0 * 6378137;
322 constexpr double INCHES_PER_METER = 39.370078;
323 const double mapWidthInches = mapExtent.
width() * METERS_PER_DEGREE * INCHES_PER_METER;
325 double scale = mapWidthInches * mapDpi /
static_cast< double >( mapSize.width() );
336 return actualMapScale;
348 auto readMatrixFromElement = [](
const QDomElement & matrixElement )->
QgsTileMatrix
351 matrix.mZoomLevel = matrixElement.attribute( QStringLiteral(
"zoomLevel" ) ).toInt();
352 matrix.mMatrixWidth = matrixElement.attribute( QStringLiteral(
"matrixWidth" ) ).toInt();
353 matrix.mMatrixHeight = matrixElement.attribute( QStringLiteral(
"matrixHeight" ) ).toInt();
355 matrixElement.attribute( QStringLiteral(
"xMin" ) ).toDouble(),
356 matrixElement.attribute( QStringLiteral(
"yMin" ) ).toDouble(),
357 matrixElement.attribute( QStringLiteral(
"xMax" ) ).toDouble(),
358 matrixElement.attribute( QStringLiteral(
"yMax" ) ).toDouble()
361 matrix.mScaleDenom = matrixElement.attribute( QStringLiteral(
"scale" ) ).toDouble();
362 matrix.mTileXSpan = matrixElement.attribute( QStringLiteral(
"tileXSpan" ) ).toDouble();
363 matrix.mTileYSpan = matrixElement.attribute( QStringLiteral(
"tileYSpan" ) ).toDouble();
364 matrix.mCrs.
readXml( matrixElement );
368 const QDomNodeList children = element.childNodes();
369 for (
int i = 0; i < children.size(); i++ )
371 const QDomElement matrixElement = children.at( i ).toElement();
372 if ( matrixElement.tagName() == QLatin1String(
"rootMatrix" ) )
375 QgsTileMatrix matrix = readMatrixFromElement( matrixElement );
382 const QDomElement rootElement = element.firstChildElement( QStringLiteral(
"rootMatrix" ) );
383 if ( !rootElement.isNull() )
385 mRootMatrix = readMatrixFromElement( rootElement );
393 QDomElement setElement = document.createElement( QStringLiteral(
"matrixSet" ) );
396 auto writeMatrixToElement = [&document](
const QgsTileMatrix & matrix, QDomElement & matrixElement )
398 matrixElement.setAttribute( QStringLiteral(
"zoomLevel" ), matrix.
zoomLevel() );
399 matrixElement.setAttribute( QStringLiteral(
"matrixWidth" ), matrix.
matrixWidth() );
400 matrixElement.setAttribute( QStringLiteral(
"matrixHeight" ), matrix.
matrixHeight() );
408 matrixElement.setAttribute( QStringLiteral(
"tileXSpan" ),
qgsDoubleToString( matrix.mTileXSpan ) );
409 matrixElement.setAttribute( QStringLiteral(
"tileYSpan" ),
qgsDoubleToString( matrix.mTileYSpan ) );
416 QDomElement matrixElement = document.createElement( QStringLiteral(
"matrix" ) );
417 writeMatrixToElement( *it, matrixElement );
418 setElement.appendChild( matrixElement );
421 QDomElement rootElement = document.createElement( QStringLiteral(
"rootMatrix" ) );
423 setElement.appendChild( rootElement );
430 QVector<QgsTileXYZ> tiles;
433 for (
int tileRow = range.
startRow(); tileRow <= range.
endRow(); ++tileRow )
437 QgsTileXYZ tile( tileColumn, tileRow, zoomLevel );
447 if ( !tiles.contains( replacement ) )
448 tiles.append( replacement );
@ Esri
No scale doubling, always rounds down when matching to available tile levels.
@ MapBox
Uses a scale doubling approach to account for hi-DPI tiles, and rounds to the nearest tile level for ...
TileAvailability
Possible availability states for a tile within a tile matrix.
@ UseLowerZoomLevelTile
Tile is not available at the requested zoom level, it should be replaced by a tile from a lower zoom ...
@ NotAvailable
Tile is not available within the matrix, e.g. there is no content for the tile.
@ AvailableNoChildren
Tile is available within the matrix, and is known to have no children (ie no higher zoom level tiles ...
@ Available
Tile is available within the matrix.
This class represents a coordinate reference system (CRS).
bool readXml(const QDomNode &node)
Restores state from the given DOM node.
bool writeXml(QDomNode &node, QDomDocument &doc) const
Stores state to the given Dom node in the given document.
Qgis::DistanceUnit mapUnits
A class to represent a 2D point.
The class is used as a container of context for various read/write operations on other objects.
A rectangle specified with double values.
double xMinimum() const
Returns the x minimum value (left side of rectangle).
double yMinimum() const
Returns the y minimum value (bottom side of rectangle).
double width() const
Returns the width of the rectangle.
double xMaximum() const
Returns the x maximum value (right side of rectangle).
double yMaximum() const
Returns the y maximum value (top side of rectangle).
double height() const
Returns the height of the rectangle.
Contains information about the context of a rendering operation.
QPainter * painter()
Returns the destination QPainter for the render operation.
double rendererScale() const
Returns the renderer map scale.
QSize outputSize() const
Returns the size of the resulting rendered image, in pixels.
QgsRectangle mapExtent() const
Returns the original extent of the map being rendered.
QgsCoordinateTransform coordinateTransform() const
Returns the current coordinate transform for the context.
virtual QDomElement writeXml(QDomDocument &document, const QgsReadWriteContext &context) const
Writes the set to an XML element.
Qgis::ScaleToTileZoomLevelMethod mScaleToTileZoomMethod
void addGoogleCrs84QuadTiles(int minimumZoom=0, int maximumZoom=14)
Adds tile matrices corresponding to the standard web mercator/GoogleCRS84Quad setup.
QgsCoordinateReferenceSystem crs() const
Returns the coordinate reference system associated with the tiles.
std::function< Qgis::TileAvailability(QgsTileXYZ id) > mTileAvailabilityFunction
double scaleForRenderContext(const QgsRenderContext &context) const
Calculates the correct scale to use for the tiles when rendered using the specified render context.
int minimumZoom() const
Returns the minimum zoom level for tiles present in the set.
double scaleToZoom(double scale) const
Calculates a fractional zoom level given a map scale denominator.
QMap< int, QgsTileMatrix > mTileMatrices
QVector< QgsTileXYZ > tilesInRange(QgsTileRange range, int zoomLevel) const
Returns a list of tiles in the given tile range.
int maximumZoom() const
Returns the maximum zoom level for tiles present in the set.
QgsTileMatrix tileMatrix(int zoom) const
Returns the tile matrix corresponding to the specified zoom.
Qgis::TileAvailability tileAvailability(QgsTileXYZ id) const
Returns the availability of the given tile in this matrix.
QgsTileMatrix rootMatrix() const
Returns the root tile matrix (usually corresponding to zoom level 0).
virtual bool readXml(const QDomElement &element, QgsReadWriteContext &context)
Reads the set from an XML element.
void addMatrix(const QgsTileMatrix &matrix)
Adds a matrix to the set.
double calculateTileScaleForMap(double actualMapScale, const QgsCoordinateReferenceSystem &mapCrs, const QgsRectangle &mapExtent, const QSize mapSize, const double mapDpi) const
Calculates the correct scale to use for the tiles when rendered using the specified map properties.
void dropMatricesOutsideZoomRange(int minimumZoom, int maximumZoom)
Deletes any existing matrices which fall outside the zoom range specified by minimumZoom to maximumZo...
bool isEmpty() const
Returns true if the matrix set is empty.
QgsTileMatrix mRootMatrix
void setRootMatrix(const QgsTileMatrix &matrix)
Sets the root tile matrix (usually corresponding to zoom level 0).
std::function< Qgis::TileAvailability(QgsTileXYZ id, QgsTileXYZ &replacement) > mTileReplacementFunction
int scaleToZoomLevel(double scale, bool clamp=true) const
Finds the best fitting (integer) zoom level given a map scale denominator.
Defines a matrix of tiles for a single zoom level: it is defined by its size (width *.
QgsRectangle tileExtent(QgsTileXYZ id) const
Returns extent of the given tile in this matrix.
QPointF mapToTileCoordinates(const QgsPointXY &mapPoint) const
Returns row/column coordinates (floating point number) from the given point in map coordinates.
QgsTileRange tileRangeFromExtent(const QgsRectangle &mExtent) const
Returns tile range that fully covers the given extent.
static QgsTileMatrix fromWebMercator(int zoomLevel)
Returns a tile matrix for the usual web mercator.
QgsRectangle extent() const
Returns extent of the tile matrix.
int matrixWidth() const
Returns number of columns of the tile matrix.
QgsCoordinateReferenceSystem crs() const
Returns the crs of the tile matrix.
double scale() const
Returns scale denominator of the tile matrix.
QgsPointXY tileCenter(QgsTileXYZ id) const
Returns center of the given tile in this matrix.
int matrixHeight() const
Returns number of rows of the tile matrix.
static QgsTileMatrix fromTileMatrix(int zoomLevel, const QgsTileMatrix &tileMatrix)
Returns a tile matrix based on another one.
int zoomLevel() const
Returns the zoom level of the tile matrix.
static QgsTileMatrix fromCustomDef(int zoomLevel, const QgsCoordinateReferenceSystem &crs, const QgsPointXY &z0TopLeftPoint, double z0Dimension, int z0MatrixWidth=1, int z0MatrixHeight=1)
Returns a tile matrix for a specific CRS, top left point, zoom level 0 dimension in CRS units.
Range of tiles in a tile matrix to be rendered.
int endColumn() const
Returns index of the last column in the range.
int endRow() const
Returns index of the last row in the range.
int startRow() const
Returns index of the first row in the range.
int startColumn() const
Returns index of the first column in the range.
Stores coordinates of a tile in a tile matrix set.
static Q_INVOKABLE double fromUnitToUnitFactor(Qgis::DistanceUnit fromUnit, Qgis::DistanceUnit toUnit)
Returns the conversion factor between the specified distance units.
T qgsEnumKeyToValue(const QString &key, const T &defaultValue, bool tryValueAsKey=true, bool *returnOk=nullptr)
Returns the value corresponding to the given key of an enum.
QString qgsDoubleToString(double a, int precision=17)
Returns a string representation of a double.
QString qgsEnumValueToKey(const T &value, bool *returnOk=nullptr)
Returns the value for the given key of an enum.
#define BUILTIN_UNREACHABLE
#define QgsDebugMsgLevel(str, level)
const QgsCoordinateReferenceSystem & crs