Процесс т.н. "обучения" представляет
из себя следующее: пользователь "обучает" программу, на какие классы
необходимо разбить изображение, самостоятельно выделив полигоны с определенной
цветовой палитрой. Например, в нашем случае, выделялись части леса, поля,
урбанизированные территории, какие-либо водные объекты. Далее программа
самостоятельно определяет, какую часть изображения занимает тот или иной класс.
Получив данные о площадях, занимаемых лесом, городом, полем, озерами, находим долю
(в процентах) от общей площади бассейна. Результаты по каждому водосбору
представлены в главе "Расчет".
3. Расчеты и анализ результатов
.1 Общее описание существующих гидрологических
характеристик
Термин "гидрографические характеристики" определяется как совокупность морфометрических и морфологических характеристик водных объектов и их водосборов, дающих достаточно полное представление о характере, форме, размерах, протяженности водных объектов и некоторых физико-географических особенностях их водосборов.
Морфометрические характеристики представляют
собой количественные показатели водных объектов и водосборов, а морфологические
- качественно-количественные показатели строения поверхности водосборов.
.2 Определение морфометрических характеристик
Площадь водосбора. В представленной ниже таблице
приведены результаты расчетов площади водосбора для каждого выбранного поста. В
первых столбцах представлена общая информация о постах, их местоположение. В
столбце 6 представлен результат расчета, описание которого приведено в разделе
2.4 главы 2. В данном столбце отображена площадь водосбора в начальных единицах
растра - пикселях, иными словами, этот столбец показывает, сколько ячеек
цифровой модели занимает тот или иной бассейн.
В нашем случае, разрешение модели составляет 90
метров, т.е. один пиксель на местности будет занимать квадратный полигон со
стороной 90 метров. Но, ввиду того, что исходный растр претерпел
перепроецирование в прямоугольную систему координат (Пулково 1942, зона 7),
начальное разрешение изменилось и составляет 74 метра. Таким образом, можно
легко вычислить, какую площадь в км2 занимает один пиксель, а именно
0,074х0,074=0,005477 км2. Теперь можно перейти от площади в пикселях
к площади в квадратных километрах, просто перемножив значения в столбцах 6 и 7.
Результат расчета площади водосбора в км2 представлен в столбце 8. В
столбце 9 приведены данные из ОГХ 1970-1975 гг., собственно с этими данными и
производится сравнение, анализ, расчет относительной погрешности.
Таблица
№
река
пост
Широта
Долгота
S
в пикселях
S
пикселя
Площадь
расчетная
Площадь
отн.
погр.
1
р.Ока
д.Вендерево
52,533
36,017
92119
0,005477
504,58
513
1,64
2
р.Ока
д.Семенково
52,572
35,951
164127
0,005477
899,00
912
1,43
3
р.Крома
сл.Черкасская
52,7
35,8
156545
0,005477
857,47
853
0,52
4
р.Рыбница
д.Любаново
52,845
36,11
125169
0,005477
685,61
709
3,30
5
р.Цоп
с.Новолуние
52,9
35,783
125944
0,005477
689,86
689
0,12
6
р.Нерунь
д.Орловка
52,933
36,917
274830
0,005477
1505,38
1500
0,36
7
р.Ока
г.Орел
52,967
36,083
893420
0,005477
4893,69
4890
0,08
8
р.Орлик
г.Орел
52,967
36,067
99717
0,005477
546,20
544
0,40
9
р.Зуша
г.Новосиль
52,967
37,083
358503
0,005477
1963,69
1970
0,32
10
р.Ока
д.Костомарово
52,983
36,083
896156
0,005477
4908,68
4900
0,18
11
р.Орлица
д.Большой
Рог
52,989
35,69
18552
0,005477
101,62
96
5,85
12
р.Раковка
д.Ратановка
53,04
37,256
67892
0,005477
371,88
364
2,16
13
р.Оптуха
с.Платово
53,1
36,233
138566
0,005477
758,99
620
22,42
14
р.Грунец
д.Гагаринские
Хутора
53,133
37,214
9391
0,005477
51,44
50
2,88
15
р.Нугрь
д.Городище
53,269
35,675
32039
0,005477
175,49
175
0,28
16
р.Вытебеть
с.Льгов
53,283
35,413
64448
0,005477
353,01
338
4,44
17
руч.Ракитня
д.Мымрино
53,297
35,489
11240
0,005477
61,57
63,3
2,74
18
р.Рыдань
д.Камынина
53,323
35,749
22168
0,005477
121,42
120
1,19
19 с.Локно
53,332
35,796
68376
0,005477
374,53
359
4,33
20
р.Шковка
д.Нижняя
Шкава
53,353
35,439
10922
0,005477
59,83
58
3,15
21
р.Снытка
д.Малая
Кутьма
53,374
36,013
26949
0,005477
147,61
151
2,24
22
р.Зуша
г.Мценск
53,417
36,533
1147339
0,005477
6284,53
6000
4,74
23
р.Нугрь
г.Болхов
53,433
36,017
184737
0,005477
1011,89
1010
0,19
24
р.Орс
д.Орс
53,434
35,947
28272
0,005477
154,86
155
0,09
25
р.Вытебеть
д.Шваново
53,523
35,59
168400
0,005477
922,41
920
0,26
26
р.Иста
с.Араны
53,585
36,477
99680
0,005477
546,00
548
0,37
27
р.Плава
г.Плавск
53,712
37,291
149199
0,005477
817,23
810
0,89
28
р.Ока
г.Белев
53,817
36,15
3205649
0,005477
17558,88
17500
0,34
29
р.Жиздра
с.Дубровка
53,883
35,083
345895
0,005477
1894,63
1900
0,28
30
р.Упа
г.Крапивна
53,946
37,157
1163547
0,005477
6373,31
6320
0,84
31
р.Упа
с.Орлово
53,95
37,15
1503471
0,005477
8235,23
8210
0,31
32
р.Солова
д.Захаровка
53,967
37,45
71567
0,005477
392,01
400
2,00
33
р.Жиздра
г.Козельск
54,033
35,8
1273786
0,005477
6977,14
6940
0,54
34
р.Ока
г.Лихвин(Чекалин)
54,096
36,251
5076580
0,005477
27806,87
27900
0,33
35
р.Упа
г.Тула
54,2
37,65
649499
0,005477
3557,62
3550
0,21
36
р.Упа
с.Маслово
54,207
37,508
737302
0,005477
4038,56
4040
0,04
37
р.Ока
с.Голодское(Перемышль)
54,29
36,188
6953884
0,005477
38089,76
38000
0,24
38
р.Осетр
д.Хрусловка
54,397
38,228
83477
0,005477
457,24
462
1,03
39
р.Ока
г.Калуга
54,5
36,267
10057890
0,005477
55091,89
54900
0,35
40
р.Безвель
д.Утешево
54,515
35,681
24903
0,005477
136,41
146
6,57
41
р.Угра
с.Всходы
54,67
34,099
353825
0,005477
1938,07
1890
2,54
42
р.Ресса
д.Лазино
54,681
35,053
441738
0,005477
2419,61
2440
0,84
43
р.Угра
пгш
Товарково
54,683
35,933
2801350
0,005477
15344,34
15300
0,29
44
р.Большое
Смедово
д.Острога
54,697
38,441
43721
0,005477
239,48
241
0,63
45
р.Таруса
д.Похвинсево
54,733
37,1
156158
0,005477
855,35
872
1,91
46
р.Угра
д.Мокрое(г.Юхнов)
54,746
35,223
1970645
0,005477
10794,17
10800
0,05
47
р.Ока
г.Кашира
54,833
38,15
12607157
0,005477
69055,45
68700
0,52
48
р.Путынка
д.Малахово
54,833
36,283
29275
0,005477
160,35
153
4,81
49
р.Лопасня
д.Кубасово
54,981
37,836
166381
0,005477
911,35
923
1,26
50
р.Протва
с.Спас-Загорье
55,05
36,6
643880
0,005477
3526,84
3640
3,11
51
р.Угра
д.Мухино(Комаровка)
55,052
34,659
888408
0,005477
4866,24
5790
15,95
52
р.Воря
свх
Ильинский
55,112
35,06
141020
0,005477
772,43
825
6,37
53
р.Северка
с.Покровское
55,2
38,317
146081
0,005477
800,16
792
1,03
54
р.Нара 55,207
36,943
171173
0,005477
937,60
971
3,44
55
р.Северка
с.Большое
Алексеевское
55,242
38,208
114096
0,005477
624,96
642
2,65
56
р.Протва
г.Верея
55,35
36,183
145474
0,005477
796,83
935
14,78
57
р.Нара
г.Наро-Фоминск
55,383
36,75
114792
0,005477
628,77
665
5,45
58
р.Пахра
д.Макарово
55,465
37,673
321699
0,005477
1762,10
1760
0,12
59
р.Пахра
Стрелковская
фабрика
55,467
37,633
309405
0,005477
1694,76
1690
0,28
60
р.Москва
пос.Марфин
Брод
55,533
35,967
249214
0,005477
1365,06
1380
1,08
61
вдхр.Можайское
Можайский
гидроузел
55,534
35,952
248729
0,005477
1362,41
1380
1,27
62
р.Воинка
д.Беззубово
55,546
35,803
8894
0,005477
48,72
64,1
24,00
63
р.Сежа
д.Баскаково
55,547
34,754
58964
0,005477
322,97
291
10,99
64
р.Гжать
г.Гагарин
55,55
34,983
80735
0,005477
442,22
406
8,92
65
р.Москва
д.Макарово
55,552
36,027
274023
0,005477
1500,96
1420
5,70
66
р.Гжать
г.Гжатск
55,553
35
82764
0,005477
453,34
418
8,45
67
р.Гжать
д.Кормино(Кормино-Шеляги)
55,583
34,987
93590
0,005477
512,64
484
5,92
68
р.Сежа
д.Новое
Село
55,589
34,673
72013
0,005477
394,45
375
5,19
69
р.Искона
д.Новинки
55,633
36,033
87636
0,005477
480,02
472
1,70
70
р.Лусянка
д.Черники
55,65
35,667
29099
0,005477
159,39
170
6,24
71
р.Москва
Перервинский
гидроузел
55,667
37,717
1654883
0,005477
9064,59
8850
2,42
72
р.Вазуза
с.Дугино
55,683
34,2
128450
0,005477
703,58
734
4,14
73
р.Москва
д.Барсуки
55,683
35,6
139816
0,005477
765,84
755
1,44
74
р.Руза
г.Руза
55,699
36,189
347021
0,005477
1900,80
1930
1,51
75
р.Москва
г.Звенигород
55,733
36,867
947372
0,005477
5189,21
5000
3,78
76
р.Москва
с.Петро-Дальнее
55,733
37,167
1367533
0,005477
7490,63
7320
2,33
77
р.Москва
Бабьегородская
плотина
55,737
37,605
1537289
0,005477
8420,47
8170
3,07
78
р.Руза
д.Покров
55,75
36,133
206760
0,005477
1132,52
1140
0,66
79
р.Яуза
Яузский
гидроузел
55,75
37,683
82087
0,005477
449,63
449
0,14
80
р.Озерна
д.Карповка
55,762
36,263
101540
0,005477
556,18
483
15,15
81
р.Касня
с.Субботники
55,802
34,426
261253
0,005477
1431,01
1420
0,78
82
р.Москва
с.Павшино
55,802
37,333
1408495
0,005477
7715,00
7550
2,19
83
р.Гжать
д.Шакуты
55,81
34,706
366828
0,005477
2009,29
1900
5,75
84
р.Истра
с.Павловская
Слобода
55,817
37,1
352723
0,005477
1932,03
1950
0,92
85
р.Озерна
д.Городище
55,85
36,35
67357
0,005477
368,95
364
1,36
86
р.Лихоборка
г.Москва,свх
Марфино
55,854
37,559
7910
0,005477
43,33
43
0,76
87
р.Малая
Истра
д.Киселево
55,883
36,767
53327
0,005477
292,10
280
4,32
88
р.Волошня
д.Чертаново
55,9
35,917
18099
0,005477
99,14
91,2
8,70
89
р.Руза
д.Красное
Село 35,518
63418
0,005477
347,37
339
2,47
90
р.Истра
с.Бужарово
55,98
36,802
187490
0,005477
1026,97
1010
1,68
91
р.Воря
с.Мишнево
55,983
38,217
179407
0,005477
982,70
947
3,77
92
р.Воря
с.Каблуково
56,025
38,17
160828
0,005477
880,93
901
2,23
93
р.Вазуза
д.Золотилово
56,043
34,493
1039531
0,005477
5694,01
5510
3,34
94
р.Осуга
д.Коротнево(Кортнево)
56,052
34,435
193779
0,005477
1061,42
1230
13,71
95
р.Осуга
д.Красное
56,089
34,012
57720
0,005477
316,16
310
1,99
96
р.Нудоль
д.Кузнецово
56,127
36,656
53217
0,005477
291,50
291
0,17
97
р.Шоша
с.Хлопово-Городище
56,303
35,426
61863
0,005477
338,85
339
0,04
Расчет относительной погрешности производился по
формуле:
На графике зависимости относительной погрешности
от площади водосбора отчетливо заметна тенденция к снижению погрешности с
ростом площади водосбора. Таким образом, принимая во внимание также и поведение
линии тренда на графике, можно сказать, что в большинстве случаев относительная
погрешность не превышает 5-7 %, исключая единичные случаи, связанные, возможно,
с неточностями самой цифровой модели рельефа, отсутствующими данными и интерполированными
участками. Также стоит учесть, что вычисление малых площадей водосбора,
используя модель с 90-метровым разрешением, будет всегда сопровождаться
большими погрешностями. Модель чисто физически не может позволить практически
безошибочно оконтурить малый водосбор, учитывая то, что сторона минимальной
единицы - пикселя, составляет 90 метров. В то же время, это мало сказывается на
расчетах больших площадей. Погрешность, конечно же, также присутствует, но на
общем фоне полученного результата они мало заметны, и составляют относительно
малую долю.
Рассматривая результаты более детально, можно
выявить следующее: относительная погрешность в подавляющем большинстве случаев
не превышает 10% для бассейнов с площадью до 1500 км2, для бассейнов
с площадью от 1,5 до 7,5 тысяч км2 колеблется в пределах 5-7%, для
более крупных бассейнов погрешность расчетов не превышает 1-2%.
Уклон. Средняя высота. В приведенной ниже
таблице, помимо гидрографических данных, приведены результаты расчета средней
высоты и уклона водосбора. Череда необходимых для расчета операций подробно
описана в разделе 2.4 главы 2.
Итак, жирным шрифтом в таблице выделены
рассчитанные данные в среде Arcgis.
По возможности, при наличии соответствующих данных из ОГХ, результаты
сравнивались и анализировались, рассчитывалась относительная погрешность
расчетов. Таким образом, расчет погрешности при определении средней высоты
водосбора производился аналогичным образом с расчетом погрешности при
определении площади водосбора. Ее значения представлены в столбце 6. Как можно
заметить, погрешность не превышает 5%, причем наблюдается тенденция к снижению
погрешности с ростом площади водосбора. Рассматривая конкретнее, можно выделить
некоторые условные интервалы, в которых колеблется значение погрешности в
зависимости от площади водосбора. Площади до 2000 км2 соответствует
погрешность 4-5%, от 2000 до 5000 тысяч - 2-4%, более крупные бассейны - менее
2%.
Что касается уклонов, то здесь не было
возможности сравнить результаты, ввиду того, что в ОГХ 1970-1975 гг. не
оказалось данных о уклонах бассейнов, были лишь средние уклоны собственно
водотоков. Все значения уклонов приведены в таблице в столбцах 2 и 3. Жирным
шрифтом выделен результат расчета среднего уклона водосбора.
.3 Определение гидрографические характеристик
Конкретно в ходе работы рассматривались
залесенность, озерность и доля урбанизированных территорий. Результаты также
представлены в таблице. Жирным шрифтом выделены данные, полученные в среде Arcgis
посредством расчетов, остальное - данные из ОГХ. Озерность представлена в
столбце 7. Как мы можем заметить, то значения практически идентичны справочным
данным. В случае с гидрографическими характеристиками не имеет смысла
рассчитывать какие-либо погрешности, ввиду того, что с течением времени (а это
практически 40 лет) они могли естественным образом измениться, и сравнивать
данные со снимков 1-2 годичной давности с данными 1975 года абсолютно
бессмысленно. Говоря конкретно об озерности, то можно заметить, что доля озер в
общей картине ландшафтов по прошествии 40 лет мало изменилась.
Что касается залесенности территории, то здесь
можно заметить, как ни странно, увеличение данного показателя в большинстве
случаев (столбец 8). Обратная ситуация только лишь с бассейном реки Воря,
расположенным на северо-востоке Московской области. Здесь залесенность
сократилась с 72% до 42%, что наверняка связано с хозяйственным развитием
территории, ростом урбанизации (кстати, доля урбанизированных территорий здесь
максимальна среди рассматриваемых бассейнов). В остальных же случаях доля леса
либо осталась практически неизменной, либо даже возросла. Например, на 2% в
бассейне Зуши, Оки до д. Костомарово, и на 19% в пределах водосбора реки Угры.
Данное явление, возможно, связано с практически полным прекращением
использования дерева в качестве энергетического ресурса и переходом на
углеводородное сырье, что сократило огромные вырубки леса.
Также рассматривалась доля урбанизированной
территории в общей картине ландшафтов (столбец 9). Но, к сожалению, ввиду того,
что справочных данных по этому показателю нет, сравнить их не с чем, можно лишь
дать некоторое пояснение к имеющимся значениям. Максимальное долю
урбанизированных территорий мы наблюдаем в пределах реки Вори (37%), как уже
было упомянуто выше. Это связано, собственно, с географическим положением
бассейна, лежащего в пределах крупнейшей в стране городской агломерации.
Наоборот же, урбанизированные территории полностью отсутствуют на водосборе
реки Угры (до пос. Всходы), что, в свою очередь, опять же связано с его
положением: до поселка Всходы бассейн Угры относительно удален от крупных
административных центров, а в его пределах находятся только небольшие деревни и
сёла Смоленской и Калужской областей.
Что касается распаханности и заболоченности
территорий, то данные расчеты не удалось произвести по одной простой причине -
на спутниковых снимках с трудом можно различить границы пашни в общей площади
полей, а также сложно выделить отдельные болота.
Заключение
Данная курсовая работа по сути является более
углубленным вариантом курсовой работы 2-го курса с более расширенным спектром
задач. В ходе работы была дана оценка точности определения морфометрических
характеристик водосбора с помощью цифровых моделей рельефа, а также оценена
возможность применения настольного ПО вкупе со спутниковыми данными для
определения гидрографических характеристик. Все поставленные задачи были
выполнены. Произведен разбор и анализ использования цифровых моделей рельефа
для определения морфометрических характеристик водосбора, в частности его
площади, уклона, средней высоты, а также произведен анализ использования
современного ПО для определения гидрографических характеристик водосбора.
Выполнив данную курсовую работу, можно смело заявить, что компьютерные возможности
вкупе со спутниковыми данными позволяют решать в частности и проблемы
гидрологии, и существенно облегчают и ускоряют процессы определения
морфометрических и гидрографических характеристик водных объектов.