GHS - حجم ساخته‌شده (Built Volume)

معرفی

لایه حجم ساخته‌شده GHS برآوردی از حجم ساختمان‌ها و زیرساخت‌ها ارائه می‌دهد. این داده‌ها برای مطالعات تراکم شهری و برنامه‌ریزی شهری بسیار مفید هستند.

ویژگی‌ها

باند

نصب پیش‌نیازها

pip install rasterio numpy

دسترسی مستقیم با پایتون

خواندن فایل GeoTIFF محلی

import rasterio
import numpy as np
from pathlib import Path

GHS_ROOT = Path("data/GHS")
BUILT_V_DIR = GHS_ROOT / "BUILT_V"

# لیست فایل‌ها
tif_files = sorted(BUILT_V_DIR.glob("GHS_*_E*_*.tif"))
print(f"تعداد فایل‌ها: {len(tif_files)}")

# خواندن یک فایل
if tif_files:
    with rasterio.open(tif_files[0]) as src:
        built_v = src.read(1)
        print(f"ابعاد: {src.width} x {src.height}")
        print(f"حجم میانگین ساخته‌شده: {np.mean(built_v):.1f} m³")

محاسبه حجم کل ساختمان‌ها

import rasterio
import numpy as np
from pathlib import Path

GHS_ROOT = Path("data/GHS")
BUILT_V_DIR = GHS_ROOT / "BUILT_V"

files_2020 = list(BUILT_V_DIR.glob("GHS_*_E2020_*.tif"))

if files_2020:
    with rasterio.open(files_2020[0]) as src:
        built_v = src.read(1)
        pixel_area = src.res[0] * src.res[1]  # مساحت هر پیکسل (m²)

        # حجم کل ساختمان‌ها
        total_volume = np.sum(built_v) * pixel_area
        print(f"حجم کل ساختمان‌ها: {total_volume:,.0f} m³")
        print(f"حجم میانگین در هر پیکسل: {np.mean(built_v):.1f} m³")

مقایسه تغییرات حجم ساخت‌وساز

import rasterio
import numpy as np
from pathlib import Path

GHS_ROOT = Path("data/GHS")
BUILT_V_DIR = GHS_ROOT / "BUILT_V"

files_1990 = list(BUILT_V_DIR.glob("GHS_*_E1990_*.tif"))
files_2020 = list(BUILT_V_DIR.glob("GHS_*_E2020_*.tif"))

if files_1990 and files_2020:
    with rasterio.open(files_1990[0]) as src:
        built_v_1990 = src.read(1)

    with rasterio.open(files_2020[0]) as src:
        built_v_2020 = src.read(1)

    # تغییرات حجم
    change = built_v_2020 - built_v_1990
    print(f"تغییرات حجم میانگین: {np.mean(change):.1f} m³")
    print(f"مناطق با افزایش حجم: {np.sum(change > 0)} پیکسل")
    print(f"مناطق با کاهش حجم: {np.sum(change < 0)} پیکسل")

لینک‌های مفید

• GHSL Data Portal