
Ю. И. Кантемиров
ВВЕДЕÐИЕ
Ð Ð°Ð´Ð°Ñ€Ð½Ð°Ñ Ñъемка выполнÑетÑÑ Ð² ультракоротковолновой (ÑверхвыÑокочаÑтотной) облаÑти радиоволн, подразделÑемой на X-, C-, L- и P-диапазоны (табл. 1). Съемка в каждом из вышеперечиÑленных диапазонов имеет Ñвои плюÑÑ‹ и минуÑÑ‹. Ð”Ð»Ñ Ð·Ð°Ð´Ð°Ñ‡ мониторинга Ñмещений земной поверхноÑти, зданий и Ñооружений по каждой конкретной территории подбираютÑÑ Ð´Ð°Ð½Ð½Ñ‹Ðµ в одном или неÑкольких из Ñтих диапазонов иÑÑ…Ð¾Ð´Ñ Ð¸Ð· типа рельефа, типа раÑтительного покрытиÑ, ожидаемых величин Ñмещений и Ñ‚. д.
Таблица 1. Диапазоны радиоволновой облаÑти Ñлектромагнитного Ñпектра
Входными данными Ð´Ð»Ñ Ð¾Ð±Ñ€Ð°Ð±Ð¾Ñ‚ÐºÐ¸ в Ñпециализированных программных комплекÑах ÑвлÑÑŽÑ‚ÑÑ Ð¸Ð½Ñ‚ÐµÑ€Ñ„ÐµÑ€Ð¾Ð¼ÐµÑ‚Ñ€Ð¸Ñ‡ÐµÑÐºÐ°Ñ Ð¿Ð°Ñ€Ð° (либо цепочка) радарных Ñнимков. Ограничением Ð´Ð»Ñ Ð²Ð¾Ð·Ð¼Ð¾Ð¶Ð½Ð¾Ñти интерферометричеÑкой обработки пары (или цепочки) радарных Ñнимков ÑвлÑÑŽÑ‚ÑÑ Ð¿Ñ€Ð¾ÑтранÑÑ‚Ð²ÐµÐ½Ð½Ð°Ñ Ð¸ Ð²Ñ€ÐµÐ¼ÐµÐ½Ð½Ð°Ñ Ð±Ð°Ð·Ñ‹ [1].
ПроÑтранÑÑ‚Ð²ÐµÐ½Ð½Ð°Ñ Ð±Ð°Ð·Ð° (или Ð±Ð°Ð·Ð¾Ð²Ð°Ñ Ð»Ð¸Ð½Ð¸Ñ) предÑтавлÑет Ñобой раÑÑтоÑние между орбитальными положениÑми радиолокатора при Ñъемке изображений, ÑоÑтавлÑющих интерферометричеÑкую пару (риÑ. 1).

