grow.region function

Get adjacent flanks from regions