sharpetronics/QGIS
1
0
Fork 0
mirror of https://github.com/qgis/QGIS.git synced 2025-02-28 00:17:30 -05:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #3056 from rouault/gdal_make_autowarpedvrt_for_rpc

Browse Source 
[GDAL provider] Use GDALAutoCreateWarpedVRT() when the dataset has only RPC
This commit is contained in:
Nyall Dawson 2016-05-11 05:08:55 +10:00
commit 21c6486a26
1 changed files with 12 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/providers/gdal/qgsgdalprovider.cpp
View File

@ -2503,7 +2503,8 @@ void QgsGdalProvider::initBaseDataset()
|| mGeoTransform[2] != 0.0
|| mGeoTransform[4] != 0.0
|| mGeoTransform[5] > 0.0 ) )
|| GDALGetGCPCount( mGdalBaseDataset ) > 0 )
|| GDALGetGCPCount( mGdalBaseDataset ) > 0
|| GDALGetMetadata( mGdalBaseDataset, "RPC" ) )
{
QgsLogger::warning( "Creating Warped VRT." );
@ -2580,7 +2581,16 @@ void QgsGdalProvider::initBaseDataset()
if ( !crsFromWkt( GDALGetProjectionRef( mGdalDataset ) ) &&
!crsFromWkt( GDALGetGCPProjection( mGdalDataset ) ) )
{
QgsDebugMsg( "No valid CRS identified" );
if ( mGdalBaseDataset != mGdalDataset &&
GDALGetMetadata( mGdalBaseDataset, "RPC" ) )
{
// Warped VRT of RPC is in EPSG:4326
mCrs.createFromOgcWmsCrs( "EPSG:4326" );
}
else
{
QgsDebugMsg( "No valid CRS identified" );
}
}
//set up the coordinat transform - in the case of raster this is mainly used to convert