Deep Zoom画像を完全復元：タイル画像からBigTIFFへの変換技術

dzc_output.xml              # メタデータ
dzc_output_files/
  ├── 0/
  │   └── 0_0.jpg          # レベル0の唯一のタイル
  ├── 1/
  │   ├── 0_0.jpg
  │   └── ...
  └── 16/                   # 最高レベル
      ├── 0_0.jpg
      ├── 0_1.jpg
      └── ...               # 数万枚のタイル

def fetch_xml_info(xml_url):
    root = ET.fromstring(response.content)

    # 複数の名前空間を試行
    namespaces = [
        '{http://schemas.microsoft.com/deepzoom/2008}Size',
        '{http://schemas.microsoft.com/deepzoom/2009}Size',
        'Size'
    ]

    for ns in namespaces:
        image_elem = root.find('.//' + ns)
        if image_elem is not None:
            break

    width = int(image_elem.attrib['Width'])
    height = int(image_elem.attrib['Height'])
    return config

def find_actual_max_level(base_url, format_ext):
    """実際に存在する最大レベルを検出"""
    for level in range(20, -1, -1):
        url = f"{base_url}{level}/0_0.{format_ext}"
        try:
            response = requests.head(url, timeout=10)
            if response.status_code == 200:
                return level
        except:
            continue
    return None

def download_tiles(tiles_list, base_url, level, format_ext):
    session = requests.Session()
    downloaded_tiles = []

    with ThreadPoolExecutor(max_workers=10) as executor:
        futures = {
            executor.submit(download_tile, base_url, level,
                          col, row, format_ext, session): (col, row)
            for col, row in tiles_list
        }

        with tqdm(total=len(futures)) as pbar:
            for future in as_completed(futures):
                result = future.result()
                if result:
                    downloaded_tiles.append(result)
                pbar.update(1)

    return downloaded_tiles

def reconstruct_image(tiles, tile_size, overlap):
    for col, row, tile_img in tiles:
        # オーバーラップを考慮した配置座標
        x = col * (tile_size - overlap)
        y = row * (tile_size - overlap)
        canvas.paste(tile_img, (x, y))

def save_bigtiff(image, output_path):
    # まずPNGで保存（ファイルサイズ制限なし）
    png_file = output_path.replace('.tif', '.png')
    image.save(png_file, format='PNG', compress_level=6)

    # tifffileライブラリでBigTIFFに変換
    import tifffile
    img_array = np.array(image)
    tifffile.imwrite(
        output_path,
        img_array,
        bigtiff=True,           # BigTIFF有効化
        compression='deflate',   # 圧縮
        tile=(256, 256)         # タイル化
    )

# 主要ライブラリ
import requests          # HTTP通信
from PIL import Image    # 画像処理
import numpy as np       # 配列操作
import tifffile          # BigTIFF対応
from concurrent.futures import ThreadPoolExecutor
from tqdm import tqdm    # プログレスバー

download_deepzoom.py           # 単一画像処理
batch_download_deepzoom.py     # バッチ処理
├─ fetch_xml_info()           # XML解析
├─ find_actual_max_level()    # レベル自動検出
├─ download_tile()            # タイルダウンロード
├─ reconstruct_image()        # 画像復元
└─ save_bigtiff()             # BigTIFF保存

session = requests.Session()
# 同一セッションでHTTP接続を再利用
response = session.get(url)

try:
    response = session.get(url, timeout=30)
    if response.status_code == 200:
        return Image.open(BytesIO(response.content))
except Exception as e:
    print(f"Failed to download tile: {e}")
    return None