function [Dout] = gazdag(Din,v,nx,nz,nt,dx,dz,dt,fmax,con)
%GAZDAG: Constant velocity migration using Gazdag method.
%
% [Dout] = gazdag(Din,v,nx,nz,nt,dx,dz,dt,fmax,con)
%
% Note that the i/o Matrices Din and Dout
% change according to the flag con
%
% IN con = 0 Migration flag
% Din(nt,nx): the poststack data
% conj = 1 Modeling flag
% Din(nz,nx): image of the subsurface
% dx,dx: image grid size in meters
% dt: sampling interval in secs
% fmax: maximum freq. in the data (Hz)
%
%
% OUT con = 1 Modeling flag
% Dout(nt,nx): data
% conj = 0 Migration flag
% Dout(nz,nx): migrated image
%
% This code needs the input variables S ad ARG that
% are generated by propagator.m. Note that propagator.m
% computes the downward and the upward operator choose
% one according to variable conj.
%
%
% Author(s): M.D.Sacchi and H. Kuehl (sacchi@phys.ualberta.ca)
% Copyright 1988-2003 SeismicLab
% Revision: 1.2 Date: Dec/2002
%
% Signal Analysis and Imaging Group (SAIG)
% Department of Physics, UofA
%
[S,ARG] = propagator(v,nx,nz,nt,dx,dz,dt,con);
if1=2;
if2=floor(fmax*(nt*dt))+1;
n1 = if1; n2 = if2;
n3 = nt-n2+2; n4=nt-n1+2;
Band = zeros(1,nt);
Band(1,n1:n2) = 1;
Band(1,n3:n4) = 1;
tc = (nt-1)*dt;
% Migration
sumv = ones(1,nt);
sumv = Band;
if con==0
FFTD=fft2(Din);
IMAGE = zeros(nz,nx);
for iz=1 : nz
FFTD = FFTD.*S;
aperture = ones(nt,nx);
t0 = (iz-1)*dz/v;
index=find(ARG <= (t0/tc)^2);
aperture(index) = 0.0;
IMAGE(iz,:) = sumv*(aperture.*FFTD/nt);
end
Dout=real(ifft(IMAGE,[],2));
end
% Modeling
spray = ones(nt,1);
spray = Band';
if con==1
IMAGE=fft(Din,[],2);
FFTD = zeros(nt,nx);
for iz=nz:-1:1
aperture = ones(nt,nx);
t0 = (iz-1)*dz/v;
index=find(ARG <= (t0/tc)^2+.1);
aperture(index) = 0.0;
FFTD = FFTD+aperture.*(spray*IMAGE(iz,:));
FFTD = FFTD.*S;
end
Dout=real(ifft2(FFTD));
end
return;
没有合适的资源?快使用搜索试试~ 我知道了~
相移法偏移MATLAB代码Gazdag相移法,经测试可完整使用.zip
共5个文件
m:3个
html:1个
asv:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 5 下载量 197 浏览量
2023-03-01
15:25:46
上传
评论
收藏 5KB ZIP 举报
温馨提示
相移法偏移MATLAB代码Gazdag相移法,经测试可完整使用.zip
资源推荐
资源详情
资源评论
收起资源包目录
相移法偏移MATLAB代码Gazdag相移法,经测试可完整使用.zip (5个子文件)
相移法偏移MATLAB代码Gazdag相移法,经测试可完整使用
MIGRATION
impulse.m 859B
gazdag.m 2KB
index.html 309B
propagator.asv 1KB
propagator.m 1KB
A
共 5 条
- 1
资源评论
- weixin_513221702023-09-08这个资源内容超赞,对我来说很有价值,很实用,感谢大佬分享~
- 灰灰是条狗2023-07-24总算找到了想要的资源,搞定遇到的大问题,赞赞赞!
- LIDAN_1232023-12-21总算找到了想要的资源,搞定遇到的大问题,赞赞赞!
- he199207212023-06-07资源中能够借鉴的内容很多,值得学习的地方也很多,大家一起进步!
- cwj_upc2023-12-13资源不错,内容挺好的,有一定的使用价值,值得借鉴,感谢分享。
GZM888888
- 粉丝: 121
- 资源: 2902
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功