最近需要處理一些GSMap的遙感影像，GSMaP是全球降水計劃（GPM）時代最主要的高分辨率降水產品，下載需要申請，然后通過ftp地址進行下載。下載好以后用ENVI打開，發現沒有頭文件，自定義頭文件之后，影像極其的奇怪，大概是下面這個樣子，實在不符合想象，如果有同學知道怎么回事，希望您能告知我，感謝。現在我打算用python先轉成tif格式。
#!/usr/bin/env python
# -*- coding:utf-8 -*-
# Name : gsmap2TIFF.py
# Author : zengsk in NanJing
# Created: 2019/8/24 23:26
'''
說明：1.該腳本是讀取gsmap小時尺度降水數據，輸出為tiff
2.運行環境需要安裝python2 需要arcpy模塊
3.使用arcgis自帶的python環境(有arcpy模塊)
4.運行結果可以直接用arcgis打開
'''
# 導入module
import os
import glob
import numpy as np
import arcpy
import warnings
warnings.simplefilter("ignore") # 忽略警告
# 原始降水數據文件夾,可根據自己本地情況修改
sPath = r'./data'
oDir = r"./tiff/"
for fileName in glob.glob(sPath + '\*.dat'):
print("Processing... {0}".format(fileName))
ds = np.fromfile(fileName, dtype=np.float32)
ds = np.resize(ds, [1200, 3600])
ds[ds 0] = -999.00
ds[np.isnan(ds)] = -999.00 # NODATA_value
# 輸出為TIFF(注意：要用到arcpy模塊)
if not os.path.exists(oDir):
os.makedirs(oDir)
TiffName = oDir + os.sep + os.path.basename(fileName)[0:-4] + '.tif' # 輸出文件名(可根據實際情況改)
# arcpy.NumPyArrayToRaster()不清楚輸入參數可以查看arcpy的官方文檔
# 矩陣轉為柵格
raster = arcpy.NumPyArrayToRaster(ds, arcpy.Point(0, -60.0),
x_cell_size=0.1, y_cell_size=0.1, value_to_nodata=-999.00)
# 添加地理坐標系 GCS_WGS_1984
spatialRef = arcpy.SpatialReference(4326)
arcpy.DefineProjection_management(raster, spatialRef)
raster.save(TiffName)
print("\n++++++ Data Processing Successfully Completed ! ++++++")
到此這篇關于Python讀取GSMap數據的問題的文章就介紹到這了,更多相關Python讀取GSMap數據內容請搜索腳本之家以前的文章或繼續瀏覽下面的相關文章希望大家以后多多支持腳本之家！
