遍历下载xyz瓦片数据

遍历并下载XYZ瓦片数据通常是指在一个地理信息系统（GIS）项目中，我们从Web地图服务（如TMS、WMTS等）获取具有层级结构（例如`z`, `x`, `y`）的地图图像数据。这个过程一般通过编程脚本完成，比如Python中的Rasterio库或JavaScript的TileJSON API。

以下是步骤概览：

1. 确定瓦片URL模式：每个瓦片都有一个唯一的地址，通常由服务提供者给出，格式通常是 `<base_url>/<z>/<x>/<y>.png` 或 `.jpg`。其中，`<z>` 是层级（zoom level），`<x>` 和 `<y>` 分别是经度和纬度坐标轴上的索引。

2. 定义范围：确定你想下载的区域范围，包括最小和最大层级（zoom range）、以及左上角和右下角的经纬度坐标 (`minZoom`, `maxZoom`, `west`, `south`, `east`, `north`）。

3. 创建循环：使用for循环遍历每一对 `(zoom, x, y)`，并在适当范围内生成瓦片 URL，然后下载文件到本地。

import requests
from PIL import Image

def download_tile(base_url, zoom, x, y, dest_dir):
    tile_url = f"{base_url}/{zoom}/{x}/{y}.png"
    response = requests.get(tile_url)
    if response.status_code == 200:
        with open(f"{dest_dir}/{zoom}_{x}_{y}.png", 'wb') as f:
            f.write(response.content)
        img = Image.open(f)
        # 可能需要对图片进行处理
        img.save(f)

# 示例范围
min_zoom = 6
max_zoom = 15
west = -90
east = 180
north = 90

# 下载指定范围内的瓦片
for zoom in range(min_zoom, max_zoom + 1):
    for x in range(int((east - west) / (2 ** zoom))):
        for y in range(int((north - south) / (2 ** zoom))):
            download_tile(your_base_url, zoom, x, y, "tiles_download_directory")