python gdal安装与简单使用

以你之姓@ 2023-10-11 11:31 24阅读 0赞

原文链接：[python gdal安装与简单使用][python gdal]

gdal安装  
方式一：在网址 https://www.lfd.uci.edu/~gohlke/pythonlibs/\#gdal 下载对应python版本的whl文件，在命令行中pip install whl文件完整路径安装（windows方式）。

方式二：  
命令行conda/pip search gdal查看版本，选择合适的版本（我的2.2.4），如果没有，使用方式一。  
命令行conda/pip install gdal=版本号，注意加上版本号，否则可能安装上老版本（windows/linux都可用）。

其他：centos 编译、配置 gdal

测试  
gdal包用于处理栅格数据，ogr用于处理矢量数据。  
以下程序为gdal处理栅格的简单应用。

1 from osgeo import gdal
     2 import numpy as np
     3 np.set_printoptions(threshold=np.inf)#使print大量数据不用符号...代替而显示所有
     4 
     5 dataset = gdal.Open("E:/RS_data/caijian1214/caijian.tif")
     6 
     7 print(dataset.GetDescription())#数据描述
     8 
     9 print(dataset.RasterCount)#波段数
    10 
    11 cols=dataset.RasterXSize#图像长度
    12 rows=(dataset.RasterYSize)#图像宽度
    13 
    14 xoffset=cols/2
    15 yoffset=rows/2
    16 
    17 band = dataset.GetRasterBand(3)#取第三波段
    18 r=band.ReadAsArray(xoffset,yoffset,1000,1000)#从数据的中心位置位置开始，取1000行1000列数据
    19 
    20 band = dataset.GetRasterBand(2)
    21 g=band.ReadAsArray(xoffset,yoffset,1000,1000)
    22 
    23 band = dataset.GetRasterBand(1)
    24 b=band.ReadAsArray(xoffset,yoffset,1000,1000)
    25 
    26 import cv2
    27 import matplotlib.pyplot as plt
    28 
    29 
    30 img2=cv2.merge([r,g,b])
    31 plt.imshow(img2)
    32 plt.xticks([]),plt.yticks([]) # 不显示坐标轴
    33 plt.show()

![987843-20190724110236844-196412220.png][]

转载于:https://www.cnblogs.com/mumu122GIS/p/11236694.html

[python gdal]: https://blog.csdn.net/nima1994/article/details/79207805/
[987843-20190724110236844-196412220.png]: https://img2018.cnblogs.com/blog/987843/201907/987843-20190724110236844-196412220.png