| | |
|---|
| | | 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() |
|---|
| | | |
|---|
| | | coord = str(x) + ", " + str(y) + ", " + str(sid) |
|---|
| | | print(coord) |
|---|
| | | |
|---|
| | | f = open(args.out.replace(".tif", "_cs.txt"), 'w') |
|---|
| | | f.write(coord) |
|---|
| | | f.close() |
|---|
| | | |
|---|
| | | |
|---|
| | | # 初始化 |
|---|
| | | def init(): |
|---|
| | | # QgsApplication.setPrefixPath("C:\Program Files\QGIS 3.16", True) |
|---|
| | |
|---|
| | | |
|---|
| | | load_layers(prj, args) |
|---|
| | | merge(prj, args) |
|---|
| | | get_coord(prj, args) |
|---|
| | | |
|---|
| | | qgs.exitQgis() |
|---|
| | | |
|---|
