QGIS API Documentation 3.32.0-Lima (311a8cb8a6)
qgsclassificationstandarddeviation.cpp
Go to the documentation of this file.
1/***************************************************************************
2 qgsclassificationstandarddeviation.h
3 ---------------------
4 begin : September 2019
5 copyright : (C) 2019 by Denis Rouzaud
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 <QObject>
17
20#include "qgsapplication.h"
21
22const QString QgsClassificationStandardDeviation::METHOD_ID = QStringLiteral( "StdDev" );
23
24
26 : QgsClassificationMethod( SymmetricModeAvailable )
27{
28
29}
30
32{
33 return QObject::tr( "Standard Deviation" );
34}
35
37{
38 return METHOD_ID;
39}
40
42{
44 copyBase( c );
45 c->mStdDev = mStdDev;
46 return c;
47}
48
50{
51 return QgsApplication::getThemeIcon( "classification_methods/mClassificationStandardDeviation.svg" );
52}
53
54
55QList<double> QgsClassificationStandardDeviation::calculateBreaks( double &minimum, double &maximum,
56 const QList<double> &values, int nclasses )
57{
58 // C++ implementation of the standard deviation class interval algorithm
59 // as implemented in the 'classInt' package available for the R statistical
60 // prgramming language.
61
62 // Returns breaks based on 'prettyBreaks' of the centred and scaled
63 // values of 'values', and may have a number of classes different from 'classes'.
64
65 // If there are no values to process: bail out
66 if ( values.isEmpty() )
67 return QList<double>();
68
69 double mean = 0.0;
70 mStdDev = 0.0;
71 const int n = values.count();
72
73 for ( int i = 0; i < n; i++ )
74 {
75 mean += values[i];
76 }
77 mean = mean / static_cast< double >( n );
78
79 double sd = 0.0;
80 for ( int i = 0; i < n; i++ )
81 {
82 sd = values[i] - mean;
83 mStdDev += sd * sd;
84 }
85 mStdDev = std::sqrt( mStdDev / n );
86
87 // if not symmetric, the symmetry point is the mean
88 mEffectiveSymmetryPoint = symmetricModeEnabled() ? symmetryPoint() : mean;
89
90 QList<double> breaks = QgsSymbolLayerUtils::prettyBreaks( ( minimum - mEffectiveSymmetryPoint ) / mStdDev, ( maximum - mEffectiveSymmetryPoint ) / mStdDev, nclasses );
91 makeBreaksSymmetric( breaks, 0.0, symmetryAstride() ); //0.0 because breaks where computed on a centered distribution
92
93 for ( int i = 0; i < breaks.count(); i++ )
94 breaks[i] = ( breaks[i] * mStdDev ) + mEffectiveSymmetryPoint;
95
96 return breaks;
97}
98
99QString QgsClassificationStandardDeviation::labelForRange( const double lowerValue, const double upperValue, QgsClassificationMethod::ClassPosition position ) const
100{
101 const QString lowerLabel = valueToLabel( lowerValue );
102 const QString upperLabel = valueToLabel( upperValue );
103
104 switch ( position )
105 {
106 case LowerBound:
107 return QStringLiteral( "< %1" ).arg( upperLabel );
108 case Inner:
109 {
110 QString label( labelFormat() );
111 label.replace( QLatin1String( "%1" ), lowerLabel ).replace( QLatin1String( "%2" ), upperLabel );
112 return label;
113 }
114 case UpperBound:
115 return QStringLiteral( "≥ %1" ).arg( lowerLabel );
116 }
117 return QString();
118}
119
120
121QString QgsClassificationStandardDeviation::valueToLabel( const double value ) const
122{
123 const double normalized = ( value - mEffectiveSymmetryPoint ) / mStdDev;
124 return QObject::tr( " %1 Std Dev" ).arg( QLocale().toString( normalized, 'f', 2 ) );
125}
126
127
128void QgsClassificationStandardDeviation::writeXml( QDomElement &element, const QgsReadWriteContext &context ) const
129{
130 Q_UNUSED( context )
131
132 element.setAttribute( QStringLiteral( "std_dev" ), QString::number( mStdDev, 'f', 16 ) );
133 element.setAttribute( QStringLiteral( "effective_symmetry_point" ), QString::number( mEffectiveSymmetryPoint, 'f', 16 ) );
134}
135
136void QgsClassificationStandardDeviation::readXml( const QDomElement &element, const QgsReadWriteContext &context )
137{
138 Q_UNUSED( context )
139
140 mStdDev = element.attribute( QStringLiteral( "std_dev" ), QStringLiteral( "1.0" ) ).toDouble();
141 mEffectiveSymmetryPoint = element.attribute( QStringLiteral( "effective_symmetry_point" ), QStringLiteral( "0.0" ) ).toDouble();
142}
static QIcon getThemeIcon(const QString &name, const QColor &fillColor=QColor(), const QColor &strokeColor=QColor())
Helper to get a theme icon.
QgsClassificationMethod is an abstract class for implementations of classification methods.
double symmetryPoint() const
Returns the symmetry point for symmetric mode.
bool symmetricModeEnabled() const
Returns if the symmetric mode is enabled.
static void makeBreaksSymmetric(QList< double > &breaks, double symmetryPoint, bool astride)
Remove the breaks that are above the existing opposite sign classes to keep colors symmetrically bala...
ClassPosition
Defines the class position.
@ LowerBound
The class is at the lower bound.
@ UpperBound
The class is at the upper bound.
@ Inner
The class is not at a bound.
bool symmetryAstride() const
Returns if the symmetric mode is astride if true, it will remove the symmetry point break so that the...
QString labelFormat() const
Returns the format of the label for the classes.
void copyBase(QgsClassificationMethod *c) const
Copy the parameters (shall be used in clone implementation)
QgsClassificationCustom is an implementation of QgsClassificationMethod based on standard deviation.
QIcon icon() const override
The icon of the method.
QString name() const override
The readable and translate name of the method.
void writeXml(QDomElement &element, const QgsReadWriteContext &context) const override
Writes extra information about the method.
QString labelForRange(double lowerValue, double upperValue, ClassPosition position) const override
Returns the label for a range.
QString id() const override
The id of the method as saved in the project, must be unique in registry.
QgsClassificationMethod * clone() const override
Returns a clone of the method.
void readXml(const QDomElement &element, const QgsReadWriteContext &context) override
Reads extra information to apply it to the method.
The class is used as a container of context for various read/write operations on other objects.
static QList< double > prettyBreaks(double minimum, double maximum, int classes)
Computes a sequence of about 'classes' equally spaced round values which cover the range of values fr...
As part of the API refactoring and improvements which landed in the Processing API was substantially reworked from the x version This was done in order to allow much of the underlying Processing framework to be ported into c