根據(jù)前面的文章python處理衛(wèi)星云圖獲取亮溫值可以從hdf文件中獲取某時(shí)刻某經(jīng)緯度的亮溫值。

然后就是把這些讀取的亮溫值可視化，繪制成TBB的衛(wèi)星云圖。

先展示一下結(jié)果：

1.讀取hdf文件

詳見上篇文章，這里不再贅述，再貼一下代碼。

主要是從全圓盤標(biāo)對(duì)稱文件經(jīng)緯度對(duì)照表獲取經(jīng)緯度和從hdf文件中獲取亮溫值。

# 從全圓盤標(biāo)對(duì)稱文件經(jīng)緯度對(duì)照表獲取經(jīng)緯度 def getLatLonFromdat():lonlatfile = 'F:/Satellite_Imagery/Code/NOM_ITG_2288_2288(0E0N)_LE.dat'with open(lonlatfile, 'rb') as f:#lon_fy = np.fromfile(f, count=2288 * 2288, dtype='float32') + 79 # 先存經(jīng)度，根據(jù)衛(wèi)星的不同加上對(duì)應(yīng)的經(jīng)度值#lat_fy = np.fromfile(f, count=2288 * 2288, dtype='float32') # 再存緯度data = np.fromfile(f, dtype='float32')data = data.reshape([2288, 2288, 2], order='F')#lon = lon_fy.reshape([2288, 2288], order='F')#lat = lat_fy.reshape([2288, 2288], order='F')lon = data[:, :, 0] + 104.5lat = data[:, :, 1]return lon, lat # 從hdf文件中獲取亮溫值 def getTBBFromhdf():hdfFile = h5py.File('F:/IR/Satellite_Imagery/IR_data/利奇馬/FY2G_FDI_ALL_NOM_20190811_1200.hdf', 'r')db1 = hdfFile['/CALChannelIR1']hw1 = hdfFile['/NOMChannelIR1']# db2 = hdfFile['/CALChannelIR2']# hw2 = hdfFile['/NOMChannelIR2']# db3 = hdfFile['/CALChannelIR3']# hw3 = hdfFile['/NOMChannelIR3']# db4 = hdfFile['/CALChannelIR4']# hw4 = hdfFile['/NOMChannelIR4']infoh = hdfFile['/NomFileInfo']# 查看衛(wèi)星的經(jīng)緯度lat_hdf = infoh[0][3]lon_hdf = infoh[0][4]# print(lat_hdf)# print(lon_hdf)hw = hw1[()]db = db1[()] # 獲取定標(biāo)表的值tb = np.zeros(shape=(2288, 2288)) # 2288*2288的圖像每個(gè)具體的亮溫值for i in range(2288):for j in range(2288):if hw[i][j] == 65535 or hw[i][j] == 65534:tb[i][j] = 0else:a = hw[i][j]tb[i][j] = db[0][a]tb = tb.Treturn tb

2. 畫圖

fig = plt.figure(figsize=(8, 6)) m = Basemap(projection='cyl', llcrnrlat=10, llcrnrlon=110, urcrnrlat=40, urcrnrlon=140) # 使用Basemap繪制地圖，這里可以讀取對(duì)應(yīng)的地圖shp文件。 m.drawcoastlines(color='black') m.drawstates(color='black') m.drawcountries(color='black') x, y = m(lon, lat) # 將lats / lons轉(zhuǎn)換為地圖投影坐標(biāo) # 繪制輪廓圖 # 這里data就是計(jì)算的亮溫值，x,y就是經(jīng)緯度投影的坐標(biāo) cf = m.contourf(x, y, data, levels=np.linspace(180, 301, 400), cmap='jet') cbar = m.colorbar(cf, location='right', size='5%', pad='2%')font = {'family': 'serif','color': 'darkred','weight': 'normal','size': 16,} cbar.set_label('Brightness Temperature ( K )', fontdict=font) m.drawmeridians(np.arange(110, 140, 5), labels=[0, 0, 0, 1]) m.drawparallels(np.arange(10, 40, 5), labels=[1, 0, 0, 0]) # 將最佳路徑集上的經(jīng)緯度映射到地圖上，再把該點(diǎn)繪制在地圖上 best_lon, best_lat = m(best_lon, best_lat) m.plot(best_lon, best_lat, 'o', color='fuchsia', ms=5) plt.xticks(fontsize=20) plt.yticks(fontsize=20) plt.title('Brightness Temperature(2019081112) ', fontdict=font) plt.savefig('test.png') plt.show()

注：讀取地圖文件時(shí)，可以使用readshapefile函數(shù)。

以2019081112時(shí)刻為例，對(duì)應(yīng)地圖10-40，110-140的區(qū)域，讀取對(duì)應(yīng)的hdf文件，繪制該點(diǎn)的衛(wèi)星云圖的亮溫值。

這是找了好久，調(diào)整過的結(jié)果，希望可以幫助看到這篇博客的人，有問題請(qǐng)?jiān)u論。

總結(jié)

以上是生活随笔為你收集整理的python绘制登陆时的卫星云图（TBB）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。