【实例简介】
【实例截图】
【核心代码】
function [arymig,tmig,xmig]=kirk(aryin,aryvel,dt,dx,params)
% KIRK: simplified Kirchhoff time migration
%
% [arymig,tmig,xmig]=kirk(aryin,aryvel,t,x,params)
%
% KIRK is a simple post stack Kirchhoff time migration routine. This is
% just simple NMO and sum. No time derivative, no cosine factor, and
% nearest neighbor interpolation.
%
% aryin ... matrix of zero offset data. One trace per column.
% aryvel ... velocity information. The are 3 possibilities:
% 1) if a scalar, then a constant velocity migration with
% velocity=aryvel is performed.
% 2) if a vector, then it must be the same length as the number
% of rows in aryin. In this case it is assumed to be an rms
% velocity function (of time) which is applied at all positions
% along the section.
% 3) if a matrix, then it must be the same size as aryin. Here it
% is assumed to give the rms velocity for each sample location.
% t ... if a scalar, this is the time sample rate in SECONDS.
% If a vector, it gives the time coordinates for the rows of aryin.
% x ... if a scalar, this is the spatial sample rate (in units
% consistent with the velocity information. If a vector, then
% it gives the x coordinates of the columns of aryin
% params ... vector of migration parameters
% params(1) ... migration aperture in physical length units
%