آموزش تداخل سنجی راداری (InSAR) برای نقشه فرونشست زمین: مفاهیم، SNAP و محدودیت‌ها

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

تداخل سنجی راداری چیست؟ از فاز تا جابه‌جایی

رادار دیافراگم مصنوعی (SAR) موج مایکروویو می‌فرستد و بازگشت آن را ثبت می‌کند. هر پیکسل تصویر SLC (Single Look Complex) دو مقدار دارد: دامنه (شدت بازگشت) و فاز (موقعیت دقیق موج در لحظه بازگشت). اگر سطح زمین بین دو تصویربرداری به اندازه d در راستای دید ماهواره (LOS) جابه‌جا شود، مسیر رفت و برگشت موج 2d تغییر می‌کند و فاز به اندازه زیر عوض می‌شود:

Δφ = (4π / λ) × d        →        d = λ × Δφ / (4π)

طول موج باند C سنتینل ۱ حدود ۵٫۵۵ سانتی‌متر است؛ پس یک چرخه کامل فاز (۲π) معادل λ/2 یعنی حدود ۲٫۷۷ سانتی‌متر جابه‌جایی LOS است. در تداخل‌نگار، این چرخه‌ها به صورت نوارهای رنگی تکرارشونده (فرینج) دیده می‌شوند. یک کاسه فرونشست با ۵ فرینج یعنی حدود ۱۴ سانتی‌متر جابه‌جایی در مرکز نسبت به حاشیه پایدار. مبانی رادار را در رادار دیافراگم مصنوعی SAR با مثال و مشخصات ماهواره را در ماهواره سنتینل ۱ بخوانید.

اجزای فاز تداخل‌نگار

فاز تداخل‌نگار فقط جابه‌جایی نیست؛ مجموع چند جزء است و هنر InSAR جدا کردن آن‌هاست:

DInSAR، SBAS و PS: کدام روش؟

این آموزش روی DInSAR در SNAP تمرکز دارد، چون پایه فهم همه روش‌هاست و برای فرونشست‌های سریع دشت‌های ایران (چند سانتی‌متر در ماه) نتیجه قابل مشاهده‌ای می‌دهد. برای گزارش نرخ سالانه در پایان‌نامه یا گزارش رسمی، حتماً به SBAS یا PS بروید.

انتخاب و دانلود تصاویر سنتینل ۱

1. محصول: Sentinel-1 IW SLC (نه GRD؛ GRD فاز ندارد و برای InSAR قابل استفاده نیست).
2. مدار یکسان: هر دو تصویر باید از یک Relative Orbit و یک جهت (Ascending یا Descending) باشند.
3. فاصله زمانی: هرچه کوتاه‌تر، همدوسی بیشتر. برای دشت‌های کشاورزی ایران ۱۲ تا ۲۴ روز معمولاً بهترین است؛ فاصله‌های طولانی‌تر در زمستان و بهار به دلیل تغییر رطوبت و کشت، همدوسی را از بین می‌برد. از دی ۱۴۰۰ با از کار افتادن سنتینل ۱B، بازدید مجدد ۱۲ روزه شد و با فعال شدن سنتینل ۱C از ۱۴۰۴ دوباره امکان جفت‌های ۶ روزه فراهم شده است.
4. خط مبنای عمودی (Perpendicular Baseline): کمتر از ۱۵۰ متر ترجیحاً؛ ابزار InSAR Stack Overview در SNAP آن را نشان می‌دهد.
5. فصل: جفت‌های تابستانه و خشک در ایران همدوسی بهتری دارند؛ از تصاویر بعد از بارندگی و برف پرهیز کنید.

دانلود از Copernicus Data Space Ecosystem یا ASF Vertex انجام می‌شود. هر تصویر SLC حدود ۴ تا ۸ گیگابایت است؛ با TOPSAR-Split فقط Subswath و Burstهای منطقه را نگه دارید تا پردازش روی کامپیوتر معمولی با ۱۶ گیگابایت رم ممکن شود.

گام‌به‌گام تداخل سنجی راداری در SNAP

