python读取GeoTIFF文件

在Python中，我们可以使用一些专业的库来读取和处理GeoTIFF文件，因为GeoTIFF是一种地理空间图像数据的标准文件格式。其中最常用的是`rasterio`和`geopandas`库。

1. `rasterio`库：它是一个高级的GDAL（Geospatial Data Abstraction Library）接口，专门用于读写多种栅格文件格式，包括GeoTIFF。你可以使用`with`语句打开文件并获取其元数据和数据：

from rasterio import open

# 打开GeoTIFF文件
with rio.open('file_path.tif') as src:
    data = src.read()  # 获取所有波段的数据
    metadata = src.meta  # 获取元数据

2. `geopandas`库：如果你需要将GeoTIFF视为GIS数据结构（如DataFrame），可以先使用`rasterio`读取，然后转换到`geopandas`：

import geopandas as gpd
from rasterio import read_as_array

# 使用rasterio读取数据
band_data = read_as_array('file_path.tif')

# 创建GeoDataFrame
data_df = gpd.GeoDataFrame(index=range(len(band_data)), 
                          columns=['Band_{}'.format(i) for i in range(band_data.shape[0])],
                          geometry=gpd.points_from_xy(*np.meshgrid(src.bounds.left, src.bounds.bottom)))

data_df['band_values'] = band_data