利用python GDAL库读写geotiff格式的遥感影像

   
   

     
     

      
      

     
     
     
     

      
      

       
       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)
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
                data[:,:,i]=dt.ReadAsArray(
       
       0,
       
       0,col,row)
      
      
     
     

     
     

      
      

     
     
     
     

      
      
     
       
       return data
      
      
     
     

     
     

      
      

     
     
     
     

      
      
 
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
       def array2raster(outpath,array,geoTransform,proj):
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           cols=array.shape[
       
       1]
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           rows=array.shape[
       
       0]
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           driver=gdal.GetDriverByName(
       
       ‘Gtiff’)
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           outRaster=driver.Create(newRasterfn,cols,rows,
       
       1,gdal.GDT_Byte)
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           outRaster.SetGeoTransform(geoTransform)
       
       #参数2,6为水平垂直分辨率，参数3,5表示图片是指北的
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           outband=outRaster.GetRasterBand(
       
       1)
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           outband.WriteArray(array)
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           outRaster.SetProjection(proj)
       
       #将几何对象的数据导出为wkt格式
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           outRaster.FlushCache()
      
      
     
     

     
     

      
      

     
     
     
     

      
      
 
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
       if _name==
       
       “_main_”:
      
      
     
     

     
     

      
      

     
     
     
     

      
      
 
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           data,geoTransform,proj=read_tiff(
       
       ‘d:/a.tif’)
      
      
     
     

     
     

      
      

     
     
     
     

      
      
 
      
      
     
     

     
     

      
      

     
     
     
     

      
      

       
           array2raster(
       
       “d:/b.tif”,np.zeros[
       
       2400,
       
       2400],geoTransform,proj)
      
      
     
     
   
   

利用python GDAL库读写geotiff格式的遥感影像
并生成与原影像具有相同地理坐标和投影坐标的geotiff格式图片



        </div>
            </div>
           转自：https://blog.csdn.net/qq_20340733/article/details/78316347