Python 有多种处理 TIFF 文件的方法，以下为你介绍一些常见库及其使用方法： ### tifffile 库 tifffile 是一个开源项目，由 Christoph Gohlke 维护，可用于在 Python 环境下读取和写入 TIFF 文件。该项目基于 Python 编程语言，利用 NumPy 数组处理图像数据，并通过多种压缩和预测方案支持多种格式的 TIFF 文件 [^1]。 读取 TIFF 文件示例代码： ```python import tifffile # 读取 TIFF 文件 image = tifffile.imread('example.tif') print(image.shape) ``` 写入 TIFF 文件示例代码： ```python import tifffile import numpy as np # 创建一个示例数组 data = np.random.randint(0, 256, (100, 100), dtype=np.uint8) # 写入 TIFF 文件 tifffile.imsave('output.tif', data) ``` ### 其他处理 TIFF 文件的操作 还涉及经纬度、行列号、投影坐标之间的转换，获取指定位置的值，使用 SHP 文件裁剪 TIFF 文件等操作 [^2]。以下给出简单示例，实际应用中可能需要根据具体需求进行调整。 #### 经纬度转投影坐标 ```python import pyproj # 定义投影坐标系 proj = pyproj.Proj(proj='utm', zone=33, ellps='WGS84') # 经纬度坐标 lon, lat = 120, 30 # 转换为投影坐标 x, y = proj(lon, lat) print(x, y) ``` #### 获取指定位置的值 ```python import rasterio # 打开 TIFF 文件 with rasterio.open('example.tif') as src: # 获取指定位置的值 row, col = 10, 20 value = src.read(1)[row, col] print(value) ``` #### 使用 SHP 文件裁剪 TIFF 文件 ```python import rasterio from rasterio.mask import mask import geopandas as gpd import json # 读取 SHP 文件 shapefile = gpd.read_file('shapefile.shp') geoms = shapefile.geometry.values # 将几何对象转换为 GeoJSON 格式 geoms = [json.loads(geoms[i].to_json()) for i in range(len(geoms))] # 打开 TIFF 文件 with rasterio.open('example.tif') as src: # 裁剪 TIFF 文件 out_image, out_transform = mask(src, geoms, crop=True) out_meta = src.meta.copy() # 更新元数据 out_meta.update({"driver": "GTiff", "height": out_image.shape[1], "width": out_image.shape[2], "transform": out_transform}) # 保存裁剪后的 TIFF 文件 with rasterio.open('cropped.tif', "w", **out_meta) as dest: dest.write(out_image) ```