fly/terragear
1
0
Fork 0
Code Activity

gdalchop: interpolate nodata pixels before processing.

Browse source 
This yields much nicer scenery tiles as terrafit does not have to work
with these missing data fields. Tested at LOWI.
This commit is contained in:
Christian Schmitt 2012-12-20 13:08:51 +01:00
parent 20751dd50b
commit 6d50ec366c
1 changed files with 5 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
src/Prep/GDALChop/gdalchop.cxx
View file

@ -247,6 +247,10 @@ void ImageInfo::GetDataChunk(int *buffer,
xformData.col_step = colstep;
xformData.row_step = rowstep;
// Interpolate nodata pixels in the source band
GDALRasterBand *mask = dataset->GetRasterBand(srcband)->GetMaskBand();
GDALFillNodata(dataset->GetRasterBand(srcband), mask, 100, 0, 0, NULL, NULL, NULL);
// TODO: check if this image can actually cover part of the chunk
/* establish the full source to target transformation */
@ -255,7 +259,6 @@ void ImageInfo::GetDataChunk(int *buffer,
int srcBandNumbers[] = { srcband };
int dstBandNumbers[] = { 1 };
//double dstNodata = (double)nodata;
double srcNodataReal;
double srcNodataImag = 0.0;
int srcHasNodataValue;
@ -476,7 +479,7 @@ int main(int argc, const char **argv)
west = std::min(west, iwest );
}
SG_LOG(SG_GENERAL, SG_INFO, "Complete bounds n=" << north << " s=" << south << " e=" << east << " w=" << west);
SG_LOG(SG_GENERAL, SG_INFO, "Bounds of all datasets: n=" << north << " s=" << south << " e=" << east << " w=" << west);
/*
* Step 2: If no tiles were specified, go through all tiles contained in