利用python画等温图

珮瑶瑶

第一次利用python画气象图，不是很完美，请大神批评指正。
绘制图形如下：


脚本见附件：

咖喱嫑洋葱

珮瑶瑶 发表于 2019-8-27 20:29
新手不太理解，请把代码贴出来吧，谢谢

1. lon,lat = np.meshgrid(lon,lat)
   
2. olon = np.linspace(lon_min,lon_max,插值格点数)
   
3. olat = np.linspace(同上)
   
4. olon,olat = np.meshgrid(olon,olat)
   
5. func = Rbf(lon, lat, value,function='linear')
   
6. value_new = func(olon, olat)

复制代码

lon，lat，value的形状要一样

天边的那边

我是来学习的，python真强大

Camb_L

还在学习Python中

咖喱嫑洋葱

请问python怎么能够平滑等值线？

seafairy

入门最难，谢谢分享。

忽如一夜春风来

咖喱嫑洋葱 发表于 2019-8-16 11:09
请问python怎么能够平滑等值线？

Scipy插值
from scipy.interpolate import Rbf

……

olon = np.linspace(lon_min,lon_max,插值格点数)
olat = np.linspace(同上)
olon,olat = np.meshgrid(olon,olat)
func = Rbf(lon, lat, value,function='linear')
value_new = func(olon, olat)


……
然后利用 lon,lat, value_new画图

咖喱嫑洋葱

忽如一夜春风来 发表于 2019-8-20 15:28
Scipy插值
from scipy.interpolate import Rbf

我用的就是1°×1°的经纬度和对应的格点数据，
我在想python上有没有类似grads中九点平滑之类的效果

忽如一夜春风来

咖喱嫑洋葱 发表于 2019-8-21 15:58
我用的就是1°×1°的经纬度和对应的格点数据，
我在想python上有没有类似grads中九点平滑之类的效果

没有找到相关直接使用的函数，这里建议的是插值成更密的格点，达到平滑的目的

咖喱嫑洋葱

忽如一夜春风来 发表于 2019-8-20 15:28
Scipy插值
from scipy.interpolate import Rbf

我用Rbf(lon, lat, value,function='linear')，报错，ValueError: All arrays must be equal length.
我的lon.shape()=(360,),lat.shape()=(181,),value.shape()=(181,360)，length肯定是不一样。
不知道问题出在哪？
这三个的栅格数必须一样吗？

忽如一夜春风来

咖喱嫑洋葱 发表于 2019-8-21 17:55
我用Rbf(lon, lat, value,function='linear')，报错，ValueError: All arrays must be equal length.
我 ...

emmmmm，怎么说呢，你需要让lon.shape和lat.shape 还有value.shape是一样的，就是360*181个点，每个点给一个经度一个纬度.就是说他们的维度都变成（181,360）#其实我当时用的时候好像二维数组不能用，所以我还用了flatten()转换成了一维的（181*360）。至于中间赋值过程，你给循环填充或者别的什么方法都行。lon[:,i]=lon[ i .] 。