QGIS API Documentation  3.0.2-Girona (307d082)
qgsruggednessfilter.cpp
Go to the documentation of this file.
1 /***************************************************************************
2  qgsruggednessfilter.cpp - description
3  -----------------------
4  begin : August 7th, 2009
5  copyright : (C) 2009 by Marco Hugentobler
6  email : marco dot hugentobler at karto dot baug dot ethz dot ch
7  ***************************************************************************/
8 
9 /***************************************************************************
10  * *
11  * This program is free software; you can redistribute it and/or modify *
12  * it under the terms of the GNU General Public License as published by *
13  * the Free Software Foundation; either version 2 of the License, or *
14  * (at your option) any later version. *
15  * *
16  ***************************************************************************/
17 
18 #include "qgsruggednessfilter.h"
19 #include <cmath>
20 
21 QgsRuggednessFilter::QgsRuggednessFilter( const QString &inputFile, const QString &outputFile, const QString &outputFormat )
22  : QgsNineCellFilter( inputFile, outputFile, outputFormat )
23 {
24 
25 }
26 
28  : QgsNineCellFilter( QString(), QString(), QString() )
29 {
30 
31 }
32 
33 float QgsRuggednessFilter::processNineCellWindow( float *x11, float *x21, float *x31,
34  float *x12, float *x22, float *x32, float *x13, float *x23, float *x33 )
35 {
36  //the formula would be that easy without nodata values...
37  /*
38  //return *x22; //test: write the raster value of the middle cell
39  float diff1 = *x11 - *x22;
40  float diff2 = *x21 - *x22;
41  float diff3 = *x31 - *x22;
42  float diff4 = *x12 - *x22;
43  float diff5 = *x32 - *x22;
44  float diff6 = *x13 - *x22;
45  float diff7 = *x23 - *x22;
46  float diff8 = *x33 - *x22;
47  return sqrt(diff1 * diff1 + diff2 * diff2 + diff3 * diff3 + diff4 * diff4 + diff5 * diff5 + diff6 * diff6 + diff7 * diff7 + diff8 * diff8);
48  */
49 
50  if ( *x22 == mInputNodataValue )
51  {
52  return mOutputNodataValue;
53  }
54 
55  double sum = 0;
56  if ( *x11 != mInputNodataValue )
57  {
58  sum += ( *x11 - *x22 ) * ( *x11 - *x22 );
59  }
60  if ( *x21 != mInputNodataValue )
61  {
62  sum += ( *x21 - *x22 ) * ( *x21 - *x22 );
63  }
64  if ( *x31 != mInputNodataValue )
65  {
66  sum += ( *x31 - *x22 ) * ( *x31 - *x22 );
67  }
68  if ( *x12 != mInputNodataValue )
69  {
70  sum += ( *x12 - *x22 ) * ( *x12 - *x22 );
71  }
72  if ( *x32 != mInputNodataValue )
73  {
74  sum += ( *x32 - *x22 ) * ( *x32 - *x22 );
75  }
76  if ( *x13 != mInputNodataValue )
77  {
78  sum += ( *x13 - *x22 ) * ( *x13 - *x22 );
79  }
80  if ( *x23 != mInputNodataValue )
81  {
82  sum += ( *x23 - *x22 ) * ( *x23 - *x22 );
83  }
84  if ( *x33 != mInputNodataValue )
85  {
86  sum += ( *x33 - *x22 ) * ( *x33 - *x22 );
87  }
88 
89  return std::sqrt( sum );
90 }
91 
float processNineCellWindow(float *x11, float *x21, float *x31, float *x12, float *x22, float *x32, float *x13, float *x23, float *x33) override
Calculates output value from nine input values.
QgsRuggednessFilter(const QString &inputFile, const QString &outputFile, const QString &outputFormat)
float mOutputNodataValue
The nodata value of the output layer.
float mInputNodataValue
The nodata value of the input layer.
Base class for raster analysis methods that work with a 3x3 cell filter and calculate the value of ea...