matlab正弦光栅条纹及调整

目录

频率均匀正弦光栅条纹

频率不均匀正弦光栅

幅度不均匀正弦光栅

---

频率均匀正弦光栅条纹

I=zeros(512,512);
for i=1:512
  for j=1:512
    I(i,j)=128+127*cos(j/32*2*pi+pi/4);
  end  
end
I1=mat2gray(I);
figure(1),imshow(I1)

mat2gray，实现矩阵归一化，使得矩阵中的每个元素都在[0,1]之间，其中，0表示黑色，1表示白色。结果如下图所示，使用mesh函数，可以清楚的发现是正弦变化的光栅曲线。其中，j/32*2*pi为函数的频率，pi/4为相位。 

频率不均匀正弦光栅

改变正弦函数的频率，能够得到一个频率不均匀的正弦光栅。这里I1尺寸太大，使用imresize()函数将原图像调整为600×600的矩阵，imresize函数有两种调用方式：

1. I1=imresize(I,scale)，返回图像为I1，I1为原图像I的scale倍，其中，原图像I可以为灰度图像，也可以为彩色图像。
2. I1=imresize(I,[num_row,num_col])，返回图像I1为num_row×num_col的矩阵。

for i=1:1000
  for j=1:1000
    f=300/log(1000-j);
    I(i,j)=127*(1-cos(f/16*2*pi-pi/4)); 
  end  
end
I1=mat2gray(I);
I2=imresize(I1,[600,600]);
subplot(121),imshow(I2),title('I2');
figure,mesh(I2),xlabel('x'),ylabel('y'),zlabel('I2');

 幅度不均匀正弦光栅

对I(i,j)整体进行处理，可以实现对正弦光栅幅值的调整，此处使用Gabor窗口函数对原函数进行处理。

s=200;
for i=1:512
  for j=1:512
    I(i,j)=127*(1-cos(j/32*2*pi+pi/4))*exp(-(i.^2+j.^2)/2/s^2);
  end  
end
I1=mat2gray(I);
figure(1),imshow(I1);
mesh(I1),xlabel('x'),ylabel('y'),zlabel('I2');