Błąd OpenStreetMap PyQGIS

Question

Mam problem z opracowaniem wtyczki Pythona, próbuję wyświetlić mapę OpenStreetMap za warstwami, tylko mapa nie jest zsynchronizowana z warstwami, nie może jej znaleźć ani nawet lokalizacji.

Położenie geograficzne mapy nie odpowiada położeniu moich warstw.

#Vlayer is defined above the code / Vlayer est définie au-dessus du code 

urlWithParams = 'type=xyz&url=http://a.tile.openstreetmap.org/%7Bz%7D/%7Bx%7D/%7By%7D.png&zmax=19&zmin=0&crs=EPSG3857'
maplayer = QgsRasterLayer(urlWithParams, 'OpenStreetMap', 'wms')
settings = QgsMapSettings()
settings.setLayers([vlayer,maplayer])
settings.setBackgroundColor(QColor(255, 255, 255))
settings.setOutputSize(QSize(800, 690))

settings.setExtent(vlayer.extent())

render = QgsMapRendererParallelJob(settings)

render.start()
render.waitForFinished()

img = render.renderedImage()
img.save(QgsApplication.qgisSettingsDirPath()+"python/plugins/nro_print/img/render.png", "png")

wynik:

Answer 1

Musisz ustawić docelowy CRS w ustawieniach modułu renderującego. Poniższe działa dla mnie.

Utwórz testową warstwę pamięci

vlayer = QgsVectorLayer('Point?crs=epsg:4647&field=ID:integer(3)&field=value:double','test',"memory")

mypts = [QgsPointXY(32500000,5700000),
  QgsPointXY(32500010,5700010),
  QgsPointXY(32500010,5700060),
  QgsPointXY(32500040,5700090)
]

for p in mypts:
    pt = QgsGeometry.fromPointXY(p)
    f = QgsFeature(vlayer.fields())
    f.setGeometry(pt)
    vlayer.dataProvider().addFeatures([f])

vlayer.updateExtents()

Dodaj warstwę OSM, zanotuj httpsi tile.nie a.tile..

urlWithParams = 'type=xyz&url=https://tile.openstreetmap.org/%7Bz%7D/%7Bx%7D/%7By%7D.png&zmax=19&zmin=0&crs=EPSG3857'
maplayer = QgsRasterLayer(urlWithParams, 'OpenStreetMap', 'wms')

Zdefiniuj ustawienia, to Twoje.

settings = QgsMapSettings()
settings.setLayers([vlayer,maplayer])
settings.setBackgroundColor(QColor(255, 255, 255))
settings.setOutputSize(QSize(800, 690))

To jest to, co musiałem ustawić, w przeciwnym razie warstwa OSM nie będzie renderowana. edycja: Najwyraźniej musi być taki sam jak CRS vlayera. Jest to również wspomniane w sekcji 10.2 książki kucharskiej PyQIS .

settings.setDestinationCrs(vlayer.crs())

settings.setExtent(vlayer.extent().buffered(30))

render = QgsMapRendererParallelJob(settings)

render.start()
render.waitForFinished()

img = render.renderedImage()
img.save("c:/temp/render.png", "png")

Wynik, który jest prawidłowo umieszczony

Pełny kod jako GIST .