第一次在微信公眾號寫的啥鳥文,感覺很low,希望慢慢進步吧,其實頭腦裡想法的,但是一寫字就難了,看來得訓練一把
在學習的過程中,我逐漸發現分享能夠產生更大的價值,也不清楚這種習慣是好是壞,但在分享的過程中我對學習的東西一定是加深了印象。因此我想把學到的東西一點一點得積累下來,希望在運用的時候儘量做到遊刃有餘。
高通濾波和低通濾波在學ENVI的時候接觸過,但那時候沒仔細看,也沒有運用過,直到後來的研究生複試我才被問到這個問題,我很想解決這個問題,因此在網上查找了一些資料和Python實現的方法。至於它的實際運用,大家可以在網上搜索到。
在Python中做高通濾波和低通濾波主要使用的是scipy庫中signal模塊,scipy的官網有所介紹,參數解釋也是比較詳細的。那就說到這裡吧,現在開始我們的任務。1,先構建一個濾波器(這裡使用的是巴特沃斯濾波器),2,再對數據做濾波處理。命令signal.butter(A,B,btype),裡面的參數分別是濾波器階數,臨界頻率和最後的類型,類型有高低通,帶通濾波和帶阻濾波。這裡選擇highpass。濾波器構建好後就開始濾吧,
FFF = signal.filtfilt(m, n, y,axis = 0)
這裡的mn是上一步構建的濾波器,y是一維時間序列數據,axis這裡是指定哪一維數據,直接出結果吧,數據是自己輸入的,還沒有使用xlrd讀取xlsx的數據來分析,先用這個小例子測試一下吧,過後再遷移運用