Instrument Response是什麼?
我們使用儀器測量物理量,再送至記錄器儲存,中間所經過的轉換會改變原始訊號的單位、以及形狀等等,為了將收到的訊號還原成真實的物理量,需要對原始資料進行儀器響應(instrument response)的修正,每一種儀器與記錄器都有自己的修正值,若修正不對就有可能獲得錯誤的觀察。就像我們用破爛耳機聽音樂,然後說音樂很難聽一樣,品質不好的資料也會影響研究結果,所以必須重視儀器響應的修正。
Poles and Zeros
這個故事要從我們的聲音說起,大家都知道平常聽自己的聲音跟別人聽你的聲音是不同的,這是因為我們聽見的其實是經過頭殼傳過來的聲音,跟經過空氣傳給別人的會不太一樣,這是因為聲音的性質在穿過頭的時候稍稍被改變了;當地震動經過地震儀變成訊號的時候也會有這樣的現象,簡單的說就是訊號的頻譜跟相位稍稍的被儀器改變了。描述儀器改變訊號的方程式叫做Transfer Function,它大概的形狀長成這樣:
其中使分子的值等於零的Zi叫做zero,使分母的值等於零的Pj叫做pole,不同儀器的Transfer Function不同,所以有獨有的poles and zeros。當我們在去除儀器的響應時,除了單位要變,也需要將儀器的扭曲修正,就需要使用到poles and zeros。
修正儀器響應可在SAC當中完成,但是首先必須擁有dataless的seed檔。
再來可以使用rdseed
$ rdseed -R -f [dataless seed] 取出RESP檔
$ rdseed -p -f [dataless seed] 取出PZ檔
