Hey Matt,<br><br>I recently encountered this exact same problem and would like to find a gdal based solution. I ended up solving the problem in Grass by creating a grayscale mask of the replacement area and used a convolution filter to feather the edge. I used this in r.mapcalc to weight the addition of the two rasters.<br>
<br>It seems this approach could be duplicated with gdal / python. The main issue would be the creation of the mask to weight the feathering. The rest should be pretty simple math.<br><br>Let me know if you've heard any other solutions, I'd be interested in pursuing this more.<br>
<br>Thanks,<br><br>Jamie Adams<br><br><div class="gmail_quote">On Feb 19, 2008 5:04 PM, Matt Wilkie <<a href="mailto:matt.wilkie@gov.yk.ca">matt.wilkie@gov.yk.ca</a>> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello,<br><br>Does anyone have any advice on how to feather rasters together so there<br>is not an abrubt transition from one to the other?<br><br>The circumstance: I have 16bit integer raster (a DEM) with some areas in<br>
it I want to fix (lakes which should be flat). I've created a second<br>raster with the desired replacement values (fixed_lakes). However if I<br>blindly just substitute fixed_lakes.tif for dem.tif there are some<br>pretty abrubt transitions which are sometimes worse than the original<br>
problem.<br><br>I've done this before using Gimp and 8bit images but gimp can't load<br>16bit images without losing half the bits. (in layer pallet stack<br>fixed.tif over original.tif > transparency to selection for fixed layer<br>
> invert selection > expand X pixels > feather > flatten image; restore<br>spatial reference with gdalcopyproj.py)<br><br>thanks in advance for your thoughts,<br><br>--<br>matt wilkie<br>--------------------------------------------<br>
Geographic Information,<br>Information Management and Technology,<br>Yukon Department of Environment<br>10 Burns Road * Whitehorse, Yukon * Y1A 4Y9<br>867-667-8133 Tel * 867-393-7003 Fax<br><a href="http://environmentyukon.gov.yk.ca/geomatics/" target="_blank">http://environmentyukon.gov.yk.ca/geomatics/</a><br>
--------------------------------------------<br><br>_______________________________________________<br>gdal-dev mailing list<br><a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br><a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div><br>