1. TOPSAR-Split: برای هر تصویر Subswath (IW1، IW2 یا IW3) و Burstهای منطقه و پلاریزاسیون VV را انتخاب کنید.
2. Apply Orbit File: مدار دقیق (Precise Orbit) را اعمال کنید؛ برای تصاویر جدید ممکن است چند هفته طول بکشد تا منتشر شود.
3. Back-Geocoding: دو تصویر را با DEM (SRTM 1Sec یا Copernicus 30m با دانلود خودکار) هم‌ثبت کنید. تصویر قدیمی‌تر Master باشد.
4. Enhanced Spectral Diversity: برای هم‌ثبتی دقیق بین Burstها در حالت TOPS ضروری است؛ بدون آن پرش فاز در مرز Burstها می‌بینید.
5. Interferogram Formation: تداخل‌نگار و باند همدوسی (Coherence) ساخته می‌شود؛ گزینه Subtract flat-earth phase فعال باشد.
6. TOPSAR-Deburst: خطوط سیاه بین Burstها را حذف کنید.
7. TopoPhase Removal: فاز توپوگرافی با DEM حذف می‌شود.
8. Multilooking: مثلاً ۴ در برد و ۱ در آزیموت، تا پیکسل تقریباً مربعی ۱۵ تا ۲۰ متری و نویز کمتر داشته باشید.
9. Goldstein Phase Filtering: با توان ۱ و اندازه پنجره ۶۴، فرینج‌ها واضح‌تر می‌شوند.
10. SNAPHU Export: خروجی را برای بازکردن فاز آماده کنید و در پوشه خروجی SNAPHU را اجرا کنید (فرمان آماده در فایل snaphu.conf هست). در نسخه‌های جدید SNAP افزونه snaphu-unwrapping این کار را داخل نرم‌افزار انجام می‌دهد.
11. SNAPHU Import: فاز بازشده را به محصول برگردانید.
12. Phase to Displacement: فاز بازشده به جابه‌جایی LOS به متر تبدیل می‌شود.
13. Range Doppler Terrain Correction: خروجی را زمین‌مرجع کنید (مثلاً WGS84 یا UTM) و به GeoTIFF صادر کنید.

این زنجیره را در Graph Builder ذخیره کنید تا برای جفت‌های بعدی تکرار شود. نرم‌افزار را از صفحه رسمی ESA SNAP بگیرید؛ انجمن STEP بهترین جا برای رفع خطاهای پردازش است. برای آشنایی گسترده‌تر با پردازش تصاویر، آموزش جامع پردازش تصاویر ماهواره‌ای را ببینید.

مثال عددی: فرونشست یک دشت در جنوب‌شرق ایران

فرض کنید جفت تصویری با فاصله ۹۶ روز (تیر تا مهر) از یک دشت کشاورزی پسته‌کاری پردازش کرده‌ایم. در مرکز کاسه فرونشست، پس از کم کردن ناحیه مرجع، جابه‌جایی LOS برابر ۶ سانتی‌متر (دور شدن از ماهواره) است.

• فاز معادل: Δφ = 4π × 0.06 / 0.0555 ≈ 13.6 rad، یعنی حدود ۲٫۲ فرینج.
• تبدیل به قائم: با فرض اینکه جابه‌جایی افقی ناچیز است و زاویه برخورد ۳۹ درجه: dv = dLOS / cos(39°) = 0.06 / 0.777 ≈ 7.7 cm.
• نرخ سالانه تقریبی: 7.7 × 365 / 96 ≈ 29 cm/yr.

این عدد در حد نرخ‌هایی است که برای برخی دشت‌های کرمان گزارش شده، اما فقط از یک جفت به دست آمده است. تا ۱ فرینج (حدود ۳ سانتی‌متر) می‌تواند اثر اتمسفر باشد؛ پس عدم قطعیت این برآورد چند ده درصد است. با ۳۰ تا ۶۰ تصویر یک سال و SBAS، این عدم قطعیت به چند میلی‌متر در سال کاهش می‌یابد.

همدوسی: معیار اعتماد به هر پیکسل

در دشت‌های کشاورزی، همدوسی مزارع آبی اغلب پایین است و نقشه فرونشست «سوراخ‌سوراخ» می‌شود. روستاها، جاده‌ها و باغ‌ها نقاط قابل اعتماد را فراهم می‌کنند و روش PS دقیقاً از همین نقاط پایدار استفاده می‌کند.

تحلیل خروجی InSAR در GIS

1. ماسک همدوسی: در ArcGIS Pro با SetNull("coh" < 0.3, "disp") پیکسل‌های نامعتبر را حذف کنید.
2. تبدیل واحد و علامت: جابه‌جایی را به سانتی‌متر ببرید و قرارداد علامت (منفی = فرونشست) را در راهنما بنویسید. در SNAP مقدار مثبت معمولاً یعنی نزدیک شدن به ماهواره.
3. پرکردن خلأها با احتیاط: برای نقشه نمایشی می‌توانید با درون‌یابی (IDW یا کریجینگ) خلأهای کوچک را پر کنید، اما در محاسبات از پیکسل‌های اصلی استفاده کنید.
4. روی‌هم‌گذاری: نقشه فرونشست را با چاه‌های بهره‌برداری، افت سطح آب زیرزمینی، ضخامت رسوبات آبرفتی، کاربری اراضی و زیرساخت‌ها (راه‌آهن، خطوط لوله) مقایسه کنید.
5. پروفیل: با ابزار Profile در ArcGIS Pro یک برش عرضی از مرکز کاسه بگیرید تا گرادیان فرونشست مشخص شود؛ گرادیان شدید محل احتمالی شکاف‌های زمین است.
6. اعتبارسنجی: با ایستگاه‌های GNSS دائم یا ترازیابی دقیق سازمان نقشه‌برداری کشور مقایسه کنید.

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

