生成多频外差的光栅图像【Matlab】

leaf_csdn

已于 2022-10-13 13:48:12 修改

阅读量537

点赞数 2

分类专栏：结构光三维重建文章标签： matlab 开发语言 c++

于 2022-10-13 11:39:42 首次发布

本文链接：https://blog.csdn.net/weixin_43956164/article/details/127297355

版权

结构光三维重建专栏收录该内容

16 篇文章 47 订阅

订阅专栏

背景介绍

在matlab中生成多频外差的光栅图像，其中3种频率的选择参考如下文献：
Liu S , Feng W , Zhang Q , et al. Three-dimensional shape measurement of small object based on tri-frequency heterodyne method[C]// 2015 International Conference on Optical Instruments and Technology: Optoelectronic Measurement Technology and Systems. International Society for Optics and Photonics, 2015.
$\lambda_1 = \frac{1}{s-\sqrt{s}},\lambda_2=\frac{1}{s},\lambda_3=\frac{1}{s+(\sqrt{s}+1)}$

代码段

% 生成三频四步条纹光栅图案，创建一个generatePic的文件夹存放生成图像
clc,clear;
width = 1280;
height = 720;

% 垂直条纹
Iv = zeros(height,width,12);
Ia = 150;
Ib = 100;
s = 36;
f = [s+sqrt(s)+1,s,s-sqrt(s)];
f = f/width;
for j=1:width
    for k=1:3
        Iv(:,j,(k-1)*4+1) = round(Ia+Ib*cos(2*pi*f(k)*(j-0.5)+0*pi/2));
        Iv(:,j,(k-1)*4+2) = round(Ia+Ib*cos(2*pi*f(k)*(j-0.5)+1*pi/2));
        Iv(:,j,(k-1)*4+3) = round(Ia+Ib*cos(2*pi*f(k)*(j-0.5)+2*pi/2));
        Iv(:,j,(k-1)*4+4) = round(Ia+Ib*cos(2*pi*f(k)*(j-0.5)+3*pi/2));
    end
end

for n=1:12
    strNum = num2str(n);
    filename = ['generatePic\',strNum,'.bmp'];
    imwrite(uint8(Iv(:,:,n)),filename);
end

% 水平条纹
Ih = zeros(height,width,12);
s = 25;
f = [s+sqrt(s)+1,s,s-sqrt(s)];
f = f/height;
for i=1:height
    for k=1:3
        Ih(i,:,(k-1)*4+1) = round(Ia+Ib*cos(2*pi*f(k)*(i-0.5)+0*pi/2));
        Ih(i,:,(k-1)*4+2) = round(Ia+Ib*cos(2*pi*f(k)*(i-0.5)+1*pi/2));
        Ih(i,:,(k-1)*4+3) = round(Ia+Ib*cos(2*pi*f(k)*(i-0.5)+2*pi/2));
        Ih(i,:,(k-1)*4+4) = round(Ia+Ib*cos(2*pi*f(k)*(i-0.5)+3*pi/2));
    end
end
for n=1:12
    strNum = num2str(n+12);
    filename = ['generatePic\',strNum,'.bmp'];
    imwrite(uint8(Ih(:,:,n)),filename);
end