On 2019-06-21 at 21:10 +02, Robert Nuske <[email protected]> wrote...
> If there is no way to only treat the islands/holes (which i don't know
> about), it will change the outer boundary as well. Fjord like
> structures will vanish. Neighboring areas might even merge if close
> enough. And since one can not control the buffer style of sharp edges
> (like in postgis, geos), everything will become "rounder".
Can you work in raster space or do you have to remain with resolution-agnostic
vectors? I've solved this issue with maintaining coastlines in raster space by
finding and saving the largest (or several largest) clumps (e.g. the ocean
surrounding my domain), filling all holes which may involve filling the sea or
growing the coastline, and then re-applying the saved clumps, thereby removing
any artificially rounded coast.
r.clump input=raster output=clumps
clump_ID_largest=$(r.stats -c clumps sort=desc | head -n1 | cut -d" " -f1)
r.mapcalc "mask_ocean = if(clumps == ${clump_ID_largest}, null(), 1)"
r.mfilter -z input=raster output=raster_overfilled filter=filter.txt
r.mapcalc "raster_filled = raster_overfilled * mask_ocean" --o
-k.
_______________________________________________
grass-user mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/grass-user
