[已解决]python作图有一行语句看不懂，求教

晋立里

用matplotlib作图，有一行语句不知道什么意思，向各位求教了。

...
no2=nc_no2[20,0,:,:]

ny=no2.shape[0];nx=no2.shape[1]
lon,lat=map.makegrid(nx,ny)

x, y = map(lon, lat)
...
a=map.contourf(x,y,no2)

问题1：
加粗的这一步不知道什么意思
但如果直接用a=map.contourf(lon,lat,no2)的话出不来图

请问x, y = map(lon, lat)这一步什么意思呢，没有找到用法

问题2：
另外还有个地方，
nc_no2是四维的，分别时间、层数、纬度、经度
上面的代码中我想去经纬度，为什么用了ny=no2.shape[0];nx=no2.shape[1]，
不应该是3,4两个维度吗，不应该是ny=no2.shape[2];nx=no2.shape[3]吗？

先谢谢各位了

weinihou

你竟然没有把  map= XXX  写上。。。。
如果是basemap（当然，我瞎猜的，不是的话无视好了），那是在做坐标投影，就是把需要的经纬度变成绘图的时候对应投影下的坐标，你可以输出看看，如果是这样的话单位是距离，好像是m


你要看你的数据读出来是怎样的，比如你输出下  no2.shape
你取用no2的时候不是已经变了么~~
“no2=nc_no2[20,0,:,:]”
取用了第一维和第二维分别取用了第21和第1个，只是单个，所以相当于你只得到个二维的啊

8828

x,y是经纬度的坐标，contourf的关键字里默认用x,y 你可以设置用经纬度

晋立里

8828 发表于 2017-11-14 21:41
x,y是经纬度的坐标，contourf的关键字里默认用x,y 你可以设置用经纬度

啊我明白了！谢谢前辈！

晋立里

weinihou 发表于 2017-11-14 22:01
你竟然没有把  map= XXX  写上。。。。
如果是basemap（当然，我瞎猜的，不是的话无视好了），那是在做坐 ...

啊我明白了！两个问题都懂了，非常非常感谢！！！

weinihou

晋立里 发表于 2017-11-15 01:18
啊我明白了！两个问题都懂了，非常非常感谢！！！