python读取水色遥感数据及可视化（截取特定区域）

冰原1

学习札记：
python在水色遥感中的应用（1）：水色遥感海表温度读取及可视化

1.python开源共享，便于和我一样的初学者学习，所以代码及数据共享啦
2.论坛里python读取NC格式的帖子也蛮多，从全球数据提取特定区域，这一块稍微有点模糊，可能初学者会有困扰，这里写出详细步骤


数据：水色遥感传感器MODIS全球数据，NC格式，
A201006.L3m_MO_SST_sst_4km.rar (28.36 MB, 下载次数: 160)

2017-4-23 23:05 上传

点击文件名下载附件



代码：
#class
import netCDF4 as nc
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
import numpy as np
import cmaps



#read nc dataset
obj=nc.Dataset('E:\\python\\A201006.L3m_MO_SST_sst_4km.nc')
print obj

#read variables
lat = obj.variables['lat'][:]
lon = obj.variables['lon'][:]
sst = obj.variables['sst'][:]

＃extract  variable in South China Sea （12-22Ｎ; 105-122E）
latindex = np.arange(len(lat))
lonindex = np.arange(len(lon))
latindex2 = np.max(latindex[lat >= 12])
latindex22 = np.min(latindex[lat <= 24])
lonindex2 = np.min(lonindex[lon >= 105])
lonindex22 = np.min(lonindex[lon >= 122])

#Extract  scope of variable
sst2 = sst[latindex22:latindex2, lonindex2:lonindex22]
sst3 = sst2[::-1, :]
lat2 = lat[latindex22:latindex2]
lat2 = lat2[::-1]
lon2 = lon[lonindex2:lonindex22]

#set grid
x, y = m(*np.meshgrid(lon2, lat2))

#basemap
m = Basemap(projection = 'cyl',llcrnrlat = 12,urcrnrlat = 24,llcrnrlon = 105,urcrnrlon = 122,resolution='l')

fig = plt.figure()
ax = fig.add_axes([0.05,0.05,0.9,0.9])
im = m.pcolormesh(x, y, sst3, vmin = 24, vmax = 33, shading='flat',cmap = cmaps.BkBlAqGrYeOrReViWh200,latlon=True)

m.fillcontinents(color='grey')
m.drawparallels(np.arange(12.,24.,4.),labels=[1,0,0,0]) # draw parallels
m.drawmeridians(np.arange(105.,122.,4.),labels=[0,0,0,1]) # draw meridians
m.drawcoastlines()
m.drawcountries()



#set sclae of colorbar
cbar = m.colorbar(im)
cbar.set_clim(24,33)
cbar.set_ticks(np.linspace(24,32,9))
cbar.set_ticklabels(('24', '26', '27','28','29','30','31','32','32'))

#add title
plt.title('SST')

#output
plt.savefig('sst.png', dpi=100, bbox_inches='tight')

plt.show()



Result





参考资料
1.http://bbs.06climate.com/forum.php?mod=viewthread&tid=43521（家园论坛）
2.http://matplotlib.org/basemap/users/examples.html
3.http://bbs.06climate.com/forum.php?mod=viewthread&tid=50262（家园论坛）

最后谢谢家园里所有写python帖子的大神，学习了。


2017/04/23

冰原1

附件里的数据，帖子功能里找了好久，不知道如何取消-5的金钱，如果版主看到，方便的话请帮忙取消。谢谢

adam_5

正在学习python，好东西，谢谢楼主

lgz

好东西，谢谢楼主

贫道敬孔

多谢楼主分享{:5_213:}

状态无

{:5_235:}{:5_235:}{:5_235:}{:5_235:}{:5_235:}{:5_235:}{:5_235:}{:5_235:}

Square

多谢楼主分享

Square

多谢楼主分享

Square

多谢楼主分享

Square

多谢楼主分享