tif数据与某省shp矢量文件有一致的坐标系

可以使用GDAL库中的Python绑定来读取tif数据并指定其坐标系为与某省.shp文件一致的坐标系。具体步骤如下:

1.导入必要的库

python 复制代码
from osgeo import gdal, osr

2.读取省.shp文件并获取其坐标系

python 复制代码
# 读取某省.shp文件
province_ds = ogr.Open('province.shp')
province_layer = province_ds.GetLayer()
python 复制代码
# 获取某省.shp文件的坐标系
province_srs = province_layer.GetSpatialRef()

3.读取tif数据并获取其坐标系

python 复制代码
# 读取tif数据
tif_ds = gdal.Open('data.tif')
python 复制代码
# 获取tif数据的坐标系
tif_srs = osr.SpatialReference()
tif_srs.ImportFromWkt(tif_ds.GetProjection())

4.将tif数据的坐标系转换为与省.shp文件一致的坐标系

python 复制代码
# 创建坐标系转换对象
transform = osr.CoordinateTransformation(tif_srs, province_srs)

进行坐标系转换

python 复制代码
transformed_tif_ds = gdal.AutoCreateWarpedVRT(tif_ds, None, tif_srs.ExportToWkt())

5.使用转换后的tif数据进行后续操作

python 复制代码
# 读取转换后的tif数据
transformed_tif_data = transformed_tif_ds.ReadAsArray()

6.进行后续操作

相关推荐
悠悠小茉莉25 分钟前
Win11 安装 Visual Studio(保姆教程 - 更新至2025.07)
c++·ide·vscode·python·visualstudio·visual studio
m0_6256865540 分钟前
day53
python
Real_man1 小时前
告别 requirements.txt,拥抱 pyproject.toml和uv的现代Python工作流
python
站大爷IP2 小时前
Python文件操作的"保险箱":with语句深度实战指南
python
运器1232 小时前
【一起来学AI大模型】算法核心:数组/哈希表/树/排序/动态规划(LeetCode精练)
开发语言·人工智能·python·算法·ai·散列表·ai编程
巴里巴气4 小时前
selenium基础知识 和 模拟登录selenium版本
爬虫·python·selenium·爬虫模拟登录
19894 小时前
【零基础学AI】第26讲:循环神经网络(RNN)与LSTM - 文本生成
人工智能·python·rnn·神经网络·机器学习·tensorflow·lstm
JavaEdge在掘金4 小时前
Redis 数据倾斜?别慌!从成因到解决方案,一文帮你搞定
python
ansurfen4 小时前
我的第一个AI项目:从零搭建RAG知识库的踩坑之旅
python·llm
前端付豪5 小时前
20、用 Python + API 打造终端天气预报工具(支持城市查询、天气图标、美化输出🧊
后端·python