ฉันต้องการสร้างสคริปต์ Python ใน QGIS ที่จะใช้ข้อมูลบางอย่าง (shp + tif ในตัวอย่างของฉัน) ใช้ tamplate Map Composer (จากไฟล์) และส่งออกองค์ประกอบที่สร้างไปยังรูปภาพ png
แทบจะไม่มีประสบการณ์ในการเขียนโปรแกรม (น้อยกว่าความรู้พื้นฐานของ Python) ฉัน google snippets โค้ดบางส่วนและพยายามทำให้มันทำงานร่วมกันได้ ฉันใช้รหัส Map Composer จากบางคำถามที่ตอบก่อนหน้านี้: บันทึกมุมมองนักแต่งเพลงพิมพ์ / แผนที่ QGIS เป็น PNG / PDF โดยใช้ Python (โดยไม่ต้องเปลี่ยนอะไรในเค้าโครงที่มองเห็น)?
ฉันจัดการโหลดข้อมูลและเทมเพลต Map Composer (มีแผนที่เดียวและรายการคำอธิบายแผนภูมิ) แต่ภาพ png ที่ส่งออกของฉันมีกรอบแผนที่ว่างเปล่า (ไม่มีข้อมูลเวกเตอร์ / แรสเตอร์ภายในเฟรม) ไอเท็มคำอธิบายดูดี
ความช่วยเหลือเกี่ยวกับการทำให้รหัสนี้ทำงานอย่างไร
from qgis.core import *
import qgis.utils
from PyQt4 import QtCore, QtGui
from qgis import core, gui
# ADD VECTOR LAYER
data_folder = "D:/QGIS/dane/"
granica = "granica_SZ VI_UTM34.shp"
granica_name = granica[0:-4]
granica = data_folder + granica
granica_style = "granica_style.qml"
granica_style = data_folder + granica_style
granica = iface.addVectorLayer(granica, granica_name, "ogr")
granica.loadNamedStyle(granica_style)
if not granica.isValid():
print "Granica failed to load or already loaded!"
qgis.utils.iface.legendInterface().setLayerVisible(granica, True)
# ADD RASTER LAYER
landsat = "SZ_VI_LC8-190-022_2015-111_LGN00_OLI_TIRS_atm.TIF"
landsat_name = landsat[0:-4]
landsat = data_folder + landsat
landsat_style = "landsat_style.qml"
landsat_style = data_folder + landsat_style
landsat = iface.addRasterLayer(landsat, landsat_name)
landsat.loadNamedStyle(landsat_style)
qgis.utils.iface.legendInterface().setLayerVisible(landsat, True)
iface.zoomFull()
# MOVE RASTER DOWN (change order)
canvas = qgis.utils.iface.mapCanvas()
layers = canvas.layers()
root = QgsProject.instance().layerTreeRoot()
landsat_old = root.findLayer(landsat.id())
landsat_move = landsat_old.clone()
parent = landsat_old.parent()
parent.insertChildNode(2, landsat_move)
parent.removeChildNode(landsat_old)
# USE MAP COMPOSER TEMPLATE TO EXPORT IMAGE
from qgis.gui import QgsMapCanvas, QgsLayerTreeMapCanvasBridge
from PyQt4.QtXml import QDomDocument
from PyQt4.QtGui import QImage
from PyQt4.QtGui import QPainter
from PyQt4.QtCore import QSize
template_path = data_folder + 'template.qpt'
template_file = file(template_path)
# Set output DPI
dpi = 300
canvas = QgsMapCanvas()
template_file = file(template_path)
template_content = template_file.read()
template_file.close()
document = QDomDocument()
document.setContent(template_content)
ms = canvas.mapSettings()
composition = QgsComposition(ms)
composition.loadFromTemplate(document, {})
# You must set the id in the template
map_item = composition.getComposerItemById('map')
map_item.setMapCanvas(canvas)
map_item.zoomToExtent(canvas.extent())
# You must set the id in the template
legend_item = composition.getComposerItemById('legend')
legend_item.updateLegend()
composition.refreshItems()
dpmm = dpi / 25.4
width = int(dpmm * composition.paperWidth())
height = int(dpmm * composition.paperHeight())
# create output image and initialize it
image = QImage(QSize(width, height), QImage.Format_ARGB32)
image.setDotsPerMeterX(dpmm * 1000)
image.setDotsPerMeterY(dpmm * 1000)
image.fill(0)
# render the composition
imagePainter = QPainter(image)
composition.renderPage(imagePainter, 0)
imagePainter.end()
image.save(data_folder + "out3.png", "png")
QgsApplication.exitQgis()