网站建设资讯

NEWS

网站建设资讯

python安装gdal的方法-创新互联

小编给大家分享一下python安装gdal的方法,相信大部分人都还不怎么了解,因此分享这篇文章给大家参考一下,希望大家阅读完这篇文章后大有收获,下面让我们一起去了解一下吧!

创新互联建站服务项目包括松江网站建设、松江网站制作、松江网页制作以及松江网络营销策划等。多年来,我们专注于互联网行业,利用自身积累的技术优势、行业经验、深度合作伙伴关系等,向广大中小型企业、政府机构等提供互联网行业的解决方案,松江网站推广取得了明显的社会效益与经济效益。目前,我们服务的客户以成都为中心已经辐射到松江省份的部分城市,未来相信会继续扩大服务区域并继续获得客户的支持与信任!

GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。

方式一:在网址 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都可用)。

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

以上是python安装gdal的方法的所有内容,感谢各位的阅读!相信大家都有了一定的了解,希望分享的内容对大家有所帮助,如果还想学习更多知识,欢迎关注创新互联行业资讯频道!


名称栏目:python安装gdal的方法-创新互联
本文链接:http://njwzjz.com/article/ceihhd.html