- 积分
- 3251
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-2-18
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 棉花猪仔 于 2021-10-29 11:13 编辑
尝试用kriging插值报错
- fnu = 'D:/Project/smfkx/webapps/fkx/data/picture/test/data/925_u.dat'
- ncol = numasciicol(fnu)
- nrow = numasciirow(fnu)
- a = asciiread(fnu,shape=(nrow,ncol))
- lon = a[:,0]
- lat = a[:,1]
- u = a[:,2]
- x = arange(109.1, 117.7,0.09)
- y = arange(20, 26, 0.09)
- prg = griddata((lon2, lat2), u, xi=(x, y), method='kriging')[0]
- #Plot
- figure(figsize=[700, 400], newfig=False)
- axesm()
- geoshow('cn_province', edgecolor='lightgray')
- layer = contourfm(x, y, prg)
- colorbar(layer)
- xlim(109.1, 117.7)
- ylim(20, 26)
复制代码 看了下包是有的
环境 Java 1.8 meteoinfo 3.2
|
|