谁抢了默默的冬瓜 2018-12-01 13:10 采纳率: 0%
浏览 5560
已结题

Python contourf函数如何只把大于某值的区域填色?

如题,如何把>0.5的区域填色,而其他区域不填色?

 names = ['num','lon','lat','XX','wind']
data = pd.read_table('D:/data/wind/out/fallingarea/2017050306.000',sep=' ',
                     skiprows=[0,1],names=names,index_col='num')
norm = Normalize()

#初始化地图
m = Basemap(projection='cyl',llcrnrlat=24.5,urcrnrlat=29.5,llcrnrlon=103,urcrnrlon=110.5)

# 将经纬度点转换为地图映射点
data['lon'], data['lat'] = m(*(data.lon.values, data.lat.values))

# 生成经纬度的栅格数据
numcols, numrows = 100, 100
xi = np.linspace(data['lon'].min(), 110, numcols)
yi = np.linspace(data['lat'].min(), 30, numcols)
xi, yi = np.meshgrid(xi, yi)

# 插值
x, y, z = data['lon'].values,data['lat'].values, data['wind'].values
#zi = gd((data[['lon','lat']]), data['wind'].values,(xi, yi),method='cubic')

func = Rbf(x,y,z,function='linear')
zi = func(xi,yi)

# 设置地图细节
m.drawcounties()
m.drawmapboundary()
m.readshapefile('D:/map/guizhoushp/Guizhou_city','states',drawbounds=True)
# 等值面图绘制
plt.scatter(x,y,z*20,marker='o',cmap='cool',alpha=0.8)
con = m.contourf(xi, yi, zi,1,alpha=0.8)
#画线的label
plt.clabel(con,inline=True,fontsize=10,colors='k')
  • 写回答

1条回答 默认 最新

  • devmiao 2018-12-01 15:55
    关注
    评论

报告相同问题?

悬赏问题

  • ¥15 有两个非常“自以为是”烦人的问题急期待大家解决!
  • ¥30 STM32 INMP441无法读取数据
  • ¥100 求汇川机器人IRCB300控制器和示教器同版本升级固件文件升级包
  • ¥15 用visualstudio2022创建vue项目后无法启动
  • ¥15 x趋于0时tanx-sinx极限可以拆开算吗
  • ¥15 pyqt信号槽连接写法
  • ¥500 把面具戴到人脸上,请大家贡献智慧,别用大模型回答,大模型的答案没啥用
  • ¥15 任意一个散点图自己下载其js脚本文件并做成独立的案例页面,不要作在线的,要离线状态。
  • ¥15 各位 帮我看看如何写代码,打出来的图形要和如下图呈现的一样,急
  • ¥30 c#打开word开启修订并实时显示批注