上周開始到這幾天，要讀取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