根據前面的文章python處理衛星云圖獲取亮溫值可以從hdf文件中獲取某時刻某經緯度的亮溫值。

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

先展示一下結果：

1.讀取hdf文件

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

主要是從全圓盤標對稱文件經緯度對照表獲取經緯度和從hdf文件中獲取亮溫值。

# 從全圓盤標對稱文件經緯度對照表獲取經緯度 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 # 先存經度，根據衛星的不同加上對應的經度值#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']# 查看衛星的經緯度lat_hdf = infoh[0][3]lon_hdf = infoh[0][4]# print(lat_hdf)# print(lon_hdf)hw = hw1[()]db = db1[()] # 獲取定標表的值tb = np.zeros(shape=(2288, 2288)) # 2288*2288的圖像每個具體的亮溫值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繪制地圖，這里可以讀取對應的地圖shp文件。 m.drawcoastlines(color='black') m.drawstates(color='black') m.drawcountries(color='black') x, y = m(lon, lat) # 將lats / lons轉換為地圖投影坐標 # 繪制輪廓圖 # 這里data就是計算的亮溫值，x,y就是經緯度投影的坐標 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]) # 將最佳路徑集上的經緯度映射到地圖上，再把該點繪制在地圖上 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()

注：讀取地圖文件時，可以使用readshapefile函數。

以2019081112時刻為例，對應地圖10-40，110-140的區域，讀取對應的hdf文件，繪制該點的衛星云圖的亮溫值。

這是找了好久，調整過的結果，希望可以幫助看到這篇博客的人，有問題請評論。

總結

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

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。