تفسیر فرونشست برای دشت‌های ایران

نقشه جابه‌جایی زمانی ارزش مدیریتی پیدا می‌کند که به علت‌ها و پیامدها وصل شود. چند الگوی تفسیری که در دشت‌های ایران بارها دیده شده است:

• کاسه‌های فرونشست منطبق بر تراکم چاه‌ها: بیشینه فرونشست معمولاً در بخش میانی دشت با ضخیم‌ترین رسوبات ریزدانه (رسی و سیلتی) رخ می‌دهد، نه لزوماً جایی که بیشترین افت آب را دارد. لایه ضخامت آبرفت از مطالعات ژئوفیزیک وزارت نیرو در این تفسیر بسیار کمک می‌کند.
• ارتباط با کاربری: روی‌هم‌گذاری با نقشه کاربری اراضی معمولاً نشان می‌دهد کاسه‌ها زیر اراضی کشاورزی آبی و باغ‌ها متمرکزند؛ و آشکارسازی تغییرات کاربری نشان می‌دهد گسترش کشت آبی پیش از شروع فرونشست بوده است یا نه.
• فرونشست فصلی و دائمی: در سری زمانی، نوسان فصلی (بالا آمدن جزئی در زمستان و پایین رفتن در تابستان) در کنار روند نزولی دائمی دیده می‌شود. بخش دائمی، تراکم غیرقابل برگشت آبخوان و از دست رفتن ظرفیت ذخیره است.
• شکاف‌ها و زیرساخت‌ها: بیشترین خطر در محل‌های تغییر شدید نرخ فرونشست است، نه لزوماً در مرکز کاسه؛ این مناطق را برای راه‌آهن، خطوط انتقال آب و گاز و ساختمان‌ها اولویت پایش قرار دهید.

برای تحلیل‌های مکانی تکمیلی مثل ارتباط فرونشست با شبکه آبراهه و حوضه‌ها، کاربرد GIS در هیدرولوژی و برای یادگیری ساختارمند سنجش از دور، دوره جامع سنجش از دور و فهرست دوره‌های GISLAND را ببینید. تحلیل نهایی و ترکیب لایه‌ها در ArcGIS Pro را در دوره ArcGIS Pro پیشرفته تمرین می‌کنید.

محدودیت‌های تداخل سنجی راداری که باید بدانید

• فقط یک مؤلفه: InSAR فقط جابه‌جایی در راستای LOS را می‌سنجد. تفکیک قائم و شرقی-غربی به ترکیب مدار صعودی و نزولی نیاز دارد و جابه‌جایی شمالی-جنوبی تقریباً دیده نمی‌شود.
• حداکثر گرادیان: اگر جابه‌جایی بین دو پیکسل مجاور بیش از حدود یک‌چهارم طول موج (حدود ۱٫۴ سانتی‌متر) باشد، فرینج‌ها در هم می‌روند و Unwrapping شکست می‌خورد. در فرونشست‌های خیلی سریع، از جفت‌های کوتاه‌تر یا باند L استفاده کنید.
• خطای Unwrapping: SNAPHU گاهی به اندازه مضربی از ۲π (یعنی ۲٫۸ سانتی‌متر) جهش می‌کند؛ ناپیوستگی‌های ناگهانی در نقشه جابه‌جایی را مشکوک بدانید.
• اثر اتمسفر: بخار آب می‌تواند الگوهایی شبیه فرونشست بسازد، به‌ویژه در کوهپایه‌ها که با ارتفاع همبسته است. الگویی که در یک جفت هست و در جفت بعدی نیست، احتمالاً اتمسفری است.
• از دست رفتن همدوسی: پوشش گیاهی، کشت، برف و آب سیگنال را از بین می‌برند.
• نسبی بودن: همه مقادیر نسبت به نقطه مرجع‌اند.
• GEE کمکی نمی‌کند: مجموعه سنتینل ۱ در گوگل ارث انجین GRD است و فاز ندارد؛ برای InSAR باید SLC را جداگانه پردازش کنید.

محصولات آماده: قبل از پردازش ببینید

پروژه COMET-LiCSAR دانشگاه لیدز تداخل‌نگارهای سنتینل ۱ را برای بخش بزرگی از مناطق فعال تکتونیکی، از جمله ایران، به صورت رایگان منتشر می‌کند و ابزار LiCSBAS از روی آن‌ها سری زمانی SBAS می‌سازد. برای شروع یک پایان‌نامه، بررسی این محصولات پیش از دانلود ده‌ها گیگابایت تصویر، هفته‌ها زمان صرفه‌جویی می‌کند.

خطاهای رایج در پردازش InSAR و راه حل

پرسش‌های متداول درباره تداخل سنجی راداری

مطالب مرتبط