Summary
This chapter covers overlay operations—intersection, union, and difference—that let you combine or compare geometric datasets by computing their overlapping, merged, or exclusive areas.

General domain of usage
geospatial analysis

**Video Description:**  
This video introduces overlay operations in geospatial analysis, focusing on the concepts of intersection, union, and difference. You will see how these operations are used to combine or compare spatial datasets, and why they are essential for spatial data analysis workflows.

オーバーレイ操作は地理空間解析の基本的な要素であり、異なる空間データセット間の比較、結合、関係抽出を可能にします。これらの操作は、「2つのフィーチャが重なっている場所はどこか」「2つの領域の合計面積はどれくらいか」「ある領域のうち、他の領域でカバーされていない部分はどこか」といった問いに答える際に特に有用です。最も一般的なオーバーレイ操作は、**intersection（交差）**、**union（和集合）**、**difference（差分）**の3つであり、それぞれ独自の分析目的があります。

**intersection（交差）**操作は、2つのデータセットが重なっている空間領域を特定するのに役立ちます。これは、環境保護区域と開発予定地が重なる領域を特定したい場合などに重要です。一方、**union（和集合）**操作は、両方のデータセットの全てのジオメトリを結合し、空間範囲と属性を統合します。これは、元のいずれかのデータセットでカバーされている全領域を表すデータセットを作成する際に有用です。**difference（差分）**操作は、一方のデータセットからもう一方の領域を差し引くことで、最初のデータセットにのみ存在する部分を明らかにします。これは、例えば土地利用図から水域を除外したい場合などに利用されます。

オーバーレイ操作は通常、**geopandas**の`overlay()`メソッドを使用して実行され、これらの複雑な空間操作を効率化し、属性データの適切な管理を保証します。

import geopandas as gpd

# Create two simple GeoDataFrames with polygons
poly1 = gpd.GeoDataFrame({'id': [1]}, geometry=gpd.GeoSeries.from_wkt(['POLYGON((0 0, 2 0, 2 2, 0 2, 0 0))']))
poly2 = gpd.GeoDataFrame({'id': [2]}, geometry=gpd.GeoSeries.from_wkt(['POLYGON((1 1, 3 1, 3 3, 1 3, 1 1))']))

# Intersection: find overlapping area
intersection = gpd.overlay(poly1, poly2, how='intersection')
print("Intersection result:")
print(intersection)

# Union: combine all areas
union = gpd.overlay(poly1, poly2, how='union')
print("\nUnion result:")
print(union)


import unittest
import user_code
import ast
import re   
import importlib
import csv
import unittest
import importlib
import geopandas as gpd
from shapely import wkt

class TestTask(unittest.TestCase):
    def setUp(self):
        # Create two polygons
        self.poly1 = gpd.GeoDataFrame({'id': [1]}, geometry=gpd.GeoSeries.from_wkt(['POLYGON((0 0, 2 0, 2 2, 0 2, 0 0))']))
        self.poly2 = gpd.GeoDataFrame({'id': [2]}, geometry=gpd.GeoSeries.from_wkt(['POLYGON((1 1, 3 1, 3 3, 1 3, 1 1))']))

    def test_intersection(self):
        import user_code
        importlib.reload(user_code)
        result = user_code.perform_overlay(self.poly1, self.poly2, "intersection")
        _dynamic_test(
            self,
            result is not None and len(result) == 1 and result.geometry.iloc[0].equals(wkt.loads('POLYGON ((2 1, 2 2, 1 2, 1 1, 2 1))')),
            "Intersection returns correct single polygon of overlap region.",
            f"Expected one intersected polygon covering (1,1)-(2,2), got: {result.geometry if result is not None else None}",
        )

    def test_union(self):
        import user_code
        importlib.reload(user_code)
        result = user_code.perform_overlay(self.poly1, self.poly2, "union")
        _dynamic_test(
            self,
            result is not None and len(result) == 3,
            "Union returns three polygons as expected.",
            f"Expected three polygons in union result, got: {len(result) if result is not None else None}",
        )
        # Check that all three expected polygons are present
        expected_geoms = [
            wkt.loads('POLYGON ((2 1, 2 2, 1 2, 1 1, 2 1))'), # overlap
            wkt.loads('POLYGON ((0 0, 2 0, 2 1, 1 1, 1 2, 0 2, 0 0))'), # left outside
            wkt.loads('POLYGON ((2 1, 3 1, 3 3, 1 3, 1 2, 2 2, 2 1))'), # right outside
        ]
        found = [False, False, False]
        for geom in result.geometry:
            for i, exp in enumerate(expected_geoms):
                if geom.equals(exp):
                    found[i] = True
        _dynamic_test(
            self,
            all(found),
            "Union contains all expected polygons.",
            f"Union missing expected polygons: {found}",
        )

    def test_invalid_operation(self):
        import user_code
        importlib.reload(user_code)
        try:
            user_code.perform_overlay(self.poly1, self.poly2, "symmetric_difference")
        except ValueError as e:
            _dynamic_test(
                self,
                "Invalid operation" in str(e),
                "ValueError raised correctly for invalid operation.",
                f"Expected ValueError with 'Invalid operation', got: {e}",
            )
        else:
            _dynamic_test(
                self,
                False,
                "Should have raised ValueError for invalid operation.",
                "Did not raise ValueError for invalid operation.",
            )

