undef("poly3_yao")
; Detrend
; x: array
function poly3_yao(yao:numeric)
local poly3, sizes, Nyear, Nlat, Nlon, x1, x2, x3, x0, x, iLat, iLon, b
begin
sizes = dimsizes(yao)
Nyear = sizes(0)
Nlat = sizes(1)
Nlon = sizes(2)
x1 = fspan(0, Nyear-1, Nyear)
x2 = x1^2.0
x3 = x1^3.0
x0 = x1
x0(:) = 1.0
x = (/x0, x1, x2, x3/)
poly3 = yao
do iLat = 0,Nlat-1
do iLon = 0,Nlon-1
b = reg_multlin(yao(:,iLat,iLon), x, True)
poly3(:,iLat,iLon) = b(0)*x0(:)+b(1)*x1(:)+b(2)*x2(:)+b(3)*x3(:)
end do
end do
poly3 = where(ismissing(yao), poly3@_FillValue, poly3)
return(poly3)
end
三次插值dtrend
最后編輯于 ：
?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明：文章內容（如有圖片或視頻亦包括在內）由作者上傳并發布，文章內容僅代表作者本人觀點，簡書系信息發布平臺，僅提供信息存儲服務。
平臺聲明：文章內容（如有圖片或視頻亦包括在內）由作者上傳并發布，文章內容僅代表作者本人觀點，簡書系信息發布平臺，僅提供信息存儲服務。
- 文/潘曉璐 我一進店門，熙熙樓的掌柜王于貴愁眉苦臉地迎上來，“玉大人，你說我怎么就攤上這事?！?“怎么了？”我有些...
- 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著，像睡著了一般。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發上，一...
- 文/蒼蘭香墨 我猛地睜開眼，長吁一口氣：“原來是場噩夢啊……” “哼！你這毒婦竟也來了？” 一聲冷哼從身側響起，我...
推薦閱讀更多精彩內容
- 在世界上各地，各種各樣的電腦運行著各自不同的操作系統為大家服務，這些電腦在表達同一種信息的時候所使用的方法是千差萬...
- TCP/IP協議中，TCP協議提供可靠的連接服務，采用三次握手建立一個連接，如下圖所示。主機A為客戶機，主機B為服...
- 答:建立連接的過程是利用客戶服務器模式，假設主機A為客戶端，主機B為服務器端。 （1）TCP的三次握手過程：主機A...
- 去年，一部紀錄片火爆全球-楊子江大鱷，記錄馬云和他的阿里巴巴一路風雨到享譽世界的鴻篇巨制！ 馬云在2014的年會上...
