=Paper=
{{Paper
|id=Vol-2033/30_paper
|storemode=property
|title=Неконтролируемая классификация мультиспектральных изображений
(Unsupervised Classification of Multispectral Images)
|pdfUrl=https://ceur-ws.org/Vol-2033/30_paper.pdf
|volume=Vol-2033
|authors=Valeria S. Sidorova
}}
==Неконтролируемая классификация мультиспектральных изображений
(Unsupervised Classification of Multispectral Images)==
<pdf width="1500px">https://ceur-ws.org/Vol-2033/30_paper.pdf</pdf>
<pre>
          UNSUPERVISED CLASSIFICATION OF MULTISPECTRAL IMAGES

                                           Valeria S. Sidorova

             Institute of Computational Mathematics and Mathematical Geophysics SB RAS,
                                          Novosibirsk, Russia

                                                 Abstract

       In the proposed hierarchical histogram algorithm, the clustering detail is found, which is different
in the subdomains of the vector space of spectral features depending on the given separation of the
clusters. The question of reducing the dimension of the space of data attributes is also being considered.
The application of the algorithm for uncontrolled classification of terrestrial cover using various satellite
remote sensing data is illustrated.

      Keywords: remote sensing, uncontrolled classification, multidimensional histogram, cluster
separation, dimensionality of space, hierarchical algorithm
     НЕКОНТРОЛИРУЕМАЯ КЛАССИФИКАЦИЯ МУЛЬТИСПЕКТРАЛЬНЫХ
                        ИЗОБРАЖЕНИЙ

                                  Сидорова В.С.
Институт вычислительной математики и математической геофизики СО РАН, Новосибирск

     В предложенном иерархическом гистограммном алгоритме отыскивается детальность кластери-
 зации, различная в подобластях векторного пространства спектральных признаков в зависимости от
 заданной отделимости кластеров. Также рассматривается вопрос о сокращении размерности про-
 странства признаков данных. Иллюстрируется применение алгоритма для неконтролируемой класси-
 фикации земного покрова по различным спутниковым данным дистанционного зондирования.
     Ключевые слова: дистанционное зондирование, неконтролируемая классификация, многомерная
 гистограмма, разделимость кластеров, размерность пространства, иерархический алгоритм.

      Введение. Предлагается неконтролируемая иерархическая классификация мультиспек-
тральных спутниковых данных с заданием порога отделимости кластеров и сокращением раз-
мерности данных.
      Внутри кластеров осуществляется кластеризация по унимодальным кластерам мето-
дом [1]. Быстрый непараметрический не итеративный алгоритм [1] разделяет векторное про-
странство признаков по унимодальным кластерам, модальные векторы которых соответ-
ствуют локальным гистограммным максимумам, а границы проходят по долинам гисто-
граммы. Он решает три задачи кластеризации одновременно, используя метод графов: нахо-
дит локальные максимумы гистограммы, объявляя их корнями деревьев-кластеров, проводит
границы между кластерами по долинам гистограммы и относит все вектора признаков к своим
кластерам. Алгоритм является жестким: каждый вектор принадлежит только одному кла-
стеру. Детальность кластеризации в алгоритме [1] регулируется заданием числа уровней пред-
варительного квантования, одинакового для всего векторного пространства. Задается отсече-
нием младших битов в байтах различных векторных направлений. Каждое такое отсечение
соответствует уменьшению векторного пространства вдвое. Алгоритм [1] был реализован ав-
тором для ЭВМ "БЭСМ-6", а затем для РС[2,3]. Для получения критерия качества классифи-
кации предложено ввести разделимость и отделимость кластеров. В [4] предложен иерархиче-
ский алгоритм, определяющий детальность в зависимости от кластерной разделимости подоб-
ластей. Но иногда бывает важно выявить предельную детальность для данной разделимости.
В другом иерархическом алгоритме [5] ставится другая цель: для подобластей пространства
векторов найти такие свои максимальные детальности, при которых еще порождаются дочер-
ние кластеры с разделимостью ниже заданной. Это позволяет с одной стороны исследовать
структуру многоспектральных данных достаточно подробно, причем на разных иерархиче-
ских уровнях, с другой стороны, регулировать детальность исследования.
      Кроме того, можно по-разному задавать закон изменения детальности. Если различные
