Обработка изображений сейчас модное направление в области разработки программного обеспечения. Начать обработку изображения можно уже на стадии фотографирования. По поводу качества такой обработки сказать ничего не могу, но такие возможности фотоаппаратов, как обнаружение лица и других объектов вполне себе работают и причем неплохо. Другое дело – это распознавание лиц, тут пока технологии отстают. Мой ноутбук никак не хочет идентифицировать меня по лицу!

Основой обработки изображения является выделение области изображения со схожими характеристиками, такое, как, например, выделение образов. Список алгоритмов для выделения образов или сегментов уже давно известен и обсуждается многими. Хотя, что значит обсуждается? Просто добавляется каждый раз новая статистика по работе алгоритмов на практике – как работает, с какой точностью.

Когда я слышу о задаче выделения областей, у меня сразу возникает несколько вопросов:

— По какому критерию производить сравнение?
— Как определить, что пиксели относятся к одной области?

Начнем с первого пункта. Тут я постараюсь описать все довольно просто.

Все изображение состоит из маленьких частиц-пикселей. На сегодняшний день пиксель может быть представлен двумя способами – это с помощью RGB или YCrCb. Это так называемое представление цветного пространства пикселя. Между двумя этими представлениями существуют формулы перехода от одного пространства к другому. На мой взгляд, цветное пространство YCbCr используется только в формате кодирования JPEG, так как позволяет значительно уменьшить хранимую информацию о пикселях. При работе с растровыми изображениями для анализа чаще всего рассматривают формат представления RGB, хотя в последнее время к нему добавили еще одну компоненту – это так называемая “яркость“. Таким образом, для рассмотрения алгоритмов в этой статье я выбираю представление пикселей в RGB – пространстве, при необходимости может добавится еще одна, четвертая компонента.

Вопрос отнесения пикселей к одной области изображения представляет собой интересную задачу. Алгоритм, который определяет отношение пикселя к заданной области, называется «Выращивание регионов, дробление-слияние». Смысл метода в том, чтобы постепенно добавлять к определённой области пиксели, удовлетворяющие характеристикам области.

Для присоединения пикселя к области существует несколько алгоритмов. Это и алгоритм “k-средних“, и “watershed” алгоритм. Смысл алгоритмов определить критерий, границу, по которой будет происходить отбор пикселей.

Для выбора этой границы произведем анализ изображения. Будем оценивать изображение по распределению пикселей и их значению. Наиболее популярный инструмент для разработки и анализа алгоритмов, на мой взгляд, это Matlab. В последних версиях Matlab появился ToolBox для обработки и анализа изображения. Тут реализованы почти все функции, которые необходимы для этого.

Одна из самых популярных функций в Matlab для анализа – это построение гистограмм. Тут надо заметить, что стандартная функция для считывания изображения в память – это imread(). Это функция на выходе формирует трехмерный массив пикселей, где число строк и столбцов равно ширине и высоте изображения соответственно. Итак, после применения этой функции получаем массивы пикселей в пространстве RGB.

Для построения гистограммы требуется одномерный массив, поэтому построим три гистограммы.

Наше изображение

Для построения гистограмм изображений предусмотрен специальный инструмент imhist, который строит распределения пикселей по цветам (256 цветов).

Гистограммы по трем матрицам.

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

Эти результаты были получены, если принять гипотезу о том, что одинаковые пиксели соответствуют одной области.

Здесь мы постарались выделить области пикселей, которых наибольшее количество в канале, далее можно использовать рекуррентную функцию для выделения подобных областей до тех пор, пока изображение не будет разбито на области по значению пикселей.

Для получения данного результата был применен следующий цикл.

for i=1:1:x
  for j=1:1:y
    if(abs( fgrey(i,j)-position ) < 2)
      fgrey(i,j)=t;
      count=count+1;
    end
  end
end

Здесь position – это значение пикселей, по которым происходит выборка областей. В этом выражении принимается поправка, что отклонение от выбранного значения может составлять 2 пункта.

Метод простой и легко реализуемый, хотя, конечно, он далек от совершенства. Вот тут я могу предложить небольшую модификацию метода, которая будет рассмотрена ниже.

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

Рассмотрим тут его упрощенный вариант.

Будем рассматривать массив пикселей из трех подряд идущих значений. В нашем случае возьмем три пикселя. Найдем их среднее значение.

среднее=(f(i-1,j,z)+f(i,j,z)+f(i+1,j,z))/3;

