[lena_bmp]=imread (' 0 _lena_horiz_lines. BMP); % to read the original image
Lena_fft=FFT (lena_bmp, [], 2); % the original image to one dimensional FFT transform
figure; Imshow (lena_fft); The title (' original image spectrum); % according to the original spectrum
Find_Line=abs ((lena_fft (:, 2) - min (min (lena_fft)))/(Max (Max (lena_fft)) - min (min (lena_fft))));
% positioning interference fringe line
For I=1:512
If Find_Line (I)==0% if is interference fringe number of rows processed
For j=1:512
If I==1% if it is the first line after use two lines to predict
Lena_fft (I, j)=lena_fft lena_fft (I + 1, j) + (I + 1, j) - lena_fft (I + 2, j);
Elseif I==512% if it is the last line in the first two lines to predict
Lena_fft (I, j)=lena_fft (I - 1, j) + lena_fft (I - 1, j) - lena_fft (I - 2, j);;
Else % to predict the rest with adjacent values
Lena_fft (I, j)=(lena_fft (I + 1, j) + lena_fft (I - 1, j))/2;
End
End
End
End

figure; Imshow (lena_fft); The title (' processed Lena spectrum diagram). % showed processed spectrum
Lena_fix=ifft (lena_fft, [], 2); % of the one dimensional FFT transformation
figure; Imshow (lena_fix, [0255]); The processed image is the title (' '); The processed image is % shows

CodePudding user response:

Do you have image

CodePudding user response:

You can go to landsat tm bands

CodePudding user response:

What's the name of this method?