Введение
Цифровая модель рельефа (ЦМР) и цифровая модель местности (ЦММ) являются важными инструментами для создания трехмерного представления местности и используются в различных областях науки и производства, таких как география, геология, экология, городское планирование и многие другие. Эти модели позволяют глубже понять характер местности, ее направления, склоны, а также методы трехмерного отображения рельефа поверхности Земли. ЦМР широко применяется в геоинформационных системах для решения различных задач – от простого анализа данных до сложного моделирования, направленного на определение оптимальных мест для строительства инженерных сооружений. Кроме того, цифровые модели рельефа находят многочисленные применения в гражданском строительстве, в частности при проектировании канализационных сетей, выделении контурных линий, анализе склонов и естественных направлений изучаемого участка. Они также полезны при проведении работ по гражданскому планированию и в оценке объемов земляных работ и обратной засыпки в строительных проектах [1-3].
Существует несколько типов цифровых моделей рельефа. Цифровая модель рельефа (ЦМР) отображает рельеф местности на основе высот, отображая его в целом. Цифровая модель рельефа на основе поверхности (digital surface models DSM) учитывает все объекты, находящиеся над землей, такие как деревья, здания и другие сооружения. Наземная цифровая модель рельефа (DTM) представляет поверхность Земли с учетом удаления различных объектов, чтобы показать только естественный рельеф [4-6].
Целью исследования является использование компьютерного моделирования для имитации аэрофотосъемки с использованием беспилотного летательного аппарата для повышения точности создания цифровой модели территории.
Материалы и методы исследования
Для достижения этой цели проводится исследование влияния увеличения количества контрольных точек, полученных в результате геодезической съемки предполагаемой территории. Данный подход позволяет достичь наилучших результатов при создании реальной цифровой карты с учетом оптимальных финансовых и временных затрат, что способствует повышению эффективности процесса картографирования и геодезических изысканий.
Существует несколько общедоступных универсальных цифровых моделей рельефа, которые отличаются друг от друга по точности (5, 10, 12,5, 30 и 90 м), таких как SRTM, ALOS AW3D30, ASTER GDEM, TANDEM, ALOS PALSAR, EMG2008 [7]. В качестве исходных данных для исследований использована цифровая модель рельефа SRTM, полученная из общедоступных интернет-источников, таких как OpenStreetMap. Для визуализации, редактирования и моделирования аэрофотосъемки на базе модели SRTM применялось программное обеспечение Blender 3D. Для этого необходимо настроить программное обеспечение и установить специальное приложение под названием GSI, которое позволяет соединяться с различными источниками общедоступной цифровой модели рельефа SRTM [8; 9]. На рисунке 1 показано, как использовать приложение GSI, а также как выбирать подходящую цифровую модель местности.
Также использовалось программное обеспечение Metashape для получения ортофотопланов и цифровой модели местности на основании данных аэрофотосъемки. В программном продукте Civil 3D построены поверхности на основе топографической съемки для каждого района исследования, а результаты экспортированы в формате LandXML [10].
Рис. 1. Функции GSI (слева), OSM (справа) в программе Blender 3D Источник: составлено автором на основе результатов исследования
Рис. 2. Порядок исследования Источник: составлено автором на основе результатов исследования
В статье представлен анализ точности общедоступной цифровой модели рельефа SRTM, а также оценка ее точности. Для этого использовались геодезические данные в виде следующих топографических планов:
1. Топографический план в масштабе 1:500 для строительства инженерных проектов жилого комплекса «I Love» (г. Москва, СВАО).
2. Топографический план в масштабе 1:500 для строительства инженерных проектов жилого комплекса «Бирюлево» (г. Москва, ЮАО).
3. Модель SRTM была применена для каждого исследуемого района и получена с использованием программного обеспечения для 3D-моделирования, такого как Blender.
Порядок исследования, представленный на рисунке 2, начинается с подготовки и планирования работы, а также сбора информации. На первом этапе осуществляется сбор топографической съемки для каждого исследуемого района.
На втором этапе проводится имитационное моделирование в программе Blender 3D, которое проходит несколько шагов. В первый шаг входят импорт и подготовка модели в Blender 3D. В рамках данного исследования авторы импортировали модель из известных источников общедоступной цифровой модели местности, таких как Get OSM и Get Elevation (SRTM). После выбора источников выбираем участок, нажимаем на клавиатуре букву E, чтобы импортировать выбранную область в программу Blender. Таким образом, получим цифровую модель местности для районов жилого комплекса «I Love» и жилого комплекса «Бирюлево» в среду программного продукта Blender 3D, что показано на рисунках 3а и 3б.
Второй шаг в имитационном моделировании в программных продуктах Blender 3D – это создание опорных геодезических точек для преобразования импортируемых моделей в местные системы координат. Для этого необходимо смоделировать опорные точки. Используем одну из фигур, которую можно добавить в программу (например, куб), разделяем её на две части и раскрашиваем в разные цвета, чтобы четко обозначить центр куба, который будет служить координатой опорной точки. Повторяем действия для других опорных точек, и в структуре проекта появятся все опорные точки. Перемещаем эти точки на высотную модель SRTM и определяем местонахождение каждой точки, которая соответствует топографическому плану района исследования [11-13], что показано на рисунках 4а и 4б.
Третий шаг – это моделирование съемки. На основании длины и ширины участка, а также максимальной и минимальной высоты модуля можно определить размер пикселей по формуле:
Pxl = Хп / Хи, (1)
где Хп – физический размер матрицы и Хи – размер изображения.
Рис. 3а. Модель ЦММ района ЖК «I Love» Источник: составлено автором на основе результатов исследования
Рис. 3б. Модель ЦММ района ЖК «Бирюлево» Источник: составлено автором на основе результатов исследования
Рис. 4а. Пример маркеров в качестве опорных точек Источник: составлено автором на основе результатов исследования
Рис. 4б. Пример маркеров опорных точек Источник: составлено автором на основе результатов исследования
Также определяем параметры камеры аэрофотосъемки, её калибровку, расчет маршрутов и значений продольного и поперечного перекрытия между снимками, а также рассчитываем высоту моделируемой съемки с учетом разрешения, необходимого для получения требуемой модели [14].
Расчетная высота фотографирования Нф [15]
Нф = 2R∙GSD∙(fк / δ), (2)
где Нф – высота фотографирования относительно заданной поверхности (м);
fк – фокусное расстояние фотосистемы (м);
GSD – проекция пикселя матрицы ЦАФА на заданной поверхности (м);
R – разрешение на матрице (цикл. /пиксель)
δ – размер пикселя матрицы (м).
Расчетное продольное перекрытие аэрофотоснимков, рх % [15]
рх = рхо + Δрh +Δрн, (3)
где рхо – минимальное продольное перекрытие = 60%,
Δрh – поправка за рельеф,
. (4)
где Δh – разность между самой высокой точкой и самой низкой точкой участка,
Δрн – навигационная поправка за ошибку самолетовождения.
Расчетное поперечное перекрытие аэрофотоснимков, ру % [15].
Для масштаба 1:10 000 и крупнее
. (5)
Продольный базис фотографирования, Вх (м)
, (6)
где MФ = НФ / f : знаменатель масштаба фотографирования
Поперечный базис фотографирования, Ву (м) [15]
, (7)
где XП, YП: физический размер матрицы (м)
Интервал фотографирования τ (сек.):
, (8)
где XП – физический размер фото матрицы вдоль линии полета (м),
W – путевая скорость БВС (м/с).
Количество аэроснимков в одном маршруте, Nх [15]:
Nx = (Lx / Bx) + N3; N3 = 4, (9)
где Lx – длина участка съёмки (м)
Количество маршрутов, Nу [15]:
Ny = (Ly / By) + 0.5∙j; j = 2, (10)
где Ly – ширина участка (м)
Таблица 1
Расчетные данные аэрофотосъемки для двух районов: «I Love» и «Бирюлево»
Название элементов |
Усл. обозначение |
ЖК «I Love» |
ЖК «Бирюлево» |
Длина участка, км |
Lx |
0.38 |
0,7 |
Ширина участка, км |
Ly |
0.4 |
0,58 |
Максимальная высота участка, м |
h max.у. |
183.00 |
178 |
Минимальная высота участка, м |
h min.y. |
170.00 |
160 |
Высота среднего уровня участка, м |
h ср.у. |
176.50 |
169 |
Физический размер матрицы |
Хп (lx) |
0.0086 |
0.0086 |
Уп (lу) |
0.0129 |
0.0129 |
|
Размер изображения |
Хи |
3648 |
3648 |
Уи |
5472 |
5472 |
|
Фокусное расстояние (м) |
f |
0.0088 |
0.0088 |
Размер пикселя (м) |
pxl |
0,00000236 |
0,00000236 |
Количество маршрутов |
5 |
5 |
|
Количество снимков |
52 |
72 |
|
Высота съемки, м |
186 |
186 |
Источник: составлено автором на основе результатов исследования.
Все результаты расчетов и параметры аэрофотосъемки для двух районов исследования представлены в таблице 1.
В результате моделирования съемки с использованием программных продуктов Blender 3D были получены данные аэрофотосъемки для двух районов исследования. Эти данные можно обрабатывать в специализированном программном обеспечении Metashape для построения цифровой модели местности и ортофотоплана. Следует отметить, что программа используется для решения задач построения трехмерных моделей на основе фотографий и является универсальной, позволяя создавать высоко детализированные 3D-модели только по фотографиям.
При работе в программе Metashape после указания параметров камеры, добавления и выравнивания снимков строится плотное облако точек и модель привязывается по четырем опорным геодезическим точкам, далее создаются 3D-модель, ЦММ и ортофотоплан [16].
После этого возвращаемся к облаку точек и увеличиваем количество опорных точек до восьми, продолжая построение цифровой модели местности и ортофотоплана. На последнем этапе увеличиваем количество опорных точек до двадцати и вновь формируем цифровую модель местности и ортофотоплан.
Повторяем эту работу для второго района исследования. В итоге авторы получали три цифровые модели местности и три ортофотоплана для каждого района.
Первая модель привязана по четырем опорным точкам, вторая – по восьми точкам, а третья, финальная модель – по максимальному количеству 20 геодезических точек (рис. 5).
На рисунках 6а, 6б представлено, как был получен ортофотоплан для каждого района работы, который был максимально привязан к опорным геодезическим точкам. В данном примере использовано 20 точек.
Экспортируем цифровую модель местности для каждого примера в формате LandXML.
На последнем этапе исследований можно провести сравнение двух поверхностей – существующей, полученной с помощью CIVIL3D, и поверхности, созданной на основе общедоступной цифровой модели местности с использованием методики имитационного моделирования в программном продукте BLENDER.
Рис. 5. Порядок построения цифровой модели местности в программе MetaShape Источник: составлено автором на основе результатов исследования
Рис. 6а. Ортофотоплан, ЖК «I Love» Источник: составлено автором на основе результатов исследования
Рис. 6б. Ортофотоплан, ЖК «Бирюлево» Источник: составлено автором на основе результатов исследования
Рис. 7а. Сравнение двух поверхностей по высоте Источник: составлено автором на основе результатов исследования
Рис. 7б. Точность модуля при увеличении количества геодезических опорных точек (ЖК «I Love») Источник: составлено автором на основе результатов исследования
В результате сравнения среднее квадратическое отклонение между существующей поверхностью и поверхностью, созданной на основе общедоступной цифровой модели местности, привязанной к четырем опорным точкам, составляет 1,6 метра. При увеличении количества опорных точек до восьми среднее квадратическое отклонение уменьшилось до 0,7 метра. Для поверхности общедоступной цифровой модели местности, привязанной к 20 опорным геодезическим точкам, средняя квадратическая ошибка составила 0,4 метра. Это вся информация по жилому комплексу «I Love». Все результаты сравнения представлены на рисунках 7а, 7б и в таблице 2.
При повторении такого сравнения для второго района исследования – жилого комплекса «Бирюлево» средняя квадратическая ошибка при сравнении между существующей поверхностью и поверхностью общедоступной цифровой модели местности, привязанной к четырем опорным точкам, составила 4,6 метра. При увеличении количества опорных точек до восьми среднее квадратическое отклонение уменьшилось до 1,6 метра. При привязке к 20 опорным геодезическим точкам оно снизилось до 1,3 метра. Результаты сравнения представлены на рисунках 8а, 8б и в таблице 3.
Таблица 2
Сравнение высот точек модели SRTM и точек топографического плана (ЖК «I Love»)
ЖК «I Love» |
4 опорные точки |
8 опорных точек |
20 опорных точек |
|||||||||
Номер точки |
Восточное положение |
Северное положение |
Отметка |
Отметка |
∆H |
Отметка |
∆H |
Отметка |
∆H |
|||
1 |
7988,050 |
16048,269 |
154,280 |
154,750 |
-0,470 |
1,487 |
154,530 |
-0,250 |
0,206 |
154,830 |
-0,550 |
0,033 |
2 |
7970,611 |
16040,365 |
155,360 |
154,470 |
0,890 |
2,541 |
154,420 |
0,940 |
2,703 |
155,000 |
0,360 |
1,191 |
3 |
8010,270 |
16041,162 |
154,000 |
155,430 |
-1,430 |
0,527 |
154,900 |
-0,900 |
0,038 |
154,750 |
-0,750 |
0,000 |
4 |
8014,127 |
16030,577 |
153,930 |
155,620 |
-1,690 |
0,972 |
154,960 |
-1,030 |
0,106 |
154,750 |
-0,820 |
0,008 |
5 |
7962,349 |
16027,380 |
155,450 |
154,220 |
1,230 |
3,740 |
154,460 |
0,990 |
2,870 |
155,110 |
0,340 |
1,148 |
6 |
8033,319 |
16018,391 |
154,060 |
156,040 |
-1,980 |
1,628 |
155,140 |
-1,080 |
0,141 |
154,810 |
-0,750 |
0,000 |
7 |
8003,732 |
16016,393 |
153,670 |
155,430 |
-1,760 |
1,115 |
154,890 |
-1,220 |
0,266 |
154,740 |
-1,070 |
0,115 |
8 |
7980,941 |
16015,794 |
154,070 |
155,030 |
-0,960 |
0,066 |
154,700 |
-0,630 |
0,005 |
155,010 |
-0,940 |
0,043 |
9 |
7954,352 |
16014,196 |
155,040 |
153,610 |
1,430 |
4,554 |
154,550 |
0,490 |
1,426 |
155,650 |
-0,610 |
0,015 |
10 |
7942,558 |
16001,184 |
155,380 |
154,420 |
0,960 |
2,769 |
154,610 |
0,770 |
2,173 |
156,180 |
-0,800 |
0,005 |
11 |
7975,062 |
16001,786 |
153,970 |
155,590 |
-1,620 |
0,839 |
154,830 |
-0,860 |
0,024 |
155,140 |
-1,170 |
0,192 |
12 |
8005,359 |
16002,187 |
153,780 |
155,720 |
-1,940 |
1,528 |
155,110 |
-1,330 |
0,392 |
154,840 |
-1,060 |
0,108 |
13 |
8032,445 |
16003,189 |
153,760 |
156,230 |
-2,470 |
3,119 |
155,390 |
-1,630 |
0,857 |
154,900 |
-1,140 |
0,167 |
14 |
8055,519 |
16003,390 |
153,930 |
156,760 |
-2,830 |
4,520 |
155,490 |
-1,560 |
0,733 |
155,120 |
-1,190 |
0,210 |
15 |
8072,625 |
15982,483 |
154,270 |
157,610 |
-3,340 |
6,948 |
155,960 |
-1,690 |
0,972 |
155,450 |
-1,180 |
0,201 |
Примечание: δ = 1,6 δ = 0,7 δ = 0,4.
Источник: составлено автором на основе результатов исследования.
Таблица 3
Сравнение высот точек модели SRTM и точек топографического плана (ЖК «Бирюлево»)
ЖК «Бирюлево» |
4 опорные точки |
8 опорных точек |
13 опорных точек |
|||||||||
Номер точки |
Восточное положение |
Северное положение |
Отметка |
Отметка |
∆H |
Отметка |
∆H |
Отметка |
∆H |
|||
1 |
-245689,23 |
-120709,838 |
697,72 |
696,33 |
1,39 |
17,25 |
700,11 |
-2,39 |
2,09 |
697,58 |
0,14 |
0,26 |
2 |
-244923,8 |
-120745,455 |
700 |
693,86 |
6,14 |
79,26 |
697,71 |
2,29 |
10,46 |
698,57 |
1,43 |
3,25 |
3 |
-245328,37 |
-120810,163 |
696,54 |
693,97 |
2,57 |
28,44 |
697,71 |
-1,17 |
0,05 |
695,54 |
1 |
1,88 |
4 |
-245691,25 |
-120805,006 |
700 |
697,65 |
2,35 |
26,14 |
700,55 |
-0,55 |
0,16 |
699,93 |
0,07 |
0,19 |
5 |
-245686,25 |
-120895,947 |
703,43 |
706,45 |
-3,02 |
0,07 |
704,45 |
-1,02 |
0,01 |
703,84 |
-0,41 |
0,00 |
6 |
-245325,39 |
-120927,733 |
695,89 |
697,21 |
-1,32 |
2,08 |
696,57 |
-0,68 |
0,07 |
694,73 |
1,16 |
2,35 |
7 |
-245331,69 |
-121017,699 |
696,88 |
700 |
-3,12 |
0,13 |
697,53 |
-0,65 |
0,09 |
697,05 |
-0,17 |
0,04 |
8 |
-245329,38 |
-121137,62 |
698,37 |
706,38 |
-8,01 |
27,53 |
701,72 |
-3,35 |
5,79 |
700,41 |
-2,04 |
2,78 |
9 |
-244895,95 |
-121171,393 |
696,19 |
699,38 |
-3,19 |
0,18 |
695,98 |
0,21 |
1,33 |
696,94 |
-0,75 |
0,14 |
10 |
-245327,24 |
-121271,315 |
700 |
707,46 |
-7,46 |
22,06 |
701,75 |
-1,75 |
0,65 |
700,99 |
-0,99 |
0,38 |
11 |
-244882,03 |
-121268,739 |
697,48 |
702,18 |
-4,7 |
3,75 |
696,22 |
1,26 |
4,86 |
698,44 |
-0,96 |
0,35 |
12 |
-244881,29 |
-121368,043 |
698,89 |
707,28 |
-8,39 |
31,66 |
700,05 |
-1,16 |
0,05 |
698,55 |
0,34 |
0,51 |
13 |
-245317,29 |
-121356,055 |
699,21 |
708,37 |
-9,16 |
40,92 |
702,53 |
-3,32 |
5,64 |
702,86 |
-3,65 |
10,75 |
Примечание: δ = 4,6 1,6 δ = 1,3.
Источник: составлено автором на основе результатов исследования.
Рис. 8а. Сравнение двух поверхностей по высоте Источник: составлено автором на основе результатов исследования
Рис. 8б. Точность модуля при увеличении количества геодезических опорных точек (ЖК «Бирюлево») Источник: составлено автором на основе результатов исследования.
Выводы
1. В исследовании выполнена оценка точности цифровой модели рельефа, которая демонстрирует, что средняя квадратическая погрешность высот понижается за счет добавления опорных точек.
2. Из исследования видно, что цифровая модель местности лучше зарекомендовала себя на «открытых» равнинных участках, чем на «закрытых» участках или на участках с уклоном.
3. Рекомендуется использовать цифровую модель местности для проектирования только после добавления минимально 12 и более опорных точек в зависимости от уровня и свойства сложности рельефа.
4. Цифровая модель рельефа SRTM является важным бесплатным исходным материалом для многих исследований, но она по-прежнему содержит значительные высотные погрешности.
Библиографическая ссылка
Насереддин Х.Х., Алексашина Е.В., Акель М.А. ПОВЫШЕНИЕ ТОЧНОСТИ ОБЩЕДОСТУПНОЙ ЦИФРОВОЙ МОДЕЛИ МЕСТНОСТИ С ИСПОЛЬЗОВАНИЕМ ИМИТАЦИОННОГО МОДЕЛИРОВАНИЯ АЭРОФОТОТОПОГРАФИЧЕСКОЙ СЬЕМКИ В ПРОГРАММНОМ ОБЕСПЕЧЕНИИ BLENDER // Успехи современного естествознания. 2025. № 6. С. 71-83;URL: https://natural-sciences.ru/ru/article/view?id=38405 (дата обращения: 03.07.2025).
DOI: https://doi.org/10.17513/use.38405