nc文件中的X Y变量的单位是米,我现在需要更加X Y坐标进行绘制海图。如何将米转换成经纬度进行绘制海图
望采纳,
参考示例
from netCDF4 import Dataset
from mpl_toolkits.basemap import Basemap
ncfile = Dataset('filename.nc', 'r')
x = ncfile.variables['x'][:]
y = ncfile.variables['y'][:]
from pyproj import Proj
inproj = Proj(init='epsg:3857')
outproj = Proj(init='epsg:4326')
lon, lat = Proj.transform(inproj, outproj, x, y)
# 绘制海图
map = Basemap(projection='merc', llcrnrlon=lon.min(), llcrnrlat=lat.min(), urcrnrlon=lon.max(), urcrnrlat=lat.max(), resolution='f')
map.drawcoastlines()
map.plot(lon, lat, 'b-', markersize=0.1)
ncfile.close()