🗺️ برنامه اجرایی و تشریحی «مرحله سوم: محاسبات آمار فضایی و ساخت آتلاس نقشه‌های حرارتی»

(از گام ۳: تلفیق داده‌ها، تحلیل مکان‌مند و انطباق بر نظام‌های ۹‌گانه)

در مراحل اول و دوم، داده‌های خام شهروندان استخراج، گمنام‌سازی، نرمال‌سازی و برچسب‌گذاری موضوعی شده و بار احساسی متون باز استخراج گردید. اما نقاط و خطوط پراکنده روی نقشه به‌تنهایی برای تصمیم‌گیری کالبدی و تدوین دستور کار توسعه فضایی شهر کافی نیستند.

مرحله سوم از گام سوم، بخش محاسبات ریاضی و آمار فضایی (Spatial Statistics) است. در این مرحله، با بکارگیری الگوریتم‌های پیشرفته تحلیل فضایی در GIS، تجمع پراکنده نظرات به کانون‌های داغ (Hotspots) و سرد (Coldspots) معنادار تبدیل می‌شوند. خروجی این مرحله، «آتلاس نقشه‌های حرارتی شهر» تفکیک‌شده بر نظام‌های ۹‌گانه مطالعاتی است که زون‌های حاد مداخلات کالبدی و اصلاحات شهری را با دقت ریاضی اثبات می‌کند.

                                    ┌─────────────────────────────────────────────────────────┐
                                    │    خط فرآیند ۵‌گانه محاسبات آمار فضایی و ساخت آتلاس      │
                                    └──────────────────────────┬──────────────────────────────┘
                                                               │
        ┌──────────────────────────────┬───────────────────────┼───────────────────────┬──────────────────────────────┐
        ▼                              ▼                       ▼                       ▼                              ▼
┌───────────────┐              ┌───────────────┐       ┌───────────────┐       ┌───────────────┐              ┌───────────────┐
│ ۱. تخمین      │              │ ۲. تحلیل      │       │ ۳. اعمال      │       │ ۴. خوشه‌بندی  │              │ ۵. رندرینگ و  │
│  تراکم هسته‌ای │              │  کانون‌های داغ │       │   ضرایب وزنی  │       │   DBSCAN      │              │   تولید آتلاس │
│  (KDE)        │              │ (Getis-Ord)   │       │  دموغرافیک    │       │  حذف نویزها   │              │   نقشه‌های حرارتی│
└───────┬───────┘              └───────┬───────┘       └───────┬───────┘       └───────┬───────┘              └───────┬───────┘
  • پهنای باند ۳۰۰ تا ۵۰۰m     • محاسبات Z-Score       • وزن ۱.۵ گروه‌های بی‌صدا • جداسازی داده‌های پرت   • طیف رنگی استاندارد
  • تابع Epanechnikov          • سطح اطمینان ۹۵٪ و ۹۹٪ • تعدیل بر تراکم جمعیت  • تعیین مرز دقیق خوشه‌ها • تفکیک بر ۹ نظام
  • رسترسازی شعاعی             • تفکیک Hot/Cold        • موازنه آراء محلات     • Concave Hulls                • خروجی GeoTIFF/GeoJSON

🛠️ بخش اول: تشریح تفصیلی گام‌به‌گام الگوریتم‌ها و محاسبات آمار فضایی

1️⃣ گام ۱: تخمین تراکم هسته‌ای (Kernel Density Estimation - KDE)

هدف: تبدیل نقاط پراکنده عارضه‌ها، خطرات و مسائل ارسالی شهروندان به یک لایه رستری پیوسته (Continuous Surface) از میزان تراکم مشکلات در سطح شهر.
اقدامات عملیاتی:
1. تعیین پهنای باند شعاع جستجو (Bandwidth / Search Radius): تنظیم شعاع اثر بر اساس مقیاس فضایی محلات (معمولاً ۳۰۰ تا ۵۰۰ متر متناسب با حوزه پیاده‌روی محله‌محور).
2. انتخاب تابع هسته (Kernel Function): بکارگیری تابع Epanechnikov یا Quartic برای کاهش وزن اثر عارضه‌ها با افزایش فاصله از مرکز نقطه.
3. محاسبه تراکم فضایی رستری: ایجاد سلول‌های رستری ۱۰×۱۰ متری که در آن ارزش هر سلول (Cell Value) نشان‌دهنده تعداد عارضه‌های ثبت‌شده در شعاع اطراف آن است.
4. تفکیک لایه‌های رستری بر اساس موضوع: اجرای مجزای الگوریتم KDE برای هر یک از موضوعات (مثلاً تراکم عارضه‌های ترافیکی، تراکم آلودگی‌ها، تراکم کمبود خدمات).

2️⃣ گام ۲: تحلیل کانون‌های داغ و سرد آماری (Getis-Ord Gi* Spatial Statistics)

