擅长:python、mysql、java
<p>我看到<code>lat_inds</code>有一个明显的问题,因为它有max-shape <code>1015*1359</code>,但是您尝试使用它作为纬度的索引,其大小为<code>1015</code>。所以在我看来,你应该先找到<code>lat_inds</code>和<code>lon_inds</code>的相似值,这两个点同时满足lon和lat的限制,然后使用这个数组来处理平坦的数据。比如:</p>
<pre><code>uni_ind=numpy.intersect1d(lat_inds,lon_inds)
T_subset=np.ravel(a.variables['T2'])[uni_ind]
</code></pre>
<p>将数组转换回二维可能还存在一些问题,因为我假设原始数据不是在柱坐标系中,因此结果子集可能不是矩形的。
这个代码没有测试,如果你共享原始数据文件,我也可以这样做。在</p>
<p>编辑:
为了正确绘图,使用遮罩比较容易,这个例子应该有足够的信息。在</p>
^{pr2}$