Python读取tif图像并提取经纬度信息

首先，需要安装rasterio和geopandas库。可以通过pip命令进行安装：

pip install rasterio geopandas

接下来，使用以下代码读取tif图像并提取经纬度信息：

import rasterio
import geopandas as gpd
from shapely.geometry import Point
# 打开tif图像文件
with rasterio.open('path/to/image.tif') as src:
# 获取地理参考信息
gcps = src.gcps
# 获取地理坐标系信息
crs = src.crs
# 获取行数和列数
rows, cols = src.shape
# 创建一个空的GeoDataFrame，用于存储经纬度数据
gdf = gpd.GeoDataFrame()
# 循环遍历每个像素点，并提取其经纬度信息
for i in range(rows):
for j in range(cols):
# 获取像素点的坐标值
x = gcps[0][j]
y = gcps[0][i]
# 将坐标值转换为经纬度数据
lon, lat = x[0], y[0]
# 创建一个点对象，并添加到GeoDataFrame中
point = Point(lon, lat)
gdf = gdf.append({'geometry': point}, ignore_index=True)
# 设置GeoDataFrame的坐标系为输入tif图像的坐标系
gdf.crs = crs

以上代码将读取tif图像中的每个像素点，并提取其经纬度信息。提取的经纬度数据将存储在一个GeoDataFrame中，该DataFrame的坐标系与输入tif图像的坐标系相同。可以根据需要将提取的经纬度数据用于其他应用场景，例如地理信息系统（GIS）分析和地图可视化等。请注意，以上代码仅适用于具有地理参考信息的tif图像文件。如果tif图像没有地理参考信息，则需要使用其他方法进行坐标转换。