fkfilter.m

这是matlab在地球物理数据处理方面的源码

function [newtraces] = fkfilter(traces,mask)

% function fkfilter - Applies an fkfilter operation
% to a set of traces using a previously developed and 
% saved mask.  Note that the number of data points
% are doubled in the Fourier domain due to padding.



[m,n] = size(traces); 
tt = ones(m,1); temp = hanning(41); tt(1:21) = temp(1:21); 
tt(m-20:m) = temp(21:41);  
[x,y] = meshgrid(hanning(n),tt);  tracemask = x.*y;
traces = tracemask.*traces;  % smooth out some rough edges for ffting later

newtraces = zeros(2*m,2*n);  
newtraces(m/2+1:3*m/2,n/2+1:3*n/2) = traces;
% imagesc(newtraces); 	

z = fft2(newtraces);  
z = fftshift(z); 

z = mask.*z;  
z = fftshift(z); newtraces = real(ifft2(z));  
newtraces = newtraces(m/2+1:3*m/2,n/2+1:3*n/2);