Example 2: Constrained Coarsening
...gridding, we will make sure that the object is represented also in the coarse model. To this end, we will partition the grid into a uniform 5x5 coarse grid and in addition preserve the edges of the object. f = flipud(double(sum(imread('data/facies1.png'),3)))'; G = computeGeometry(cartGrid(size(f)) ); ps = ...