КачеÑтво результатов интерферометричеÑкой обработки напрÑмую завиÑит от величины перпендикулÑрной ÑоÑтавлÑющей базовой линии. С одной Ñтороны, качеÑтво получаемой интерферометричеÑким методом карты Ñмещений земной поверхноÑти возраÑтает Ñ ÑƒÐ¼ÐµÐ½ÑŒÑˆÐµÐ½Ð¸ÐµÐ¼ длины перпендикулÑрной базовой линии. Ð’ Ñлучае базовой линии, равной нулю, интерферограмма, раÑÑÑ‡Ð¸Ñ‚Ð°Ð½Ð½Ð°Ñ Ð¿Ð¾ такой паре Ñнимков, вообще Ñодержит только фазу Ñмещений. С другой Ñтороны, при превышении некоторого критичеÑкого Ð·Ð½Ð°Ñ‡ÐµÐ½Ð¸Ñ Ð±Ð°Ð·Ð¾Ð²Ð¾Ð¹ линии интерферометричеÑÐºÐ°Ñ Ð¾Ð±Ñ€Ð°Ð±Ð¾Ñ‚ÐºÐ° ÑтановитÑÑ Ð² принципе невозможной вÑледÑтвие проÑтранÑтвенной декоррелÑции.
КритичеÑкое значение перпендикулÑрной проÑтранÑтвенной базы Ð´Ð»Ñ ÐºÐ°Ð¶Ð´Ð¾Ð¹ пары Ñнимков может быть вычиÑлено по Ñледующей формуле:
где: Ð’n, cr — критичеÑÐºÐ°Ñ Ð±Ð°Ð·Ð¾Ð²Ð°Ñ Ð»Ð¸Ð½Ð¸Ñ;
λ — длина зондирующей волны радиолокатора;
Rr — проÑтранÑтвенное разрешение в направлении наклонной дальноÑти.
ОÑтальные аргументы ÑоответÑтвуют обозначениÑм к риÑ. 1: H1 и H2 — выÑота орбиты радиолокатора; R1 и R2 – дальноÑть (путь зондирующей волны); θ1 и θ2 – углы между дальноÑтью и выÑотой; Bâ”´ и Bâ•‘ — перпендикулÑÑ€Ð½Ð°Ñ Ð¸ Ð¿Ð°Ñ€Ð°Ð»Ð»ÐµÐ»ÑŒÐ½Ð°Ñ ÑоÑтавлÑющие базовой линии; Ð’ — Ñ€ÐµÐ·ÑƒÐ»ÑŒÑ‚Ð¸Ñ€ÑƒÑŽÑ‰Ð°Ñ Ð±Ð°Ð·Ð¾Ð²Ð°Ñ Ð»Ð¸Ð½Ð¸Ñ.
По Ñтой формуле, критичеÑÐºÐ°Ñ Ð¿ÐµÑ€Ð¿ÐµÐ½Ð´Ð¸ÐºÑƒÐ»ÑÑ€Ð½Ð°Ñ Ð±Ð°Ð·Ð¾Ð²Ð°Ñ Ð»Ð¸Ð½Ð¸Ñ Ð´Ð»Ñ Ð¸Ñпользуемых данных ENVISAT/ ASAR Image Mode ÑоÑтавлÑет порÑдка 900 — 1500 м, а Ð´Ð»Ñ Ð´Ð°Ð½Ð½Ñ‹Ñ… ALOS/PALSAR — 6500 м в режиме FBS и 13 000 м в режиме FBD (6500 м при перекреÑтной обработке режимов FBS и FBD). ÐžÐ¿Ñ‚Ð¸Ð¼Ð°Ð»ÑŒÐ½Ð°Ñ Ð´Ð»Ñ Ñ€Ð°Ñчета Ñмещений перпендикулÑÑ€Ð½Ð°Ñ Ð±Ð°Ð·Ð° в Ñлучае обоих вышеуказанных Ñпутников колеблетÑÑ Ð² диапазоне от 0 до 30 % от критичеÑкой базы.
Временной базой называют промежуток времени, прошедший между Ñъемкой изображений, ÑоÑтавлÑющих интерферометричеÑкую пару. ПонÑтие временной базы напрÑмую ÑвÑзано Ñ Ñ‚Ð°ÐºÐ¾Ð¹ важной проблемой, как Ð²Ñ€ÐµÐ¼ÐµÐ½Ð½Ð°Ñ Ð´ÐµÐºÐ¾Ñ€Ñ€ÐµÐ»ÑциÑ, Ð²Ð¾Ð·Ð½Ð¸ÐºÐ°ÑŽÑ‰Ð°Ñ Ð·Ð° Ñчет изменений рельефа, раÑтительноÑти, влажноÑти, шероховатоÑти и других ÑвойÑтв отражающей радарный луч поверхноÑти, произошедших за период между Ñъемками. Проблема временной декоррелÑции может быть решена увеличением длины зондирующей волны (что повышает «проÑвечивающую» ÑпоÑобноÑть радиоволн) либо Ñокращением временной базы (Ñ‚. е. промежутка между Ñъемками).
Ð’ ходе мониторинга Ñмещений на меÑторождении Тенгиз был применен вышеуказанный подход интерферометричеÑкой обработки пар Ñнимков, отобранных из вÑех возможных пар по принципу наименьших базовых линий. При Ñтом результат обработки в итоге был предÑтавлен в виде карт Ñмещений за периоды между первой Ñъемкой и поочередно каждой из поÑледующих Ñъемок. То еÑть результат показывает развитие Ñмещений во времени, а также итоговую карту Ñмещений за веÑÑŒ период наблюдений.
Помимо проÑтранÑтвенной и временной баз, важным параметром, определÑющим возможноÑть или невозможноÑть интерферометричеÑкой обработки, ÑвлÑетÑÑ Ñ€Ð°Ð·Ð½Ð¾Ñть положений ДоплеровÑкого центроида Ð²Ñ€Ð°Ñ‰ÐµÐ½Ð¸Ñ Ð—ÐµÐ¼Ð»Ð¸ Ð´Ð»Ñ Ñнимков пары. Ð’Ñе обработанные пары характеризуютÑÑ Ð½Ð¸Ð·ÐºÐ¸Ð¼Ð¸ либо нормальными значениÑми Ñтого параметра, что положительно ÑказалоÑÑŒ на обработке.
ЕÑли Ð·Ð½Ð°Ñ‡ÐµÐ½Ð¸Ñ Ð¿Ñ€Ð¾ÑтранÑтвенной и временной баз, а также разницы положений ДоплеровÑкого центроида позволÑÑŽÑ‚ выполнить интерферометричеÑкую обработку, то ÑтановитÑÑ Ð²Ð¾Ð·Ð¼Ð¾Ð¶Ð½Ñ‹Ð¼ раÑÑчитать интерферограмму.
Каждый радарный Ñнимок интерферометричеÑкой пары Ñодержит в Ñебе амплитудный и фазовый Ñлой. Ðмплитудный Ñлой более пригоден Ð´Ð»Ñ Ð²Ð¸Ð·ÑƒÐ°Ð»ÑŒÐ½Ð¾Ð³Ð¾ анализа. Ð ÐµÐ·ÑƒÐ»ÑŒÑ‚Ð¸Ñ€ÑƒÑŽÑ‰Ð°Ñ Ñ„Ð°Ð·Ð° Ф, Ð¿Ð¾Ð»ÑƒÑ‡ÐµÐ½Ð½Ð°Ñ Ð² ходе интерферометричеÑкой обработки фазовых Ñлоев Ñнимков интерферометричеÑкой пары, в общем Ñлучае ÑоÑтоит из Ñледующих компонентов:
Ф = Фtopo + Фdef + Фatm + Фn,
где:
Фtopo – фазовый набег за Ñчет обзора топографии под Ð´Ð²ÑƒÐ¼Ñ Ñ€Ð°Ð·Ð½Ñ‹Ð¼Ð¸ углами;
Фdef – фазовый набег за Ñчет ÑÐ¼ÐµÑ‰ÐµÐ½Ð¸Ñ Ð¾Ñ‚Ñ€Ð°Ð¶Ð°ÑŽÑ‰ÐµÐ¹ поверхноÑти в период между Ñъемками;
Фatm – фазовый набег за Ñчет Ñ€Ð°Ð·Ð»Ð¸Ñ‡Ð¸Ñ Ð´Ð»Ð¸Ð½ оптичеÑких путей из-за Ð¿Ñ€ÐµÐ»Ð¾Ð¼Ð»ÐµÐ½Ð¸Ñ Ð² Ñреде раÑпроÑÑ‚Ñ€Ð°Ð½ÐµÐ½Ð¸Ñ Ñигнала;
Фn – вариации фазы в результате Ñлектромагнитного шума.
ÐепоÑредÑтвенно интерферометричеÑÐºÐ°Ñ Ð¾Ð±Ñ€Ð°Ð±Ð¾Ñ‚ÐºÐ° пары Ñнимков в общем Ñлучае ÑоÑтоит из неÑкольких базовых шагов:
- Совмещение оÑновного и вÑпомогательного радарных изображений интерферометричеÑкой пары (в автоматичеÑком режиме либо Ñ Ñ€ÑƒÑ‡Ð½Ñ‹Ð¼ вводом контрольных точек).
- Ð“ÐµÐ½ÐµÑ€Ð°Ñ†Ð¸Ñ Ð¸Ð½Ñ‚ÐµÑ€Ñ„ÐµÑ€Ð¾Ð³Ñ€Ð°Ð¼Ð¼Ñ‹, ÑвлÑющейÑÑ Ñ€ÐµÐ·ÑƒÐ»ÑŒÑ‚Ð°Ñ‚Ð¾Ð¼ комплекÑного поÑлементного Ð¿ÐµÑ€ÐµÐ¼Ð½Ð¾Ð¶ÐµÐ½Ð¸Ñ Ð¾Ñновного Ð¸Ð·Ð¾Ð±Ñ€Ð°Ð¶ÐµÐ½Ð¸Ñ Ð¸ вÑпомогательного изображениÑ, геометричеÑки Ñовмещенного Ñ Ð¾Ñновным.
- Разделение компонентов фазы Фdef и Фtopo за Ñчет Ñинтеза фазы рельефа Ñ Ð¿Ð¾Ð¼Ð¾Ñ‰ÑŒÑŽ имеющейÑÑ Ñ†Ð¸Ñ„Ñ€Ð¾Ð²Ð¾Ð¹ модели рельефа (ЦМР) или цифровой модели меÑтноÑти (ЦММ), либо модели ÑллипÑоида, либо путем Ð·Ð°Ð´Ð°Ð½Ð¸Ñ Ñредней по площади Ñнимка выÑотной отметки (Ð´Ð»Ñ Ñ€Ð°Ð²Ð½Ð¸Ð½Ð½Ñ‹Ñ… облаÑтей) Ñ Ð¿Ð¾Ð»ÑƒÑ‡ÐµÐ½Ð¸ÐµÐ¼ на выходе дифференциальной интерферограммы, на которой один полный Ñпектр цветов (интерференционный цикл) ÑоответÑтвует ÑмещениÑм земной поверхноÑти, равным половине длины волны радиолокатора.
- Ð¤Ð¸Ð»ÑŒÑ‚Ñ€Ð°Ñ†Ð¸Ñ Ð¸Ð½Ñ‚ÐµÑ€Ñ„ÐµÑ€Ð¾Ð³Ñ€Ð°Ð¼Ð¼Ñ‹, позволÑÑŽÑ‰Ð°Ñ Ð² определенной Ñтепени уменьшить фазовый шум (помехи) за Ñчет Ð·Ð°Ð³Ñ€ÑƒÐ±Ð»ÐµÐ½Ð¸Ñ Ð²Ñ‹Ñ…Ð¾Ð´Ð½Ð¾Ð¹ карты Ñмещений земной поверхноÑти.
- Получение файла когерентноÑти Ð´Ð»Ñ Ð¾Ð±Ð»Ð°Ñти Ð¿ÐµÑ€ÐµÐºÑ€Ñ‹Ñ‚Ð¸Ñ Ð´Ð²ÑƒÑ… Ñнимков, ÑоÑтавлÑющих интерферометричеÑкую пару, в значениÑÑ… от 0 до 1 Ð´Ð»Ñ ÐºÐ°Ð¶Ð´Ð¾Ð¹ пары ÑоответÑтвующих друг другу пикÑелей.
- Развертка фазы (процедура перехода от отноÑительных значений фазы к абÑолютным) и уÑтранение разрывов фазы.
- ÐšÐ¾Ñ€Ñ€ÐµÐºÑ†Ð¸Ñ Ð·Ð½Ð°Ñ‡ÐµÐ½Ð¸Ñ Ð±Ð°Ð·Ð¾Ð²Ð¾Ð¹ линии по наземным контрольным точкам.
- Преобразование абÑолютных значений фазы в ÑÐ¼ÐµÑ‰ÐµÐ½Ð¸Ñ Ð·ÐµÐ¼Ð½Ð¾Ð¹ поверхноÑти в миллиметрах.
ÐŸÐ¾Ð»Ð½Ð°Ñ Ñ†ÐµÐ¿Ð¾Ñ‡ÐºÐ° интерферометричеÑкой обработки, а также многие другие дополнительные функции реализованы в иÑпользованном авторами программном комплекÑе SARscape (разработчик SARMAP, ШвейцариÑ), ÑвлÑющемÑÑ Ð´Ð¾Ð¿Ð¾Ð»Ð½Ð¸Ñ‚ÐµÐ»ÑŒÐ½Ñ‹Ð¼ модулем программы ENVI (Exelis VIS, СШÐ).
Ð”Ð»Ñ Ð°Ð½Ð°Ð»Ð¸Ð·Ð° не пар, а многопроходных цепочек интерферометричеÑких радарных Ñнимков в SARscape реализованы модификации радарной интерферометрии: технологии интеферометрии поÑтоÑнных раÑÑеивателей (Persistent Scatterers) и интерферометрии малых базовых линий (SBas).
ПОСЛЕДОВÐТЕЛЬÐÐЯ ИÐТЕРФЕРОМЕТРИЯ ПОСТОЯÐÐЫХ Ð ÐССЕИВÐТЕЛЕЙ Ð ÐДÐÐ ÐОГО СИГÐÐЛР(PS)
Ðтот вариант радарной интерферометрии характеризуетÑÑ Ð¼Ð°ÐºÑимально возможной точноÑтью оценки Ñмещений (2–4 мм по выÑоте). Входными данными Ð´Ð»Ñ Ð¾Ð±Ñ€Ð°Ð±Ð¾Ñ‚ÐºÐ¸ должны ÑвлÑтьÑÑ Ð½Ðµ менее 30 Ñнимков одной и той же территории за разные даты, Ñделанных в одной и той же геометрии Ñъемки Ñпутникового радиолокатора [2–4].
Программой автоматичеÑки выбираетÑÑ Ð¾Ñновное изображение, на которое автоматичеÑки Ñ Ñ‚Ð¾Ñ‡Ð½Ð¾Ñтью до 1/100 пикÑÐµÐ»Ñ ÐºÐ¾Ñ€ÐµÐ³Ð¸ÑтрируютÑÑ Ð¾Ñтальные Ñнимки интерферометричеÑкой цепочки. Далее программа Ñтроит так называемые интерферограммы (комплекÑно поÑлементно перемноженные фазовые Ñлои радарных Ñнимков) по каждой паре Ñнимков. Затем Ð´Ð»Ñ ÐºÐ°Ð¶Ð´Ð¾Ð¹ пары оцениваютÑÑ Ð²ÐµÐ»Ð¸Ñ‡Ð¸Ð½Ñ‹ когерентноÑти (меры коррелÑции фаз радарных Ñнимков). Также Ð´Ð»Ñ ÐºÐ°Ð¶Ð´Ð¾Ð¹ пары ÑтроÑÑ‚ÑÑ ÐºÐ°Ñ€Ñ‚Ñ‹ величин Ñтандартных отклонений амплитуд Ñнимков.
Затем программой определÑÑŽÑ‚ÑÑ Ñ‚Ð¾Ñ‡ÐºÐ¸ — поÑтоÑнные (или уÑтойчивые) раÑÑеиватели радарного Ñигнала. Ð”Ð»Ñ Ð²Ñ‹Ð±Ð¾Ñ€Ð° точек иÑпользуетÑÑ Ð½ÐµÑколько порогов (порог коррелÑции амплитуд, порог когерентноÑти, порог проÑтранÑтвенного и временного отклонений величин Ñмещений первой итерации и Ñ‚. д.). ПоÑле того как поÑтоÑнные раÑÑеиватели определены, Ð´Ð»Ñ Ð½Ð¸Ñ… выполнÑетÑÑ Ð¿Ñ€Ð¾Ñ†ÐµÐ´ÑƒÑ€Ð° оценки фазовых разноÑтей и мультивременной развертки фазы Ð´Ð»Ñ Ñ‚Ð¾Ñ‡ÐµÑ‡Ð½Ñ‹Ñ… целей. Именно в разноÑти фаз каждого Ñнимка «зашита» величина Ñмещений за период между Ñъемками Ñтих Ñнимков.
Таким образом, Ð´Ð»Ñ ÐºÐ°Ð¶Ð´Ð¾Ð¹ из выбранных точек воÑÑтанавливаетÑÑ Ñ…Ñ€Ð¾Ð½Ð¾Ð»Ð¾Ð³Ð¸Ñ Ð¸Ð·Ð¼ÐµÐ½ÐµÐ½Ð¸Ñ Ñ„Ð°Ð·Ñ‹ во времени, ÐºÐ¾Ñ‚Ð¾Ñ€Ð°Ñ Ð·Ð°Ñ‚ÐµÐ¼ математичеÑки переÑчитываетÑÑ Ð² ÑÐ¼ÐµÑ‰ÐµÐ½Ð¸Ñ Ð² миллиметрах. Дополнительно в процеÑÑе обработки применÑетÑÑ Ñпециальный фильтр, удалÑющий возможное влиÑние атмоÑферы на интерферометричеÑкую фазу.
Результатом обработки ÑвлÑетÑÑ Ð²ÐµÐºÑ‚Ð¾Ñ€Ð½Ñ‹Ð¹ файл точек, в атрибутах которых запиÑаны:
- ÑÐ¼ÐµÑ‰ÐµÐ½Ð¸Ñ Ð½Ð° каждую дату Ñъемки;
- ÑÑ€ÐµÐ´Ð½ÐµÐ³Ð¾Ð´Ð¾Ð²Ð°Ñ ÑкороÑть Ñмещений;
- ÑÑƒÐ¼Ð¼Ð°Ñ€Ð½Ð°Ñ Ð²ÐµÐ»Ð¸Ñ‡Ð¸Ð½Ð° Ñмещений;
- когерентноÑть;
- выÑота над ÑллипÑоидом WGS-84.
Пример результата обработки по технологии PS приведен на риÑ. 2. ОÑновным недоÑтатком данной технологии ÑвлÑетÑÑ ÐµÐµ применимоÑть только Ð´Ð»Ñ Ð·Ð°Ñтроенных территорий, а также Ð´Ð»Ñ Ð¾Ñ‚Ð´ÐµÐ»ÑŒÐ½Ñ‹Ñ… зданий и Ñооружений при Ñъемке в выÑоком разрешении.

