pip install osgeo - Shell-Bash

简介

osgeo 是一个与开源地理空间软件相关的Python模块。它提供了许多功能，例如读取各种地理空间格式、操作地理空间数据、投影转换等。该模块旨在为Python程序员和GIS开发者提供一个方便的方式来处理地理空间数据。

安装

使用 pip 工具可方便地进行 osgeo 模块的安装，执行以下命令：

pip install osgeo

使用

导入 osgeo 模块：

import osgeo

使用 osgeo 模块中的 gdal 库读取一个GeoTIFF格式的图像：

from osgeo import gdal

ds = gdal.Open('path/to/image.tif')

# 输出图像行列数
print(ds.RasterXSize, ds.RasterYSize)

# 输出投影信息
print(ds.GetProjection())

# 输出地理参考信息
print(ds.GetGeoTransform())

示例

以下是一个使用 osgeo 模块的示例，该示例读取一个GeoTIFF格式的红外遥感图像，计算图像的植被指数（Normalized Difference Vegetation Index，简称NDVI），并将NDVI值保存到一个GeoTIFF格式的新图像中。

import numpy as np
from osgeo import gdal

# 打开红外图像
ir_ds = gdal.Open('path/to/infrared.tif')
ir_band = ir_ds.GetRasterBand(1)
ir_data = ir_band.ReadAsArray()

# 打开可见光图像
vis_ds = gdal.Open('path/to/visible.tif')
vis_band = vis_ds.GetRasterBand(1)
vis_data = vis_band.ReadAsArray()

# 计算NDVI
ndvi_data = (ir_data - vis_data) / (ir_data + vis_data)

# 创建新图像
driver = gdal.GetDriverByName('GTiff')
ndvi_ds = driver.Create('path/to/ndvi.tif', ir_ds.RasterXSize, ir_ds.RasterYSize, 1, gdal.GDT_Float32)
ndvi_ds.SetProjection(ir_ds.GetProjection())
ndvi_ds.SetGeoTransform(ir_ds.GetGeoTransform())

# 写入数据
ndvi_band = ndvi_ds.GetRasterBand(1)
ndvi_band.WriteArray(ndvi_data)

# 释放资源
ndvi_band = None
ndvi_ds = None
ir_band = None
ir_ds = None
vis_band = None
vis_ds = None

整合其他库

与其他库结合使用，可以更加方便地处理地理空间数据。例如 geopandas 库可以用来读取和分析各种地理空间格式数据文件。

import geopandas as gpd
from osgeo import gdal

# 读取Shapefile文件
gdf = gpd.read_file('path/to/shapefile.shp')

# 读取GeoTIFF格式图像
ds = gdal.Open('path/to/image.tif')

# 投影转换
gdf = gdf.to_crs(ds.GetProjection())