地震wiggle图显示

地震剖面图通常在一屏上显示多道数据，每道数据用一条wiggle曲线描述，wiggle曲线是由许多小线段组成，线段顶点横坐标、纵坐标分别与地震波振幅、时间深度一一对应，把这些顶点两两相连，便可实现wiggle曲线的绘制

代码链接
下面是网上找的代码，对于line的警告可以删除erasmode和background这两个属性
function wigb (a,scal,x,z,amx)

%WIGB: Plot seismic data using wiggles.

% WIGB(a,scal,x,z,amx)

% IN a: seismic data

% scale: multiple data by scale

% x: x-axis;

% z: vertical axis (time or depth)

% If only ‘a’ is enter, ‘scal,x,z,amn,amx’ are decided automatically;

% otherwise, ‘scal’ is a scalar; ‘x, z’ are vectors for annotation in

% offset and time, amx are the amplitude range.

% Author(s): Xingong Li (Exxon-Mobil)

% Revision: 1.2 Date: Dec/2002

% Signal Analysis and Imaging Group (SAIG)

% Department of Physics, UofA

if nargin == 0, nx=10;nz=10; a = rand(nz,nx)-0.5; end;

[nz,nx]=size(a);

trmx= max(abs(a));

if (nargin <= 4); amx=mean(trmx); end;

if (nargin <= 2); x=[1:nx]; z=[1:nz]; end;

if (nargin <= 1); scal =1; end;

if nx <= 1 ; disp(’ ERR: PlotWig: nx has to be more than 1’);return;end;

% take the average as dx

dx1 = abs(x(2:nx)-x(1:nx-1));

dx = median(dx1);

dz=z(2)-z(1);

xmx=max(max(a)); xmn=min(min(a));

if scal == 0; scal=1; end;

a = a * dx /amx;

a = a * scal;

fprintf(’ PlotWig: data range [%f, %f], plotted max %f \n’,xmn,xmx,amx);

% set display range

x1=min(x)-2.0dx; x2=max(x)+2.0dx;

z1=min(z)-dz; z2=max(z)+dz;

set(gca,‘NextPlot’,‘add’,‘Box’,‘on’, …

‘XLim’, [x1 x2], …

‘YDir’,‘reverse’, …

‘YLim’,[z1 z2]);

fillcolor = [0 0 0];

linecolor = [0 0 0];

linewidth = 0.1;

z=z’; % input as row vector

zstart=z(1);

zend =z(nz);

for i=1:nx,

if trmx(i) ~= 0; % skip the zero traces

tr=a(:,i); % — one scale for all section

s = sign(tr) ;

i1= find( s(1:nz-1) ~= s(2:nz) ); % zero crossing points

npos = length(i1);

%12/7/97

zadd = i1 + tr(i1) ./ (tr(i1) - tr(i1+1)); %locations with 0 amplitudes

aadd = zeros(size(zadd));

[zpos,vpos] = find(tr >0);

[zz,iz] = sort([zpos; zadd]); % indices of zero point plus positives

aa = [tr(zpos); aadd];

aa = aa(iz);

% be careful at the ends

if tr(1)>0, a0=0; z0=1.00;

else, a0=0; z0=zadd(1);

end;

if tr(nz)>0, a1=0; z1=nz;

else, a1=0; z1=max(zadd);

end;

zz = [z0; zz; z1; z0];

aa = [a0; aa; a1; a0];

zzz = zstart + zz*dz -dz;

patch( aa+x(i) , zzz, fillcolor);
%对于此处的警告可以删除erasmode和background
line( ‘Color’,[1 1 1],‘EraseMode’,‘background’, …

‘Xdata’, x(i)+[0 0], ‘Ydata’,[zstart zend]); % remove zero line

%‘LineWidth’,linewidth, …

%12/7/97 ‘Xdata’, x(i)+[0 0], ‘Ydata’,[z0 z1]*dz); % remove zero line

line( ‘Color’,linecolor,‘EraseMode’,‘background’, …

‘LineWidth’,linewidth, …

‘Xdata’, tr+x(i), ‘Ydata’,z); % negatives line

else % zeros trace

line( ‘Color’,linecolor,‘EraseMode’,‘background’, …

‘LineWidth’,linewidth, …

‘Xdata’, [x(i) x(i)], ‘Ydata’,[zstart zend]);

end;

地震wiggle图显示相关推荐

一维地震子波合成记录c语言,地震子波波形显示及一维地震合成记录
<地震子波波形显示及一维地震合成记录>由会员分享,可在线阅读,更多相关<地震子波波形显示及一维地震合成记录(22页珍藏版)>请在人人文库网上搜索. 1.地震勘探原理实验一地震子 ...
R语言使用aov函数进行双因素方差分析（Two-way factorial ANOVA）、使用HH包中的interaction2wt函数为任何阶的双因素方差分析可视化主效应和交互作用图、箱图显示主效应
R语言使用aov函数进行双因素方差分析(Two-way factorial ANOVA).使用HH包中的interaction2wt函数为任何阶的双因素方差分析可视化主效应和交互作用图(Main ef ...
%matplotlib inline是jupyter notebook里的命令, 意思是将那些用matplotlib绘制的图显示在页面里而不是弹出一个窗口
%matplotlib inline是jupyter notebook里的命令, 意思是将那些用matplotlib绘制的图显示在页面里而不是弹出一个窗口终端输入jupyter notebook, ...
《Splunk智能运维实战》——3.11　制作折线图显示项目浏览量和购买量随时间的变化...
本节书摘来自华章计算机<Splunk智能运维实战>一书中的第3章,第3.11节,作者 [美]乔史·戴昆(Josh Diakun),保罗R.约翰逊(Paul R. Johnson),德莱克· ...
MATLAB可视化实战系列（二十四）-三维可视化如何利用圆锥图显示向量场？
前言在物理学中,我们学习了矢量,知道了磁力线等矢量图.如何画出优美的矢量场图呢?科学计算软件matlab可以很容易地绘制出矢量场图.下面跟随博主一起来看下matlab高级可视化内容,如何利用圆锥图显 ...
入门话题1. 在Web中控制图的显示外观？把一张500*800 的图, 显示成180*110 的小图....
问题: 在web中,如何把一张500*800 的图, 显示成180*110 的小图? 方法: a. 利用脚本控制. 在onload 里加个函数. 1<img src ="Imag ...
Linux命令之pstree - 以树状图显示进程间的关系
2019独角兽企业重金招聘Python工程师标准>>> pstree命令以树状图显示进程间的关系(display a tree of processes).ps命令可以显示当前正在运 ...
Linux 命令之 pstree -- 以树状图显示进程/查看进程信息
文章目录介绍常用选项命令示例介绍 Linux 系统中 pstree 命令的英文全称是"process tree",即将所有进程以树状图显示,树状图将会以 pid (如果有指 ...
爆料图显示iPhone 14 Pro及Max机身更厚摄像头凸起也更多
近日,据国外媒体报道,在苹果的春季新品发布会举行之后,分析师和研究机构的重点也就转移到了苹果秋季新品发布会将推出的新品上,尤其是下半年将推出的iPhone 14系列智能手机,有关这一系列智能手机的消息 ...

地震wiggle图显示

地震wiggle图显示相关推荐

最新文章

热门文章