1

13693261870

2024-09-04 dbed0354a00fde9db87862b9f03db0fbf0d3aabc

 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()

@@ -68,6 +68,12 @@
            print("authid: " + args.authid)
            prj.setCrs(layer.crs())

        # for j in range(1, layer.bandCount() + 1):
        #     try:
        #         layer.dataProvider().setNoDataValue(j, 0)
        #     except Exception as e:
        #         print(e)

        prj.addMapLayer(layer)


@@ -81,8 +87,8 @@
        'INPUT': layers,
        'DATA_TYPE': 5,
        'OUTPUT': args.out,
        'NODATA_INPUT': 0,
        'NODATA_OUTPUT': 0,
        # 'NODATA_INPUT': 0,
        # 'NODATA_OUTPUT': 0,
        'OPTIONS': ''
    }
    print(ops)
@@ -96,6 +102,31 @@

    ops = get_ops(prj, args)
    processing.run("gdal:merge", ops)


# 获取中心点坐标
def get_coord(prj, args):
    layer = QgsRasterLayer(args.out, "layer_tif")
    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)

    f = open(args.out.replace(".tif", "_cs.txt"), 'w')
    f.write(coord)
    f.close()


# 初始化
@@ -113,6 +144,7 @@

    load_layers(prj, args)
    merge(prj, args)
    get_coord(prj, args)

    qgs.exitQgis()