هدف: متمایز ساختن تراکم‌های تصادفی از کانون‌های داغ واقعی و معنادار به لحاظ آماری (با حذف فاکتور شانس).
اقدامات عملیاتی:
1. محاسبه شاخص Gi* برای هر سلول: مقایسه تراکم محلی عارضه‌ها نسبت به کل سطح شهر با فرمول استاندارد Getis-Ord: \[\text{Gi}^* = \frac{\sum{j=1}^{n} w{ij}xj - \bar{X}\sum{j=1}^{n} w{ij}}{S \sqrt{\frac{n\sum{j=1}^{n} w{ij}^2 - (\sum{j=1}^{n} w_{ij})^2}{n-1}}}\]
2. محاسبه شاخص‌های Z-Score و P-Value:
  - Z-Score > +1.96 (P-Value < 0.05): کانون داغ معنادار آماری با اطمینان ۹۵٪ (Hotspot - تراکم بسیار حاد مشکلات یا احساسات منفی).
  - Z-Score > +2.58 (P-Value < 0.01): کانون داغ بحرانی با اطمینان ۹۹٪ (Critical Hotspot).
  - Z-Score < -1.96: کانون سرد معنادار آماری (Coldspot - پهنه‌های فاقد مسئله یا واجد احساس مثبت و رضایت).

3️⃣ گام ۳: اعمال ضرایب وزنی حمایتی و تعدیل جمعیتی (Demographic & Population Weighting)

هدف: جلوگیری از خطای تمرکز مطلق داده‌ها در محلات پرجمعیت و اعمال عدالت توزیعی در نقشه‌ها.
اقدامات عملیاتی:
1. تعدیل بر اساس نرخ جمعیت محله (Population-Normalized Density): تقسیم نمره تراکم عارضه‌ها بر تراکم جمعیت محله تا محلات کم‌جمعیت‌تر اما واجد مشکلات حاد کالبدی ندیده گرفته نشوند.
2. اعمال ضریب وزنی حمایتی به گروه‌های بی‌صدا: ضرب کردن نمره شدت عارضه‌ها و احساسات منفی ارسالی از سوی اقشار محروم (زنان سرپرست خانوار، معلولان، سالمندان و ساکنان حاشیه‌ای) در ضریب ۱.۵ تا ۲.۰.

4️⃣ گام ۴: خوشه‌بندی فضایی مبتنی بر تراکم و حذف نویزها (DBSCAN / HDBSCAN)

هدف: حذف داده‌های پرت (Outliers) و مرزگذاری دقیق جغرافیایی پیرامون کانون‌های بحران.
اقدامات عملیاتی:
1. تنظیم پارامترهای الگوریتم DBSCAN:
  - Eps (شعاع همسایگی): ۲۰۰ متر.
  - MinPts (حداقل نقاط تشکیل خوشه): ۵ نقطه.
2. تفکیک داده‌های پرت: حذف نقاط منفردی که خارج از خوشه‌های اصلی قرار دارند به عنوان نویز یا نظر غیرتکرارپذیر.
3. مرزگذاری محدوده خوشه‌ها (Cluster Boundary Generation): رسم چندضلعی‌های بسته محصورکننده (Concave Hulls) حول خوشه‌های شناسایی‌شده برای تعیین دقیق پهنه مداخله کالبدی در طرح جامع.

5️⃣ گام ۵: رندرینگ استاندارد و ساخت آتلاس نقشه‌های حرارتی (Heatmap Atlas Production)

هدف: تدوین مجلد و لایه‌های استاندارد بصری جهت درج در گزارش تشخیص و ارائه به کارفرمایان و شورای راهبری.
اقدامات عملیاتی:
1. استانداردسازی طیف رنگی (Color Ramp Standard):
  - طیف قرمز/بنفش (Red/Purple): تراکم حاد مشکلات (Z-Score > 2.58).
  - طیف نارنجی/زرد (Orange/Yellow): تراکم متوسط مشکلات محله‌ای.
  - طیف سبز/آبی (Green/Blue): پهنه‌های واجد کیفیت مطلوب، امنیت فضایی و سرمایه اجتماعی.
2. تفکیک آتلاس بر اساس نظام‌های ۹‌گانه: خروجی‌گیری مجزا در قالب ۹ شیت نقشه حرارتی موضوعی (شامل: ۱. حرارتی آلودگی‌ها/محیط‌زیست، ۲. حرارتی فرسودگی مسکن، ۳. حرارتی معابر ناامن ترافیکی، ۴. حرارتی کمبود خدمات محله‌ای، و...).
3. خروجی‌گیری چندگانه: تولید فایلهای رستری سنگین (GeoTIFF) برای پردازش‌های بعدی و فایلهای سبک برداری (GeoJSON) برای پیشخوان برخط شهری.

📅 بخش دوم: برنامه اجرایی، زمان‌بندی (SOP) و جدول نقشه‌ها (RACI)

این مرحله ظرف مدت ۷ روز کاری طبق جدول زمان‌بندی و تفکیک مسئولیت‌های زیر اجرا می‌شود:

