pythongdal安装与简单使用-创新互联

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都可用)。

其他:https://www.jb51.net/article/87675.htm

测试


gdal包用于处理栅格数据,ogr用于处理矢量数据。

以下程序为gdal处理栅格的简单应用。

from osgeo import gdal
import numpy as np
np.set_printoptions(threshold=np.inf)#使print大量数据不用符号...代替而显示所有
dataset = gdal.Open("E:/RS_data/caijian1214/caijian.tif")
print(dataset.GetDescription())#数据描述
print(dataset.RasterCount)#波段数
cols=dataset.RasterXSize#图像长度
rows=(dataset.RasterYSize)#图像宽度
xoffset=cols/2
yoffset=rows/2
band = dataset.GetRasterBand(3)#取第三波段
r=band.ReadAsArray(xoffset,yoffset,1000,1000)#从数据的中心位置位置开始,取1000行1000列数据
band = dataset.GetRasterBand(2)
g=band.ReadAsArray(xoffset,yoffset,1000,1000)
band = dataset.GetRasterBand(1)
b=band.ReadAsArray(xoffset,yoffset,1000,1000)
import cv2
import matplotlib.pyplot as plt
img2=cv2.merge([r,g,b])
plt.imshow(img2)
plt.xticks([]),plt.yticks([]) # 不显示坐标轴
plt.show()

网站标题:pythongdal安装与简单使用-创新互联
转载来源:http://azwzsj.com/article/ddjdso.html