Grassfire transform
In image processing, the grassfire transform is the computation of the distance from a pixel to the border of a region. It can be described as "setting fire" to the borders of an image region to yield descriptors such as the region's skeleton or medial axis. Harry Blum introduced the concept in 1967.[1]
Motivation
A region's skeleton can be a useful descriptor, because it describes things such as the symmetry of the region as well as subparts, depressions and protrusions.[2] It also provides a way of relating the interior of a region to the shape of the boundary. In the grassfire transform, the skeleton forms at the points in the region where the "fires" meet. In the literature this is described as the locus of meeting waveforms.[2]
Another advantage of using the outcome of the grassfire transform as a descriptor is that it is invertible. Assuming information about when the medial axis or skeleton is created by meeting waveforms is kept, then the skeleton can be restored by radiating outward.[1]
Example Algorithm
The algorithm below is a simple two pass method for computing the Chessboard distance from the border of a region. Of course there are several other algorithms for performing the grassfire transform.
for each row in image left to right
for each column in image top to bottom
if(pixel is in region){
set pixel to 1 + minimum value of the north and west neighbours
}else{
set pixel to zero
}
}
}
for each row right to left
for each column bottom to top
if(pixel is in region){
set pixel to min(value of the pixel,1 + minimum value of the south and east neighbours)
}else{
set pixel to zero
}
}
}
Below is the result of this transform. It is important to note that the most intense lines make up the skeleton.
 Source image |
 Result image |
Applications
The grassfire transform can be abstracted to suit a variety of computing problems. It has been shown that it can be extended beyond the context of images to arbitrary functions.[3] This includes applications in energy minimization problems such as the those handled by the Viterbi algorithm, max-product belief propagation, resource allocation, and in optimal control methods.[3]
It can also be used to compute the distance between regions by setting the background to be as a region.
See also
References
- 1 2 Blum, Harry. A Transformation for extracting new descriptors of shape, 1967,"http://pageperso.lif.univ-mrs.fr/~edouard.thiel/rech/1967-blum.pdf",6/8/2012
- 1 2 Leymarie, F. and Levine, M.,Simulating the Grassfire Transform Using an Active contour Model,IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol.14, NO 1, January 1992,"http://portal.acm.org/citation.cfm?id=123438"
- 1 2 Felzenszwalb, P. and Huttenlocher, D. , Distance Transforms of Sampled Functions ,"http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.88.1647&rep=rep1&type=pdf"