Branch data Line data Source code
1 : : /***************************************************************************
2 : : qgsrasterdataprovider.cpp - DataProvider Interface for raster layers
3 : : --------------------------------------
4 : : Date : Mar 11, 2005
5 : : Copyright : (C) 2005 by Brendan Morley
6 : : email : morb at ozemail dot com dot au
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 "qgsproviderregistry.h"
19 : : #include "qgsrasterdataprovider.h"
20 : : #include "qgsrasteridentifyresult.h"
21 : : #include "qgsprovidermetadata.h"
22 : : #include "qgsrasterprojector.h"
23 : : #include "qgslogger.h"
24 : : #include "qgsmessagelog.h"
25 : : #include "qgsapplication.h"
26 : :
27 : : #include <QTime>
28 : : #include <QMap>
29 : : #include <QByteArray>
30 : : #include <QVariant>
31 : :
32 : : #define ERR(message) QgsError(message, "Raster provider")
33 : :
34 : 0 : void QgsRasterDataProvider::setUseSourceNoDataValue( int bandNo, bool use )
35 : : {
36 : 0 : if ( mUseSrcNoDataValue.size() < bandNo )
37 : : {
38 : 0 : for ( int i = mUseSrcNoDataValue.size(); i < bandNo; i++ )
39 : : {
40 : 0 : mUseSrcNoDataValue.append( false );
41 : 0 : }
42 : 0 : }
43 : 0 : mUseSrcNoDataValue[bandNo - 1] = use;
44 : 0 : }
45 : :
46 : 0 : QgsRasterBlock *QgsRasterDataProvider::block( int bandNo, QgsRectangle const &boundingBox, int width, int height, QgsRasterBlockFeedback *feedback )
47 : : {
48 : 0 : QgsDebugMsgLevel( QStringLiteral( "bandNo = %1 width = %2 height = %3" ).arg( bandNo ).arg( width ).arg( height ), 4 );
49 : 0 : QgsDebugMsgLevel( QStringLiteral( "boundingBox = %1" ).arg( boundingBox.toString() ), 4 );
50 : :
51 : 0 : std::unique_ptr< QgsRasterBlock > block = std::make_unique< QgsRasterBlock >( dataType( bandNo ), width, height );
52 : 0 : if ( sourceHasNoDataValue( bandNo ) && useSourceNoDataValue( bandNo ) )
53 : : {
54 : 0 : block->setNoDataValue( sourceNoDataValue( bandNo ) );
55 : 0 : }
56 : :
57 : 0 : if ( block->isEmpty() )
58 : : {
59 : 0 : QgsDebugMsg( QStringLiteral( "Couldn't create raster block" ) );
60 : 0 : block->setError( { tr( "Couldn't create raster block." ), QStringLiteral( "Raster" ) } );
61 : 0 : block->setValid( false );
62 : 0 : return block.release();
63 : : }
64 : :
65 : : // Read necessary extent only
66 : 0 : QgsRectangle tmpExtent = boundingBox;
67 : :
68 : 0 : if ( tmpExtent.isEmpty() )
69 : : {
70 : 0 : QgsDebugMsg( QStringLiteral( "Extent outside provider extent" ) );
71 : 0 : block->setError( { tr( "Extent outside provider extent." ), QStringLiteral( "Raster" ) } );
72 : 0 : block->setValid( false );
73 : 0 : block->setIsNoData();
74 : 0 : return block.release();
75 : : }
76 : :
77 : 0 : double xRes = boundingBox.width() / width;
78 : 0 : double yRes = boundingBox.height() / height;
79 : : double tmpXRes, tmpYRes;
80 : 0 : double providerXRes = 0;
81 : 0 : double providerYRes = 0;
82 : 0 : if ( capabilities() & Size )
83 : : {
84 : 0 : providerXRes = extent().width() / xSize();
85 : 0 : providerYRes = extent().height() / ySize();
86 : 0 : tmpXRes = std::max( providerXRes, xRes );
87 : 0 : tmpYRes = std::max( providerYRes, yRes );
88 : 0 : if ( qgsDoubleNear( tmpXRes, xRes ) ) tmpXRes = xRes;
89 : 0 : if ( qgsDoubleNear( tmpYRes, yRes ) ) tmpYRes = yRes;
90 : 0 : }
91 : : else
92 : : {
93 : 0 : tmpXRes = xRes;
94 : 0 : tmpYRes = yRes;
95 : : }
96 : :
97 : 0 : if ( tmpExtent != boundingBox ||
98 : 0 : tmpXRes > xRes || tmpYRes > yRes )
99 : : {
100 : : // Read smaller extent or lower resolution
101 : :
102 : 0 : if ( !extent().contains( boundingBox ) )
103 : : {
104 : 0 : QRect subRect = QgsRasterBlock::subRect( boundingBox, width, height, extent() );
105 : 0 : block->setIsNoDataExcept( subRect );
106 : 0 : }
107 : :
108 : : // Calculate row/col limits (before tmpExtent is aligned)
109 : 0 : int fromRow = std::round( ( boundingBox.yMaximum() - tmpExtent.yMaximum() ) / yRes );
110 : 0 : int toRow = std::round( ( boundingBox.yMaximum() - tmpExtent.yMinimum() ) / yRes ) - 1;
111 : 0 : int fromCol = std::round( ( tmpExtent.xMinimum() - boundingBox.xMinimum() ) / xRes );
112 : 0 : int toCol = std::round( ( tmpExtent.xMaximum() - boundingBox.xMinimum() ) / xRes ) - 1;
113 : :
114 : 0 : QgsDebugMsgLevel( QStringLiteral( "fromRow = %1 toRow = %2 fromCol = %3 toCol = %4" ).arg( fromRow ).arg( toRow ).arg( fromCol ).arg( toCol ), 4 );
115 : :
116 : 0 : if ( fromRow < 0 || fromRow >= height || toRow < 0 || toRow >= height ||
117 : 0 : fromCol < 0 || fromCol >= width || toCol < 0 || toCol >= width )
118 : : {
119 : : // Should not happen
120 : 0 : QgsDebugMsg( QStringLiteral( "Row or column limits out of range" ) );
121 : 0 : block->setError( { tr( "Row or column limits out of range" ), QStringLiteral( "Raster" ) } );
122 : 0 : block->setValid( false );
123 : 0 : return block.release();
124 : : }
125 : :
126 : : // If lower source resolution is used, the extent must be aligned to original
127 : : // resolution to avoid possible shift due to resampling
128 : 0 : if ( tmpXRes > xRes )
129 : : {
130 : 0 : int col = std::floor( ( tmpExtent.xMinimum() - extent().xMinimum() ) / providerXRes );
131 : 0 : tmpExtent.setXMinimum( extent().xMinimum() + col * providerXRes );
132 : 0 : col = std::ceil( ( tmpExtent.xMaximum() - extent().xMinimum() ) / providerXRes );
133 : 0 : tmpExtent.setXMaximum( extent().xMinimum() + col * providerXRes );
134 : 0 : }
135 : 0 : if ( tmpYRes > yRes )
136 : : {
137 : 0 : int row = std::floor( ( extent().yMaximum() - tmpExtent.yMaximum() ) / providerYRes );
138 : 0 : tmpExtent.setYMaximum( extent().yMaximum() - row * providerYRes );
139 : 0 : row = std::ceil( ( extent().yMaximum() - tmpExtent.yMinimum() ) / providerYRes );
140 : 0 : tmpExtent.setYMinimum( extent().yMaximum() - row * providerYRes );
141 : 0 : }
142 : 0 : int tmpWidth = std::round( tmpExtent.width() / tmpXRes );
143 : 0 : int tmpHeight = std::round( tmpExtent.height() / tmpYRes );
144 : 0 : tmpXRes = tmpExtent.width() / tmpWidth;
145 : 0 : tmpYRes = tmpExtent.height() / tmpHeight;
146 : :
147 : 0 : QgsDebugMsgLevel( QStringLiteral( "Reading smaller block tmpWidth = %1 height = %2" ).arg( tmpWidth ).arg( tmpHeight ), 4 );
148 : 0 : QgsDebugMsgLevel( QStringLiteral( "tmpExtent = %1" ).arg( tmpExtent.toString() ), 4 );
149 : :
150 : 0 : std::unique_ptr< QgsRasterBlock > tmpBlock = std::make_unique< QgsRasterBlock >( dataType( bandNo ), tmpWidth, tmpHeight );
151 : 0 : if ( sourceHasNoDataValue( bandNo ) && useSourceNoDataValue( bandNo ) )
152 : : {
153 : 0 : tmpBlock->setNoDataValue( sourceNoDataValue( bandNo ) );
154 : 0 : }
155 : :
156 : 0 : if ( !readBlock( bandNo, tmpExtent, tmpWidth, tmpHeight, tmpBlock->bits(), feedback ) )
157 : : {
158 : 0 : QgsDebugMsg( QStringLiteral( "Error occurred while reading block" ) );
159 : 0 : block->setError( { tr( "Error occurred while reading block." ), QStringLiteral( "Raster" ) } );
160 : 0 : block->setValid( false );
161 : 0 : block->setIsNoData();
162 : 0 : return block.release();
163 : : }
164 : :
165 : 0 : int pixelSize = dataTypeSize( bandNo );
166 : :
167 : 0 : double xMin = boundingBox.xMinimum();
168 : 0 : double yMax = boundingBox.yMaximum();
169 : 0 : double tmpXMin = tmpExtent.xMinimum();
170 : 0 : double tmpYMax = tmpExtent.yMaximum();
171 : :
172 : 0 : for ( int row = fromRow; row <= toRow; row++ )
173 : : {
174 : 0 : double y = yMax - ( row + 0.5 ) * yRes;
175 : 0 : int tmpRow = std::floor( ( tmpYMax - y ) / tmpYRes );
176 : :
177 : 0 : for ( int col = fromCol; col <= toCol; col++ )
178 : : {
179 : 0 : double x = xMin + ( col + 0.5 ) * xRes;
180 : 0 : int tmpCol = std::floor( ( x - tmpXMin ) / tmpXRes );
181 : :
182 : 0 : if ( tmpRow < 0 || tmpRow >= tmpHeight || tmpCol < 0 || tmpCol >= tmpWidth )
183 : : {
184 : 0 : QgsDebugMsg( QStringLiteral( "Source row or column limits out of range" ) );
185 : 0 : block->setIsNoData(); // so that the problem becomes obvious and fixed
186 : 0 : block->setError( { tr( "Source row or column limits out of range." ), QStringLiteral( "Raster" ) } );
187 : 0 : block->setValid( false );
188 : 0 : return block.release();
189 : : }
190 : :
191 : 0 : qgssize tmpIndex = static_cast< qgssize >( tmpRow ) * static_cast< qgssize >( tmpWidth ) + tmpCol;
192 : 0 : qgssize index = row * static_cast< qgssize >( width ) + col;
193 : :
194 : 0 : char *tmpBits = tmpBlock->bits( tmpIndex );
195 : 0 : char *bits = block->bits( index );
196 : 0 : if ( !tmpBits )
197 : : {
198 : 0 : QgsDebugMsg( QStringLiteral( "Cannot get input block data tmpRow = %1 tmpCol = %2 tmpIndex = %3." ).arg( tmpRow ).arg( tmpCol ).arg( tmpIndex ) );
199 : 0 : continue;
200 : : }
201 : 0 : if ( !bits )
202 : : {
203 : 0 : QgsDebugMsg( QStringLiteral( "Cannot set output block data." ) );
204 : 0 : continue;
205 : : }
206 : 0 : memcpy( bits, tmpBits, pixelSize );
207 : 0 : }
208 : 0 : }
209 : 0 : }
210 : : else
211 : : {
212 : 0 : if ( !readBlock( bandNo, boundingBox, width, height, block->bits(), feedback ) )
213 : : {
214 : 0 : QgsDebugMsg( QStringLiteral( "Error occurred while reading block" ) );
215 : 0 : block->setIsNoData();
216 : 0 : block->setError( { tr( "Error occurred while reading block." ), QStringLiteral( "Raster" ) } );
217 : 0 : block->setValid( false );
218 : 0 : return block.release();
219 : : }
220 : : }
221 : :
222 : : // apply scale and offset
223 : 0 : block->applyScaleOffset( bandScale( bandNo ), bandOffset( bandNo ) );
224 : : // apply user no data values
225 : 0 : block->applyNoDataValues( userNoDataValues( bandNo ) );
226 : 0 : return block.release();
227 : 0 : }
228 : :
229 : 0 : QgsRasterDataProvider::QgsRasterDataProvider()
230 : 0 : : QgsDataProvider( QString(), QgsDataProvider::ProviderOptions(), QgsDataProvider::ReadFlags() )
231 : 0 : , QgsRasterInterface( nullptr )
232 : 0 : , mTemporalCapabilities( std::make_unique< QgsRasterDataProviderTemporalCapabilities >() )
233 : 0 : {
234 : :
235 : 0 : }
236 : :
237 : 0 : QgsRasterDataProvider::QgsRasterDataProvider( const QString &uri, const ProviderOptions &options,
238 : : QgsDataProvider::ReadFlags flags )
239 : 0 : : QgsDataProvider( uri, options, flags )
240 : 0 : , QgsRasterInterface( nullptr )
241 : 0 : , mTemporalCapabilities( std::make_unique< QgsRasterDataProviderTemporalCapabilities >() )
242 : 0 : {
243 : 0 : }
244 : :
245 : 0 : QgsRasterDataProvider::ProviderCapabilities QgsRasterDataProvider::providerCapabilities() const
246 : : {
247 : 0 : return QgsRasterDataProvider::NoProviderCapabilities;
248 : : }
249 : :
250 : 0 : int QgsRasterDataProvider::colorInterpretation( int bandNo ) const
251 : : {
252 : : Q_UNUSED( bandNo )
253 : 0 : return QgsRaster::UndefinedColorInterpretation;
254 : : }
255 : :
256 : : //
257 : : //Random Static convenience function
258 : : //
259 : : /////////////////////////////////////////////////////////
260 : :
261 : : // TODO
262 : : // (WMS) IdentifyFormatFeature is not consistent with QgsRaster::IdentifyFormatValue.
263 : : // IdentifyFormatHtml: better error reporting
264 : 0 : QgsRasterIdentifyResult QgsRasterDataProvider::identify( const QgsPointXY &point, QgsRaster::IdentifyFormat format, const QgsRectangle &boundingBox, int width, int height, int /*dpi*/ )
265 : : {
266 : 0 : QgsDebugMsgLevel( QStringLiteral( "Entered" ), 4 );
267 : 0 : QMap<int, QVariant> results;
268 : :
269 : 0 : if ( format != QgsRaster::IdentifyFormatValue || !( capabilities() & IdentifyValue ) )
270 : : {
271 : 0 : QgsDebugMsg( QStringLiteral( "Format not supported" ) );
272 : 0 : return QgsRasterIdentifyResult( ERR( tr( "Format not supported" ) ) );
273 : : }
274 : :
275 : 0 : if ( !extent().contains( point ) )
276 : : {
277 : : // Outside the raster
278 : 0 : for ( int bandNo = 1; bandNo <= bandCount(); bandNo++ )
279 : : {
280 : 0 : results.insert( bandNo, QVariant() );
281 : 0 : }
282 : 0 : return QgsRasterIdentifyResult( QgsRaster::IdentifyFormatValue, results );
283 : : }
284 : :
285 : 0 : QgsRectangle finalExtent = boundingBox;
286 : 0 : if ( finalExtent.isEmpty() )
287 : 0 : finalExtent = extent();
288 : :
289 : 0 : if ( width == 0 )
290 : : {
291 : 0 : width = capabilities() & Size ? xSize() : 1000;
292 : 0 : }
293 : 0 : if ( height == 0 )
294 : : {
295 : 0 : height = capabilities() & Size ? ySize() : 1000;
296 : 0 : }
297 : :
298 : : // Calculate the row / column where the point falls
299 : 0 : double xres = ( finalExtent.width() ) / width;
300 : 0 : double yres = ( finalExtent.height() ) / height;
301 : :
302 : 0 : int col = static_cast< int >( std::floor( ( point.x() - finalExtent.xMinimum() ) / xres ) );
303 : 0 : int row = static_cast< int >( std::floor( ( finalExtent.yMaximum() - point.y() ) / yres ) );
304 : :
305 : 0 : double xMin = finalExtent.xMinimum() + col * xres;
306 : 0 : double xMax = xMin + xres;
307 : 0 : double yMax = finalExtent.yMaximum() - row * yres;
308 : 0 : double yMin = yMax - yres;
309 : 0 : QgsRectangle pixelExtent( xMin, yMin, xMax, yMax );
310 : :
311 : 0 : for ( int i = 1; i <= bandCount(); i++ )
312 : : {
313 : 0 : std::unique_ptr< QgsRasterBlock > bandBlock( block( i, pixelExtent, 1, 1 ) );
314 : :
315 : 0 : if ( bandBlock )
316 : : {
317 : 0 : double value = bandBlock->value( 0 );
318 : :
319 : 0 : results.insert( i, value );
320 : 0 : }
321 : : else
322 : : {
323 : 0 : results.insert( i, QVariant() );
324 : : }
325 : 0 : }
326 : 0 : return QgsRasterIdentifyResult( QgsRaster::IdentifyFormatValue, results );
327 : 0 : }
328 : :
329 : 0 : double QgsRasterDataProvider::sample( const QgsPointXY &point, int band,
330 : : bool *ok, const QgsRectangle &boundingBox, int width, int height, int dpi )
331 : : {
332 : 0 : if ( ok )
333 : 0 : *ok = false;
334 : :
335 : 0 : const auto res = identify( point, QgsRaster::IdentifyFormatValue, boundingBox, width, height, dpi );
336 : 0 : const QVariant value = res.results().value( band );
337 : :
338 : 0 : if ( !value.isValid() )
339 : 0 : return std::numeric_limits<double>::quiet_NaN();
340 : :
341 : 0 : if ( ok )
342 : 0 : *ok = true;
343 : :
344 : 0 : return value.toDouble( ok );
345 : 0 : }
346 : :
347 : 0 : QString QgsRasterDataProvider::lastErrorFormat()
348 : : {
349 : 0 : return QStringLiteral( "text/plain" );
350 : : }
351 : :
352 : 0 : bool QgsRasterDataProvider::writeBlock( QgsRasterBlock *block, int band, int xOffset, int yOffset )
353 : : {
354 : 0 : if ( !block )
355 : 0 : return false;
356 : 0 : if ( !isEditable() )
357 : : {
358 : 0 : QgsDebugMsg( QStringLiteral( "writeBlock() called on read-only provider." ) );
359 : 0 : return false;
360 : : }
361 : 0 : return write( block->bits(), band, block->width(), block->height(), xOffset, yOffset );
362 : 0 : }
363 : :
364 : : // typedef QList<QPair<QString, QString> > *pyramidResamplingMethods_t();
365 : 0 : QList<QPair<QString, QString> > QgsRasterDataProvider::pyramidResamplingMethods( const QString &providerKey )
366 : : {
367 : 0 : QList<QPair<QString, QString> > methods = QgsProviderRegistry::instance()->pyramidResamplingMethods( providerKey );
368 : 0 : if ( methods.isEmpty() )
369 : : {
370 : 0 : QgsDebugMsg( QStringLiteral( "provider pyramidResamplingMethods returned no methods" ) );
371 : 0 : }
372 : 0 : return methods;
373 : 0 : }
374 : :
375 : 0 : bool QgsRasterDataProvider::hasPyramids()
376 : : {
377 : 0 : const QList<QgsRasterPyramid> pyramidList = buildPyramidList();
378 : 0 : return std::any_of( pyramidList.constBegin(), pyramidList.constEnd(), []( QgsRasterPyramid pyramid ) { return pyramid.getExists(); } );
379 : 0 : }
380 : :
381 : 0 : void QgsRasterDataProvider::setUserNoDataValue( int bandNo, const QgsRasterRangeList &noData )
382 : : {
383 : 0 : if ( bandNo >= mUserNoDataValue.size() )
384 : : {
385 : 0 : for ( int i = mUserNoDataValue.size(); i < bandNo; i++ )
386 : : {
387 : 0 : mUserNoDataValue.append( QgsRasterRangeList() );
388 : 0 : }
389 : 0 : }
390 : 0 : QgsDebugMsgLevel( QStringLiteral( "set %1 band %1 no data ranges" ).arg( noData.size() ), 4 );
391 : :
392 : 0 : if ( mUserNoDataValue[bandNo - 1] != noData )
393 : : {
394 : : // Clear statistics
395 : 0 : int i = 0;
396 : 0 : while ( i < mStatistics.size() )
397 : : {
398 : 0 : if ( mStatistics.value( i ).bandNumber == bandNo )
399 : : {
400 : 0 : mStatistics.removeAt( i );
401 : 0 : mHistograms.removeAt( i );
402 : 0 : }
403 : : else
404 : : {
405 : 0 : i++;
406 : : }
407 : : }
408 : 0 : mUserNoDataValue[bandNo - 1] = noData;
409 : 0 : }
410 : 0 : }
411 : :
412 : 0 : QgsRasterDataProviderTemporalCapabilities *QgsRasterDataProvider::temporalCapabilities()
413 : : {
414 : 0 : return mTemporalCapabilities.get();
415 : : }
416 : :
417 : 0 : const QgsRasterDataProviderTemporalCapabilities *QgsRasterDataProvider::temporalCapabilities() const
418 : : {
419 : 0 : return mTemporalCapabilities.get();
420 : : }
421 : :
422 : 0 : QgsRasterDataProvider *QgsRasterDataProvider::create( const QString &providerKey,
423 : : const QString &uri,
424 : : const QString &format, int nBands,
425 : : Qgis::DataType type,
426 : : int width, int height, double *geoTransform,
427 : : const QgsCoordinateReferenceSystem &crs,
428 : : const QStringList &createOptions )
429 : : {
430 : 0 : QgsRasterDataProvider *ret = QgsProviderRegistry::instance()->createRasterDataProvider(
431 : 0 : providerKey,
432 : 0 : uri, format,
433 : 0 : nBands, type, width,
434 : 0 : height, geoTransform, crs, createOptions );
435 : 0 : if ( !ret )
436 : : {
437 : 0 : QgsDebugMsg( "Cannot resolve 'createRasterDataProviderFunction' function in " + providerKey + " provider" );
438 : 0 : }
439 : :
440 : : // TODO: it would be good to return invalid QgsRasterDataProvider
441 : : // with QgsError set, but QgsRasterDataProvider has pure virtual methods
442 : :
443 : 0 : return ret;
444 : 0 : }
445 : :
446 : 0 : QString QgsRasterDataProvider::identifyFormatName( QgsRaster::IdentifyFormat format )
447 : : {
448 : 0 : switch ( format )
449 : : {
450 : : case QgsRaster::IdentifyFormatValue:
451 : 0 : return QStringLiteral( "Value" );
452 : : case QgsRaster::IdentifyFormatText:
453 : 0 : return QStringLiteral( "Text" );
454 : : case QgsRaster::IdentifyFormatHtml:
455 : 0 : return QStringLiteral( "Html" );
456 : : case QgsRaster::IdentifyFormatFeature:
457 : 0 : return QStringLiteral( "Feature" );
458 : : default:
459 : 0 : return QStringLiteral( "Undefined" );
460 : : }
461 : 0 : }
462 : :
463 : 0 : QString QgsRasterDataProvider::identifyFormatLabel( QgsRaster::IdentifyFormat format )
464 : : {
465 : 0 : switch ( format )
466 : : {
467 : : case QgsRaster::IdentifyFormatValue:
468 : 0 : return tr( "Value" );
469 : : case QgsRaster::IdentifyFormatText:
470 : 0 : return tr( "Text" );
471 : : case QgsRaster::IdentifyFormatHtml:
472 : 0 : return tr( "Html" );
473 : : case QgsRaster::IdentifyFormatFeature:
474 : 0 : return tr( "Feature" );
475 : : default:
476 : 0 : return QStringLiteral( "Undefined" );
477 : : }
478 : 0 : }
479 : :
480 : 0 : QgsRaster::IdentifyFormat QgsRasterDataProvider::identifyFormatFromName( const QString &formatName )
481 : : {
482 : 0 : if ( formatName == QLatin1String( "Value" ) ) return QgsRaster::IdentifyFormatValue;
483 : 0 : if ( formatName == QLatin1String( "Text" ) ) return QgsRaster::IdentifyFormatText;
484 : 0 : if ( formatName == QLatin1String( "Html" ) ) return QgsRaster::IdentifyFormatHtml;
485 : 0 : if ( formatName == QLatin1String( "Feature" ) ) return QgsRaster::IdentifyFormatFeature;
486 : 0 : return QgsRaster::IdentifyFormatUndefined;
487 : 0 : }
488 : :
489 : 0 : QgsRasterInterface::Capability QgsRasterDataProvider::identifyFormatToCapability( QgsRaster::IdentifyFormat format )
490 : : {
491 : 0 : switch ( format )
492 : : {
493 : : case QgsRaster::IdentifyFormatValue:
494 : 0 : return IdentifyValue;
495 : : case QgsRaster::IdentifyFormatText:
496 : 0 : return IdentifyText;
497 : : case QgsRaster::IdentifyFormatHtml:
498 : 0 : return IdentifyHtml;
499 : : case QgsRaster::IdentifyFormatFeature:
500 : 0 : return IdentifyFeature;
501 : : default:
502 : 0 : return NoCapabilities;
503 : : }
504 : 0 : }
505 : :
506 : 0 : QList<double> QgsRasterDataProvider::nativeResolutions() const
507 : : {
508 : 0 : return QList< double >();
509 : : }
510 : :
511 : 0 : bool QgsRasterDataProvider::ignoreExtents() const
512 : : {
513 : 0 : return false;
514 : : }
515 : :
516 : 0 : QgsPoint QgsRasterDataProvider::transformCoordinates( const QgsPoint &point, QgsRasterDataProvider::TransformType type )
517 : : {
518 : 0 : Q_UNUSED( point )
519 : : Q_UNUSED( type )
520 : 0 : return QgsPoint();
521 : : }
522 : :
523 : 0 : bool QgsRasterDataProvider::userNoDataValuesContains( int bandNo, double value ) const
524 : : {
525 : 0 : QgsRasterRangeList rangeList = mUserNoDataValue.value( bandNo - 1 );
526 : 0 : return QgsRasterRange::contains( value, rangeList );
527 : 0 : }
528 : :
529 : 0 : void QgsRasterDataProvider::copyBaseSettings( const QgsRasterDataProvider &other )
530 : : {
531 : 0 : mDpi = other.mDpi;
532 : 0 : mSrcNoDataValue = other.mSrcNoDataValue;
533 : 0 : mSrcHasNoDataValue = other.mSrcHasNoDataValue;
534 : 0 : mUseSrcNoDataValue = other.mUseSrcNoDataValue;
535 : 0 : mUserNoDataValue = other.mUserNoDataValue;
536 : 0 : mExtent = other.mExtent;
537 : 0 : mProviderResamplingEnabled = other.mProviderResamplingEnabled;
538 : 0 : mZoomedInResamplingMethod = other.mZoomedInResamplingMethod;
539 : 0 : mZoomedOutResamplingMethod = other.mZoomedOutResamplingMethod;
540 : 0 : mMaxOversampling = other.mMaxOversampling;
541 : :
542 : : // copy temporal properties
543 : 0 : if ( mTemporalCapabilities && other.mTemporalCapabilities )
544 : : {
545 : 0 : *mTemporalCapabilities = *other.mTemporalCapabilities;
546 : 0 : }
547 : 0 : }
548 : :
549 : 0 : static QgsRasterDataProvider::ResamplingMethod resamplingMethodFromString( const QString &str )
550 : : {
551 : 0 : if ( str == QLatin1String( "bilinear" ) )
552 : : {
553 : 0 : return QgsRasterDataProvider::ResamplingMethod::Bilinear;
554 : : }
555 : 0 : else if ( str == QLatin1String( "cubic" ) )
556 : : {
557 : 0 : return QgsRasterDataProvider::ResamplingMethod::Cubic;
558 : : }
559 : 0 : else if ( str == QLatin1String( "cubicSpline" ) )
560 : : {
561 : 0 : return QgsRasterDataProvider::ResamplingMethod::CubicSpline;
562 : : }
563 : 0 : else if ( str == QLatin1String( "lanczos" ) )
564 : : {
565 : 0 : return QgsRasterDataProvider::ResamplingMethod::Lanczos;
566 : : }
567 : 0 : else if ( str == QLatin1String( "average" ) )
568 : : {
569 : 0 : return QgsRasterDataProvider::ResamplingMethod::Average;
570 : : }
571 : 0 : else if ( str == QLatin1String( "mode" ) )
572 : : {
573 : 0 : return QgsRasterDataProvider::ResamplingMethod::Mode;
574 : : }
575 : 0 : else if ( str == QLatin1String( "gauss" ) )
576 : : {
577 : 0 : return QgsRasterDataProvider::ResamplingMethod::Gauss;
578 : : }
579 : 0 : return QgsRasterDataProvider::ResamplingMethod::Nearest;
580 : 0 : }
581 : :
582 : 0 : void QgsRasterDataProvider::readXml( const QDomElement &filterElem )
583 : : {
584 : 0 : if ( filterElem.isNull() )
585 : : {
586 : 0 : return;
587 : : }
588 : :
589 : 0 : QDomElement resamplingElement = filterElem.firstChildElement( QStringLiteral( "resampling" ) );
590 : 0 : if ( !resamplingElement.isNull() )
591 : : {
592 : 0 : setMaxOversampling( resamplingElement.attribute( QStringLiteral( "maxOversampling" ), QStringLiteral( "2.0" ) ).toDouble() );
593 : 0 : setZoomedInResamplingMethod( resamplingMethodFromString( resamplingElement.attribute( QStringLiteral( "zoomedInResamplingMethod" ) ) ) );
594 : 0 : setZoomedOutResamplingMethod( resamplingMethodFromString( resamplingElement.attribute( QStringLiteral( "zoomedOutResamplingMethod" ) ) ) );
595 : 0 : enableProviderResampling( resamplingElement.attribute( QStringLiteral( "enabled" ) ) == QLatin1String( "true" ) );
596 : 0 : }
597 : 0 : }
598 : :
599 : 0 : static QString resamplingMethodToString( QgsRasterDataProvider::ResamplingMethod method )
600 : : {
601 : 0 : switch ( method )
602 : : {
603 : : case QgsRasterDataProvider::ResamplingMethod::Nearest:
604 : 0 : return QStringLiteral( "nearestNeighbour" );
605 : : case QgsRasterDataProvider::ResamplingMethod::Bilinear:
606 : 0 : return QStringLiteral( "bilinear" );
607 : : case QgsRasterDataProvider::ResamplingMethod::Cubic:
608 : 0 : return QStringLiteral( "cubic" );
609 : : case QgsRasterDataProvider::ResamplingMethod::CubicSpline:
610 : 0 : return QStringLiteral( "cubicSpline" );
611 : : case QgsRasterDataProvider::ResamplingMethod::Lanczos:
612 : 0 : return QStringLiteral( "lanczos" );
613 : : case QgsRasterDataProvider::ResamplingMethod::Average:
614 : 0 : return QStringLiteral( "average" );
615 : : case QgsRasterDataProvider::ResamplingMethod::Mode:
616 : 0 : return QStringLiteral( "mode" );
617 : : case QgsRasterDataProvider::ResamplingMethod::Gauss:
618 : 0 : return QStringLiteral( "gauss" );
619 : : }
620 : : // should not happen
621 : 0 : return QStringLiteral( "nearestNeighbour" );
622 : 0 : }
623 : :
624 : 0 : void QgsRasterDataProvider::writeXml( QDomDocument &doc, QDomElement &parentElem ) const
625 : : {
626 : 0 : QDomElement providerElement = doc.createElement( QStringLiteral( "provider" ) );
627 : 0 : parentElem.appendChild( providerElement );
628 : :
629 : 0 : QDomElement resamplingElement = doc.createElement( QStringLiteral( "resampling" ) );
630 : 0 : providerElement.appendChild( resamplingElement );
631 : :
632 : 0 : resamplingElement.setAttribute( QStringLiteral( "enabled" ),
633 : 0 : mProviderResamplingEnabled ? QStringLiteral( "true" ) : QStringLiteral( "false" ) );
634 : :
635 : 0 : resamplingElement.setAttribute( QStringLiteral( "zoomedInResamplingMethod" ),
636 : 0 : resamplingMethodToString( mZoomedInResamplingMethod ) );
637 : :
638 : 0 : resamplingElement.setAttribute( QStringLiteral( "zoomedOutResamplingMethod" ),
639 : 0 : resamplingMethodToString( mZoomedOutResamplingMethod ) );
640 : :
641 : 0 : resamplingElement.setAttribute( QStringLiteral( "maxOversampling" ),
642 : 0 : QString::number( mMaxOversampling ) );
643 : 0 : }
644 : :
645 : 0 : QString QgsRasterDataProvider::colorInterpretationName( int bandNo ) const
646 : : {
647 : 0 : return colorName( colorInterpretation( bandNo ) );
648 : : }
|
