From 3be2cdf37ce7765880b5351cde33bfd393cae3df Mon Sep 17 00:00:00 2001
From: 13693261870 <252740454@qq.com>
Date: 星期三, 28 六月 2023 19:57:58 +0800
Subject: [PATCH] 1
---
ExportMap/Sources/merge.py | 11 ++++++++++-
1 files changed, 10 insertions(+), 1 deletions(-)
diff --git a/ExportMap/Sources/merge.py b/ExportMap/Sources/merge.py
index 3171edf..b0552f5 100644
--- a/ExportMap/Sources/merge.py
+++ b/ExportMap/Sources/merge.py
@@ -29,7 +29,7 @@
parser.add_argument("-src", type=str, default=get_full_path(), required=False)
parser.add_argument("-qgz", type=str, default=r"xyz.qgz", required=False)
parser.add_argument("-file", type=str, default=r"D:\xyz\dem\dem.txt", required=False)
- parser.add_argument("-out", type=str, default=r"D:\xyz\dem\dem\astgtm2.tif", required=False)
+ parser.add_argument("-out", type=str, default=r"D:\xyz\dem\out.tif", required=False)
return parser.parse_args()
@@ -104,7 +104,16 @@
rect = layer.extent()
x = (rect.xMinimum() + rect.xMaximum()) / 2
y = (rect.yMinimum() + rect.yMaximum()) / 2
+
sid = layer.crs().authid()
+ if sid != "EPSG:4490" and sid != "EPSG:4326":
+ sid = "EPSG:4490" if layer.crs().toWkt().find("CGCS2000") > -1 else "EPSG:4326"
+ transform = QgsCoordinateTransform(QgsCoordinateReferenceSystem(layer.crs().authid()), QgsCoordinateReferenceSystem(sid), prj)
+
+ p = QgsPoint(x, y)
+ p.transform(transform)
+ x = p.x()
+ y = p.y()
coord = str(x) + ", " + str(y) + ", " + str(sid)
print(coord)
--
Gitblit v1.9.3