Maison  >  Article  >  développement back-end  >  Comment installer gdal en python

Comment installer gdal en python

(*-*)浩
(*-*)浩original
2019-06-24 15:23:038233parcourir

GDAL (Geospatial Data Abstraction Library) est une bibliothèque open source de conversion de données spatiales raster sous licence X/MIT. Il utilise un modèle de données abstrait pour exprimer les différents formats de fichiers pris en charge. Il dispose également d'une gamme d'outils de ligne de commande pour la conversion et le traitement des données.

Comment installer gdal en python

Méthode 1 : Téléchargez le fichier whl correspondant à la version python sur le site https://www.lfd.uci.edu/~gohlke/pythonlibs/ #gdal. Dans la ligne de commande pip install whl fichier d'installation du chemin complet (mode Windows). (Apprentissage recommandé : Tutoriel vidéo Python)

Méthode 2 :

Ligne de commande conda/pip recherche gdal pour afficher la version, sélectionnez la version appropriée (la mienne est 2.2.4), sinon, utilisez la première méthode.

Ligne de commande conda/pip install gdal=numéro de version, attention à bien ajouter le numéro de version, sinon une ancienne version risque d'être installée (disponible pour windows/linux).

Le package gdal est utilisé pour traiter les données raster et ogr est utilisé pour traiter les données vectorielles.
Le programme suivant est une simple application de gdal pour traiter les rasters.

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

Pour plus d'articles techniques liés à Python, veuillez visiter la colonne Tutoriel Python pour apprendre !

Ce qui précède est le contenu détaillé de. pour plus d'informations, suivez d'autres articles connexes sur le site Web de PHP en chinois!

Déclaration:
Le contenu de cet article est volontairement contribué par les internautes et les droits d'auteur appartiennent à l'auteur original. Ce site n'assume aucune responsabilité légale correspondante. Si vous trouvez un contenu suspecté de plagiat ou de contrefaçon, veuillez contacter admin@php.cn