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()

posted @ 2019-07-24 11:03  mumu122  阅读(29730)  评论(0编辑  收藏  举报