спектральные направления не эквивалентны, то и менять их можно по-разному. И это может
привести к сокращению размерности векторного пространства признаков.
      Этапы иерархического алгоритма. Параметром детальности кластеризации является
число уровней квантования векторного пространства, обозначим его n. На каждом этапе пред-
ложенного иерархического алгоритма находим такое число n (и соответствующие новые век-
торы g), при котором распределение по унимодальным кластерам, полученное методом [1],
дает абсолютный минимум мере (2), в диапазоне изменения 255>n> n1.
                                                                               j
     В [6] были предложены: мера отделимости унимодального кластера m (n) (1), и мера
средней разделимости K (n) кластеров m(n) (2):
                                                     B (n)j
                                              1
                                  m (n)  j
                                    j
                                                  j   
                                         B (n)* H (n) i 1
                                                           hi j (n),                        (1)


                                                  156
                                               1 K(n) j
                                     m(n)          m (n) ,
                                              K(n) j 1
                                                                                         (2)

                                                                    j
где hij (n) значение гистограммы в i-той точке границы кластера j, B (n) число точек границы
              j
кластера, H (n) максимальное значение гистограммы.
      В [6] показано, что (2) удовлетворяет требованиям, предъявляемым к мере разделимости
кластеров или качества кластеризации [7,8]. Минимумы (2) соответствуют лучшим классифи-
кациям. Всегда m j (n)  1 и m (n)  1 . Традиционные меры разделимости оперируют сред-
неквадратичным отклонением и расстоянием, которые взаимозависимы при жесткой класте-
ризации. Введенные меры позволят сравнивать изолированность распределений с тесно рас-
положенными кластерами. Более подробно о рассмотренных мерах в [6].
      На новом этапе иерархии алгоритм для каждого кластера увеличивает число уровней
квантования (детальность), полученное на предыдущем этапе, и в новом интервале находит
свое новое число и соответствующее наилучшее кластерное распределение в смысле меры(2).
Новый алгоритм рассматривает каждый дочерний кластер как отдельную область только то-
гда, когда его разделимость удовлетворяет условию (3):
                                              j
                                           m (n)  ε ,                                  (3)
где ε заданная точность отделимости кластера.
     Плохо разделенные дочерние подкластеры (не удовлетворяющие (3)) рассматриваются
как одна область. Детальность квантования увеличивается, и деление продолжается, пока
n<n1 в каждом дочернем кластере. Затем результаты анализируются снизу-вверх для плохих
кластеров, когда условие (3) нарушено. Процесс заканчивается, когда в иерархическом дереве
попадется хороший в смысле (3) предок или хорошая сестра. В первом случае данные плохого
кластера возвращаются на уровень хорошего предка, во втором кластер считается плохой вет-
вью и может быть отнесен к общему фону кластеров, не удовлетворяющих (3). Рассмотрим
изображение земной поверхности (рис.1); получен со спутника NOAA_17 17 апреля 2003г.
Изображение пятиспектральное. Его объем около 1,7 мегабайт, размер 1480*1124 пикселей,
разрешение около км* км / пиксель.
     Верхнюю левую часть изображения занимают тающие снега тайги Сибири, внизу отта-
явшая поверхность Казахстана. Правую часть снимка покрывают сплошные и полупрозрачные
облака. Пять спектральных каналов позволяют различать объекты, не различимые лишь в ви-
димом диапазоне. Например, снег и облака визуально не различимы по цвету и яркости, но в
диапазоне инфракрасного 3,55 – 3,93 мкм снег сильно поглощает излучение от Земли и выгля-
дит даже черным, а облака белые. В диапазонах 10,3 – 11,3 мкм и 11,5 - 12,5 мкм (рис.1 d и
рис.1 e) могут быть различены типы облаков (кучевые, перистые и другие, имеющие разную
высоту). Для построения многомерной гистограммы был применен алгоритм с использова-
нием чередования таблиц хеширования и сортировки [9].
     Задано ε = 0.07. Было пройдено семь этапов иерархического алгоритма, дальнейшее
