管道基础大数据平台系统开发-【CS】-ExportMap
1
13693261870
2023-06-28 3be2cdf37ce7765880b5351cde33bfd393cae3df
1
已修改1个文件
11 ■■■■■ 文件已修改
ExportMap/Sources/merge.py 11 ●●●●● 补丁 | 查看 | 原始文档 | blame | 历史
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)