clc
clear
close?all
%%?預處理
I=imread（‘0.bmp‘）;

I=rgb2gray（I）;
I0=I;
%?[m0n0]=size（I0）;

f0=0.17;
%?h1=figure;
%?imshow（I[]）;
%?title（‘左上到右下選取ROI‘）;
%?[pxpy]=ginput（2）;
%?px=round（px）;
%?py=round（py）;

px=[109;903];
py=[181;741];

I=I（py（1）:py（2）px（1）:px（2））;
I=double（I）;
figureimshow（I[]）;
title（‘未校準圖像‘）;
[mn]=size（I）;

%%?提取峰值點
%?digits（64）;
fengzhi=zeros（size（I））;
fengzhishow=zeros（size（I））;
for?aaa=1:size（I1）
????[~locs]=findpeaks（I（aaa3:end-3）‘minpeakheight‘80‘minpeakdistance‘5）;
????locs=locs+2;
????fengzhishow（aaalocs）=1;%記錄峰值點像素位置
????locs0=locs;
????for?ii=1:length（locs）
????????XX=locs（ii）-2:locs（ii）+2;
????????XX0?=?XX（1）-1;
????????XX?=?XX-XX0;
????????XX=XX‘;
????????
????????X=[XX.^2XXones（51）];
????????Y=I（aaalocs（ii）-2:locs（ii）+2）;
????????Y=Y‘;
????????C=（X‘*X）\X‘*Y;
????????C（2）?=?C（2）-2*C（1）*XX0;
????????locs（ii）=-C（2）/（2*C（1））;%以鄰域最小二乘估計的頂點位置作為精確坐標
????????%修補空缺
????????if?aaa==1
????????????continue
????????end
????????xt00?=?fengzhi（aaa-1locs0（ii）-1:locs0（ii）+1）;
????????if?sum（xt00）==0
????????????XX=locs0（ii）-2:locs0（ii）+2;
????????????XX0?=?XX（1）-1;
????????????XX?=?XX-XX0;
????????????XX=XX‘;
????????????X=[XX.^2XXones（51）];
????????????Y=I（aaa-1locs0（ii）-2:locs0（ii）+2）;
????????????Y=Y‘;
????????????Cplus=（X‘*X）\X‘*Y;
????????????Cplus（2）?=?Cplus（2）-2*Cplus（1）*XX0;
????????????locsplus?=?-Cplus（2）/（2*Cplus（1））;
????????????fengzhi（aaa-1round（locsplus））=locsplus;
????????????fengzhishow（aaa-1round（locsplus））=1;
????????end
????end
????fengzhi（aaalocs0）=locs;%將峰值點的精確x坐標寫入峰值位置
end

figureimshow（fengzhishow[]）;%fengzhishow==fengzhi>0;
%%?規定坐標原點
%提取圖像中角點坐標
Image=imread（‘原點.bmp‘）;
Image=double（rgb2gray（Image））;
Image=Image（py（1）:py（2）px（1）:px（2））;
C=corner（Image2）;%取圖中最明顯的兩個角點
h2=figure;
imshow（Image[]）;
hold?on
plot（C（:1）?C（:2）?‘rx‘）;
hold?off
C=mean（C）;
%?刪除原點左側的光柵峰值
fengzhi（:1:floor（C（1）））=0;
fengzhishow（:1:floor（C（1）））=0;
%%?計算梯形校正參數
L=bwlabel（fengzhishow8）;%檢測8鄰域連通域
num=max（L（:））;
for?a=1:num
????l=L==a;
????%刪除長度不足的連通域
????if?sum（l（:））????????fengzhishow（l）=0;
????????fengzhi（l）=0;
????????L（l）=0;
????end
end
figureimshow（fengzhishow[]）;%?不允許有光柵未提取成功
title（‘篩選后的峰值圖像‘）;
[~col]=find（L（1:）>0）;
W=（length（col）-1）./f0;
l1=L（1min（col））;
l2=L（1max（col））;
[row1col1]=find（L==l1）;%光柵最靠左的一列極值點
[row2col2]=find（L==l2）;%光柵最靠右的一列極值點
l1=zeros（1m）;
l2=zeros（1m）;
for?i=1:m
????l1（1i）=fengzhi（row1（i）col1（i））;
????l2（1i）=fengzhi（row2（i）col2（i））;
end
A=[row1ones（m1）];
b=l1‘;
p1=（A‘*A）^-1*A‘*b;
A=[row2ones（m1）];
b=l2‘;
p2=（A‘*A）^-1*A‘*b;

