I wanted to try buffering a mask raster, without converting to a shapefile and then back, so I used the SciPy

Multi-dimensional Image Module's binary_dilation function. I then used the following function to generate a circular binary structure array of the desired radius:

import numpy as np

def circular_structure(radius):
size = radius*2+1
i,j = np.mgrid[0:size, 0:size]
i -= (size/2)
j -= (size/2)
return np.sqrt(i**2+j**2) <= radius

## No comments:

## Post a Comment