QGIS使用python代码导出给定坐标图片

代码基于https://blog.youkuaiyun.com/x572722344/article/details/108121230进行修改,代码在QGIS内部编译器运行

# -*- coding: utf-8 -*-
from osgeo import ogr

# 像素[高, 宽]
px_geosize = [2.645859085290482, 2.6458015267176016]


# 待裁剪影像的坐标范围[min_x, min_y, max_x, max_y]
clip_ext = [15580385, 4335805
### WGS84 坐标经纬度查询工具 对于WGS84坐标系下经纬度的查询,可以通过多种方式实现。以下是几种常见的方式及其下载和使用方法: #### 1. **GDAL 工具** GDAL 是一个开源地理空间数据处理库,其中 `gdallocationinfo` 可用于提取栅格图像上的特定位置值,并支持不同坐标系统的转换。 - **安装 GDAL**: 用户可以从官方站点 (https://gdal.org/) 下载并安装适合操作系统的版本。 - **使用方法**: - 执行命令行工具时,指定 `-geoloc` 参数来输入投影坐标,或者通过 `-wgs84` 输入 WGS84 地理坐标[^1]。 ```bash gdallocationinfo raster_file.tif -wgs84 longitude latitude ``` 上述命令会返回对应于给定经度和纬度的位置像素值。 --- #### 2. **在线工具** 如果不想安装本地软件,可以选择一些成熟的在线服务来进行 WGS84 的查询与转换工作。 - **推荐网站**: 提供 WGS84 转换功能的一些知名平台包括但不限于: - EPSG.io: 支持基于 EPSG 编码的标准坐标系统之间的相互转换。 - MyGeodata Converter: 实现 GPS 数据与其他主流地图服务商坐标的互转[^3]。 这些工具通常提供直观界面,只需简单填写原始坐标即可获得目标格式的结果。 --- #### 3. **编程接口调用** 利用 Python 或其他语言编写脚本也可以完成自动化查询任务。下面是一个简单的例子展示如何借助 Pyproj 库执行此类运算: ```python from pyproj import Transformer def wgs84_to_custom(lon, lat): transformer = Transformer.from_crs("epsg:4326", "epsg:your_target_epsg_code") # 替换为目标EPSG编号 easting, northing = transformer.transform(lat, lon) return easting,northing result=wgs84_to_custom(119.496145 ,31.426896 ) # 示例中的北斗芯片获取的数据 [^2] print(result) ``` 注意需调整 `"epsg:your_target_epsg_code"` 部分为你实际需求的目标参考框架编码。 --- #### 4. **桌面应用程序** 除了线上资源外,还有专门设计用来管理复杂地理信息项目的独立程序可供选用。比如 QGIS 就内置强大的插件体系允许用户轻松导入导出以及变换各类矢量/栅格资料形式。 加载项目后,在属性表里找到相应字段应用表达式编辑器或几何计算器达成目的;另外 ArcMap/ArcGIS Pro 同样具备相似的功能集只是可能更加商业化一点。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值