Выявление устойчивых закономерностей распределения земельных участков в охранной зоне малой реки Ировка Республики Марий Эл по берегам от истока до устья было выполнено по классам земельного покрова ООН [1].
Целью статьи является установление по спутниковым снимкам, к какому классу земельного покрова относятся земельные участки, находящиеся по обеим берегам русла малой реки Ировка в характерных точках от истока до устья. Затем проводили вейвлет анализ [2] и сравнение обоих берегов методом факторного анализа [3] по закономерностям влияния ранга 290 характерных точек на стрежне малой реки на изменение рангов классов земельного покрова ООН.
Река Ировка протекает в Республике Марий Эл. Исток Ировка берет 3 км. на северо-запад от деревни Ирмучаш, находящейся в 5 км. от п. Параньга (рис. 1). Устье реки расположено в 136 км. по правому берегу реки Илеть. Длина реки Ировка составляет 69 км, а площадь водосборного бассейна равна 917 км2.
Рис. 1. Малая река Ировка на территории Республики Марий Эл
Вдоль реки от истока до устья выбираются определенные характерные точки. Как правило, ими становятся места в горизонтальных изломах реки, где резко изменяется кривизна русла в плане. Высота на снимке показывает средний рельеф обоих берегов.
Земельный кадастр РФ отличается по классификации с ООН по методу GAEZ [1, 4-5]. По данным [1] в [5] относительно 11 классов земельного покрова были установлены следующие ранги, соответствующие ухудшению качества поверхности земель (чем больше ранг, тем хуже): 0 – травяной покров (сенокосы + пастбища);
1 – древесно-кустарниковая растительность, ДКР (многолетние насаждения + лесные насаждения, не включающиеся в лесной фонд);
2 – лес (лесные участки);
3 – мозаика растительности с включением сельхозугодий;
4 – мозаика растительного покрова с доминированием сельхозугодий; 5 – земли для растениеводства, пашни;
6 – мокрые, увлажненные земли (болота);
7 – бесплодная земля;
8 – вода, побережье (земли под водой);
9 – лед, холодная пустыня;
10 – городская территория.
В таблице 1 приведен фрагмент расставленных нами рангов классов земельного покрова по правой и левой сторонам малой реки Ировка.
Таблица 1. Распределений рангов классов земельного покрова реки Ировка
Ранг точки стрежня реки |
Ранг земельного покрова (справа) |
Примечание |
Ранг земельного покрова (слева) |
Примечание |
0 |
0 |
Травяной покров |
0 |
Травяной покров |
1 |
0 |
Травяной покров |
3 |
Мозаика с сельхозугодиями |
2 |
0 |
Травяной покров |
3 |
Мозаика с сельхозугодиями |
… |
… |
… |
… |
… |
287 |
2 |
Лесной участок |
2 |
Лесной участок |
288 |
0 |
Травяной покров |
2 |
Лесной участок |
289 |
0 |
Травяной покров |
0 |
Травяной покров |
290 |
0 |
Травяной покров |
0 |
Травяной покров |
Целью метода Галеркина при моделировании статистических данных является получение системы линейных уравнений. для этого внасчале строятся системы дифференциально-интегральных уравнений. В методе идентификации [2] мы не строим исходные дифференциальные уравнения, а сразу же работаем с явно нелинейными уравнениями неизвестных дифференциальных и интегральных уравнений.
Каждый параметр конечного решения имеет физический смысл, поэтому колебания (вейвлет сигналы) записываются волновой формулой [2] вида
,
,
, (1)
где – показатель (зависимый фактор),
– номер составляющей модели (1),
– количество членов в модели (1),
– объясняющая переменная (влияющий фактор),
– параметры модели (1), принимающие числовые значения в ходе структурно-параметрической идентификации (табл. 2) в программной среде CurveExpert-1.40 (URL: http://www.curveexpert.net/),
– амплитуда (половина) вейвлета (ось
),
– полупериод колебания (ось
). По формуле (1) реализуется метод Декарта о наличии одного единственного алгебраического уравнения, членами которого становятся асимметричные вейвлет сигналы как кирпичики Гильберта.
Таблица 2. Параметры распределений рангов классов земельного покрова
Номер |
Вейвлет |
Коэф. корр. |
|||||||
Амплитуда (половина) колебания |
Полупериод колебания |
Сдвиг |
|||||||
|
|
|
|
|
|
|
|
||
Асимметричные вейвлет сигналы вдоль левого берега реки Ировка |
|||||||||
1 |
0.21053 |
0 |
-0.74350 |
0.16305 |
0 |
0 |
0 |
0 |
0.3879 |
2 |
5.06149e-34 |
22.49839 |
0.0034692 |
2.06191 |
0 |
0 |
0 |
0 |
|
3 |
-1.23182e-6 |
0 |
-12.70689 |
0.0077343 |
37.76789 |
-1.42205 |
0.33750 |
-2.09993 |
|
4 |
-0.0015381 |
0 |
-3.53006 |
0.10505 |
57.66170 |
2.55529e-5 |
1.97040 |
0.98736 |
|
… |
… |
… |
… |
… |
… |
… |
… |
… |
… |
Асимметричные вейвлет сигналы вдоль правого берега реки Ировка |
|||||||||
1 |
0.52580 |
0 |
0.028292 |
0 |
0 |
0 |
0 |
0 |
0.3212 |
2 |
0.21124 |
0.43034 |
1.35536e-4 |
1.49916 |
0 |
0 |
0 |
0 |
|
3 |
-0.54333 |
0 |
-0.00232518 |
0 |
19.09961 |
2.29669e-4 |
1.66878 |
1.12722 |
|
4 |
-0.417958 |
0 |
-0.00215509 |
0 |
8.48308 |
6.30657e-4 |
1.33931 |
4.88727 |
|
… |
… |
… |
… |
… |
… |
… |
… |
… |
… |
Наилучшую закономерность (рис. 2 и рис. 3) с теснотой связи по коэффициенту корреляции 0,2465 имеют земли вдоль левого берега. Она не отличается от правого берега, здесь коэффициент корреляции составляет 0,2406.
|
|
Двухчленный тренд |
Первое колебание |
|
|
Второе колебание |
Тренд и два колебания |
Рис. 1. Графики моделей распределений рангов классов земельного покрова характерных точек от истока до устья слева реки Ировка
Всего было получено для обоих берегов более 60 составляющих. Поэтому адекватность моделей повышается и достигает среднего уровня.
Растительный покров по берегам вдоль реки Ировка показал, что слева реки земельный покров больше составляют лесные участки, а справа преобладает травяной покров, но лесные площадки также присутствуют.
|
|
Двухчленный тренд |
Первое колебание |
|
|
Второе колебание |
Тренд и два колебания |
Рис. 2. Графики моделей распределений рангов классов земельного покрова характерных точек от истока до устья справа реки Ировка
Сущностью метода факторного анализа состоит в выделении из множества изучаемых факторов, влияющих на изучаемый объект, меньшего их числа, но выявляющий более существенные свойства исследуемого явления.
Корреляционный анализ включает в себя изучение связей не двух, а множества переменных, измеренных в количественной шкале на одной выборке. В этом случае вычисляются корреляции для каждой пары из этого множества переменных. Вычисления обычно проводятся на компьютере, а результатом является корреляционная матрица.
Оценим возможность применения факторного анализа и технику построения корреляционной матрицы для рангов классов земельного покрова по берегам реки Ировка справа и слева. Результаты зависимостей представлены в таблице 3.
Таблица 3. Корреляционная матрица факторного анализа и рейтинг рангов классов земельного покрова по берегам малой реки Ировка
Влияющие факторы (параметры |
Зависимые факторы (показатели |
Сумма |
Место |
|
Ранги КЗП справа |
Ранги КЗП слева |
|||
Ранги КЗП справа |
0.3212 |
0.1961 |
0.5173 |
2 |
Ранги КЗП слева |
0.2584 |
0.3879 |
0.6453 |
1 |
Сумма коэффиц. корр. |
0.5796 |
0.5840 |
1.1626 |
- |
Место |
2 |
1 |
- |
0.2907 |
В строках таблицы представлены ранги классов земельного покрова справа и слева реки Ировка, а в столбцах - эти же факторы как зависимые параметры.
В нашем случае класс земельного покрова как влияющая переменная получает коэффициент корреляции справа и слева реки соответственно 0,5173 и 0,6453. А как зависимый показатель первое место занимает также ранг класса земельного покрова ООН слева.
На основании влияющих и зависимых факторов был рассчитан коэффициент коррелятивной вариации, который применяется для сравнения разных объектов исследования. Он определяется по соотношению 1.1626 / 22 = 0.2907.
Проведенные измерения и их математический анализ показали, разные классы земельного покрова ООН вдоль реки Ировка по обоим её берегам имеют четкие закономерности по асимметричным вейвлет сигналам.