def _dynamic_test(test_case, condition, success_message, failure_message):
    if condition:
        test_case._testMethodName = success_message
        test_case.assertTrue(True, success_message)
    else:
        test_case._testMethodName = failure_message
        test_case.fail(failure_message)

def normalize_text(text):
    text = text.lower()
    text = re.sub(r"\\s{2,}", " ", text)
    text = re.sub(r"\\s*([,:?])\\s*", r"\\1 ", text)
    return text.strip()

def change_var(code: str, var_name: str, value: str):
    tree = ast.parse(code)
    lines = code.splitlines()
    changed = False
    # Collect all assignment nodes to modify
    assign_nodes = [
        (i, node)
        for i, node in enumerate(tree.body)
        if isinstance(node, ast.Assign)
        and any(isinstance(target, ast.Name) and target.id == var_name for target in node.targets)
    ]

    # If nothing to change, return unmodified code
    if not assign_nodes:
        return code

    # Perform replacements for all matching assignments (from last to first to not break line offsets)
    for i, node in reversed(assign_nodes):
        start_line = node.lineno - 1
        line = lines[start_line]
        indent = ' ' * (len(line) - len(line.lstrip()))
        lines[start_line] = f"{indent}{var_name} = {value}"
        next_line = len(lines)
        for next_node in tree.body[i+1:]:
            if hasattr(next_node, 'lineno'):
                next_line = next_node.lineno - 1
                break
        if next_line > start_line + 1:
            lines[start_line+1:next_line] = []
        changed = True

    return '\\n'.join(lines) if changed else code

if __name__ == "__main__":
    unittest.main()


test_main.py

Pythonを使用した地理空間データ解析の基礎を学びます。この中級コースでは、地理空間解析の主要な概念と実践的な応用について紹介し、実際のデータセットとハンズオンプロジェクトに焦点を当てます。Pythonと主要なライブラリを用いて空間データの操作、可視化、解析を行う方法を学び、特にgeopandasの活用に重点を置いています。インタラクティブなグリッド、動画、クイズを取り入れ、学習効果を高める構成となっています。

地理空間データの基本、座標系、Pythonによる空間解析のための主要なツールについて学びます。本セクションでは、空間データセットの取り扱いの基礎を築き、geopandasの詳細、主要なデータ構造、読み書き機能、pandasやmatplotlibとの統合について紹介します。

geopandasや可視化ライブラリを用いた、空間結合、オーバーレイ、集約などの主要な空間解析手法を深く学びます。本セクションでは、空間データセットの結合、比較、要約のための高度なgeopandasメソッドや、空間データ可視化のベストプラクティスについて解説します。

実際の地理空間データセットおよびシナリオについての確かな理解を身につけ、実践的な分析および可視化の課題に焦点を当てます。本セクションでは、プロジェクトベースの学習と地理空間分析技術の統合を重視します。

オーバーレイ操作

解答