将0-360°的海温资料转换成-180-180°

66kqkq

最近在学习python编程绘图，读取海温资料的时候，发现python没有类似于ncl中lonFlip的函数，虽然可以加载basemap获取里面的类似函数，但是basemap已经停更了，安装也不能conda傻瓜式安装。在论坛看到了一个转换经度的帖子http://bbs.06climate.com/forum.p ... &extra=page%3D3。

参考这个帖子，加了一步来赋值经度坐标属性。

#将西经(-180到0)+360
ind=np.where(np.array(lon>=180))          #如果是-180到180转换成0到360，这里应该是<=0
ind=np.array(ind).reshape(np.size(ind,1))
nlon=np.array(lon)
nlon[ind]=nlon[ind]-360                          #如果是-180到180转换成0到360，这里应该是+360

#后利用numpy.argsort按照新的经度重新排序
sorted_indices=np.argsort(nlon,axis=0)

new_lon = np.arange(-180,180,2).astype(np.float32)
#sst_lonFlip = sst_D[:,:,sorted_indices]

time = sst_D.time
lat  = sst_D.lat
sst_lonFlip = xr.DataArray(sst_D[:,:,sorted_indices],coords=[time,lat,new_lon],dims=['time','lat','lon'])