int QgsImageWarper::warpFile( const QString& input,
const QString& output,
const QgsGeorefTransform &georefTransform,
ResamplingMethod resampling,
bool useZeroAsTrans,
const QString& compression,
const QString &projection,
double destResX, double destResY )
{
if ( !georefTransform.parametersInitialized() )
return false;
CPLErr eErr;
GDALDatasetH hSrcDS, hDstDS;
GDALWarpOptions *psWarpOptions;
if ( !openSrcDSAndGetWarpOpt( input, resampling, georefTransform.GDALTransformer(), hSrcDS, psWarpOptions ) )
{
// TODO: be verbose about failures
return false;
}
double adfGeoTransform[6];
int destPixels, destLines;
eErr = GDALSuggestedWarpOutput( hSrcDS, georefTransform.GDALTransformer(),
georefTransform.GDALTransformerArgs(),
adfGeoTransform, &destPixels, &destLines );
if ( eErr != CE_None )
{
GDALClose( hSrcDS );
GDALDestroyWarpOptions( psWarpOptions );
return false;
}
// If specified, override the suggested resolution with user values
if ( destResX != 0.0 || destResY != 0.0 )
{
// If only one scale has been specified, fill in the other from the GDAL suggestion
if ( destResX == 0.0 )
destResX = adfGeoTransform[1];
if ( destResY == 0.0 )
destResY = adfGeoTransform[5];
// Make sure user-specified coordinate system has canonical orientation
if ( destResX < 0.0 )
destResX = -destResX;
if ( destResY > 0.0 )
destResY = -destResY;
// Assert that the north-up convention is fullfiled by GDALSuggestedWarpOutput (should always be the case)
assert( adfGeoTransform[0] > 0.0 );
assert( adfGeoTransform[5] < 0.0 );
// Find suggested output image extent (in georeferenced units)
double minX = adfGeoTransform[0];
double maxX = adfGeoTransform[0] + adfGeoTransform[1] * destPixels;
double maxY = adfGeoTransform[3];
double minY = adfGeoTransform[3] + adfGeoTransform[5] * destLines;
// Update line and pixel count to match extent at user-specified resolution
destPixels = ( int )((( maxX - minX ) / destResX ) + 0.5 );
destLines = ( int )((( minY - maxY ) / destResY ) + 0.5 );
adfGeoTransform[0] = minX;
adfGeoTransform[3] = maxY;
adfGeoTransform[1] = destResX;
adfGeoTransform[5] = destResY;
}
if ( !createDestinationDataset( output, hSrcDS, hDstDS, destPixels, destLines,
adfGeoTransform, useZeroAsTrans, compression,
projection ) )
{
GDALClose( hSrcDS );
GDALDestroyWarpOptions( psWarpOptions );
return false;
}
// Create a QT progress dialog
QProgressDialog *progressDialog = new QProgressDialog( mParent );
progressDialog->setWindowTitle( tr( "Progress indication" ) );
progressDialog->setRange( 0, 100 );
progressDialog->setAutoClose( true );
progressDialog->setModal( true );
progressDialog->setMinimumDuration( 0 );
// Set GDAL callbacks for the progress dialog
psWarpOptions->pProgressArg = createWarpProgressArg( progressDialog );
psWarpOptions->pfnProgress = updateWarpProgress;
psWarpOptions->hSrcDS = hSrcDS;
psWarpOptions->hDstDS = hDstDS;
// Create a transformer which transforms from source to destination pixels (and vice versa)
psWarpOptions->pfnTransformer = GeoToPixelTransform;
psWarpOptions->pTransformerArg = addGeoToPixelTransform( georefTransform.GDALTransformer(),
georefTransform.GDALTransformerArgs(),
adfGeoTransform );
// Initialize and execute the warp operation.
GDALWarpOperation oOperation;
oOperation.Initialize( psWarpOptions );
progressDialog->show();
progressDialog->raise();
progressDialog->activateWindow();
eErr = oOperation.ChunkAndWarpImage( 0, 0, destPixels, destLines );
// eErr = oOperation.ChunkAndWarpMulti(0, 0, destPixels, destLines);
destroyGeoToPixelTransform( psWarpOptions->pTransformerArg );
GDALDestroyWarpOptions( psWarpOptions );
delete progressDialog;
GDALClose( hSrcDS );
GDALClose( hDstDS );
return mWarpCanceled ? -1 : eErr == CE_None ? 1 : 0;
}
