上传者: 38594687 
                                    |
                                    上传时间: 2021-11-07 17:30:21
                                    |
                                    文件大小: 33KB
                                    |
                                    文件类型: -
                                
                            
                                
                            
                                如下所示:
from osgeo import gdal
import numpy as np
def read_tiff(inpath):
  ds=gdal.Open(inpath)
  row=ds.RasterXSize
  col=ds.RasterYSize
  band=ds.RasterCount
  geoTransform=ds.GetTransform()
  proj=ds.GetTransform()
  data=np.zeros([row,col,band])
  for i in range(band):
   dt=ds.GetRasterBand(1)