┌─────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│                             جدول زمان‌بندی و مراحل SOP محاسبات آمار فضایی (۷ روز)                        │
├──────────────────────────────┬──────────────────────────────┬───────────────────────────┬───────────────┤
│ گام A (روز ۱ و ۲)            │ گام B (روز ۳ و ۴)            │ گام C (روز ۵)             │گام D (روز ۶ و ۷)│
├──────────────────────────────┼──────────────────────────────┼───────────────────────────┼───────────────┤
│ • اجرای الگوریتم‌های KDE     │ • محاسبات Getis-Ord Gi*      │ • اجرای خوشه‌بندی DBSCAN  │ • رندرینگ طیف‌های رنگی│
│ • تنظیم پهنای باند ۳۰۰-۵۰۰m  │ • محاسبه Z-Score و P-Value   │ • رسم مرز Concave Hulls   │ • خروجی‌گیری آتلاس ۹ نظام│
└──────────────────────────────┴──────────────────────────────┴───────────────────────────┴───────────────┘

ماتریس مسئولیت‌های اجرایی (RACI Matrix):

زیرفرآیند عملیاتی مرحله سوم	متخصص ارشد GIS	تحلیل‌گر آمار فضایی	سرپرستان ۹ نظام	ناظر راه و شهرسازی
۱. تنظیم پارامترها و اجرای محاسبات رستری KDE	A / R	C	C	I
*۲. محاسبات Getis-Ord Gi و تعیین Z-Score**	C	A / R	I	I
۳. اعمال ضرایب وزنی حمایتی و تعدیل جمعیتی	C	A / R	C	I
۴. خوشه‌بندی DBSCAN و استخراج مرز خوشه‌ها	A / R	C	I	I
۵. رندرینگ و تولید آتلاس نقشه‌های حرارتی ۹ نظام	A / R	C	C	A

🔍 بخش سوم: سنجه‌ها و شاخص‌های کنترل کیفیت محاسبات آمار فضایی (Quality KPIs)

ناظر فنی اداره کل راه و شهرسازی استان قبل از ابلاغ آتلاس، شاخص‌های زیر را ممیزی و تأیید می‌کند:

سطح اطمینان آماری (Statistical Confidence Guard): تمامی کانون‌های داغ درج‌شده در آتلاس باید واجد Z-Score بالای ۱.۹۶ (سطح اطمینان ۹۵٪ به بالا) باشند؛ تراکم‌های ضعیف و تصادفی نباید به عنوان کانون داغ نمایش داده شوند.
شاخص تعدیل جمعیتی (Population Bias Control): صحت عدم سوگیری نقشه‌ها به نفع محلات پرجمعیت باید از طریق تطبیق لایه KDE تعدیل‌شده با آمار جمعیت بلوک‌ها اثبات شود.
قابلیت انطباق بر لایه‌های کالبدی: مرز خوشه‌های استخراج‌شده (Concave Hulls) باید به‌طور کامل بر معابر، بلوک‌های مسکونی یا پهنه‌های کاربری GIS قابل نگاشت و انطباق باشد.

💾 بخش چهارم: ساختار داده‌ای استاندارد خروجی کانون‌های داغ (Hotspot GeoJSON Schema)

خروجی محاسبات آمار فضایی جهت ذخیره در Geodatabase و نمایش روی پیشخوان برخط با فرمت ساختاریافته زیر ذخیره می‌شود:

{
  "type": "FeatureCollection",
  "features": [
    {
      "type": "Feature",
      "geometry": {
        "type": "Polygon",
        "coordinates": [[[45.7580, 36.7640], [45.7600, 36.7640], [45.7600, 36.7660], [45.7580, 36.7660], [45.7580, 36.7640]]]
      },
      "properties": {
        "HotspotID": "HOT-MHB-1405-0921",
        "TargetSystem": "System6_Infrastructure_and_Traffic",
        "KdeDensityValue": 84.2,
        "GetisOrd_ZScore": 2.84,
        "GetisOrd_PValue": 0.0045,
        "ConfidenceLevelPercent": 99,
        "ClusterClassification": "Critical_Traffic_and_Safety_Hotspot",
        "DemographicWeightFactor": 1.5,
        "AffectedPopulationEstimate": 3400,
        "RecommendedInterventionType": "Immediate_Infrastructure_Redesign",
        "SpatialBoundMethod": "DBSCAN_Concave_Hull"
      }
    }
  ]
}

💡 نتیجه‌گیری عملیاتی: با پیاده‌سازی این مرحله، نظرات و عارضه‌های پراکنده شهروندان به نقشه‌های حرارتی فوق‌العاده دقیق، واجد پشتوانه ریاضی و اثبات‌پذیر به لحاظ آماری تبدیل می‌شوند که کانون‌های واقعی بحران شهر را شفاف ساخته و پایه اصلی تدوین دستور کار توسعه فضایی در گزارش تشخیص خواهند بود.