- 积分
- 653
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-4-21
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2018-6-19 07:10:54
|
显示全部楼层
王老师好,可能是我讲得不够清楚,我的数据是把nc变量导出后经过批量数据处理得到的csv文档,里面是360*720格的降水次数,不带有经纬度,只有一个变量就是降水次数,所以在导入文件的过程中必须要把每个降水次数的经纬度再赋值上去。我参照了几个您的脚本之后总结出了大概这样的脚本。已知经纬度是0.25-359.75,-89.75-89.75,精度为0.5.但是我的脚本一直处理不对,卡在这里很久了,所以想请教一下您。
fn = ('/Users/v/Downloads/freq.csv')
lon = [0.25, 359.75, 0.5]
lat = [-89.75, 89.75, 0.5]
axesm()
mlayer = shaperead('/Users/v/Downloads/Meteoinfo/map/country..shp')
geoshow(mlayer, facecolor=(220,220,220), edgecolor=None)
layer = scatterm(lon, lat, fn, edge=False, size=2)
xlim(0.25, 359.75)
ylim(-89.75, 89.75)
title('Flood plain area')
ylabel('Latitude')
xlabel('Longitude')
colorbar(layer) |
|