[lena_bmp]=imread('0_lena_horiz_lines.bmp'); %讀取原始影像
lena_fft = fft(lena_bmp, [], 2); %原始影像進行一維FFT變換
figure;imshow(lena_fft);title('原始影像頻譜圖'); %顯示原始頻譜圖
Find_Line=abs((lena_fft(:,2)-min(min(lena_fft)))/(max(max(lena_fft))-min(min(lena_fft))));
%定位干擾條紋所在的行
for i = 1 : 512
if Find_Line(i)==0 %如果是干擾條紋所在行數則進行處理
for j = 1 : 512
if i==1 %如果是第一行用后兩行預測出
lena_fft(i,j)=lena_fft(i+1,j)+lena_fft(i+1,j)-lena_fft(i+2,j);
elseif i==512 %如果是最后一行用前兩行預測出
lena_fft(i,j)=lena_fft(i-1,j)+lena_fft(i-1,j)-lena_fft(i-2,j);;
else %其余的用鄰近值預測出
lena_fft(i,j)=(lena_fft(i+1,j)+lena_fft(i-1,j))/2;
end
end
end
end
figure;imshow(lena_fft);title('處理后的lena頻譜圖'); %顯示處理后的頻譜圖
lena_fix=ifft(lena_fft,[],2); %進行一維FFT反變換
figure;imshow(lena_fix,[0,255]);title('處理后的影像'); %顯示處理后的影像
uj5u.com熱心網友回復：
請問有影像嗎uj5u.com熱心網友回復：
可以去landsat tm 條帶嗎uj5u.com熱心網友回復：
這個方法叫什么名字啊？轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/57243.html
標籤：圖象工具使用