sharpetronics/QGIS
1
0
Fork 0
mirror of https://github.com/qgis/QGIS.git synced 2025-10-04 00:04:03 -04:00
Code Issues Packages Projects Releases Wiki Activity

Tests for issue #63171

Browse Source
This commit is contained in:
Alessandro Pasotti 2025-09-19 09:06:25 +02:00
parent 9074cfbcfc
commit 02946bfa56
1 changed files with 28 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
tests/src/core/testqgsgdalutils.cpp
View File

@ -44,6 +44,7 @@ class TestQgsGdalUtils : public QObject
void testPathIsCheapToOpen(); void testPathIsCheapToOpen();
void testVrtMatchesLayerType(); void testVrtMatchesLayerType();
void testMultilayerExtensions(); void testMultilayerExtensions();
void testResampleSingleBandRasterNoData();
private: private:
double identify( GDALDatasetH dataset, int band, int px, int py ); double identify( GDALDatasetH dataset, int band, int px, int py );
@ -322,5 +323,32 @@ double TestQgsGdalUtils::identify( GDALDatasetH dataset, int band, int px, int p
return value; return value;
} }
void TestQgsGdalUtils::testResampleSingleBandRasterNoData()
{
// Create 2x2 integer raster with values 1, 65535, 65535, 65535 and nodata=65535
QTemporaryDir tempDir;
const QString tempPath = tempDir.path();
const QString inputFilename = tempPath + "/input.tif";
gdal::dataset_unique_ptr ds1 = QgsGdalUtils::createSingleBandTiffDataset( inputFilename, GDT_UInt16, QgsRectangle( 1, 1, 3, 3 ), 2, 2, QgsCoordinateReferenceSystem( "EPSG:4326" ) );
QVERIFY( ds1 );
GDALRasterBandH band = GDALGetRasterBand( ds1.get(), 1 );
GDALSetRasterNoDataValue( band, 65535 );
GUInt16 buffer[4] = { 1, 65535, 65535, 65535 };
const CPLErr err = GDALRasterIO( band, GF_Write, 0, 0, 2, 2, buffer, 2, 2, GDT_UInt16, 0, 0 );
QCOMPARE( err, CE_None );
// Create destination GDAL DS
std::unique_ptr<QgsRasterBlock> dstBlock = std::make_unique<QgsRasterBlock>( Qgis::DataType::UInt16, 4, 4 );
dstBlock->setNoDataValue( 65535.0 );
gdal::dataset_unique_ptr gdalDsOutput = QgsGdalUtils::blockToSingleBandMemoryDataset( QgsRectangle( 1, 1, 3, 3 ), dstBlock.get() );
QVERIFY( QgsGdalUtils::resampleSingleBandRaster( ds1.get(), gdalDsOutput.get(), GRA_NearestNeighbour, nullptr ) );
QCOMPARE( identify( gdalDsOutput.get(), 1, 0, 0 ), 1.0 );
QCOMPARE( identify( gdalDsOutput.get(), 1, 2, 2 ), 65535.0 );
}
QGSTEST_MAIN( TestQgsGdalUtils ) QGSTEST_MAIN( TestQgsGdalUtils )
#include "testqgsgdalutils.moc" #include "testqgsgdalutils.moc"