zoukankan      html  css  js  c++  java
  • mseed读取

    上周开始到这几天,要读取miniseed格式的数据,困难很多。

    1 rdmseed.m rdmseedfast.m

    matlab2017运行rdmseed遇到有gap的数据,或者数据不完整的,会报错。暂时不懂怎么处理。sanhua说matlab版本低的可以运行,不会报错。

    rdmseedfast.m读取快速,但是读取的时间和rdmseed读取的不一致,不知道哪个是对的。

    2 没办法,试试mseed2sac。

     mseed2sac QT-5314_4-20170809062553 QT-5314_4-20170809073939

    安装与使用参考http://www.cnblogs.com/seisjun/p/6225189.html

    Building mseed2sac :
     1) $ tar xvf mseed2sac-2.0.tar.gz
     2) $ cd mseed2sac-2.0/
     3) $ make
     ... and the executable 'mseed2sac' is created.
     4) Installing : $ sudo cp mseed2sac /usr/local/bin/
    mseed2sac nc_stations.mseed
    它与rdseed相比较的弊端在于:
    但解压出来的 SAC 文件中只有台站名和台网名,与用 rdSEED 解压出来的 SAC 数据相比, 少了 cmpaz、cmpinc、stla、stlo、stel、stdp、kcmpnm 信息,因而需要手动添加台站信息。
    
    

    http://seisman.info/convert-seed-to-sac.html

    在使用 mseed2sac 时,还可以加上额外的 metafile 和 selectfile 文件,以及地震事件信息。比如:

    mseed2sac -m metafile -I selectfile -E 2006,123,15:27:08.7/-20.33/-174.03/65.5/Tonga
    其中 metafile 提供了台站的基本信息,selectfile 中列出了要提取哪些台站的数据。具体文件格式参考官方文档。
    metafile 可以自己根据格式生成,也可以直接使用 IRIS 提供的 FetchMetadata 脚本获取。


    批量处理
    #!/bin/bash
    clear all
    path="/home/ling/software/mseed2sac-master/20171109"
    cd $path
    filelist=$(ls $dir)
    for filename in $filelist
     do
       mseed2sac $filename
       done

    convert mseed to sac  

    http://blog.163.com/weijia_sun/blog/static/16659249720154763955973/

    #!/bin/sh

    # step 1: convert miniSEED to SAC
    for file in `ls`
    do

    if [ -d "${file}" ]; then
    echo "$file is a directory"
    mkdir ${file}/SAC
    rm ${file}/SAC/*
    cd ${file}/t/
    rm *SAC
    mseed2sac -O *BHE *BHZ *BHN
    cd ../../
    mv ${file}/t/*SAC ${file}/SAC/
    fi
    done
    # remove blank in the filename after the conversion
    ls -R */SAC/* |
    awk '{printf "mv "%s"",$0;gsub(/[[:blank:]]/,"");print " " $0}' > blank.list
    ls -R */SAC/*BH.* | awk '{print "rm " $0}' >> blank.list
    chmod u+x blank.list
    ./blank.list
     
  • 相关阅读:
    CF 461B Appleman and Tree
    POJ 1821 Fence
    NOIP 2012 开车旅行
    CF 494B Obsessive String
    BZOJ2337 XOR和路径
    CF 24D Broken robot
    POJ 1952 BUY LOW, BUY LOWER
    SPOJ NAPTIME Naptime
    POJ 3585
    CF 453B Little Pony and Harmony Chest
  • 原文地址:https://www.cnblogs.com/gisalameda/p/8039155.html
Copyright © 2011-2022 走看看