基于python与GDAL获取遥感影像

最近做项目，发现网上关于利用python与GDAL读取遥感影像数据的资料比较少。因此整理了一些，以方便自己学习。

     1.首先确认自己电脑上安装好python版的GDAL库。这个会以后有时间慢慢整理。

      2. 导入gdal库

         from osgeo import gdal

      3. 注册GDAL库

         gdal.AllRegister()

      4. 启动GDAL所支持的文件格式驱动，这里的format为'tiff','jpg'等GDAL支持的栅格数据

          Driver = gdal.GetDriverByName(format)

      5.  打开文件，如果文件打开成功，返回1.如果文件不能打开，则返回零。

          fh = gdal.Open( filename )

      6.  fh为文件句柄，包含了文件的行列号等信息。

          bands = fh.RasterCount  //文件波段数

          xsize = fh.RasterXSize //文件行数

          ysize   = ysize = fh.RasterYSize //文件列数

          band_type = fh.GetRasterBand(1).DataType //数据类型

          projection = fh.GetProjection()  //投影参数信息

         geotransform = fh.GetGeoTransform() //文件几何信息

         文件几何信息是一维六元数组，主要是：

              ulx = geotransform[0] //左上x坐标

              resolutionx = geotransform[1] //x方向分辨率

             rotationx =  geotransform[2] //x方向旋转

             uly = geotransform[3] //左上y坐标

             rotationy =  geotransform[4] //y方向旋转

            resolutiony = geotransform[5] //y方向分辨率

        GDAL读取文件时把遥感数据倒置了，因此y方向分辨率为赋值。

        ct = GetRasterBand(1).GetRasterColorTable() //数据颜色表

     7. 读取影像数据

          band = fh.GetRasterBand(1) //获取波段

          data = band.ReadAsArray( xoff, yoff, xsize, ysize, t_xsize, t_ysize ) //获取影像数据

          xoffset ：x方向偏移

          yoffset： y方向偏移

          xsize： x方向数据大小

          ysize： y方向数据大小 

          t_xsize : 将图像转换为的x大小 

          t_ysize : 将图像转换为的y大小

          如果您要读取一整景影像，则如下：

                data = band.ReadAsArray( 0, 0, xsize, ysize)  

          如果您分块读取影像数据操作

                data = band.ReadAsArray( 0, 0, xsize, ysize)