Example:
MyMatrix = rand(19,42) < 0.5; %some data for example purposes
D = bwdist(MyMatrix, 'cityblock');
3) Look at the bottom of the above link. It says
For cityblock, chessboard, and quasi-Euclidean distance transforms, bwdist uses the two-pass, sequential scanning algorithm described in
[2] Rosenfeld, Azriel and John Pfaltz, "Sequential operations in digital picture processing," Journal of the Association for Computing Machinery, Vol. 13, No. 4, 1966, pp. 471-494.
