построение SMAP спутниковых данных Swath (L2) - PullRequest
Новая
       36

построение SMAP спутниковых данных Swath (L2)

0 голосов
/ 15 мая 2019

Output Image looks like this.

Я пытаюсь построить продукт данных Swath Salinity для спутниковой системы SMAP (файл .nc), используя базовую карту. Я также пытался использовать картопий. 

import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
import numpy as np
import numpy.ma as ma
from netCDF4 import Dataset as nn

filename=r'D:\temp_2017_40km_v2\304\RSS_smap_sss_L2C_r14681.nc'
nc=nn(filename)
latitude_smap=nc.variables['cellat'][:,:,0]
latitude_smap=ma.masked_where(latitude_smap==0.0,latitude_smap)

longitude_smap=nc.variables['cellon'][:,:,0]
longitude_smap=ma.masked_where(longitude_smap==0.0,longitude_smap)
longitude_smap=longitude_smap-180
sss_smap=nc.variables['sss_smap'][:,:,0]

m = Basemap(projection = 'cyl',llcrnrlat=-90.0,llcrnrlon=-180.0,urcrnrlat=90.0,urcrnrlon=180.0,resolution='l') 
plt.figure(figsize=(8,6))
m.scatter(longitude_smap, latitude_smap,c=sss_smap,cmap=plt.cm.jet, s=50,linewidth=0)
m.drawcoastlines(linewidth=0.5,color = 'k')
m.drawparallels(np.arange(-90, 91, 45), labels=[True,False,False,True],color = 'k', fontsize=8)
m.drawmeridians(np.arange(0, 360, 45), labels=[True,False,False,True],color = 'k', fontsize=8)       
m.scatter(longitude_smap, latitude_smap, c=data, s=1, cmap=plt.cm.jet,edgecolors=None, linewidth=0,vmin=30,vmax=40)
plt.show()
plt.close()
...