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的年会上...