По результатам 30-проходной Ñъемки выÑвлены поÑтоÑнные раÑÑеиватели радарного Ñигнала, Ð´Ð»Ñ ÐºÐ°Ð¶Ð´Ð¾Ð³Ð¾ из которых извеÑтны ÑÐ¼ÐµÑ‰ÐµÐ½Ð¸Ñ Ð½Ð° каждую дату Ñъемки и ÑÑ€ÐµÐ´Ð½ÐµÐ³Ð¾Ð´Ð¾Ð²Ð°Ñ ÑкороÑть Ñмещений. Точки покрашены по атрибуту «ÑÑ€ÐµÐ´Ð½ÐµÐ³Ð¾Ð´Ð¾Ð²Ð°Ñ ÑкороÑть Ñмещений» (зеленые – Ñтабильные, желтые – оÑеданиÑ, краÑные — интенÑивные оÑеданиÑ). График Ñмещений за период 30 Ñъемок Ð´Ð»Ñ Ð¾Ð´Ð½Ð¾Ð¹ из точек — в правом верхнем углу риÑ. 2.
ПЕРЕКРЕСТÐÐЯ ИÐТЕРФЕРОМЕТРИЯ МÐЛЫХ БÐЗОВЫХ ЛИÐИЙ (SBAS)
Ð˜Ð½Ñ‚ÐµÑ€Ñ„ÐµÑ€Ð¾Ð¼ÐµÑ‚Ñ€Ð¸Ñ Ð¼Ð°Ð»Ñ‹Ñ… базовых линий, в отличие от интерферометрии поÑтоÑнных раÑÑеивателей, ÑвлÑетÑÑ Ð¼ÐµÐ½ÐµÐµ автоматизированным методом требующим большей квалификации иÑполнителÑ. Ð’ Ñтом Ñлучае уÑиливаетÑÑ Ð²ÐºÐ»Ð°Ð´ ÑтатиÑтики в финальный результат за Ñчет перекреÑтной обработки очень большого количеÑтва интерферометричеÑких пар при том же Ñамом количеÑтве Ñнимков [5]. Ð”Ð»Ñ Ð¾Ð±Ñ€Ð°Ð±Ð¾Ñ‚ÐºÐ¸ по Ñтому методу не обÑзательно наличие 30 Ñнимков, обработка возможна и при меньшем их количеÑтве.
Ðапример, в Ñлучае 15-проходной цепочки общее возможное количеÑтво пар Ñнимков доÑтигает 105. Из них по величине наименьшей проÑтранÑтвенной базы выбираютÑÑ, к примеру, 30–40 пар. Пары могут быть перекреÑтными (первый проход Ñо вторым, второй Ñ Ñ‚Ñ€ÐµÑ‚ÑŒÐ¸Ð¼, первый Ñ Ñ‚Ñ€ÐµÑ‚ÑŒÐ¸Ð¼, второй Ñ Ñ‡ÐµÑ‚Ð²ÐµÑ€Ñ‚Ñ‹Ð¼ и Ñ‚. п.). Обработка каждой пары ведетÑÑ Ð² полуавтоматичеÑком режиме Ñ Ð²Ñ‹Ð¿Ð¾Ð»Ð½ÐµÐ½Ð¸ÐµÐ¼ нижеуказанных четырех шагов:
- ÐвтоматичеÑÐºÐ°Ñ ÐºÐ¾Ñ€ÐµÐ³Ð¸ÑтрациÑ, раÑчет интерферограммы, Ñинтез фазы рельефа, вычитание фазы рельефа из интерферограммы, Ñ„Ð¸Ð»ÑŒÑ‚Ñ€Ð°Ñ†Ð¸Ñ Ð´Ð¸Ñ„Ñ„ÐµÑ€ÐµÐ½Ñ†Ð¸Ð°Ð»ÑŒÐ½Ð¾Ð¹ интерферограммы, раÑчет когерентноÑти, развертка фазы.
- Ðабор точек Ñ Ð¸Ð·Ð²ÐµÑтными координатами и выÑотами Ð´Ð»Ñ ÐºÐ¾Ñ€Ñ€ÐµÐºÑ†Ð¸Ð¸ орбитальных параметров.
- РаÑчет Ñкорректированных дифференциальных интерферограмм и развернутых фаз.
- ИнверÑÐ¸Ñ Ð¿Ð¾Ð»ÑƒÑ‡ÐµÐ½Ð½Ñ‹Ñ… перекреÑтных во времени развернутых фаз по методике SBas Ñ Ð²Ð¾ÑÑтановлением поÑледовательной во времени иÑтории Ñмещений. При Ñтом в итоге также воÑÑтанавливаетÑÑ Ñ…Ñ€Ð¾Ð½Ð¾Ð»Ð¾Ð³Ð¸Ñ Ñмещений от первого Ñнимка цепочки до поÑледнего (пример — на риÑ. 3).

