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

Test for transformCoordinates()

Browse Source
This commit is contained in:
Martin Dobias 2020-05-17 20:48:17 +02:00 committed by Nyall Dawson
parent 68bbf46ebb
commit 6c5cb2ee77
1 changed files with 30 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
tests/src/core/testqgsgdalprovider.cpp
View File

@ -58,6 +58,7 @@ class TestQgsGdalProvider : public QObject
void bandNameWithDescription(); // test band name for when description available (#16047)
void interactionBetweenRasterChangeAndCache(); // test that updading a raster invalidates the GDAL dataset cache (#20104)
void scale0(); //test when data has scale 0 (#20493)
void transformCoordinates();
private:
QString mTestDataDir;
@ -369,5 +370,34 @@ void TestQgsGdalProvider::scale0()
delete provider;
}
void TestQgsGdalProvider::transformCoordinates()
{
// Test implementation of QgsRasterDataProvider::transformCoordinates()
QString raster = QStringLiteral( TEST_DATA_DIR ) + "/float1-16.tif";
QgsDataProvider *provider = QgsProviderRegistry::instance()->createProvider( QStringLiteral( "gdal" ), raster, QgsDataProvider::ProviderOptions() );
QgsRasterDataProvider *rp = dynamic_cast< QgsRasterDataProvider * >( provider );
QVERIFY( rp );
QVERIFY( rp->isValid() );
// forward transform - image coordinates to georeferenced coordinates
QgsPoint pt1Layer = rp->transformCoordinates( QgsPoint( 0, 0 ), QgsRasterDataProvider::TransformImageToLayer ); // bottom-left corner
QgsPoint pt2Layer = rp->transformCoordinates( QgsPoint( 4, 0 ), QgsRasterDataProvider::TransformImageToLayer );
QgsPoint pt3Layer = rp->transformCoordinates( QgsPoint( 4, 4 ), QgsRasterDataProvider::TransformImageToLayer ); // top-right corner
QCOMPARE( pt1Layer, QgsPoint( 106.0, -7.0, 0 ) );
QCOMPARE( pt2Layer, QgsPoint( 106.8, -7.0, 0 ) );
QCOMPARE( pt3Layer, QgsPoint( 106.8, -6.2, 0 ) );
// inverse transform - georeferenced coordinates to image coordinates
QgsPoint pt1Image = rp->transformCoordinates( QgsPoint( 106.0, -7.0 ), QgsRasterDataProvider::TransformLayerToImage ); // bottom-left corner
QgsPoint pt2Image = rp->transformCoordinates( QgsPoint( 106.8, -7.0 ), QgsRasterDataProvider::TransformLayerToImage );
QgsPoint pt3Image = rp->transformCoordinates( QgsPoint( 106.8, -6.2 ), QgsRasterDataProvider::TransformLayerToImage ); // top-right corner
QCOMPARE( pt1Image, QgsPoint( 0, 0, 0 ) );
QCOMPARE( pt2Image, QgsPoint( 4, 0, 0 ) );
QCOMPARE( pt3Image, QgsPoint( 4, 4, 0 ) );
}
QGSTEST_MAIN( TestQgsGdalProvider )
#include "testqgsgdalprovider.moc"