maltab中对查找表函数makelut和applylut的理解

查找表（LUT）是预先 计算出每个可能邻域形状的像素值，然后把这些值储存到一个向量中。

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代码，知

[plain]  view plain copy print ?

1. if (n == 2)
2. lut = zeros(16,1);
3. for k = 1:16
4. a = reshape(fliplr(dec2bin(k-1,4) == '1'), 2, 2);
5. lut(k) = feval_r(fun, a, params{:});
6. 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使用查找表来处理二值图像.

跟踪知：

[plain]  view plain copy print ?

1. 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。