安装gdal库
sudo add-apt-repository ppa:ubuntugis/ppa
sudo apt-get update
sudo apt-get install gdal-bin
sudo apt-get install libgdal-dev
export CPLUS_INCLUDE_PATH=/usr/include/gdal
export C_INCLUDE_PATH=/usr/include/gdal
gdal-config --version #(得到version)
pip install GDAL==version
代码实现
from osgeo import gdal
# 影像文件路径
image_file = '你的tif图像地址'
# 打开影像文件
dataset = gdal.Open(image_file)
# 检查是否成功打开影像
if not dataset:
print(f"无法打开影像文件:{image_file}")
else:
print("影像成功打开!")
# 打印影像基本信息
print(f"影像大小:{dataset.RasterXSize} x {dataset.RasterYSize}")
print(f"波段数量:{dataset. RasterCount}")
print(f"坐标投影:{dataset. GetProjection()}")
# 获取地理变换信息
geotransform = dataset.GetGeoTransform()
if geotransform:
print(f"左上角坐标:({geotransform[0]}, {geotransform[3]})")
print(f"像素分辨率:({geotransform[1]},{geotransform[5]})")
# 读取第一波段数据
band = dataset.GetRasterBand(1)
band_data = band.ReadAsArray()
# 打印波段统计信息
print (f"第一波段最小值:{band.GetMinimum()}")
print (f"第一波段最大值:{band.GetMaximum()}")