python qgis 加载xyztile

在Python中，QGIS是一个强大的地理信息系统库，用于处理空间数据和进行地图制图。如果你想在QGIS中加载XYZ瓦片（也称为Web tiles），可以使用其内置的在线层功能或者通过GDAL（Geospatial Data Abstraction Library）插件。

以下是使用QGIS加载XYZ瓦片的基本步骤：

1. 安装所需的库：首先确保已安装了`gdal`和`qgis.core`库。如果没有，你可以通过pip安装：

   pip install gdal qgis

2. 使用`gdal.AddRasterBand()`加载瓦片：创建一个虚拟文件，并逐级添加瓦片URL作为每个级别（zoom level）的数据来源。例如，对于TMS（Tile Map Service）瓦片：

   from osgeo import gdal

   url_template = "http://example.com/tiles/{z}/{x}/{y}.png"
   driver = gdal.GetDriverByName('GTiff')
   dataset = driver.Create('/vsimem/tiled_raster.tif', None, None, 0)
   for zoom in range(max_zoom_level):
       band = dataset.GetRasterBand(zoom + 1) # QGIS从1开始计数
       options = ['TILED=YES', 'SPARSE_OK=TRUE']
       band.SetNoDataValue(0)  # 如果有透明度设置为None
       url = url_template.format(z=zoom, x=(2**zoom)*i+x, y=(2**zoom)*j+y)
       band.WriteUrl(0, url, options)

3. 创建QGIS在线层：将刚刚生成的临时文件添加到QGIS项目中，作为WMS服务或栅格数据源：

   from qgis.core import QgsRasterLayer

   raster_layer = QgsRasterLayer("/vsimem/tiled_raster.tif", "Tiled Layer")
   if raster_layer.isValid():
       QgsProject.instance().addMapLayer(raster_layer)