1、makelut
makelut函数可以提供给用户构造一个查找表。
lut = makelut(f, n)
其中f是函数句柄,n可以是2或3。对于2*2邻域总共有16种不同,对3*3邻域有512种变化。
以下我将举例说明:
例1、lut = makelut('sum(x(:)) == 4', 2);
跟踪makelut代码,知
- if (n == 2)
- lut = zeros(16,1);
- for k = 1:16
- a = reshape(fliplr(dec2bin(k-1,4) == '1'), 2, 2);
- lut(k) = feval_r(fun, a, params{:});
- end
对于2*2邻域,总共有16种模版,对每一种模版求出他的模版矩阵a,再对矩阵a用fun(在此为'sum(x(:)) == 4')计算值,将计算得到的值作为lut的元素。
当k=4时,a=[1 0;1 0]。sun(a(:))=2,所以此时的lut(4)=0;
当k=16, a=[1 1;1 1]。 sun(a(:))=4,所以此时的lut(4)=1;
这样我们就能构造出一个查找表。
2、applylut
函数applylut使用查找表来处理二值图像.
跟踪知:
- B = applylutc(A,lut);
而applylutc是在applylutc.mexw32实现。
applylutc.c,函数applylutc的实现,代码参考http://code.google.com/p/mirone/source/browse/trunk/mex/applylutc.c?spec=svn1874&r=1874
函数入口为mexFunction,主要完成2*2或3*3模版像素值的计算,
对3*3模版个像素的权重为
1 8 64
2 16 128
4 32 256
再根据像素值找到对应的查找表,第1个元素对应编号为0的模板,第512个元素对应编号为511的模板。第i个元素值为1,表示可以去掉对应模板的中间的像素,将其像素值从1改为0。