您所在位置：主题：跪求关于震相数据的提取程序
萌萌的臭猴子
[专家分：0] 发布于 2012-09-09 17:29:00
南北带子网震相数据
(台网代码：NX)
日期 时间 纬度 经度 深度 震级 台站 参考地点
(北京时) (°) (°) (km) 数目
2008/12/28 20:44:01.9 38.93 106.45 5 ML2.4 5 宁夏平罗
台站 震中距 清晰度 震相 到时 震级 周期 振幅 振幅
(km) 及初动 (s) 名称 (μm)
YCI 0.10 SME 3.240
YCI 154.4 SG 20:44:48.7 ML2.5 0.10 SMN 3.000
YCI EC PG 20:44:29.8
TXN 0.20 SME 1.530
TXN 181.8 SG 20:44:56.1 ML2.4 0.20 SMN 0.970
TXN ED PG 20:44:34.8
SZS 0.20 SME 13.790
SZS 42.9 SG 20:44:14.9 ML2.7 0.20 SMN 24.080
SZS ED PG 20:44:10.0
ZHW 0.40 SME 1.630
ZHW 183.4 SG 20:44:54.6 ML2.8 0.30 SMN 1.570
共计10条震相数据。
这样的一个样本数据，提取其中2008/12/28 20:44:01.9 38.93 106.45 5 ML2.4 5 宁夏平罗这一行的年，月，日，时，分，秒，纬度，经度，深度，震级和台站数目，为了好识别部分数据之间用0.00隔开，输出为
0812280.002044 01.900.0038.930.00106.450.00 50.00 2.40.00 50.00
从“YCI 0.10 SME 3.240” 开始数据读取其中的太站名，初动，和到时，例如下一行数据
“YCI 154.4 SG 20:44:48.7 ML2.5 0.10 SMN 3.000” 如果有到时数据20:44:48.7，内提取并输出“YCI SG204448.70 ”
向第一行没有到时数据则跳过，继续下一行。整个文件读取之后结果如下
0812280.002044 01.900.0038.930.00106.450.00 50.00 2.40.00 50.00
YCI SG204448.70YCI PG204429.80TXN SG204456.10TXN PG204434.80SZS SG204414.90SZS PG204410.00六组数据一行
ZHW SG204454.600.000000000000000000000000000000000000000000000000000000000000000000000不够六组的用0.00补齐