деление кластеров приводило к нарушению условия разделимости (3) для большой части под-
кластеров.
     Пространственное расположение кластеров на карте хорошо соответствуют земным объ-
ектам. Интересно заметить, что три кластера, полученные автоматически новым алгоритмом
соответствуют районам добычи полезных ископаемых открытым способом. Они отмечены
черным и желтым тонами на карте и локализуются в районах Экибастуза, Баянаула - угледо-
быча, в Майкаине золото открытым способом; южнее Семипалатинска также уголь и золото.
Зоны отечественных разработок, в частности, Кузбасс находятся в заснеженной области под
полупрозрачными облаками и не выделились в отдельный кластер. Для всего изображения по-
лучено 120 кластеров с отделимостью меньшей ε  0.07 . Для сравнения: при той же макси-


                                                  157
мальной детальности n = 50, одинаковой для всего пятиспектрального изображения алгорит-
мом [1] получено 2686 кластеров. Время работы нового алгоритма – несколько минут на од-
ноядерном компьютере.




                        Рис.1. Кластерная карта седьмого этапа иерархии.

      Другим важным аспектом является уменьшение размерности данных. Квантование про-
странства признаков может производиться по разным правилам. До сих пор в каждом спектральном
направлении число уровней квантования сохранялось одинаковым. Однако, в общем случае, данные
вытянуты вдоль какого-то направления, и правило квантования, обеспечивающее наименьшую потерю
информации, требует различного подхода в различных направлениях, а именно: квантование должно
сохранять ячейку квантования в форме гиперкуба (а не гиперпараллелепипеда). Это условие будет вы-
полнено, если число уровней квантования вдоль каждой оси собственного пространства пропорцио-
нально квадратному корню из соответствующего собственного числа. (Собственное число характери-
зует разброс вдоль оси), а именно:
                                    N e1 N e2       N
                                                 ek ,                                     (1)
                                    S e1 S e 2      S ek
где N e1 , N e 2 ,  , N ek числа уровней квантования вдоль для соответствующих собственных
                                               2         2     2
векторов по k ортонормированным осям, а S e1 , S e 2 ,  , S ek собственные числа.
     Для перевода пространства в пространство собственных векторов применялся метод
Якоби [10]. Зададим максимальное число уровней квантования в собственном пространстве
равным N em =255, таково обычное число уровней серого для данных дистанционного зонди-
рования по каждому измерению. Тогда, в соответствии с пропорциями (1) может быть найдено
число уровней квантования и по другим осям собственного пространства. Для задач кластери-
зации это число должно быть больше или равно 2, иначе эта компонента одинакова для всех
векторов и никакой роли в кластеризации не играет. Таким образом, если отношение             S em /
S ex < 2, то соответствующая ось x может не рассматриваться, и мы получаем сокращение раз-
мерности пространства признаков.

                                                   158
     Алгоритм был апробирован на семиканальном (в видимом и инфракрасном диапазонах)
изображении Омской области ИСЗ «Landsat-8» 08.02.2014 (размер 3 161×2 590 пикселей, раз-
решение 15 м), предоставленном сибирским центром ФГБУ «НИЦ «Планета»». Алгоритм кла-
стеризации предварительно осуществляет сокращение размерности векторного пространства
спектральных признаков с семи до трех. Детальность, различная по полученным кластерам,
определяется делимым иерархическим гистограммным алгоритмом для предельной отдели-
мости кластеров d = 0,15 (0 < d < 1). Иерархичность получающихся кластерных карт отражает
иерархичность реальных объектов на космоснимках. Выбор числа этапов осуществляется сов-
местно с экспертом. Например, на первом этапе иерархии получено всего 6 кластеров, два из
которых (красный и черный) соответствуют дымам ТЭЦ Омской области. Для десяти этапов
иерархии получено 27 унимодальных кластеров (Рис. 2); специалисты НИЦ «Планета» уста-
новили, что расположение ярко розового и фиолетового кластеров соответствует загрязнению
территории Омской области.
     Можно также рассматривать не все пространство векторов при сокращении размерности,
а лишь часть его, относящуюся к каждому кластеру. Поскольку алгоритм иерархический, и
каждый кластер делится на подкластеры на каждом этапе иерархии в соответствии с измене-
нием детальности представления данных, то размерность каждого подкластера может изме-
нится. При решении задачи внутри кластера (построении ковариационной матрицы) исполь-
зуется уже построенная ранее гистограмма признаков в виде определенным образом органи-
зованного списка. Рассмотрим пример




Рис.2. Кластерная карта, полученная делимым иерархическим гистограммным алгоритмом. 15 этапов
  иерархии. d= 0, 015. Получено 54 кластера (включая маленькие). Загрязнение: лиловые оттенки.

                                                159
Рис. 3.a) Кластеризация по пяти спектральным каналам иерархическим гистограммным алгоритмом;
   5 этапов иерархии; задано d=0,1. Получено 22 кластера. б) Кластеризация иерархическим гисто-
граммным алгоритмом с поиском размерности по кластерам; 15 этапов иерархии; задано d=0,1; полу-
                                        чено 40 кластеров.

      Анализируется изображение поверхности Земли со спутника NOAA 17 от 7.04.2003,
полный кадр (1328x624) пикселей представлен в пяти спектральных каналах (один в видимой
части спектра, остальные в инфракрасной), 4 мегабайт. В нижней части снимка формирование
вихря, озера; в верхней в основном – тающие снега, тайга Сибири. У левой границы-Омск, у
правой-Абакан. Хорошо прослеживается железная дорога: Омск - Чаны - Барабинск - Новоси-
бирск. На рис.3a кластерная карта при глобальном сокращении размерности. На рис. 3b при
покластерном. При большей детальности кластеризации (Рис. 3b.) кластер облаков делятся на
унимодальные подкластеры с заданной отделимостью 0,1, и некоторые из них требуют пятис-
пектрального рассмотрения. Это полупрозрачные облака. Время вычислений оказалось в три
раза меньше, чем для полностью пятиспектрального варианта и составило несколько минут на
одноядерном компьютере.
      Работа выполнена при финансовой поддержке РФФИ (проект № 16-07-00066) и Про-
граммы I.33П фундаментальных исследований Президиума РАН (проект № 0315-2015-0012)).

                                            ЛИТЕРАТУРА

[1]  Narendra P.M. and Goldberg M. A non-parametric clustering scheme for LANDSAT // Pattern Recog-
     nition. 1977. Т 9. P. 207 -215.
[2]   Cидорова В.С. Кластеризация многоспектральных изображений с помощью анализа многомер-
     ной гистограммы // Новосибирск. Сб.: Математические и технические проблемы обработки изоб-
     ражений. СО АН СССР. 1986. С. 52-57.
[3]   Сидорова В.С. Классификация многоспектральных космических изображений поверхности
     Земли с помощью разделения многомерной гистограммы по унимодальным кластерам // Ж. Вест-
     ник КазНУ., сер. географическая. 2004. N 2(19). С. 206-210.
[4] V.S. Sidorova. Automatic Hierarchical Clustering Algorithm for Remote Sensing Data. // Pattern
     Recognition and Image Analysis. 2011. Vol.21. No.2. P. 328 - 331.
[5] V.S. Sidorova. Detecting Clusters of Specified Separability for Multispectral Data on Various      Hi-
     erarchical Levels // Pattern Recognition and Image Analysis. 2014. Vol. 24. No. 1. P. 151-155.
[6]   Сидорова В.С. Оценка качества классификации многоспектральных изображений гистограм -
     мным методом // Автометрия. 2007. Том 43. №1. С. 37- 43.
[7]   M. Halkidi, Y. Batistakis and M. Vazirgiannis. On clustering validation techniques // Journal of In-
     telligent Information Systems. 2001. No.17 (2-3), P.107-132.
[8] Keinosuke Fukunaga. Introduction to Statistical Pattern Recognition. Academic Press, New York and
     London, 1972.
[9] Сидорова В.С. Многомерная гистограмма и разделение векторного пространства признаков по
     унимодальным кластерам. Труды конференции GraphiCon’2005. 2005. C. 267-274.
[10] Калиткин Н.Н. Численные методы. Москва. “Наука”. 1978. С. 512




                                                      160

</pre>