Совзонд - Геоматика

36• коррекция орбит по наземным опорнымточкам (GCP);• пересчет фазы в высоты;• геокодирование ЦМР в систему WGS-4, проекцииUTM-40 (шаг сетки 50 м).Создание интерферограмм выполнялось в трехвариантах: без входной ЦМР и с входными ЦМРGTOPO30 и SRTM.Опорными точками для коррекции орбит (GCP)служили 25 наземных пунктов, координаты и высотыкоторых были определены из GPS-наблюдений.В общей сложности было построено 24 радарныхЦМР, из которых вырезались прямоугольные фрагментыразмером 25х30 км на исследуемую территориюг. Перми и преобразовывались в матрицывысот с шагом сетки 50х50 м.Оценка точности выполнялась по разности высотΔH радарных ЦМР H DEM и топографическойповерхности H Topo на контрольных точках:ΔH = H DEM ― H TopoРис. 1. Выбранная для построения ЦМР территориягорода Перми с контрольными точкамиПОСТРОЕНИЕ ЦМР ПО РАДАРНЫМ СНИМКАМДля построения ЦМР использовалась многопроходнаясерия радарных снимков ALOS/PALSAR,выполненных с 23.01.2008 по 16.12.2010 г. Из 12снимков с поляризацией HH, образующих 66 пар,было выбрано 8 пар с временными базисами 46, 92,138 и 368 дней с отношением пространственногобазиса B к его критическому значению B êð , в интервале~ 0.05 – 0.15, для которых значения разностивысот ΔH 2π , соответствующие изменению фазы на2π, близки к перепаду высот на выбранной территории(табл. 1).Обработка снимков проводилась в программномпродукте ENVI SARscape по стандартнойметодике [2]:• оценка базисных линий;• создание интерферограмм;• адаптивная фильтрация (метод Goldstain) и расчеткогерентности;• развертка фазы (метод MCF);Координаты и высоты H Topo 90 контрольных точекбыли определены из GPS-наблюдений (на рис. 1показаны красным цветом). Значения H Topo вычислялисьпо высотам в узлах сетки интерполированиемкубическими полиномами.В качестве основных показателей точностивысот радарных ЦМР были приняты следующиепараметры:• среднее значение разности высот ΔH (систематическаяошибка);• средняя квадратическая ошибка RMSE ΔH ;• средняя абсолютная ошибка MAE ΔH ;• вероятная линейная ошибка LE90 ΔH ;• минимальное ΔH min и максимальное ΔH maxзначения ошибок.ОЦЕНКА ТОЧНОСТИ ВЫСОТ ЦМР ПОКОНТРОЛЬНЫМ GPS-ТОЧКАМПолученные значения оценок параметров точностирадарных ЦМР, без входной ЦМР и для двухвидов входных ЦМР GTOPO30 и SRTM, приведеныв табл. 2–4.ГЕОМАТИКА №3'2012