tana1=p1（1）;
tana2=p2（1）;
a1=atan（tana1）;%左直線傾斜角
a2=atan（tana2）;
jiaodu=（a1+a2）./2;
tanaa=tan（（a1-a2）./2）;

%%?生成反向補償光柵投影
m=?574;
n=?574;
x=1:n;
y=1:m;
y=m-y+1;
[xy]=meshgrid（xy）;

S?=?1.12;

xt=x;
xt=-y*tanaa+xt.*（W+2*y*tanaa）./W;
xt=xt*cos（jiaodu）-y*sin（jiaodu）;
%?I0=0.5+0.5*cos（2*pi*f0*xt）;
xy=5;
jg=2*pi/x

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-04-04?16:05??條紋校準\
?????文件?????2359350??2017-05-26?10:10??條紋校準\0.bmp
?????文件?????2359350??2017-05-26?10:19??條紋校準\1.bmp
?????文件?????2156186??2018-04-04?15:53??條紋校準\C01.mat
?????文件????????3345??2018-04-04?15:53??條紋校準\Correction.m
?????文件????????2083??2018-04-04?15:58??條紋校準\jieguofenxi_PMP1.m
?????文件????????2091??2018-04-04?15:58??條紋校準\jieguofenxi_PMP2.m
?????文件???????17396??2016-04-06?14:55??條紋校準\Miguel_2D_unwrapper.cpp
?????文件???????11264??2016-04-15?18:33??條紋校準\Miguel_2D_unwrapper.mexw64
?????文件?????????427??2018-04-04?15:52??條紋校準\Original.m
?????目錄???????????0??2018-04-04?15:53??條紋校準\PMPafter1\
?????文件?????2359350??2017-05-26?10:22??條紋校準\PMPafter1\01.bmp
?????文件?????2359350??2017-05-26?10:22??條紋校準\PMPafter1\02.bmp
?????文件?????2359350??2017-05-26?10:22??條紋校準\PMPafter1\03.bmp
?????文件?????2359350??2017-05-26?10:23??條紋校準\PMPafter1\04.bmp
?????文件?????2359350??2017-05-26?10:23??條紋校準\PMPafter1\05.bmp
?????文件?????2359350??2017-05-26?10:26??條紋校準\PMPafter1\06.bmp
?????文件?????2359350??2017-05-26?10:26??條紋校準\PMPafter1\07.bmp
?????文件?????2359350??2017-05-26?10:26??條紋校準\PMPafter1\08.bmp
?????文件?????2359350??2017-05-26?10:27??條紋校準\PMPafter1\09.bmp
?????文件?????2359350??2017-05-26?10:27??條紋校準\PMPafter1\10.bmp
?????目錄???????????0??2018-04-04?15:53??條紋校準\PMPafter2\
?????文件?????2359350??2017-05-26?10:22??條紋校準\PMPafter2\01.bmp
?????文件?????2359350??2017-05-26?10:22??條紋校準\PMPafter2\02.bmp
?????文件?????2359350??2017-05-26?10:22??條紋校準\PMPafter2\03.bmp
?????文件?????2359350??2017-05-26?10:23??條紋校準\PMPafter2\04.bmp
?????文件?????2359350??2017-05-26?10:23??條紋校準\PMPafter2\05.bmp
?????文件?????2359350??2017-05-26?10:28??條紋校準\PMPafter2\11.bmp
?????文件?????2359350??2017-05-26?10:29??條紋校準\PMPafter2\12.bmp
?????文件?????2359350??2017-05-26?10:29??條紋校準\PMPafter2\13.bmp
?????文件?????2359350??2017-05-26?10:29??條紋校準\PMPafter2\14.bmp
............此處省略26個文件信息