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
     
  • 相关阅读:
    springboot 无法访问静态资源
    webrtc源码阅读理解一
    c++ include的顺序重要吗?
    简说yuv
    i420 转 nv21
    ffmpeg mp4 转 yuv、 y4m转yuv、mp4转y4m && ffplay 播放y4m、yuv
    mysql 查询json字符串条件
    Webpack4.X中sourcemap的配置 详解
    webpack如何打包生成的map文件不生效?
    vue中使用setInterval,移入暂停,移出继续
  • 原文地址:https://www.cnblogs.com/gisalameda/p/8039155.html
Copyright © 2011-2022 走看看