baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時(shí)
蟲號: 505962
注冊: 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

[交流] 分享用matlab顯示地震記錄的波形變面積圖已有9人參與

wigb函數(shù)是網(wǎng)上找的，自己添了個(gè)主程序方便使用。

主程序
function PlotRecord_Wigb(filename,nx,nt,dx,dt,zy)
%波形變面積顯示地震記錄
% clear; clf;clc;
% nx=200; %道數(shù)
% nt=2001; %采樣點(diǎn)數(shù)
% dx=10.;    %道間距
% dt=0.001; %時(shí)間采樣間隔
% zy=1.0;    %增益
% filename='r200-2001-1ms.dat';
clf;
fid=fopen(filename,'r');
[A,COUNT]=fread(fid,[nt,nx],'float32');
fclose(fid);

x=0:dx: (nx-1)*dx;
z=0:dt: (nt-1)*dt;
wigb(A,zy,x,z);

end

子函數(shù)
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)
%  Copyright 1998-2003 Xingong
%  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.0*dx; x2=max(x)+2.0*dx;
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);

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;
end;

[ Last edited by baobiao007 on 2012-4-9 at 10:41 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

地球物理資源集錦

小胖的專輯

» 本帖已獲得的紅花（最新10朵）

逄碩

» 猜你喜歡

求助裂變徑跡已經(jīng)有0人回復(fù)
澳門大學(xué)機(jī)器人與自主系統(tǒng)碩士項(xiàng)目已經(jīng)有2人回復(fù)
地球化學(xué)論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有192人回復(fù)
南京大學(xué)能源與資源學(xué)院馬朝陽課題組博士后招聘已經(jīng)有30人回復(fù)
調(diào)劑西南科技大學(xué) 環(huán)資學(xué)院已經(jīng)有4人回復(fù)
西南科技大學(xué) 調(diào)劑綿陽已經(jīng)有10人回復(fù)
散金幣，求基金中標(biāo)��！已經(jīng)有39人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

如何用MATLAB做出地震正演模型已經(jīng)有5人回復(fù)
地震面波英文已經(jīng)有7人回復(fù)
如何用matlab產(chǎn)生markov波形已經(jīng)有3人回復(fù)
【原創(chuàng)】從皇極數(shù)經(jīng)談地震預(yù)報(bào) 已經(jīng)有6人回復(fù)

我同意叔本華的觀點(diǎn)，人們投身藝術(shù)和科學(xué)領(lǐng)域的強(qiáng)烈愿望之一就是逃離痛苦、殘酷和枯燥無味的現(xiàn)實(shí)生活，逃離自己飄忽不定的七情六欲的桎梏。--愛因斯坦

1樓 2012-04-09 10:35:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖置頂 ( 共有2個(gè) )

smithjohn

金蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 1334.9
散金: 372
紅花: 4
帖子: 673
在線: 148.8小時(shí)
蟲號: 1030430
注冊: 2010-05-28

0.8

板塊做地震的不多
樓主試過實(shí)際地震資料么？

贊一下(1人)

回復(fù)此樓

2樓2012-04-09 18:13:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時(shí)
蟲號: 505962
注冊: 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

沒有，實(shí)際地震資料要根據(jù)具體的數(shù)據(jù)格式來弄，勘探地震有自己行業(yè)標(biāo)準(zhǔn)的數(shù)據(jù)格式。天然地震的就不清楚了

贊一下(1人)

回復(fù)此樓

3樓2012-04-09 18:21:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時(shí)
蟲號: 505962
注冊: 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

引用回帖:

9樓: Originally posted by 逄碩 at 2013-05-05 19:25:49
這個(gè)導(dǎo)入的數(shù)據(jù)應(yīng)該是什么參數(shù)啊

普通二進(jìn)制地震記錄數(shù)據(jù)

贊一下(1人)

回復(fù)此樓

10樓2013-05-05 19:29:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cj000211

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 4
帖子: 13
在線: 5.7小時(shí)
蟲號: 804526
注冊: 2009-07-06
專業(yè): 應(yīng)用地球物理學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

多謝了~~~剛下載完，回去試試看

贊一下(1人)

回復(fù)此樓

11樓2013-05-07 11:20:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

lzk55882

謝謝分享

4樓2012-04-11 10:29:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

i星辰

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1811.8
帖子: 403
在線: 230小時(shí)
蟲號: 1142640
注冊: 2010-11-09
性別: GG
專業(yè): 應(yīng)用地球物理學(xué)

犀利啊，拜讀了

回復(fù)此樓

5樓2012-04-11 14:59:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sapyren

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 341.5
散金: 427
紅花: 3
帖子: 110
在線: 31.9小時(shí)
蟲號: 1405389
注冊: 2011-09-17
性別: GG
專業(yè): 大地測量學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

fname附件貼出來啊，這樣看好沒勁！

贊一下

回復(fù)此樓

做純粹的自己，享受責(zé)任與自由�。�！

6樓2012-04-26 20:07:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

應(yīng)助: 0 (幼兒園)
金幣: 5
帖子: 1
在線: 28分鐘
蟲號: 2135987
注冊: 2012-11-19
專業(yè): 應(yīng)用地球物理學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

送紅花一朵

這個(gè)導(dǎo)入的數(shù)據(jù)應(yīng)該是什么參數(shù)啊

贊一下

回復(fù)此樓

9樓2013-05-05 19:25:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 北科281學(xué)碩材料求調(diào)劑 +5	tcxiaoxx 2026-03-20	5/250	2026-03-20 21:35 by laoshidan
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen