CC BY
14
УДК 528.852
ФИЛЬТРАЦИЯ И СЖАТИЕ ГЕОИНФОРМАЦИИ
С ИСПОЛЬЗОВАНИЕМ ТЕОРИИ ВЕЙВЛЕТОВ
А. С. ЯРМОЛЕНКО
Новгородский государственный университет им. Ярослава Мудрого, г. Великий Новгород, Россия
О. В. СКОБЕНКО
УО «Белорусская государственная сельскохозяйственная академия», г. Горки, Беларусь, 213407 e-mail: [email protected]
(Поступила в редакцию 23.03.2018)
Построены ортонормированные базисы вейвлет-преобразований и фильтрации информации. Рассмотрен порядок создания вейвлет-фильтров последовательностью сверток, применением КМА-анализа для построения ортонормированного базиса вейвлет-преобразования. Разработана методика определения значения фильтра в зависимости от точности исходной геоинформации. Реализована практическая возможность вейвлет-фильтрации на основе составленных конкретных программ моделирования полей данных геинформации и изображений, сжатия данных и их фильтрации.
Фильтрация геоинформации необходима при определении отклонений показателей определенного явления от их средних или заданных (нормативных) значений при мониторинге измерений различных явлений. Сжатие геоинформации и ее фильтрация является актуальной проблемой теории математической обработки как геодезической информации (измерений), так и изображений. При этом необходимо максимально использовать полученную информацию и получать конечные результаты с достаточной точностью и минимальными затратами при ее хранении, что связано со сжатием информации. Это важно при работе с геоинформацией в геодезии, землеустройствеt природообустройстве, мониторинге земель, при ведении точного сельского хозяйства [1]. При этом алгоритмы должны обладать простотой и точностью вычислений.
Ключевые слова: базис вейвлет-преобразования, свертка, КМ- анализ, фильтрация, сжатие, точность, программа, изображение, поле данных.
We have constructed orthonormal bases of wavelet transforms and information filtration. We have examined the order of creation of wavelet filters by convolutional sequence, application of MRA-analysis for construction of orthonormal basis of wavelet transform. We have developed a technique for determining the filter value depending on the accuracy of the initial geo-information. We have realized the practical possibility of wavelet filtration on the basis of compiled specific programs for modeling the geoinformation data fields and images, data compression and theirfiltering.
The filtration of geo-information is necessary for determining the deviations of indicators of a certain phenomenon from their mean or specified (normative) values when monitoring the measurements of various phenomena. Compression of geo-information and its filtration is an actual problem in the theory of mathematical processing of both geodetic information (measurements) and images. At the same time, it is necessary to make maximum use of the obtained information and to get final results with sufficient accuracy and minimum costs when storing it, which is associated with information compression. This is important when working with geo-information in geodesy, land management, environmental monitoring, land monitoring, and while maintaining accurate agriculture. In this case, algorithms must have the simplicity and accuracy of calculations.
Key words: basis of wavelet transform, convolution, MRA analysis, filtration, compression, accuracy, program, image, data field.
Введение
К настоящему времени для моделирования объектов, сжатия и фильтрации информации широко применяются преобразования Фурье [2-7]. Однако даже в своем быстром варианте (быстрое преобразование Фурье - БПФ) оно сопряжено с большим количеством вычислений. В свою очередь в работах [6], [2] отмечается, что в отличие от преобразований Фурье возможны преобразования с другими базисами, восстанавливающие дискретные и непрерывные функции, но значительно сокращающие вычисления. Одним из таких базисов является вейвлетный. И к настоящему времени он постепенно находит практическое применение [8-19]. Следует при этом отметить, что работы [8, 10,11,15,16] носят лишь ознакомительный характер, в основательной работе [9] речь о применении вейвлетов не ведется, но детально описывается большинство фильтров подавления шумов, которые могут применяться при моделировании процессов на основе вейвлетного базиса. В трудах [12-14] приводятся основные теоретические положения по теории вейвлетов, базирующиеся на основе зарубежных исследований. И для создания технологий обработки геоинформации на основе этих работ необходимы дополнительные исследования. Работами [17-19] выполнены значительные исследования по описанию вейв-летами гравитационного поля Земли. В [20] вейвлет-преобразовании используются зарубежные пакеты при неизвестном алгоритме. Перечисленные труды имеют практическое и теоретическое значение при обработке геоинформации, но в них теория вейвлетов не доведена до инженерного уровня их
применения. На основе отмеченного в данной статье с целью разработки детальной и доступной технологии применения вейвлетов в обработке геоинформации ставятся следующие задачи; 1) построение ортонормированных базисов вейвлет-преобразований и фильтрация; 2) порядок создания вейвлет-фильтров; 3) построение фильтров последовательностью сверток; 4) применение КМА-анализа для построения ортонормированного базиса вейвлет-преобразования; 5) показать практическую возможность вейвлет-фильтрации на основе составленных конкретных программ моделирования полей данных геоинформации, сжатия данных и их фильтрации.
Основная часть
В соответствии с [2] дадим следующие определения.
Оператором циклического сдвига последовательности Z [2, стр. 125] на k позиций вправо является оператор Rk , осуществляющий создание новой последовательности RkZ по формуле:
{RkZ\n) = Z(n-к) . (1)
где n - номер элемента в создаваемой последовательности.
Согласно определения 3.7 и теоремы 3.8 [2] при некоторых заданных векторах U и v, принадлежащих тому же пространству элементов, что и вектор Z, возможно построение ортонормированного базиса вида:
B = R^M-,1 U {R^Ut:* = V, Ry> R4V,..., Rn-2V,U,RU, ru,...,Rn_2U}. (2)
где u, согласно [2] можно назвать отцовским вейвлетом, а v-материнским, символ ^ логического объединения множеств, M = N /2. В данной статье составляющие базиса, построенные по вектору U будем называть отцовскими, а по вектору v-материнскими.
Ортонормированность (2) возможна тогда лишь [2], когда система:
( U(n) V (n) )
W W (3)
аАпП V2 |U(n+M) V(n+M)
для п=0,1,2,...М-1 - унитарна.
Под унитарной понимается такая матрица [2], для которой
А-1 = А* (4)
где А-1 обратная к А матрица, а А - сопряженная к А матрица, получаемая взятием комплексно-сопряженных значений от всех элементов матрицы Ат, транспонированной к А.
Построение ортонормированных базисов вейвлет-преобразований и фильтрация. В качестве примера построения ортонормированного базиса примем векторы Хаара [2, 6].
и = 1Д/2Д/л/2Д...,о)Г (5)
V = (/-/2-¡/-Яд...^ (6)
Поскольку в [2] установлено, что матрицы А(п) (3) для векторов и и V унитарны, то пользуясь правилом (2) построим ортонормированный базис для вейвлет-преобразования вектора 2. Здесь N=4 и М=М2=2. Тогда в соответствии с (2),(1) на основании (5),(6) можно составить ортонормированный вейвлет-базис
(1 о 1 оЛ
B=(V ,R2V ,u , R2U Н
л/2
-10 10 0 10 1 0 -10 1
(7)
Тогда вейвлет-преобразование на данном этапе будет
2 = В * 2 . (8)
а обратное преобразование, следуя (8), имеет вид
2 = В12 = В * 2 = Вт2. (9)
Теперь осуществим фильтрацию сигнала, то есть его разложение по отцовским и материнским составляющим вейвлет базиса. Согласно [ 2, 6] в пространстве с комплексным скалярным произведением < , > и ортонормированным базисом К = {"1=«2..... "Л для любого у этого пр°странства справедлив°:
у=Т=Лу-и1>1. (10)
Здесь у соответствует вектору данных 2, а базисные векторы "} - всем столбцам матрицы (7). После подстановки их в (1 0) получим:
г=ф)+р(г ). (11)
Следуя [2] составляющую высоких частот сигнала представим в виде формулы:
) = (1Л2кУ)Л2кУ АгЛ0У)К У + ((гЛ2У)К У + - + ((12-У)^2-У . (12)
а составляющую низких частот в виде следующей формулы.
Р(г) = Е«^2,К2ки)К2ки = ({2,Лои))л0и + ({г,Л2и))Л2и + - + ((1Лу2-1и)*м/2-1и) . (13)
Примеры фильтров на основе базисов Шеннона с вещественными и коммплексными числами. Если воспользоваться составляющими векторов и и У (преобразований Фурье) в [2] то при N=4 с использованием обратных преобразований Фурье можно по аналогии с (10) - (13) определить высоко - и низкочастотные составляющие для данных примеров и в этих базисах (табл.1). Таблица 1. Разложение вектора ^ по отцовскому и материнскому вейвлетам в различных базисах
Исходный вектор 2 Базисы
Хаара Шеннона 1 Шеннона 2
Р(2)
1 1 3 1 ---г 4 4 1
2 2
0 1 1 1 ---г 4 4 0
2
0 1 1 1 — + — г 4 4 1
2 2
1 1 3 1 ---г 44 1
2
от
1 1 1 1 — +—г 4 4 1
2 2
0 1 1 1 ---+— г 4 4 0
2
0 1 1 1 ----г 4 4 1
2 2
1 1 2 1 1 — + — г 4 4 0
Из сравнения базисов, приведенных в табл. 1, а также базисов Мейера, Баттла-Лемарье, Добеши [6, стр. 207] преимущество следует отдавать базису Хаара по следующим причинам: базис Хаара прост в вычислении; понятна на его основании и фильтрация сигнала; так низкочастотная часть на первом этапе равна нулевому коэффициенту преобразования Фурье, высокочастотная - соответствует отклонениям сигнала от его середины; в базисах, отличных от Хаара, задаются дополнительные требования к числу n. Например в базисах Шеннона оно должно быть кратным 4, а в базисе Хаара достаточно лишь его четности; числом N определяется и вид базиса Добеши [2].
Построение фильтров последовательностью сверток. В вейвлет-разложениях (12),( 13) замечено, что коэффициенты при Л2кУ, Л2ки являются соответствующими скалярными произведениями,
которые в виде сверток можно записать так [2,6,7]: {I, Р2У) = 2 *~(2к) , (7., Л2ки) = 7 *и(2к). Настоящие выражения являются теоретической основой быстрого вейвлет-преобразования. Тогда фильтр на основе такой свертки строится следующим образом [2], [6]:
1) образуются векторы и,У; 2) осуществляются свертки I*и,I*У ; 3) вводится оператор децимации - удаления составляющих с нечетными номерами р(г *и) р(г *у) ; 4) реализуется оператор разрежающей выборки - удваивания размера вектора вставкой нуля между двумя смежными значениями и (р(г *и),и (р(г *у)); 5) Осуществляется фильтрация - строится вектор низкочастотной составляющей в виде свертки. р(г)=и*и(р(г*и)) и высокочастотной 0(г) = V*и(р(г*у)); 6) восстанавливается сигнал I = р(г)+о(г).
На этом заканчивается первый этап разложения сигнала на высоко- и низкочастотные составляющие. Число всех этапов определяется по формуле Р = 1оё2 N . Каждьш последующий этап тостотт го фазы анализа и фазы синтеза. В фазе анализа на этапе п осуществляется: 1) ввод вектора гТ-1*ип-1 и его децимация 2п = р(1т-1 *и~п-1); 2) ввод векторов ип,Уп нормированного базиса размерности N = N2п-1; 3) осуществляются децимация и разреживание сверток 1Тп *ип, 1тп *Уп .
UDUn = uDzT„Un Judv* = UD(zTvn); 4) вычисляются высокочастотная и низкочастотная составляющие вектора zn: Zn Q(zn)=Vn *u{D{zn *vn)) p(zn)=un *u(d(z* *un)) . В фазе синтеза осуществляется: 1) разреживание векторов Q(zn), p(zn) (при этом разреживание выполняется числом n-1 раз до достижения размерности исходного сигнала). U(Q(zn)),U(p(zn)); 2) операциями свертки получают высокочастотную и низкочастотную составляющие сигнала на этапе n : Qn(z) = U(Q(zn))*U1 , Pn(z)=U(p(zn))*U.
p i
Окончательный результат анализа синтеза будет: z = P1(z0QP-i(z). Настоящий алгоритм реализован в специально составленной авторами программе Sub Вейвлет Анализ Синтез (). Сигнал Z можно представить в виде разложения:
z = c0b0 + C\h + ■■■ + Cv-1bv-1 (14)
по базису B , который представляет собой совокупность ортонормированных базисов [2] B = {b, JV-1
. Очевидно, что в этом случае будет c, =(zbT). Задача заключается в том, чтобы определить ортонор-мированные векторы в разложении (14). В теории вейвлетов [2] вместо (14) принимается следующая запись:
^(v/2P-1) ^(v/ 2P-1) ^(v/ 2P1 -l) ^(v/ 2p(p-l)-l)
z = Ak=0 C1'kV- P'k + Ak=0 C2W-p,k + ^k=0 C3,kW- p,k + ■■■ +^k=0 CNl2,kW-p,k (15)
Если принять что p = log2N, то (15) перепишется в виде:
z = c1,0.-p,0 + C2,0W-p, 0 + c2,0^-p, 1 + ■■■ + CN 2,0^-1,0 + CNj 2, 1^-1,1 + ■■■ + CNj 2, Nj 2-lV-\Nj2-1 . (16)
Здесь все векторы базиса В k записаны слева на право по степени детализации вектора z. В [2] они записаны наоборот справа налево. Построение базисных векторов у/-j к осуществляется в следующем порядке [2]: 1) применяется последовательность вейвлет-фильтров U1,V1; U2,V2; ... Up,Vp;
при этом Ui, V e i2(zn2), например исходя из (5), (6), при l=2, будет j^2j = Ц ° ;
2) Каждый из базисов jk строится по формуле [2] w~jk = Rjufj, где f = gl-1Ul-1(Vi), а gl = gl-1Ul-1(Ul) при начальных значениях, равных f1 = V1, g1 = U1. В приведенных выражениях оператор U1 -1() означает l-1 - кратное разреживание вектора Vl или Ul. Элементы сверток f1 (n), gl (n) для n = 0,1V„, N-1 вычисляются по формулам:
fl (n) = Y!'\ g-1 (m)Ul-1(Vi X* - m)), gi (n) = VN-1 {gi-1 (m)U1 -1(Ui )(n - m)). (17)
КМА - анализ для построения ортонормированного базиса вейвлет-преобразования. В фундаментальных работах [2,6,7] вейвлет - разложение предполагает наличие масштабирующей и уточняющей функций. В [2] масштабирующая функция еще называется отцовским вейвлетом. В [6, 7], термин «отцовский вейвлет» в определении масштабирующей функции не применяется. В тех же работах [2, 7] уточняющая функция (материнский вейвлет) просто названа вейвлетом. В системе Хаара в соответствии с [2] отцовскиу вейвлет записан в виде:
.Ы1'0 * Х < 1 . (18)
[0, иначе
материнский вейвлет строится по формуле [2]
y(x)=q(2x -1)-^2x) (19)
и, что легко показать, имеет вид
-1,0 < x < 12
^(x) = i 1,1/2 < x < 1 . (20)
0, иначе
Отцовский и материнский вейвлеты позволяют построить ортонормированный базис вейвлет-преобразования для представления дискретно заданной геоинформации с целью ее обработки. Для этого в вейвлет - теории разработан так называемый кратномасштабный анализ (КМА - анализ). На его основе создается удобный в использовании алгоритм построения ортонормированного базиса вейвлет-преобразования. КМА-анализ базируется на функциях:
j = 2-1 ^x - k); (21)
Wjk = 2-j 2w{l-jx - k). (22)
Эти формулы приведены в [2, 6, 7, 10]. При этом в приведенных работах показателям степени в (20), (21) приписывается как положительный знак, так и отрицательный, как в нашем случае. В случае отрицательной степени график функции растягивается по оси x, а в случае положительной - сжимается. В работе нас интересует растяжение по оси x, поэтому принята запись степени с отрицательным знаком. При такой записи осуществляется уточнение значений функций в зависимости от числа ортонормированных векторов базиса вейвлет-разложения, т. е. осуществляется увеличение деталей анализируемой информации или увеличение разрешения. поэтому в зарубежной литературе [10] КМА-анализ справедливо называется многоразрешающим анализом (multiresolution analysis). Построение вейвлет-базисов в системе Хаара будем вести на основе (18) - (22) в следующем порядке:
1) построение базисного вектора нулевого приближения Ф;
2) построение последующих уточняющих базисных векторов.
1. Построение базисного вектора нулевого приближения
В основу построения всех базисных векторов положим вейвлет - базис вида [3], применяемый нами в разложении (16). Хотя в теории вейвлетов [2-5] допускается возможность построения нескольких базисных векторов нулевого приближения, в данной работе мы ограничимся лишь одним в соответствии с данным разложением. При этом такое разложение наиболее часто применимо в различных рядах, в том числе и при разложении в ряд Фурье.
В соответствии с разложением (16), принятым здесь за основу, базисный вектор pk, единст-
,2 р
вен, так как к пробегает значения от 0 до II ,Р I 11 [2]. Поскольку N=2', то к=0. Следовательно в (21)
и
} = р, к = 0 . в соответствии с (18) 0 < 2 Рх - к < 1, или к < * <(* + кV2 Р . При к = 0, будет 0 < х < 1/2 Р составляющие этого вектора определятся по формуле:
о-р/2 11,0 < х < 2р
9},к =2 ' Г . (23)
[0, иначе
Таким образом на основе (23) получится вектор:
9р,0 = 2- р 2 (1,1.....1) (24)
с количеством одинаковых членов 2р.
2. Построение последующих уточняющих базисных векторов
Последующие уточняющие базисные векторы вычисляются по (22) с учетом (20). Порядок множества непересекающихся базисных векторов [2] и^-},к }кМ()2 )-1 определяется значением]. Это значение изменяется от р-1 до 0 с шагом -1. Тогда число базисных уточняющих векторов порядка ] составит величину N2} . А значение к будет изменяться от 0 до м/2} -1. Например для ]=р будем иметь следующий уточняющий вектор (единственный) 0 . Для ]=1 уточняющими векторами порядка I будут <//-10; ; ...; . Так, при N=8 р=3 уточняющими векторами порядка 1=1 при верхнем
пределе к равном м/2 -1 = 2р /2 -1 = 3 будут ^-10; ; у/-12; v-\3. Каждый из этих векторов при определенном к также определяется по (22) с учетом материнского вейвлета (20). Тогда можно записать, что
-1,0 < 2-'х - к < 1/2 1 [-1, к21 < х <(1/2 + к)2!
1,1/2 < 2-1 и к =
21 (25)
0,иначе 0, иначе
_1_
Vi k = 2 2
_ 1
В качестве примера возьмем p=3, l=1, k=2. Тогда Vi,2 = 2
-1,4 < х < 5 1,5 < х < 6 . 0, иначе
В табл. 2 дан пример вейвлет разложения как методом сверток, так и КМА-методом.
Можно отметить, что каждая составляющая вычислялась по формуле: для низкочастотной части
, р- } ^М/2р-7 „
р,} = 12р У Р,}-1 и высокочастотной 0^} = Рг-,}-1 - Р^} .
Исходный вектор Составляющие
c1,0 • Ф-3,0 c2,0 • У-3,0 c3,0 • У-3,0 c3,1 У-2,1 c4,0 У-1,0 c4,1 • У-1,1 c4,2 •у-1,2 c4,3 • У-1,3
4 7,25 -3,25 -1 0 1 0 0 0
2 7,25 -3,25 -1 0 -1 0 0 0
3 7,25 -3,25 1 0 0 -2 0 0
7 7,25 -3,25 1 0 0 2 0 0
10 7,25 3,25 0 -1,5 0 0 1 0
8 7,25 3,25 0 -1,5 0 0 -1 0
10 7,25 3,25 0 1,5 0 0 0 -2
14 7,25 3,25 0 1,5 0 0 0 2
Здесь j - номер этапа, i-номер составляющей сигнала в группе (группа состоит соответственно из 2, 4,...2p, элементов в зависимости от порядкового номера этапа), p-число этапов. Окончательный
результат будет Z = P1 + ^Р JQp-i где P1 - составляющая нулевой частоты (вектор у-30 в примере табл.1), Qp-i - уточняющие составляющие. В приведенных преобразованиях на каждом этапе при
нецеломp возможен остаток элементов, число которых меньше числа 21р при ip=1, 2,...,2p. По этому остатку находится также среднее значение TZ, которое записывается в старшей строке pip на место этого остатка элементов. В последней строке, следующей сразу после строки с номером целой части p, обозначаемой через pf, находится среднее по всем элементам предыдущей строки Ppf. Это среднее является одним и тем же для строки Ppf+i- Подстрока qjf+i вычисляется в общем порядке qpf+\ = Ppf - Ppf+1.
Исследование эффективности сжатия геоинформации и фильтрации шумов помощью вейв-летов. В основу исследований положены специально составленные авторами программы на языках VISUAL BASIC Excel (VBE) и IDL системы ENVI по сжатию и фильтрации геоинформации. В качестве первого объекта исследований принята модель рельефа, приведенная в [8]. В VBE-программе Sub МакросВЕЙВЛсИстОш() истинные высоты точек представлены массивом сс1(), а высоты, отягощенные случайными ошибками - массивом cm(). Также как и в [8] связь настоящих массивов определяется формулой cm(i) = cc\(i) + delta, где i изменяется от 0 до N - 1, а delta= Randbetwen(-t. +t)*Std-функция языка VBE генерирования случайного числа в интервале значений квантиля от -t до +t, Std -задаваемый стандарт случайных ошибок (шумов). Именно высоты массива cm() подвержены вейвлет-разложению по частотам. При вейвлет-сжатии и соответственно фильтрации оставались самые большие по амплитуде члены разложения на всех частотах. Фильтрационный вектор представлен массивом Filt= Array(1, 0.8, 0.6, 0.5, 0.25, 0.15, 0.1, 0). Например, при значении фильтра Filter=Filt(0)=0.8 на всех частотах разложения оставались значения, большие 0.8, остальные значения обнулялись. На выходе формировался отфильтрованный массив высот Tw(). Среднее квадратическое
ccl(i))2
i=l
. В табл. 3 приведе-
ошибка отфильтрованных высот определялось по формуле Sko = |
ны значения Sko в зависимости от значения фильтра, Filter и точности отмоделированных высот, величины Std. Оставшаяся информация в процентах вычислена как процентное отношени числа наибольших оставшихся коэффициентов к их общему числу до сжатия, равному 15.
Таблица 3. Средние квадратические отклонения Sko по каждому фильтру в метрах в зависимости от стандартов Std распределения случайных ошибок высот
Filter Оставшаяся информация в процентах Значения Std
0,05 0,1 0,3 0,4 0,5
1 25 0,66 0,88 0,71 0,93 0,94
0,8 33 0,55 0,35 0,63 0,86 0,88
0,6 40 0,23 0,26 0,63 0,86 0,94
0,5 48 0,23 0,26 0,64 0,86 0,94
0,25 51 0,19 0,24 0,63 0,86 0,99
0,15 66 0,15 0,20 0,64 0,86 0,99
0,1 70 0,11 0,17 0,63 0,86 0,99
0 100 0,09 0,18 0,63 0,86 0,99
Для сохранения максимально-высокой точности исходных высот при их вейвлет-разложении сжатие информации недопустимо; даже при отсутствии сжатия (Filter=0) и наличии случайных ошибок среднее квадратическое отклонение результирующих высот больше их стандартного отклонения на
входе; сжатие информации можно допускать, но при этом следует учитывать порог понижения точности изображения рельефа; например, при трехкратном сжатии (Filter=0,8) и Std=0,1 среднее квад-ратическое отклонение высот на выходе Sko = 0,35 больше Std в 3,5 раза, в остальных случаях оно больше в 2 раза.
Одновременно авторами составлена программа pro oroi_data_corr24bitW^VE на алгоритмическом языке IDL системы ENVI вейвлет-фильтрации (сжатия) изображений.
При большем сжатии различных несовпадений больше и они представляют собой уже шумы, при которых изображение не пригодно к использованию. Таким образом, при сжатии изображения рекомендуется также подобрать вначале фильтр в соответствии с пунктами 1 -3 в случае рельефа, или при котором отсутствуют видимые несовпадения. И лишь при нем осуществлять вейвлет-сжатие.
Заключение
Методика расчета значения фильтра для сжатия информации о высотах рельефа может быть принята такой.
1. Определяется по известной методике оценка стандарта высот рельефа Std в виде средней квад-ратической ошибки съемки рельефа.
2. Для данного объекта определяются средние квадратические отклонения (СКО) по каждому фильтру.
3. Если значение средней ошибки рельефа, полученное как 0,8 Std, меньше трети высоты сечения рельефа, то по полученной экспериментально можно выбрать порог (величину СКО) понижения точности сжатой информации о рельефе при определенном значении фильтра (Filter). Если средняя ошибка рельефа больше трети высоты сечения рельефа и приближается к величине СКО, получаемой при 25-30 процентном сжатии информации (то есть трех и четырехкратном сжатии), то в таком случае значение порога (величины СКО) понижения точности сжатой информации о рельефе не имеет значения.
ЛИТЕРАТУРА
1. Шпаара, Д. Точное сельское хозяйство (Precision agriculture) / Д. Шпаара, А. Захаренко, В. Якушева ; под ред. Д. Шпаара. - СПб. - Пушкин, 2009. - 398 с.
2. Фрейзер, М. Введение в вейвлеты в свете линейной алгебры / М. Фрейзер ; пер. с англ. - М. :БИНОМ, Лаборатория знаний. 2008. - 487 с.
3. Пискунов, Н. С. Дифференциальные и интегральные исчисления для втузов т.2 / Н. С. Пискунов. - М. : Наука. -1978. - 575 с.
4. Мазурова, Е. М. Алгоритмы быстрого преобразования Фурье / Е. М. Мазурова // Известия вузов. Геодезия и аэрофотосъемка. - 2004. - №3. - С. 18-35.
5. Мазурова, Е. М. Двумерное и матричное представление быстрого преобразования Фурье / Е. М. Мазурова // Известия вузов. Геодезия и аэрофотосъемка. - 2004. - №4. - С. 3-12.
6. Малла, С. Вейвлеты в обработке сигналов / С. Малла; пер. с англ. - М. : Мир, 2005. - 671 с.
7. Чуи, К. Введение в вейвлеты / К. Чуи ; пер. с англ. - М. :Мир, 2001. - 412 с.
8. Ярмоленко, А. С. Фильтрация геоинформации в рядах Фурье / А. С. Ярмоленко, О. В. Скобенко // Известия вузов. Геодезия и аэрофотосъемка. - 2016. - №1. - С. 107-113.
9. Шовенгердт, Р. А. Дистанционное зондирование, модели и методы обработки изображений / Р. А. Шовенгерд. -М.: Техносфера, 2010. - 560 с.
10. Журкин, И. Г. Автоматизированная обработка данных дистанционного зондирования: учеб. для вузов / И. Г. Журкин, Н. К. Шавенько ; под общ. ред. И.Г. Журкина. - М.: ООО «Диона», 2013. - 456 с.
11. Дьяконов, В. П. Вейвлеты. От теории к практике / В. П. Дьяконов. - М.: Солон - Р, 2002. - 448 с.
12. Уэлстид, С. Фракталы и вейвлеты для сжатия изображений в действии / С. Уэлстид. - М. : Издательство Триумф, 2003. - 320 с.
13. Воробъев, В. И. Теория и практика вейвлет преобразования / В. И. Воробъев, В. Г. Грибунин. - СПб. : ВУС, 1999. - 204 с.
14. Яковлев, А. Н. Введение в вейвлет-преобразование : учеб. пособие / А. Н. Яковлев. - Новосибирск: Издательство НГТУ, 2003. - 104с.
15. Красильников, Н. П. Цифровая обработка 2D- и 3D-изображений : учеб. пособие / Н. П. Красильников. - СПб БХВ - Петербург, 2011. - 608 с.
16. Малинников, В. А. Анализ методов формирования мультифрактальной меры, основанных на вейвлет-обработке экспериментальных данных / В. А. Малинников, Д. В. Учаев // Известия вузов. Геодезия и аэрофотосъемка. - 2007. - №6. -С. 57-61.
17. Большаков, В. Д. Теория математической обработки геодезических измерений / В. Д. Большаков, П. А. Гайдаев -Изд. 2, перераб. и доп. - М. : «Недра», 1977. - 367 с.
