依据左超论文建立查找表

生成条纹级次可能项

-floor(width1) <= K1 <=floor (width1)

-floor(width2) <= K2 <= floor(width2)

按照条纹宽度赋当前k值

按照公式

生成查找表LUT = [ k1 k2 k2p1-k1p2 ]

去除重复值获得查找表LUT。

fun_LUT.m

%生成查找表

function [ lut ] = fun_LUT( width1,width2 )

%!!!注意:width1 > width2

%

show = 1;

Lcm = lcm(width1,width2);

number1 = Lcm/width1;

number2 = Lcm/width2;

for i = -floor(number1/2):floor(number1/2)

n = i - (-floor(number1/2))+1;

start = width1*(n-1)+1;

en_d = width1*n;

A(1,start:en_d) = i;

end

for i = -floor(number2/2):floor(number2/2)

n = i - (-floor(number2/2))+1;

start = width2*(n-1)+1;

en_d = width2*(n);

B(1,start:en_d) = i;

end

if rem(width1,2) == 0

[nouse,length_A] = size(A);

AA = zeros(1,2*length_A);

for n = 1:length_A

AA(2*n-1:2*n) = A(1,n);

end

[nouse,length_B] = size(B);

BB = zeros(1,2*length_B);

for n = 1:length_B

BB(2*n-1:2*n) = B(1,n);

end

dif = length_B-length_A;%(2*length_B-2*length_A)/2;

BB = BB(dif+1:end-dif);

lut = [AA;BB;BB(1,:)*number1-AA(1,:)*number2];

if show == 1

figure;

juxing(1:2*length_A,lut(3,:));

end

%去除重复值

lut = unique(lut','rows');

lut = sortrows(lut,3);

else

lut = [A;B;B(1,:)*number1-A(1,:)*number2];

if show == 1

juxing(1:Lcm,lut(3,:));

end

%去除重复值

lut = unique(lut','rows');

lut = sortrows(lut,3);

end

end