Это значение примем за основу. Далее будем наращивать области, также используя три подряд идущих пикселя. Если разница значений удовлетворяет заданному критерию, то мы присоединяем эту область к исходной области, если нет, отбрасываем. Для трех матриц получаем следующий результат.

То есть мы видим, что изображение распадается на области с одинаковыми цветовыми значениями.

Функция для нахождения среднего по трем пикселям.

for z=1:1:3
  for i=2:3:x-1
    for j=1:1:y 
      dev=(f(i-1,j,z)+f(i,j,z)+f(i+1,j,z))/3;
      if(dev~=0)
	break;
      end
    end
  end
end

Функция для выделения областей.

for z=1:1:3
  for i=2:3:x-1
    for j=1:1:y
      devnew=(f(i-1,j,z)+f(i,j,z)+f(i+1,j,z))/3;
      if((dev-devnew)<0.1)
	dev=devnew;
	f(i-1,j,z)=0;
	f(i,j,z)=0;
	f(i+1,j,z)=0;
      end
    end
  end
end

Погрешность для трех пикселей принята равным 0.1, данный параметр можно изменить, но для этого изображения параметр 0.1 является наилучшим.

Таким образом, после применения всех фильтров мы получаем следующее изображение.

Визуальная оценка показывает, что выделение области вторым методом, дает более плавные переходы из одной области в другую.

Дальнейшее улучшение результатов даст возможность выделить области. Для этого надо использовать метод выделения пограничных контуров.

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

Например, выбирается точка из заданной области point(i,j), далее ищется точка со схожими значения в окружающем пространстве.

for j=point_j:1:y
  if(F(point_i,j)~=0)
    [point_i,j]=Point(F,point_i-1,j);
  else
    for z=1:1:3
      f(i-2,j,z)=0;
      f(i-1,j,z)=0;
      f(i,j,z)=0;
    end
  end
end

Здесь рассмотрение стоящих рядом точек идет только в одном направлении, сдвигаясь вправо по столбам.

Применение этого метода к уже частично обработанному изображению приближает нас к цели. Конечно, тут не получен идеальный вариант, но был рассмотрен математический аппарат для выделения областей со схожими характеристиками пикселей.

Следующим улучшением может быть создание более интеллектуальной системы, которая на основании некоторых изображений создает базу коэффициентов для отбора пикселей. Такие системы называются обучающимися.

Теперь рассмотрим реализацию этого примера на языках программирования высокого уровня.

В .NET существует класс Bitmap. Данный класс предоставляет большие возможности для обработки изображений.

Данный класс, также как и Matlab, представляет изображение в виде массива пикселей. Есть лишь одно отличие – вместо размерности матрицы три, в классе Bitmap представляется изображение в виде матрицы размерностью четыре.

Применим ту же фильтрацию, что и для Matlab, сделав только одно усовершенствование. Построим для каждого цвета что-то вроде гистограммы, но только будем объединять рядом стоящие значения пикселей. В нашем случае расстояние между схожими пикселями будет равно 4.

public void FindMean(BitmapIm, intlength, charcanal)
{
  int middlepixel = 0;
  middlepixel=Im.GetPixel(460,256).G;
  int count = 0;
  int countofcolor=0;
  countofcolor=Math.Abs(256/length);
  mas_color=new int[countofcolor];
 
  for (int i = 0; i<countofcolor; i++)
    mas_color[i] = 0;
 
  if (canal=='B')
  {
    for (int i = 0; i<Im.Width; i++)
      for (int j = 0; j <Im.Height; j++)
	if (Im.GetPixel(i,j).B!=Color.Black.B)
        {
	  mas_color[Math.Abs(Im.GetPixel(i, j).B / length)]++;
	  count++;
	}
  }
 
  if (canal == 'G')
  {
    for (int i = 0; i<Im.Width; i++)
      for (int j = 0; j <Im.Height; j++)
	if (Im.GetPixel(i, j).G != Color.Black.G)
	{
	  mas_color[Math.Abs(Im.GetPixel(i, j).G / length)]++;
	  count++;
	}
  }
 
  if (canal == 'R')
  {
    for (int i = 0; i<Im.Width; i++)
      for (int j = 0; j <Im.Height; j++)
	if (Im.GetPixel(i, j).R != Color.Black.R)
	{
	  mas_color[Math.Abs(Im.GetPixel(i, j).R / length)]++;
	  count++;
	}
  }
 
  if (canal == 'A')
  {
    for (int i = 0; i<Im.Width; i++)
      for (int j = 0; j <Im.Height; j++)
	if (Im.GetPixel(i, j).A != Color.Black.A)
	{
	  mas_color[Math.Abs(Im.GetPixel(i, j).A / length)]++;
	  count++;
	}
  }
}