ЗÐКЛЮЧЕÐИЕ
Вышеприведенные данные можно обобщить Ñледующим образом:
- по результатам интерферометричеÑкой обработки пар радарных Ñнимков возможно получать карты Ñмещений земной поверхноÑти, однако в карте Ñмещений, раÑÑчитанной только лишь по одной паре Ñнимков, могут приÑутÑтвовать атмоÑферные артефакты, неточноÑти орбитальных параметров, фазовый шум и Ñ‚.д.;
- обработка многопроходных интерферометричеÑких цепочек, в отличие от обработки отдельных пар Ñнимков, позволÑет иÑключить влиÑние атмоÑферы, значительно уменьшить орбитальную погрешноÑть и вообще получить гораздо более ÑтатиÑтичеÑки выверенный результат;
- Ñреди методов интерферометричеÑкой обработки многопроходных цепочек наиболее извеÑтными ÑвлÑÑŽÑ‚ÑÑ Ð¼ÐµÑ‚Ð¾Ð´Ñ‹ PS и SBas (оба реализованы в программном комплекÑе SARscape);
- полноÑтью автоматизированный метод PS позволÑет добитьÑÑ Ð¼Ð°ÐºÑимальных точноÑтей замера Ñмещений и деформаций зданий и Ñооружений (2–4 мм по выÑоте), а также позволÑет изучать ÑÐ¼ÐµÑ‰ÐµÐ½Ð¸Ñ Ð¸ деформации в динамике, однако он применим в большей Ñтепени к заÑтроенным территориÑм и требует в качеÑтве иÑходных данных не менее 30 Ñнимков;
- полуавтоматизированный метод SBas позволÑет также изучать ÑÐ¼ÐµÑ‰ÐµÐ½Ð¸Ñ Ð¸ деформации, характеризуетÑÑ Ñ‡ÑƒÑ‚ÑŒ меньшей точноÑтью, чем PS, и большими требованиÑми к пользователю, но зато он может работать не только по заÑтроенным территориÑм, а также он не требует обÑзательных 30 Ñнимков (метод будет работать и по 12–15 Ñнимкам, однако Ð´Ð»Ñ Ð¿Ð¾Ð»ÑƒÑ‡ÐµÐ½Ð¸Ñ Ð³Ð°Ñ€Ð°Ð½Ñ‚Ð¸Ñ€Ð¾Ð²Ð°Ð½Ð½Ð¾Ð³Ð¾ результата вÑе же надо ÑтремитьÑÑ Ðº 30 проходам).
СПИСОК ЛИТЕРÐТУРЫ:
- M. Richards: «A Beginner’s Guide to Interferometric SAR Concepts and Signal Processing». IEEE Aerospace and Electronic, Vol. 22, No. 9, September 2007.
- Ferretti, C. Prati and F. Rocca: «Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry». Geoscience and Remote Sensing, IEEE Transactions on, vol. 38, no. 5, Part 1, Sept. 2000, pp. 2202 – 2212.
- Ferretti, C. Prati and F. Rocca: «Permanent scatterers in SAR interferometry». Geoscience and Remote Sensing, IEEE Transactions on, vol. 39, no. 1, Jan. 2001, pp. 8 – 20.
- Hooper, H. Zebker, P. Segall, and B. Kampes: «A new method for measuring deformation on volcanoes and other non-urban areas using InSAR persistent scatterers». Geophysical Research Letters, vol. 31, December 2004.
- P. Berardino, G. Fornaro, R. Lanari, E. Sansosti: «A new algorithm for surface deformation monitoring based on Small Baseline differential SAR Interferometry». IEEE Aerospace and Electronic, Vol. 40, No. 11, November 2002.