مبانی SAR و InSAR؛ از سیگنال راداری تا جابهجایی قابل اعتماد
راهنمای علمی و کاربردی دامنه و فاز، هندسه دید، SLC و GRD، همدوسی، بازکردن فاز، جابهجایی LOS و نقش بازتابنده راداری در پایش زمین و سازه.
فهرست این صفحه
۱. SAR چگونه تصویر میسازد؟
رادار دهانه مصنوعی (Synthetic Aperture Radar) یک سنجنده فعال است: پالس مایکروویو را به زمین میفرستد و بازتاب آن را ثبت میکند. زمان رفتوبرگشت موج، اطلاعات فاصله را فراهم میکند؛ پردازش همدوس مشاهدات هنگام حرکت ماهواره، دهانه مصنوعی را میسازد و تفکیک در امتداد مسیر پرواز را بهبود میدهد. بنابراین SAR صرفاً یک دوربین با نور متفاوت نیست.
تصویربرداری به نور خورشید وابسته نیست و در باندهای رایج ماهوارهای معمولاً از میان ابر انجام میشود؛ بااینحال، جو، بارش و تغییر رطوبت همچنان میتوانند بر اندازهگیری اثر بگذارند. روشنایی تصویر تابع زبری نسبت به طول موج، خواص دیالکتریک و رطوبت، هندسه سطح، قطبش و زاویه برخورد است. روشنتر بودن لزوماً به معنی ارتفاع بیشتر یا حرکت بیشتر نیست.
۲. دامنه، شدت و فاز؛ سه مفهوم متفاوت
هر نمونه مختلط SLC از دو مؤلفه حقیقی I و موهومی Q تشکیل میشود. دامنه A اندازه سیگنال، شدت کمیتی متناسب با مربع دامنه، و فاز موقعیت زاویهای سیگنال مختلط است. نماد j واحد موهومی است. فاز یک تصویر بهتنهایی نقشه ارتفاع یا جابهجایی نیست.
برای سطحی با پراکنشگرهای متعدد، جمع همدوس بازتابها الگویی دانهدانه به نام اسپکل ایجاد میکند. چندنگری و فیلترگذاری میتوانند نوسان شدت را کاهش دهند، اما معمولاً با تغییر تفکیک مؤثر همراهاند. در تحلیل هدف نقطهای باید اثر پردازش بر قله و پاسخ ضربه کنترل شود.
۳. هندسه دید؛ برد، آزیموت و اعوجاجها
برد مایل (Slant range) فاصله سنجنده تا هدف و آزیموت راستای حرکت ماهواره است. برد زمینی تصویر فاصله روی سطح مرجع است. زاویه برخورد نسبت به عمود سطح مرجع تعریف میشود؛ زاویه برخورد محلی با شیب زمین تغییر میکند. اختلاط این دو زاویه در تفسیر شیب یا پسپراکنش خطا ایجاد میکند.
دید جانبی باعث میشود تصویر SAR با نقشه قائم تفاوت داشته باشد. تصحیح هندسی با DEM برای مکانگذاری لازم است، اما اطلاعات ازدسترفته در سایه راداری یا اختلاط اهداف در رویهمافتادگی را بازسازی نمیکند.
| پدیده | معنی و پیامد |
|---|---|
| کوتاهشدگی (Foreshortening) | دامنه رو به رادار در تصویر فشرده میشود؛ فاصله تصویری با طول روی زمین برابر نیست. |
| رویهمافتادگی (Layover) | بازتاب بخش بالاتر و پایینتر میتواند با ترتیب معکوس یا در سلولهای مشترک ثبت شود؛ جداسازی هدف دشوار میشود. |
| سایه راداری (Shadow) | بخشی از زمین در دید مستقیم موج نیست؛ نبود سیگنال به معنی نبود تغییرشکل نیست. |
| تفکیک در برابر فاصله پیکسل | نمونهبرداری ریزتر یا بزرگنمایی تصویر الزاماً دو هدف نزدیک را از هم تفکیک نمیکند. |
۴. انتخاب داده؛ SLC یا GRD؟
برای تداخلسنجی متعارف Sentinel-1 باید اطلاعات مختلط حفظ شده باشد. محصولات SLC دامنه و فاز را در هندسه برد مایل نگه میدارند. در GRD آشکارسازی و چندنگری انجام شده و فاز لازم برای تشکیل تداخلنگاشت در دسترس نیست. از یک فایل GRD نمیتوان فاز SLC را بازیابی کرد. [۱]
Sentinel-1 در باند C با فرکانس مرکزی حدود ۵٫۴۰۵ گیگاهرتز کار میکند؛ طول موج آن تقریباً ۵٫۵۵ سانتیمتر است. حالت IW از روش TOPS و سه زیرنوار استفاده میکند. برای هر محل باید پوشش Burst و تاریخهای واقعاً موجود بررسی شود؛ فاصله زمانی ثابت را بدون کنترل آرشیو فرض نکنید. [۲]
| ویژگی | SLC | GRD |
|---|---|---|
| اطلاعات اصلی | نمونه مختلط؛ دامنه و فاز | شدت آشکارسازیشده؛ بدون فاز تداخلسنجی |
| هندسه محصول | برد مایل و آزیموت | برد زمینی روی بیضوی؛ لزوماً تصحیح توپوگرافی نشده |
| کاربرد این طرح | استخراج پاسخ هدف و تحلیل فاز | بازبینی زمینه و شدت؛ نه جایگزین SLC |
| نکته پردازشی | حفظ فاز و هممرجعسازی دقیق ضروری است | کالیبراسیون و تصحیح هندسی متناسب با کاربرد لازم است |
۵. سازگاری برداشتها و فراداده
هر جفت تصویر صرفاً به دلیل پوشش یک منطقه برای InSAR مناسب نیست. هندسه و قطبش باید سازگار باشند. مجموعه زمانی یا Stack معمولاً برای یک مدار نسبی و جهت گذر مشخص ساخته میشود؛ دادههای صعودی و نزولی جداگانه پردازش و در مرحله تفسیر ترکیب میشوند.
خط مبنای زمانی فاصله زمانی دو برداشت و خط مبنای عمود مؤلفه جدایی مدارها عمود بر دید است. فاصله زمانی بیشتر میتواند تغییر پراکنشگرها را افزایش دهد؛ خط مبنای عمود بزرگتر حساسیت به توپوگرافی و برخی افتهای همدوسی را تغییر میدهد. آستانه مناسب به هدف و هندسه وابسته است.
- ثبت شناسه محصول، زمان UTC، مدار نسبی و جهت گذر؛ فراداده را کنار تصویر نگه دارید.
- کنترل حالت تصویربرداری، قطبش یکسان، زیرنوار و Burst مشترک حاوی هدف.
- ثبت مدار مورد استفاده، DEM و مبنای ارتفاعی آن، نسخه پردازشگر و تنظیمات.
- ثبت فاصله نمونهها و تفکیک مکانی بهصورت جداگانه؛ همچنین کالیبراسیون و اطلاعات نویز.
۶. تداخلنگاشت چگونه تشکیل میشود؟
پس از هممرجعسازی، نمونههای دو تصویر باید به یک ناحیه از زمین مربوط باشند. ضرب تصویر دوم در مزدوج مختلط تصویر اول، تداخلنگاشت میسازد. فاز این حاصلضرب اختلاف فاز دو برداشت است. در TOPS، دقت هممرجعسازی آزیموت و کنترل ناپیوستگی در مرز Burst اهمیت ویژه دارد. [۳]
در روابط این صفحه، تصویر ۱ برداشت نخست و تصویر ۲ برداشت بعدی است؛ قرارداد حاصلضرب صریحاً بهصورت زیر انتخاب شده است. تغییر ترتیب ضرب، علامت اختلاف فاز را برعکس میکند. رنگهای تداخلنگاشت نمایش چرخهای فازند و بدون راهنما و پردازش، جهت حرکت را نشان نمیدهند.
۷. چه چیزهایی در اختلاف فاز حضور دارند؟
اختلاف فاز مشاهدهشده فقط ناشی از حرکت زمین نیست. برای رسیدن به جابهجایی باید سهم هندسه زمین، توپوگرافی، جو، خطای مدار و نویز بررسی شود. رابطه زیر بودجه مفهومی فاز پیش از پیچیدهشدن است؛ در خروجی واقعی هر پردازشگر، بعضی مؤلفهها ممکن است از قبل حذف یا در هم ادغام شده باشند.
| مؤلفه | تفسیر و کنترل |
|---|---|
| زمین مرجع و توپوگرافی | با هندسه مدار و DEM مدل میشوند؛ خطای DEM میتواند فاز باقیمانده وابسته به خط مبنا بسازد. |
| حرکت (deformation) | تغییر فاصله هدف در راستای دید؛ مؤلفه مورد نظر پایش است. |
| جو | تغییر تأخیر تروپوسفر و، بسته به شرایط، یونوسفر؛ میتواند شبیه حرکت دیده شود. |
| مدار و پردازش | مدار نامناسب، خطای هممرجعسازی یا مدل هندسی؛ باید با کنترل مستقل و باقیماندهها بررسی شود. |
| نویز و تغییر پراکنش | نویز حرارتی، تغییر رطوبت و پوشش، یا تغییر مرکز فاز؛ صرفاً با تبدیل واحد حذف نمیشوند. |
۸. همدوسی؛ شاخص کیفیت، نه گواه بیحرکتی
همدوسی (Coherence) میزان ارتباط مختلط دو مشاهده را اندازه میگیرد و مقدار برآوردشده آن بین صفر و یک است. براکتها در رابطه زیر نشاندهنده میانگینگیری روی نمونههای مناسباند. اندازه پنجره و ناهمگنی سطح بر برآورد اثر دارند؛ همدوسی یک عدد مستقل از روش محاسبه نیست.
گیاهان، آب، تغییر رطوبت یا عملیات زمینی ممکن است همدوسی را کاهش دهند. در مقابل، یک هدف متحرک با پاسخ پایدار میتواند همدوسی بالایی داشته باشد. همدوسی بالا نه پایداری زمین را اثبات میکند و نه تأخیر جو را حذف میکند. برای بازتابنده منفرد، SCR و کیفیت سری زمانی هدف نیز باید بررسی شوند.
۹. بازکردن فاز و ابهام چرخهها
فاز ثبتشده معمولاً در بازه منفی π تا مثبت π نمایش داده میشود؛ چند مقدار متفاوت فاز واقعی میتوانند یک رنگ داشته باشند. بازکردن فاز (Phase unwrapping) تعداد صحیح چرخههای ۲π را برآورد میکند. این کار به پیوستگی مکانی یا زمانی، کیفیت داده و مدل مورد استفاده وابسته است.
افت همدوسی، فاصله زیاد نقاط، جهش حرکت و گسست زمانی میتوانند باعث انتخاب چرخه اشتباه شوند. همواربودن نقشه بهتنهایی نشانه درستی نیست؛ اتصال مؤلفهها، باقیماندههای شبکه، سازگاری حلقهها و داده مستقل باید کنترل شوند. بازکردن فاز، تصحیح جو یا حذف توپوگرافی نیست.
۱۰. تبدیل فاز حرکت به جابهجایی؛ با قرارداد روشن
برای روشنبودن علامت، فرض میکنیم فاز انتشار سیگنال برابر منفی ۴πR تقسیم بر λ باشد، تداخلنگاشت با ترتیب بخش ۶ ساخته شود، و d_LOS مثبت به معنی حرکت هدف به سمت ماهواره باشد. در این قرارداد، افزایش فاصله ΔR مثبت است و d_LOS برابر منفی ΔR خواهد بود. λ طول موج و Δφ_def فقط فاز بازشده ناشی از حرکت پس از کنترل مؤلفههای مزاحم است.
مثال آموزشی: با طول موج تقریبی ۵۵٫۵ میلیمتر و فاز حرکت مثبت π رادیان، جابهجایی حدود مثبت ۱۳٫۹ میلیمتر به سمت ماهواره به دست میآید. یک چرخه کامل ۲π معادل ۲۷٫۸ میلیمتر تغییر فاصله است. این اعداد نتیجه پایلوت یا حد دقت دستگاه نیستند.
خطای یک چرخه در بازکردن فاز نیز حدود نصف طول موج خطای جابهجایی ایجاد میکند. در خروجی نرمافزارهای دیگر ممکن است ترتیب ضرب یا تعریف علامت متفاوت باشد؛ قرارداد را از مستندات و یک مثال معلوم کنترل کنید.
۱۱. LOS با نشست قائم برابر نیست
InSAR معمولاً تصویر جابهجایی سهبعدی روی راستای دید را اندازه میگیرد. در رابطه زیر d_E، d_N و d_U مؤلفههای شرق، شمال و بالا و مؤلفههای ℓ متعلق به بردار واحد از زمین به ماهوارهاند. حرکت عمود بر راستای دید ممکن است در این مشاهده دیده نشود.
فقط اگر حرکت افقی قابل صرفنظر باشد، میتوان جابهجایی قائم را از تقسیم LOS بر کسینوس زاویه برخورد θ برآورد کرد. برای مثال آموزشی، LOS برابر منفی ۱۰ میلیمتر در زاویه ۳۵ درجه، با فرض حرکت کاملاً قائم، حدود منفی ۱۲٫۲ میلیمتر جابهجایی قائم یا ۱۲٫۲ میلیمتر نشست میدهد.
ترکیب گذر صعودی و نزولی معمولاً برای تفکیک قائم و شرق–غرب با یک فرض اضافی درباره شمال–جنوب به کار میرود. دو مشاهده مستقل برای حل آزاد سه مؤلفه کافی نیست؛ حساسیت مدارهای نزدیکقطبی به شمال–جنوب نیز معمولاً ضعیفتر است. مقایسه با GNSS باید پس از تصویرکردن بردار GNSS روی LOS و یکسانسازی مرجع انجام شود.
۱۲. DInSAR، تحلیل سری زمانی و مرجع
DInSAR با حذف سهم مدلشده توپوگرافی از تداخلنگاشت، بررسی تغییرشکل بین برداشتها را ممکن میکند. تحلیل سری زمانی از مشاهدات متعدد و یک شبکه زمانی برای برآورد روند، تغییرات فصلی یا رویدادها استفاده میکند. افزایش تعداد تصاویر بهتنهایی تضمین حذف خطا نیست؛ توزیع زمانی و هندسه شبکه اهمیت دارند.
در PS-InSAR تمرکز بر اهدافی با رفتار فازی پایدار است؛ روشهای SBAS از شبکه جفتهای دارای خط مبنای مناسب استفاده میکنند. اینها خانواده روشاند و تعداد تصویر، آستانه و مدل ثابت و یکسانی برای همه پروژهها ندارند. بازتابنده مصنوعی میتواند یک هدف قابل ردیابی به شبکه اضافه کند.
نتیجه معمولاً نسبت به نقطه یا ناحیه مرجع و تاریخ مرجع گزارش میشود. اگر مرجع حرکت کند، آن حرکت به تمام سریها منتقل میشود. یک بازتابنده روشن الزاماً مرجع پایدار نیست. پایداری مرجع باید با شرایط زمینشناختی، سری زمانی و در صورت امکان اندازهگیری مستقل بررسی شود.
۱۳. بازتابنده راداری؛ RCS، SCR و پایداری فاز
کرنر رفلکتور با هندسه و جهتدهی مناسب یک هدف مصنوعی قوی ایجاد میکند. سطح مقطع راداری (RCS) قدرت پراکندگی هدف را به سطح معادل بیان میکند و با مساحت هندسی صفحات یا مقدار یک پیکسل یکسان نیست. RCS خطی بر حسب مترمربع و بیان لگاریتمی آن نسبت به یک مترمربع است؛ نماد dBm² را با dBm توان اشتباه نگیرید.
SCR نسبت توان پاسخ هدف به کلاتر محیط و SNR نسبت سیگنال به نویز است؛ این دو یکسان نیستند. شدت خام یا پسپراکنش کالیبره یک پیکسل نیز جای برآورد RCS هدف نقطهای را نمیگیرد. روش انتگرالگیری پاسخ، کالیبراسیون، زمینه و واحدها باید مشخص باشند. [۵]
افزایش SCR معمولاً اثر کلاتر بر فاز هدف را کاهش میدهد، اما خطای جو، مرجع متحرک یا ناپایداری پایه را حذف نمیکند. تغییر جهت، انبساط حرارتی، آبگرفتگی و تغییر مرکز فاز میتوانند اندازهگیری را تحت تأثیر قرار دهند. بازتابنده باید حرکت کمیت مورد نظر را منتقل کند؛ حرکت اتصال یا پایه لزوماً حرکت سازه اصلی نیست.
۱۴. مکانیابی هدف با اندازهگیری حرکت فرق دارد
درونیابی FFT یا Oversampling نمونهبرداری پاسخ را متراکم میکند تا مکان قله زیرپیکسلی بهتر برآورد شود؛ تفکیک فیزیکی سنجنده یا اطلاعات مستقل جدید ایجاد نمیکند. اختلاف مکان مشاهدهشده و پیشبینیشده در برد و آزیموت، مبنای ارزیابی مکانیابی مانند ALE است.
دقت مکانیابی مطلق و دقت جابهجایی نسبی فازی دو کمیت متفاوتاند. هدفی که بهخوبی مکانیابی شده ممکن است سری زمانی آلوده به جو داشته باشد؛ همچنین اندازهگیری فازی زیرپیکسلی به معنی تفکیک دو هدف با همان فاصله نیست. تصحیح تأخیر جو در مکانیابی را با برآورد صفحه فاز جوی در InSAR یکسان نگیرید.
۱۵. زنجیره پردازش و نقاط کنترل
ترتیب جزئی مراحل با روش و نرمافزار تغییر میکند؛ زنجیره زیر نقشه مفهومی کنترل کیفیت است. اجرای موفق یک فرمان یا تولید نقشه رنگی، تأیید اعتبار اندازهگیری محسوب نمیشود.
- تعریف کمیت مورد اندازهگیری، محل هدف، هندسه دید، مرجع و داده اعتبارسنجی.
- انتخاب SLCهای سازگار؛ کنترل پوشش، قطبش، مدار و فاصله زمانی.
- آمادهسازی مدار و DEM، استخراج محدوده و هممرجعسازی دقیق با کنترل TOPS.
- تشکیل تداخلنگاشت و مدلکردن سهم هندسه و توپوگرافی؛ ارزیابی همدوسی و پاسخ هدف.
- فیلترگذاری متناسب، ساخت شبکه و بازکردن فاز با کنترل چرخهها و گسستها.
- برآورد مؤلفههای مزاحم و سری زمانی با مرجع مکانی و زمانی صریح.
- تبدیل فاز حرکت به LOS، مکانگذاری خروجی و مقایسه با مشاهده مستقل هممرجع.
- ثبت عدمقطعیت، دادههای حذفشده، نسخهها و محدودیت تفسیر همراه نتیجه.
۱۶. چه زمانی میتوان به نتیجه استناد کرد؟
حساسیت فاز به جابهجایی کوچک با دقت نهایی محصول برابر نیست. نسبت طول موج به ۴π فقط ضریب تبدیل فاز به فاصله است. دقت میلیمتری یا سانتیمتری باید برای همان هدف، بازه زمانی و زنجیره پردازش با شواهد ارزیابی شود.
عدمقطعیت باید اثر نویز، جو، خطای DEM، مدل حرکت، انتخاب مرجع و خطای اندازهگیری مستقل را در نظر بگیرد. خطاهای زمانی ممکن است همبسته باشند؛ بنابراین شمار بیشتر تصاویر لزوماً عدمقطعیت را با یک قاعده ساده کاهش نمیدهد.
| برداشت نادرست | تفسیر درست |
|---|---|
| قله روشن یعنی هدف پایدار است | روشنایی، پایداری فاز و پایداری مکانیکی جداگانه ارزیابی میشوند. |
| همدوسی بالا یعنی حرکت صفر | همدوسی ارتباط سیگنالهاست؛ هدف میتواند همدوس و متحرک باشد. |
| هر چرخه رنگی یعنی نشست | چرخه ممکن است از توپوگرافی، جو یا حرکت LOS ناشی شود. |
| LOS مستقیماً نشست قائم است | تبدیل قائم به فرض یا اطلاعات اضافی درباره حرکت افقی نیاز دارد. |
| بازکردن فاز، خطای جو را حذف میکند | رفع ابهام چرخه و جداسازی جو دو مسئله متفاوتاند. |
| نقشه صاف یا نمودار خطی یعنی نتیجه صحیح | سازگاری شبکه، مرجع و اعتبارسنجی مستقل نیز ضروری است. |
۱۷. واژهنامه و مسیر مطالعه
- برای انتخاب شکل و جهتدهی هدف، بخش «طراحی و نصب» را بخوانید.
- برای مراحل اجرایی، بخشهای دریافت داده، پیشپردازش و CR-InSAR را دنبال کنید.
- برای معیارهای پذیرش و درج نتایج واقعی، به اعتبارسنجی و پرونده پایلوت مراجعه کنید.
| اصطلاح | تعریف کاربردی |
|---|---|
| SAR / InSAR | تصویربرداری راداری دهانه مصنوعی / تداخلسنجی تصاویر راداری |
| SLC / GRD | محصول مختلط دارای فاز / محصول آشکارسازیشده |
| Range / Azimuth | راستای برد / راستای مسیر پرواز |
| Track / Stack | مسیر مداری / مجموعه زمانی تصاویر سازگار |
| IW / TOPS / Burst | حالت تصویربرداری پهننوار تداخلسنجی / روش جاروب پرتو / قطعه زمانی تصویربرداری |
| Baseline / DEM | خط مبنای مدارها / مدل رقومی ارتفاع |
| Coherence / Unwrapping | همدوسی / تعیین تعداد صحیح چرخههای فاز |
| PS / DS / CR | پراکنشگر پایدار / پراکنشگر توزیعشده / کرنر رفلکتور |
| RCS / SCR / SNR | سطح مقطع راداری / نسبت سیگنال به کلاتر / نسبت سیگنال به نویز |
| LOS / APS | راستای دید / صفحه فاز جوی |
| ALE / TRF / Epoch | خطای مکانیابی مطلق / چارچوب مرجع زمینی / زمان مرجع مختصات |
۱۸. منابع و دامنه استفاده
منابع زیر برای تفکیک محصول، هندسه برداشت، پردازش و پایش هدف معرفی شدهاند. روابط با قرارداد علامت صریح و مثالهای عددی آموزشی در این صفحه تنظیم شدهاند؛ مثالها داده مشاهدهشده شرکت نیستند.
[۱] مشخصات محصولات Sentinel-1 در SentiWiki؛ [۲] مشخصات مأموریت و سنجنده؛ [۳] پردازش Sentinel-1 و TOPS؛ [۴] راهنمای مبانی InSAR آژانس فضایی اروپا؛ [۵] مستندات پایش بازتابنده GECORIS؛ [۶] معرفی تداخلسنجی NASA. برای تصمیم اجرایی، نسخه نرمافزار و فراداده همان محصول ملاک نهایی است.
برای گفتوگو درباره محل پایلوت، داده راداری یا اندازهگیری مستقل با تیم پژوهش در ارتباط باشید.
ارتباط با تیم پژوهش ←