Здесь происходит распределение значений пикселей по цветовым компонентам. А именно, пиксели с близко стоящими значениями объединяются в группу. Расстояние между такими группами в нашем примере равно 4. В дальнейшем, благодаря этому массиву, можно выделять значения пикселей, которых в изображении большинство.

При помощи разделения каналов было установлено, что наилучший результат получается при выделение областей из голубого канала, наихудший результат получен при выделении из зеленого.

Как видно, выделилась область голубого неба.

С помощью наложения некоторых фильтров, сначала голубого, а потом красного, получается следующее изображение.

Пример кода, для обработки голубого канала.

public void MarkBlue(BitmapIm)
{
  int metka = -1, count = 0, count_opposite = 0;
  // Color pixel = MarkIm.GetPixel(0, 0);
  for (int i = 0; i<Im.Width; i++)
    for (int j = 0; j <Im.Height; j++)
      if (Im.GetPixel(i, j) == Color.Black)
	//    count++;
 
  count = 0;
  FindMean(Im, 4,'B');
  int maxColor = mas_color[0];
  int maxIndex = 0;
  for (int i = 1; i< 256 / 4; i++)
    if (maxColor<mas_color[i])
    {
      maxColor = mas_color[i];
      maxIndex = i;
    }
 
  //  MessageBox.Show(Convert.ToString(mas_color[20]));
  for (int i = 0; i<Im.Width; i++)
    for (int j = 0; j <Im.Height; j++)
      if (maxIndex> (Im.GetPixel(i, j).B) / 4)
      {
	//  count++;
	Im.SetPixel(i, j, Color.Black);
      }
 
  for (int i = 0; i<Im.Width; i++)
    for (int j = 0; j <Im.Height; j++)
      if (Im.GetPixel(i, j) == Color.Black)
	count++;
}

Фактически, мы получили простой математический аппарат для выделения областей. Тут за основу принимается распределение пикселей по значению, геометрические и пространственные характеристики никак не учитываются.

В дальнейшем после применения фильтров для цвета пикселей, на мой взгляд, должен быть поставлен вопрос о правильном выделении пограничных контуров. Тут может возникнуть проблема, которая отображена на изображениях, приведенных выше. Что, если область принадлежит области “лошади”, но при этом цветовые характеристики совпадают с областью “неба”?

Думаю, для решения этой проблемы потребуется визуальная оценка результатов, сделанная человеком.

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

Я постараюсь максимально объективно сравнить два языка. В конечном итоге, выбор платформы, стиля и языка – не определяющий фактор успеха. Главное, чтобы руки разработчика не повторяли форму интеграла. Ведь язык лишь инструмент, а инструмент, как известно, продолжение рук.

Начнем.

WEB развивается уже второй десяток лет. Время статичных HTML-страниц кануло в истории, на смену которому пришел абстрактный стандарт WEB2, главное требование которого – взаимодействие с пользователем. Спорить, правильно ли это, конечно, можно, только в таком случае появляется вполне резонный вопрос: зачем читать эту статью, если она о средствах создания подобных сервисов?

Порог вхождения.

Традиционно для сравнительного анализа языков в первую очередь используют удобство синтаксиса, скорость работы и возможности.

Такой подход логичен для выбора собственного. А вот с точки зрения работодателя, качества самого программиста значительно важнее.

PHP считается простым языком. Это не удивительно: от открытия учебника до первого “HelloWorld-a” проходит не больше 5 минут. С одной стороны это хорошо. Но подводных камней не мало.

Потратив неделю на изучение PHP, можно научиться базовым приемам ввода-вывода, примитивной выборке из БД и парочке других алгоритмов обработки информации. С усложнением задач появляется соблазн изобретать “велосипеды” и “костыли” на основе старых знаний, делая код нечитабельным, уничтожая возможность дальнейшей поддержки на корню.

Работа с таким “кодером” в команде, или еще хуже, поддержка проекта, написанного им, превращается в будни школьного учителя. Естественно, это утверждение относится не ко всем, но такая проблема существует.

