wangliang1989 / oh-my-cap Goto Github PK
View Code? Open in Web Editor NEWHome Page: http://seiswave.cn/oh-my-cap/
License: GNU General Public License v3.0
Home Page: http://seiswave.cn/oh-my-cap/
License: GNU General Public License v3.0
我在教程中看到:数据必须去除仪器响应得到位移或速度记录,单位为 cm 或 cm/s
我在使用seed文件中提取的PZ文件去除仪器响应后,(得到的结果据说单位是m)将结果乘以100,然后按照流程进行得到的反演结果看起来不太对。很多波形图超出边界,且计算出的震级与地震目录中的震级相差较大。
请问通过PZ文件去除仪器响应后是否有不一样的处理步骤?
计划如下修改:
考虑加强权重设置的功能
数据预处理:
DT CAP里的-H参数,重采样的周期
FREQ 去仪器相应的四个频率
反演:
DEPTH 反演的深度
GREEN -G参数
MODEL 格林函数model
PLOT cap的-p参数,绘图
SHIFT cap -s
P_WINDOW 反演P窗口
S_WINDOW S窗口
WEIGHTFILE 权重文件
P_FILTER CAP -c参数里的滤波
S_FILTER 同上
以下三个我想不出合适名字:
-D
-W: 1
-X: 10
以上问题需要检查
目前 hk 中的四列数据分别为:厚度、vs、波速比、Qs,这不是一个常见的速度模型格式,建议修改为更通用的格式:
thickness vs vp rho Qs Qp
在反演时,是权重文件中的台名+分量来搜索sac文件的,对此要考虑有loc的时候,文件名协调的问题
配置参数化后更新文件的工作之前没有做完,做完后注意修改历史记录!
共三个地方涉及到了 dt:
cap.pl
的 -H
选项需要做的是:
#!/bin/bash
psmeca -P -JX0/10 -R-5/5/-5/5 -Sm8 -G100 -Y9.5 -X-0.7 > 1.ps <<EOF
0 0 0 -0.3 -1 -0.650 0 0 0 17
1 1 0 -0.3 -1 -0.649 0 0 0 17
EOF
ps2raster -E1080 -P -Tf 1.ps
使用以上内容作图,可以看到两个机制解其中一个黑色,然而这两个机制解应该均为白色。我这边作图结果是1.pdf
老师您好,我是用里面的gcap时一直有问题
#报错:can‘t exec 'cap.pl':NO such...
请问是gcc要降到3.4版本的吗
另外,请问咱们是否有qq群这样的交流平台大家可以方便交流
考虑在生成权重文件的脚本中加一段判断:如果 t9 不是 -12345,这个数据权重为0。
foreach (@layer) {
my $layer = trim($_);
my ($dep) = split m/\s+/, $layer;
foreach (@depth) {
if ($_ == $dep) {
print "设置错误:震源在 $_ km 深度的界面上\n";
$err++;
}
}
}
请教一个问题:计算单力源和爆炸型震源时,程序报错。非常感谢!
命令行为:$ fk.pl -Mhk/15/k -N1024/0.1 -S0 -U0 10。(计算爆炸型震源时)。 程序默认的-S2(双力偶型震源)
报错:Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Backtrace for this error:
#0 0x7f8eac02532a
#1 0x7f8eac024503
#2 0x7f8eab8b9fcf
#3 0x559541734f9f
#4 0x559541734922
#5 0x5595417329db
#6 0x55954173380b
#7 0x7f8eab89cb96
#8 0x5595417313e9
#9 0xffffffffffffffff
hk_15/10.grn.0
Unable to open hk_15/10.grn.0
hk_15/10.grn.5
Unable to open hk_15/10.grn.5
marktime.pl
中对是否安装了 TauP 进行检测。若未安装,则直接跳过标记理论到时get_depth.pl 调用 Zhu 的脚本,没有考虑到路径下可能有多个model的结果
请教一下老师,设置不同时窗长度的P和S波,在输出图件中长度为什么几乎一样呢?在参数设置文件中,更改-P的值似乎也只能整体的改变其大小,如何才可以设置最终绘制的图件中各分量波形长度与其窗口长度是成比例的呢?谢谢。
现在脚本是按一律算爆炸源来处理,应该增加是否计算的选项
现在只实现了原depth.pl的做图功能,输出深度反演的功能还未实现
输出含义也未完全弄清
可能会更换模型,修改定义模型的方式,便于修改模型。另外需要确认给 fk 计算的模型文件一定会正确。
add rdseed
老师您好!请问gCAP需要使用的商业软件 Numerical Recipes在哪里可以购买
反演的矩张量结果保存在.out的结果文件中# MomentTensor那一行。
官方脚本是这样传入gmt绘图的:
$plt3 = "| psmeca -JX1/1 -R-1/1/-1/1 -Sm8 -G100 -Y9.5 -X-0.7 -O -K >> $outps" if $tensor[1] eq "MomentTensor";
open(PLT, $plt3);
printf PLT "0 0 0 @tensor[9,4,7,6] %f %f 17\n",-$tensor[8],-$tensor[5];
因为GMT的psmeca中,所用的是RTP/RTF/USE坐标系。所以得到对应关系(MOPAD默认用NED):
9 4 7 6 -8 -5 = mrr, mtt, mff, mrt, mrf, mtf = dd nn ee dn -de -ne
或者说:
4 5 6 7 8 9 = nn ne dn ee de dd
打算重写波形拟合图,拟合误差深度图和新写台站位置示意图
marktime.pl 未考虑不同分量长度不同的情况而造成 cut 数据出错
Ubuntu22.04 编译安装GMT6.1.1,运行示例perl gmt6depth.pl 20080418093658时报错
gmt [ERROR]: Shared GMT module not found: basemap
ERROR: No module named basemap was found. This could mean one of four cases:
gmt [ERROR]: Shared GMT module not found: plot
ERROR: No module named plot was found. This could mean one of four cases:
meca [ERROR]: Internal Failure = GMT_MAP_BAD_LAT_MAX
meca (GMT_psmeca): North is outside -90 to +90 degree range
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.