MOD16 ET tif格式 数据处理

南极的七色鹿

MOD16数据处理 - 简书 (jianshu.com)
## read the GTIFF datadataset = gdal.Open(r'Q:\ET_ET_1km.tif')# get the x and y sizexsize = dataset.RasterXSizeysize = dataset.RasterYSize# get the information of orginal coordinate (upper left point) and range of each grid# coordinate[0] -----> longitude# coordinate[1] -----> longitude range of one grid# coordinate[3] -----> latitute# coordinate[5] -----> latitute range of one gridcoordinate = dataset.GetGeoTransform()# read the ET datadata = dataset.ReadAsArray()# create the grid mesh using numpyxrange = range(0, xsize)yrange = range(0, ysize)x, y = np.meshgrid(xrange, yrange)lon = coordinate[0] + x*coordinate[1] + y*coordinate[2]lat = coordinate[3] + x*coordinate[4] + y*coordinate[5]