Что касается Java, порог вхождения очень высокий. Вместо “Установка LAMP под Ubuntu” и “Скачать самоучитель PHP“, необходимо вникнуть в сервлеты, взаимодействие между ними, jsp, научиться работе с каким-нибудь фреймворком.

Плюс ООП, который в Java обязателен. Если проводить аналогию, начало работы сравнимо с высшей математикой в ВУЗ-е: всех неспособных отправить на службу родине.

Подобный подход к оценке языка не панацея. Достаточно много серьезных проектов, написанных на PHP. Более того, на данный момент в сфере WEB-разработок он доминирует.

Синтаксис и код.

Одна из обсуждаемых проблем Java – длиннейшие названия классов и стандартных методов. Чего только стоит “System.out.println()” оного вместо “echo” на PHP. Конечно, это легко решается “import-ом” в начале документа или выбором современной IDE, в которой достаточно написать первую букву метода или поля, как появится возможность выбора нужного из списка. Тем не менее, отрицать впечатляющий объем кода на Java – предвзято.

PHP тоже не идеален. Взять к примеру обязательный знак доллара перед переменной. Конечно, рациональное зерно различия между обозначениями переменных и методами(функциями) присутствует. Лично у меня это вызывает жуткую головную боль. В первую очередь от созерцания самого знака. Плюс, я так и не сумел привить у себя
привычки ставить доллар перед каждой переменной. Слишком меркантильно, что-ли…

Что касается кода, я думаю мало кто станет спорить, что серьезный проект == серьезная работа в команде. А работа в команде требует годного, понятного кода, качественную документацию, и возможность эффективно распределить работу.

Здесь принудительное ООП Java и Javadoc для документации дают ощутимый прирост в копилке аргументов “За Java”.

В PHP эти возможности так же есть. Сравнительно недавно я узнал о аналогах Javadoc в PHP- phpdoc/PhpDocumentor.

С известных версий PHP поддерживает объектный подход(ООП). Правда он смешанный, и соблазн написать простую функцию без обертки очень большой.

Техническая площадка.

С кодом мы разобрались. А вот что касается хостинга, Java в проигрыше. PHP-хостинг распространен и он дешевле.

Правда под вопросом, стоит ли экономить на хостинге, если проект серьезный? Тем не менее, то, что заменить техническую основу проекта на PHP проще – факт.

Другой вопрос, если проект используется в интрасети. Стоимость хостинга отпадает. И программирование на стороне клиента дает плюс в сторону Java. Интеграция апплетов на нем в сети интернет потерпела крах, так как для запуска они требуют имплементатор на локалхосте. Но когда клиенты доступны физически, проблем нет.

Единая платформа, без связки разных языков – очевидный плюс.

Вместо заключения.

Я постарался не затрагивать сравнения технических характеристик языков. Такой информации в интернете полно. Так же считаю, что качество разработки значительно важнее. Язык создан для человека, а не для процессора. Для второго есть ассемблер.

Претензий на объективность у меня нет. Я попытался, насколько успешно- не знаю. Буду рад аргументированной критике.

Под конец повторюсь: все в руках разработчика, язык лишь инструмент. Профессионал и пальцем болт закрутит.

Нашей целью является объединить людей этой профессии для оказания помощи друг другу в прохождении собеседования, составлении резюме, накоплении знаний.

В разделе “Вопросы на собеседовании” мы будем выкладывать самые распространенные вопросы, которые задает работодатель соискателю, откликнувшемуся на его вакансию программиста. Раздел будет постоянно пополняться, сейчас в нем можно найти вопросы по PHP, C++, .NET и базам данных.

В разделе “Вакансии” вы сможете найти самые свежие объявления от работодателей для программистов. Раздел обновляется несколько раз в день.

Мы также предлагаем вам проверить свои знания, посетив раздел “Тесты”. Зарегистрированные пользователи при желании могут разместить свой тест, который им предлагали пройти на собеседовании, а затем в личном кабинете следить за статистикой прохождения теста другими пользователями сайта.

В разделе “Книги” собраны книги для программистов. Здесь можно найти литературу по самым разным IT-тематикам: устройство компьютера и работа с “железом”, языки программирования и алгоритмы, компьютерные сети, а также многое другое.

Если вы опытный программист:

Если вы начинающий программист:

Итак,

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

Присоединяйтесь к нашей группе Вконтакте.

Следите за нашими вакансиями в Twitter  

Желаем удачи,
Команда Job-Interview.ru