VLMで地図を読み取る

def basic_map_analysis(image_path):
    prompt = """
    この地図画像を分析してください：

    1. 地図の種類（地形図、都市計画図等）
    2. 表示されている地域の特徴
    3. 主要な地物（道路、河川、建物等）
    4. 特徴的なポイント
    """
    return analyze_image(image_path, prompt)

def analyze_road_network(image_path):
    prompt = """
    この地図から道路網を分析してください：

    1. 主要道路（国道、県道等）
    2. 道路の走行方向
    3. 交差点・結節点
    4. 道路の特徴（幅員、種類等）
    """
    return analyze_image(image_path, prompt)

def analyze_land_use(image_path):
    prompt = """
    この地図から土地利用状況を分析してください：

    1. 住宅地の分布
    2. 商業地域
    3. 工業地域
    4. 農地・緑地
    5. 土地利用の傾向・特徴
    """
    return analyze_image(image_path, prompt)

def compare_maps(image_path1, image_path2):
    # 2つの地図を比較分析
    analysis1 = basic_map_analysis(image_path1)
    analysis2 = basic_map_analysis(image_path2)

    prompt = f"""
    以下の2つの地図分析結果を比較してください：

    【地図1】
    {analysis1}

    【地図2】
    {analysis2}

    共通点と相違点を整理してください。
    """
    return call_llm(prompt)

import requests
import base64
from pathlib import Path

class MapAnalyzer:
    def __init__(self, model="llava"):
        self.model = model
        self.base_url = "http://localhost:11434/api/generate"

    def encode_image(self, image_path):
        """画像をBase64エンコード"""
        with open(image_path, 'rb') as f:
            return base64.b64encode(f.read()).decode('utf-8')

    def analyze(self, image_path, prompt):
        """画像を分析"""
        image_data = self.encode_image(image_path)

        response = requests.post(
            self.base_url,
            json={
                'model': self.model,
                'prompt': prompt,
                'images': [image_data],
                'stream': False
            }
        )

        return response.json()['response']

    def basic_analysis(self, image_path):
        """基本分析"""
        prompt = """
        この地図画像を分析してください：
        1. 地図の種類
        2. 表示地域
        3. 主要な地物
        4. 特徴的なポイント
        """
        return self.analyze(image_path, prompt)

    def extract_features(self, image_path, feature_type):
        """特定の地物を抽出"""
        prompt = f"""
        この地図から{feature_type}を特定してください。
        見つかった{feature_type}をリストで回答してください。
        """
        return self.analyze(image_path, prompt)

# インスタンス作成
analyzer = MapAnalyzer()

# 基本分析
result = analyzer.basic_analysis("map.png")
print(result)

# 道路の抽出
roads = analyzer.extract_features("map.png", "道路")
print(roads)

# 建物の抽出
buildings = analyzer.extract_features("map.png", "建物")
print(buildings)

def batch_analyze_maps(folder_path, output_file):
    """フォルダ内の地図を一括分析"""
    analyzer = MapAnalyzer()
    results = []

    folder = Path(folder_path)
    for image_path in folder.glob("*.png"):
        print(f"Analyzing: {image_path.name}")

        result = {
            'filename': image_path.name,
            'analysis': analyzer.basic_analysis(image_path)
        }
        results.append(result)

    # 結果を保存
    import json
    with open(output_file, 'w', encoding='utf-8') as f:
        json.dump(results, f, ensure_ascii=False, indent=2)

    return results

# 実行
batch_analyze_maps("maps/", "analysis_results.json")

この地図から現場調査に必要な情報を抽出してください：
- アクセス道路
- 周辺施設
- 地形的な注意点
- 駐車可能な場所の推測

この都市計画図から土地評価に関連する情報を抽出してください：
- 用途地域
- 道路状況
- 周辺の土地利用
- 開発可能性

この地形図から災害リスクに関連する地形的特徴を分析してください：
- 河川との位置関係
- 傾斜地の有無
- 低地・窪地の有無
- 避難経路の確認

def thorough_analysis(image_path):
    """多角的な分析"""
    analyzer = MapAnalyzer()

    results = {
        '概要': analyzer.analyze(image_path,
            "この地図の概要を説明してください"),

        '道路': analyzer.analyze(image_path,
            "道路網について詳しく分析してください"),

        '地形': analyzer.analyze(image_path,
            "地形的特徴を分析してください"),

        '土地利用': analyzer.analyze(image_path,
            "土地利用状況を分析してください"),
    }

    return results

# QGISから画像をエクスポートしてVLMで分析
def analyze_qgis_export(qgis_project, extent, output_path):
    # 1. QGISで地図画像をエクスポート
    export_map_image(qgis_project, extent, output_path)

    # 2. VLMで分析
    analyzer = MapAnalyzer()
    result = analyzer.basic_analysis(output_path)

    # 3. 結果を返す
    return result

def create_report(image_path, output_path):
    """分析結果からレポートを作成"""
    analyzer = MapAnalyzer()

    # 詳細分析
    analysis = thorough_analysis(image_path)

    # レポート生成
    report = f"""
# 地図分析レポート

## 1. 概要
{analysis['概要']}

## 2. 道路網
{analysis['道路']}

## 3. 地形
{analysis['地形']}

## 4. 土地利用
{analysis['土地利用']}

---
*自動生成レポート*
"""

    with open(output_path, 'w', encoding='utf-8') as f:
        f.write(report)