QGIS API Documentation 3.41.0-Master (3440c17df1d)
Loading...
Searching...
No Matches
qgstolerance.cpp
Go to the documentation of this file.
1/***************************************************************************
2 qgstolerance.cpp - wrapper for tolerance handling
3 ----------------------
4 begin : March 2009
5 copyright : (C) 2009 by Richard Kostecky
6 email : csf.kostej at gmail dot com
7 ***************************************************************************
8 * *
9 * This program is free software; you can redistribute it and/or modify *
10 * it under the terms of the GNU General Public License as published by *
11 * the Free Software Foundation; either version 2 of the License, or *
12 * (at your option) any later version. *
13 * *
14 ***************************************************************************/
15
16#include "qgstolerance.h"
17#include "moc_qgstolerance.cpp"
18#include "qgsmapsettings.h"
22#include "qgspointxy.h"
23
24#include <QPoint>
25#include <cmath>
26
27
28// return ratio [mu/lu] between map units and layer units
29// this is of course only an approximation
30double _ratioMU2LU( const QgsMapSettings &mapSettings, QgsMapLayer *layer )
31{
32 const double distMU = mapSettings.mapUnitsPerPixel();
33 const QgsPointXY ptMapCenterMU = mapSettings.visibleExtent().center();
34 const QgsPointXY ptMapCenterRightMU( ptMapCenterMU.x() + distMU, ptMapCenterMU.y() );
35 const QgsPointXY ptMapCenterLU = mapSettings.mapToLayerCoordinates( layer, ptMapCenterMU );
36 const QgsPointXY ptMapCenterRightLU = mapSettings.mapToLayerCoordinates( layer, ptMapCenterRightMU );
37 const double distLU = std::sqrt( ptMapCenterLU.sqrDist( ptMapCenterRightLU ) );
38 const double ratio = distMU / distLU;
39 return ratio;
40}
41
42double QgsTolerance::toleranceInProjectUnits( double tolerance, QgsMapLayer *layer, const QgsMapSettings &mapSettings, Qgis::MapToolUnit units )
43{
44 // converts to map units
45 if ( units == Qgis::MapToolUnit::Project )
46 return tolerance;
47 else if ( units == Qgis::MapToolUnit::Pixels )
48 return tolerance * mapSettings.mapUnitsPerPixel();
49 else // units == LayerUnits
50 {
51 // [mu] = [lu] * [mu/lu]
52 return tolerance * _ratioMU2LU( mapSettings, layer );
53 }
54}
55
56
57double QgsTolerance::toleranceInMapUnits( double tolerance, QgsMapLayer *layer, const QgsMapSettings &mapSettings, Qgis::MapToolUnit units )
58{
59 // converts to layer units
60 if ( units == Qgis::MapToolUnit::Layer )
61 {
62 return tolerance;
63 }
64 else if ( units == Qgis::MapToolUnit::Pixels )
65 {
66 const double layerUnitsPerPixel = computeMapUnitPerPixel( layer, mapSettings );
67 return tolerance * layerUnitsPerPixel;
68 }
69 else // ProjectUnits
70 {
71 // [lu] = [mu] / [mu/lu]
72 return tolerance / _ratioMU2LU( mapSettings, layer );
73 }
74}
75
84
91
98
99
100double QgsTolerance::computeMapUnitPerPixel( QgsMapLayer *layer, const QgsMapSettings &mapSettings )
101{
102 // the layer is projected. Find out how many pixels are in one map unit - either horizontal and vertical direction
103 // this check might not work correctly in some cases
104 // (on a large area the pixels projected around "0,0" can have different properties from the actual point)
105 const QgsPointXY p1 = toLayerCoordinates( layer, mapSettings, QPoint( 0, 1 ) );
106 const QgsPointXY p2 = toLayerCoordinates( layer, mapSettings, QPoint( 0, 2 ) );
107 const QgsPointXY p3 = toLayerCoordinates( layer, mapSettings, QPoint( 1, 0 ) );
108 const QgsPointXY p4 = toLayerCoordinates( layer, mapSettings, QPoint( 2, 0 ) );
109 const double x = p1.sqrDist( p2 );
110 const double y = p3.sqrDist( p4 );
111 if ( x > y )
112 {
113 return std::sqrt( x );
114 }
115 else
116 {
117 return std::sqrt( y );
118 }
119}
120
121
122QgsPointXY QgsTolerance::toLayerCoordinates( QgsMapLayer *layer, const QgsMapSettings &mapSettings, QPoint point )
123{
124 const QgsPointXY pt = mapSettings.mapToPixel().toMapCoordinates( point );
125 return mapSettings.mapToLayerCoordinates( layer, pt );
126}
MapToolUnit
Type of unit of tolerance value from settings.
Definition qgis.h:4645
@ Layer
Layer unit value.
@ Project
Map (project) units.
@ Pixels
Pixels unit of tolerance.
Base class for all map layer types.
Definition qgsmaplayer.h:76
The QgsMapSettings class contains configuration for rendering of the map.
const QgsMapToPixel & mapToPixel() const
double mapUnitsPerPixel() const
Returns the distance in geographical coordinates that equals to one pixel in the map.
QgsRectangle visibleExtent() const
Returns the actual extent derived from requested extent that takes output image size into account.
QgsPointXY mapToLayerCoordinates(const QgsMapLayer *layer, QgsPointXY point) const
transform point coordinates from output CRS to layer's CRS
QgsPointXY toMapCoordinates(int x, int y) const
Transforms device coordinates to map (world) coordinates.
A class to represent a 2D point.
Definition qgspointxy.h:60
double sqrDist(double x, double y) const
Returns the squared distance between this point a specified x, y coordinate.
Definition qgspointxy.h:186
double y
Definition qgspointxy.h:64
double x
Definition qgspointxy.h:63
QgsPointXY center
T value(const QString &dynamicKeyPart=QString()) const
Returns settings value.
static const QgsSettingsEntryDouble * settingsDigitizingSearchRadiusVertexEdit
Settings entry digitizing search radius vertex edit.
static const QgsSettingsEntryEnumFlag< Qgis::MapToolUnit > * settingsDigitizingDefaultSnappingToleranceUnit
Settings entry digitizing default snapping tolerance unit.
static const QgsSettingsEntryEnumFlag< Qgis::MapToolUnit > * settingsDigitizingSearchRadiusVertexEditUnit
Settings entry digitizing search radius vertex edit unit.
static const QgsSettingsEntryDouble * settingsDigitizingDefaultSnappingTolerance
Settings entry digitizing default snapping tolerance.
static double toleranceInMapUnits(double tolerance, QgsMapLayer *layer, const QgsMapSettings &mapSettings, Qgis::MapToolUnit units=Qgis::MapToolUnit::Layer)
Static function to translate tolerance value into layer units.
static double vertexSearchRadius(const QgsMapSettings &mapSettings)
Static function to get vertex tolerance value.
static double defaultTolerance(QgsMapLayer *layer, const QgsMapSettings &mapSettings)
Static function to get default tolerance value for a layer.
static double toleranceInProjectUnits(double tolerance, QgsMapLayer *layer, const QgsMapSettings &mapSettings, Qgis::MapToolUnit units)
Static function to translate tolerance value into map units.
double _ratioMU2LU(const QgsMapSettings &mapSettings, QgsMapLayer *layer)