IMPROVING ALGORITHMS OF VIDEO SEQUENCE DATA RECOGNITION FOR IDENTIFICATION OF TRANSITION PROCESSES IN A FLOATATION MACHINE OF POTASSIUM ORE
Journal: VESTNIK OF ASTRAKHAN STATE TECHNICAL UNIVERSITY. SERIES: MANAGEMENT, COMPUTER SCIENCE AND INFORMATICS
№ 1
, 2019
Rubrics: MANAGEMENT, MODELING, AUTOMATION
Authors:
1.
Perm National Research Polytechnic University
2.
Berezniki branch of Perm National Research Polytechnic University
Type:
Article
Pages:
from 26
to 39
Status:
Published
Received:
05.01.2019
Accepted:
25.01.2019
Published:
25.01.2019
Subject area:
VAC 05.12.2013 Системы, сети и устройства телекоммуникаций
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 658.562.3:004.9
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 658.562.3:004.9
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
Language:
Russian
Keywords:
potassium ore, flotation, floatation machine, foam, computer vision, binarization, transient process, algorithm
Abstract (English):
Potash fertilizers are important for the Russian national agriculture and have become an export item. This fact results in increasing potash fertilizer production and improving potassium procession management. The object of research is floatation processes of potassium in the example of “Uralkaliy”, PJSC (Berezniki, Perm region). The aim of the research is improving algorithms of bubble recognizing in the video stream and using them to identify transient processes and situations in a flotation machine. Methods of researches include the system analysis, mathematical modeling, regression analysis, elements of automatic control theory and object identification. Algorithms for recognizing foam in the sylvinic floatation machine have been modified, which significantly increased the speed of recognizing bubbles in images of the low-quality video stream. Experiments were carried out on laboratory and industrial flotation machines, the results showing the possibility of using modified algorithms both in laboratory and industrial conditions. Video sequences of such quality were obtained and processed on the industrial floatation machine and could be used on the industrial flotation machine to identify situations and to control the machine operation. Using modified algorithms in experimental data processing allowed to identify the transient process and to clarify the time of the transient process. It has been shown that the obtained values are comparable with the data of other researchers and are not at variance with the experimental data. The error of bubble recognition has been estimated. The ways of using the data obtained for the decision support systems of the floater or of the automated control systems of the floatation machine have been shown.
Potash fertilizers are important for the Russian national agriculture and have become an export item. This fact results in increasing potash fertilizer production and improving potassium procession management. The object of research is floatation processes of potassium in the example of “Uralkaliy”, PJSC (Berezniki, Perm region). The aim of the research is improving algorithms of bubble recognizing in the video stream and using them to identify transient processes and situations in a flotation machine. Methods of researches include the system analysis, mathematical modeling, regression analysis, elements of automatic control theory and object identification. Algorithms for recognizing foam in the sylvinic floatation machine have been modified, which significantly increased the speed of recognizing bubbles in images of the low-quality video stream. Experiments were carried out on laboratory and industrial flotation machines, the results showing the possibility of using modified algorithms both in laboratory and industrial conditions. Video sequences of such quality were obtained and processed on the industrial floatation machine and could be used on the industrial flotation machine to identify situations and to control the machine operation. Using modified algorithms in experimental data processing allowed to identify the transient process and to clarify the time of the transient process. It has been shown that the obtained values are comparable with the data of other researchers and are not at variance with the experimental data. The error of bubble recognition has been estimated. The ways of using the data obtained for the decision support systems of the floater or of the automated control systems of the floatation machine have been shown.
Keywords:
potassium ore, flotation, floatation machine, foam, computer vision, binarization, transient process, algorithm
potassium ore, flotation, floatation machine, foam, computer vision, binarization, transient process, algorithm
Калийные удобрения, помимо увеличения урожайности, улучшают качественные характеристики выращиваемой продукции: способствуют повышению сопротивляемости растений к заболеваниям, повышению лежкости плодов при хранении и стойкости при транспортировке, а также улучшению их вкусовых и эстетических качеств. Для России производство калийных удобрений является важным источником валютных поступлений. Общемировую важность этой продукции подтверждает то, что российская индустрия минеральных удобрений, в отличие от нефтегазодобывающей отрасли, практически не ощутила влияния западных санкций. По оценкам аналитиков, сельскохозяйственное производство до 2050 г. должно вырасти на 70 %. В связи с этим необходимо повышать продуктивность и эффективность использования земельных угодий, в том числе за счет роста мирового ежегодного потребления калийных удобрений на 2-3 %. В настоящее время в России калийные удобрения производит ПАО «Уралкалий», заканчиваются пусконаладочные работы в ООО «Еврохим - Усольский калийный комбинат». Но существенная проблема последнего - низкое качество исходной руды, а именно большой разброс содержания KCl и нерастворимого остатка в руде. В работе [1] предложено решение проблемы путем специфической рассыпки и выемки руды со склада; в проект предприятия заложены также и другие способы бункерного усреднения руды, но окажутся ли они достаточно эффективными, покажет время. В частности, переменный состав руды может привести к тому, что в пульпе, поступающей на обогащение, будут наблюдаться колебания состава, и это потребует специальных методов идентификации ситуаций и управления процессами. Состояние проблемы Для обогащения калийных руд на Верхнекамских калийных комбинатах применяют пенную флотацию [2], в том числе используются шламовая флотация (отделение хлористых солей от нерастворимого остатка) и сильвиновая флотация (разделение хлорида калия и хлорида натрия) (рис. 1). Рис. 1. Краткая схема стадий флотационного процесса на примере флотационной обогатительной фабрики Второго Березниковского калийного производственного рудоуправления ПАО «Уралкалий», г. Березники Пермского края [3] Процесс сильвиновой флотации является головным в цепи технологических операций обогащения калийного сырья, поэтому основная задача при управлении этим процессом - обеспечение заданной производительности и поддержание необходимого качества концентрата на выходе как отдельной флотомашины (ФМ), так и всей технологической секции. Моделирование ФМ необходимо как для улучшения установок процесса при изменении пульпы (в том числе состава, как на «Еврохиме»), так и для улучшения управления при колебаниях в пределах нормы, в том числе в случае, когда существующие средства выравнивания руды переменного состава не справляются с поставленной задачей. На производстве процессы в ФМ наблюдаются визуально. В статье [4] на основании анализа видеографической информации с ФМ уточнена стохастическая модель ячейки ФМ. Кроме того, флотатор производит зрительную оценку пены, принимает решения на основании увиденного, что может вызвать определенные проблемы [5]. Вид пены (цвет, плотность, количество и размер пузырьков) имеет значение и может послужить основой для принятия решений не только в ФМ, но и, например, при производстве битумов [6]. Поэтому задача уменьшения влияния человеческого фактора при контроле вида пены в разных производствах представляется практически важной и актуальной. Применение систем технического зрения (СТЗ) позволяет этих проблем избежать, в частности, вовремя распознавать переходные процессы, начинающиеся в ФМ из-за изменения состава пульпы. Далее с помощью модели, встроенной в систему поддержки принятия решений, можно управлять процессом [7]. Конечно, параметры модели также необходимо уточнить с помощью технического зрения. Это также послужит доказательством работоспособности самой СТЗ. Системы технического зрения используются при управлении процессами флотации, но чаще в ФМ при флотации угля или полиметаллических руд [7-9]. В статье [10] определяются параметры пены не с поверхности пенного слоя, а в трехмерном измерении (с учетом глубины), что создает ряд проблем вычислительного характера. Пожалуй, одной из первых статей, в которой описываются непреодолимые вычислительные сложности при обработке видеопотока с пены, является [11], где делается вывод, что реконструкционные алгоритмы (восстанавливающие границы пузырька) успешно могут применяться только к статическим фотографиям. Работа написана по результатам исследования изображений, полученных при угольной флотации. В источнике [12] описывается определение параметров пены при грубой флотации медной руды каскадным методом. По видеоряду идентифицируют тип руды, в зависимости от которого по цветовым характеристикам определяют параметры процесса. В [13] c помощью машинного зрения контролируют флотацию угля (зависимость между флотационным концентратом и содержанием золы), однако приведенные в работе алгоритмы не могут быть применены непосредственно к флотации калийной руды из-за разницы в RGB-распределениях в изображениях. Автором [14] предложен метод определения размеров пузырьков, позволяющий идентифицировать пузырьковые кластеры, мелкие и большие пузырьки, а также неполные пузыри в рамке изображения. Следует отметить, что многие описываемые в литературе подходы к применению видеоанализа для контроля параметров флотации и использования их в целях автоматического управления приводятся скорее на уровне описания идеи, чем успешного внедрения, хотя и авторы вышеперечисленных работ, и автор [15] описывают успешную проверку своего алгоритма на идентификации целого и суб-избражения пены в промышленной ФМ. Пена ФМ полиметаллических руд плотная, глянцевая, границы пузырей хорошо распознаются. Пена в ФМ сильвиновой флотации неконтрастная, плотная, с большим разбросом диаметров пузырей. Это существенно затрудняет распознавание границ пузырей, используемое во всех вышеперечисленных работах, приводит к большим вычислительным затратам. В работе [16] предложено перейти от распознавания границ пузыря к сбору статистики по размерам бликов и расстояниям между ними. Этот подход развит в [4] до пробного программного обеспечения, использованного для распознавания фотографий пены, сделанных со вспышкой, а также искусственно созданной и подкрашенной пены. Вычислительные проблемы, связанные с уточнением распознавания цветных пен, хорошо известны. В [17] для компенсации некачественного изображения авторы предлагают матрицы вспомогательного пересчета серых и цветных изображений. К сведению, в этой статье в действительности используются очень и очень качественные изображения, полученные только на лабораторной ФМ. В [18] для повышения распознавания качества полиэдральных пузырей, каким-то образом характеризующим угольную флотацию, авторы также создали специальную лабораторную ФМ с многоисточниковым освещением и смогли добиться только пятисекундного периода обработки изображений. Само по себе это значение весьма удовлетворительное для промышленной ФМ, но вот создать такие идеальные условия съемки на ней точно не удастся. В [19] авторы для повышения качества распознавания тех же статических изображений, что и в [17], предлагают элементы искусственного интеллекта. В целом проблемы с распознаванием изображений пены, даже статических, полученных в лабораторных условиях и на высококонтрастной пене, широко описаны в научных источниках. При переходе к обработке видеопотока реального времени возникает ряд дополнительных проблем, сформулированных и частично решенных в [20, 21], где эксперименты производились на лабораторной ФМ. Дальнейшее ускорение расчетов, связанных с идентификацией бликов, и алгоритм исключения пеногона предложены в [22]. Для последней работы съемки проводились на промышленной ФМ ПАО «Уралкалий». Таким образом, целью настоящей работы является дальнейшее улучшение алгоритмов распознавания пузырей в видеопотоке и использование их для идентификации переходных процессов в ФМ калийной руды. Эксперимент и его обработка Важным вопросом при настройке алгоритма распознавания бликов является выбор порога бинаризации. Все пиксели изображения в градациях серого, имеющие яркость ниже этого порога, считаются черными, остальные - белыми. Последние как раз и создают блик, а по разным статистикам бликов определяется начало и окончание переходного процесса в ФМ. В [21] показано, что для лабораторной ФМ зависимость оптимального порога бинаризации, при котором распознается наибольшее количество бликов в кадре, от освещенности кадра может быть приближенно выражена линейно, и при этом алгоритмы обработки распознают переходный процесс в камере. Данный метод использован для видеопотока с промышленной ФМ, место которой в технологическом процессе описано выше. Обычно ФМ моделируют как последовательность ячеек идеального перемешивания [23]. При шламовой флотации с пеной уносится нерастворимый остаток, при сильвиновой - полезный продукт (рис. 2). Рис. 2. Схема флотационной камеры: q - объемный расход пульпы; p - доля полезного компонента; n - доля загрязняющего компонента в пульпе; Qp - объемный расход полиакриламида. Индексы 1-3 обозначают входной и два выходных потока флотомашины соответственно Съемка производилась на 1-й и 3-й камерах ФМ, с искусственным точечным источником света и без него (при естественном цеховом освещении). Флотационные агенты подаются именно в 1-ю и 3-ю камеры. Конкретные режимные параметры процесса не могут быть опубликованы в открытой печати, однако в 1-ю камеру их расход в несколько раз больше, чем в 3-ю. В остальных камерах ФМ происходит дофлотация. Использовалась широкоугольная видеокамера высокого разрешения 1 980 × 1 200 точек с частотой кадров 50 fps (рис. 3). Рис. 3. Исходное изображение 3-й камеры флотомашины Из полученного кадра «вырезалось» окно между осью пеногона и станиной камеры, в котором производилось выделение бликов. В терминах использованной библиотеки OpenCV [24] подобное окно называется Region Of Interest (ROI), для обработки изображений в ROI в библиотеке имеется множество функций. При съемке 1-й камеры ФМ с искусственным источником света расход реагентов был уменьшен на 1/3. После того, как наблюдаемый визуально переходный процесс завершался и небольшой дополнительной выдержки расход реагентов восстанавливался, начинался обратный переходный процесс. При съемке 3-й камеры подача реагентов полностью прекращалась на сопоставимое время, затем восстанавливалась. Переходный процесс при этом визуально не наблюдался. Далее производилась съемка 5-й камеры с источником света, 1-й и 3-й камер при естественном освещении. Эти фрагменты использовались для тестирования алгоритма. Каждое видео обрабатывалось несколько раз. При этом решались следующие задачи: 1. «Подгонка» частоты выборки кадров под частоту вращения пеногона (в данном случае она оказалась не кратной частоте кадров). В результате получали настройки алгоритма поиска «перевала освещенности», описанного в [22], при которых вертикальное положение вращающегося пеногона определялось уверенно и безошибочно. Настройки, полученные при обработке видео 1-й и 3-й камер с точечным источником, оказались приемлемыми для всех видеофрагментов. 2. Определение оптимального порога бинаризации, при котором алгоритм распознает в ROI наибольшее количество бликов, в зависимости от средней освещенности окна. В результате получали линейную регрессию, связывающую порог и освещенность. 3. При найденном оптимальном пороге бинаризации в зависимости от освещенности окна определялось количество пузырей в кадре в зависимости от номера кадра. Таким образом, для первых двух фрагментов определяли параметры переходного процесса, а на остальных убеждались в его отсутствии. При установленном (а в задаче 2 меняющемся) пороге бинаризации каждый кадр подвергался следующей обработке: 1. Вырезалось изображение в заданном окне. 2. Изображение переводилось в градации серого функцией cvtColor (здесь и далее речь идет о функциях библиотеки OpenCV). 3. По освещенности окна и другим настройкам алгоритма определения положения пеногона определялось, расположен ли он вертикально. Если нет, дальнейшая обработка не производилась, происходил переход к следующему кадру. 4. Производилась бинаризация изображения функцией cvThreshold. 5. Изображение копировалось в бинарную матрицу функцией cvarrToMat. 6. На бинаризованном изображении алгоритмом движущегося угла «нумеровались» пузыри. 7. Исключались все пузыри слишком большого и слишком малого размера (области засветки, возникшие из-за случайной волны или слишком низкого порога бинаризации, и шум). 8. Оставшиеся пузыри учитывались в статистике (количество, размер, расстояния между пузырями и т. д.). Результаты исследований и их обсуждение При обработке видео с разных камер промышленной ФМ возникли новые проблемы, не отмечавшиеся ранее. Например, потребовалась модификация метода распознавания момента параллельного к оси камеры положения пеногона, т. к. на одной из камер за полуоборот освещенность несколько раз достигала максимального значения (рис. 4). Рис. 4. Зависимость освещенности на 1-й и 3-й камерах флотомашины от номера кадра: А - «острые» пики освещенности; Б - соответствуют искомому положению пеногона «Острые» пики освещенности, отмеченные «А», соответствуют блику на пеногоне, когда он стоит почти перпендикулярно оси камеры и отражает источник света в объектив. Соответственно, в алгоритме, приведенном ранее в [22], пришлось добавить условие ограничения максимального изменения освещенности в процессе приближения к максимуму, и он стал уверенно распознавать точки «Б». Для всех камер был определен оптимальный порог бинаризации, обеспечивающий распознавание наибольшего количества бликов, в зависимости от освещения кадра. Типичный график зависимости оптимального уровня бинаризации (ось ординат) от освещенности приведен на рис. 5 (для его построения полученные экспериментальные результаты были упорядочены по росту освещенности кадра). Рис. 5. Зависимость оптимального уровня бинаризации от освещенности кадра на 1-й камере с источником освещения Очевидно, что для 1-й камеры общая тенденция зависимости оптимального порога бинаризации от средней освещенности кадра близка к линейной. Коэффициент корреляции для тренда на рис. 5 равен 0,65, его невысокое значение можно объяснить ошибкой определения освещенности и целочисленностью порога бинаризации. Процедура определения оптимального порога достаточно трудоемкая, т. к. ROI обрабатывается последовательно при уровнях бинаризации от 110 до 252, т. е. оптимизация происходит методом полного перебора. На компьютере, на котором производилась обработка, одно оптимальное значение получается за 70-300 с, что недопустимо при обработке реального видеопотока. Поэтому для данных, полученных при обработке каждого эксперимента, методом наименьших квадратов были рассчитаны коэффициенты a и b зависимости оптимального порога бинаризации p от средней освещенности кадра s: . (1) Возможность замены оптимального порога бинаризации на расчетный проверена следующим образом. При расчетном пороге бинаризации для каждой камеры повторно было произведено определение количества пузырей, что на один ROI занимало 0,7-1,1 с. Это время меньше времени полуоборота пеногона, следовательно, в дальнейшем ускорении обработки нет необходимости. Результаты в целом соответствуют полученным при оптимальном пороге (рис. 6). Рис. 6. Фрагмент тренда зависимости количества пузырей в кадре при оптимальном и расчетном порогах бинаризации Полученные коэффициенты линейной зависимости (1) и погрешность определения количества пузырей в ROI приведены в таблице. Экспериментальные данные и расчет коэффициентов линейной зависимости между средней освещенностью кадра и оптимальным порогом бинаризации Параметры Лабораторная ФМ Камера 1 Камера 3 Камера 5 Камера 1 Камера 3 Наличие источника освещения да да да да нет нет Средняя освещенность s 212,24 119,72 129,82 138,79 105,6 108,86 Количество кадров нет данных 236 204 67 313 312 Коэффициент a -0,204 0,8828 1,7992 1,0017 1,1536 1,6064 Коэффициент b 273,9 23,9217 -79,6112 19,8192 -8,7088 -37,1921 Погрешность распознавания пузырей, % нет данных 8,28 9,51 26,18 7,78 7,79 Наибольшая погрешность получена при съемках 5-й камеры, что объясняется краткостью видеофрагмента (67 полуоборотов пеногона). Очевидно, что при росте количества кадров свыше 200 погрешность определения количества пузырей приходит к асимптотическому значению порядка 7-8 %, что вполне сопоставимо с шумом исходных данных, а потому допустимо (рис. 7). Рис. 7. Зависимость относительной погрешности от количества кадров в выборке Так, из рис. 6 очевидно, что в отсутствие переходного процесса количество пузырей в кадре меняется примерно на ±50, т. е. на ±25 % от среднего значения. Зависимость коэффициентов (1) от освещенности показана на рис. 8 (значения для графика упорядочены по росту освещенности кадра). Рис. 8. Зависимость линейных коэффициентов от средней освещенности кадра Коэффициент корреляции между освещенностью и коэффициентом a равен 0,8336, между освещенностью и коэффициентом b - 0,9090. Однако если, в свою очередь, рассчитать уже зависимости самих коэффициентов от освещенности в виде a(s) = A + Bs или a(s) = A + Bs + Сs2 и использовать расчетные значения для определения порога бинаризации, погрешность определения количества пузырей составляет уже десятки процентов. С использованием расчетного порога бинаризации были обработаны все видеозаписи. Как принято при проведении идентификационных экспериментов в автоматизации технологических процессов, полученные данные были нормированы и сглажены по пяти точкам. В результате на 1-й и 3-й камерах получили очевидно выраженные переходные процессы (рис. 9), вызванные, как изложено выше, изменением подачи реагентов (по оси абсцисс - номер кадра, полученного при правильном положении пеногона; частота вращения пеногона порядка 1,75 с на полуоборот). Рис. 9. Переходные процессы на 1-й и 3-й камерах флотомашины Если рассматривать ячейки как объекты первого порядка и исключить запаздывание, то при идентификации объекта получим его уравнение , где T - постоянная времени: T ≈ 8,3 с для 1-й камеры и T ≈ 8,6 с для 3-й камеры. В [4] использовалось значение T = 7,6 с для любой камеры ФМ, полученное ранее в [23, с. 72]. В этой же монографии в другой главе используется также значение T = 10 с, полученное в других экспериментальных условиях [23, с. 95]. Следовательно, можно констатировать, что полученные путем анализа видеоряда данные как минимум не противоречат предыдущим имеющимся данным. Как и ожидалось, при отсутствии изменений на входе в ФМ среднее количество пузырей остается неизменным (рис. 10). Рис. 10. Фрагмент тренда зависимости количества пузырей в кадре с 3-й камеры флотомашины без дополнительного источника освещения Во всяком случае, явно выраженный переход от одного среднего значения к другому отсутствует. На рис. 10 представлен тренд с 3-й камеры без дополнительного источника освещения. С использованием источника освещения среднее отклонение при условии отсутствия переходного процесса примерно вдвое меньше, чем без источника. Сравнение рис. 6 и рис. 10 является подтверждением этого факта. Таким образом, выбор источника нужной яркости является важной задачей при дальнейшей доработке метода. Было также произведено исследование влияния переходного процесса на среднюю интенсивность красной (R), зеленой (G) и голубой (B) компоненты в кадре и на статистические характеристики расстояний между пузырями. Для расщепления на R-, G- и B-составляющие кадр до перевода в градации серого рассматривался как совокупность трехкомпонентных векторов, описывающих каждый пиксель. В терминах библиотеки OpenCV это модель BGR, т. е. первая компонента вектора соответствует интенсивности голубого цвета и т. д. Более сложным является составление статистических характеристик расстояний между пузырями. Так как пузырей в кадре достаточно много (порядка 300), все расстояния между всеми пузырями дадут массив длиной 3002 = 90 000. Собственно, все они и не нужны, т. к. нас интересует расстояние только между ближайшими пузырями, или центрами бликов. Поиск ближайших бликов осуществлялся по следующему алгоритму: 1. Для каждого перенумерованного блика находился его центр как середина охватывающего прямоугольника, как это предложено в [21]. 2. Для каждого центра блика составлялся массив расстояний до всех остальных центров бликов. Из этого массива выбирались несколько наименьших расстояний, после ряда экспериментов пришли к выводу, что переходные процессы стабильно наблюдаются при 8 наименьших расстояниях. Чтобы 300 раз не сортировать массив длиной 300 элементов для выбора наименьших значений, использовали следующее: - при добавлении каждого расстояния в массив контролировалось и изменялось при необходимости минимальное и максимальное расстояние до бликов; - в выборку добавлялось минимальное значение, счетчик выборки увеличивался на 1; - после заполнения массива в цикле от минимального плюс 1 до максимального расстояния массив перебирался, производился поиск расстояния, равного счетчику цикла. Если такое расстояние найдено, оно добавлялось в выборку, счетчик выборки увеличивался на 1, в исходном массиве значение расстояния обнулялось, цикл начинался сначала; - при каждом добавлении контролировалось значение счетчика выборки; как только оно достигало 8, выборка наименьших расстояний прекращалась. 3. Затем вся выборка из 8 наименьших расстояний копировалась в глобальный массив расстояний и п. 2 повторялся до окончания всех бликов. 4. Глобальный массив сортировался для расчета среднемедианного расстояния, заодно вычислялось среднее арифметическое значение всех его элементов. Обнаружено, что значения RGB-компонент могут в некоторых случаях служить индикаторами начала переходного процесса. Так, компонента R (red, красная) на рис. 11 скачкообразно изменяется в момент, когда восстанавливается прежняя подача реагентов. Рис. 11. Переходный процесс Все три компоненты меняются скачкообразно одновременно, поэтому контролировать можно любую из них, хотя наибольшее изменение наблюдалось именно у красной компоненты. Среднее арифметическое расстояние между пузырями, полученное вышеописанным способом, хорошо характеризует переходный процесс, как видно на рис. 11. Среднемедианное расстояние также изменяется в ходе переходного процесса, но, как и RGB-компоненты, не так отчетливо и не в такой степени. Поэтому из вышеописанного алгоритма можно исключить п. 4 и ограничиться контролем цветовых компонент и среднего арифметического расстояния между пузырями как вспомогательными параметрами, характеризующими наличие переходного процесса в флотационной камере. Заключение Таким образом, цель, поставленная в данной статье, достигнута. Произведены модификации алгоритмов распознавания пены в машине сильвиновой флотации. Произведены эксперименты на лабораторной и промышленной флотомашинах. Обработка данных эксперимента с использованием модифицированных алгоритмов позволила идентифицировать сам переходный процесс и определить время переходного процесса значением, сопоставимым с данными других авторов. Оценка погрешности распознавания и дальнейшее уточнение параметров флотомашины как объекта управления возможны с использованием большего количества видеоданных.
1. Grigalashvili A. S. Realizaciya modeli zapolneniya i razgruzki usrednitel'nogo sklada rudy // Sovremennaya nauka: aktual'nye problemy teorii i praktiki. Ser.: Estestvennye i tehnicheskie nauki. 2016. № 1. S. 44-48.
2. Teterina N. N., Sabirov R. H., Skvirskiy L. Ya., Kirichenko L. N. Tehnologiya flotacionnogo obogascheniya kaliynyh rud. Perm': Perm. gos. tehn. un-t, 2002. 484 s.
3. Kurmaev R. H. Flotacionnyy sposob polucheniya hlorida kaliya iz sil'vinita. Perm': Perm. gos. tehn. un-t, 1993. 83 s.
4. Zatonskiy A. V., Varlamova S. A., Malysheva A. V., Myasnikov A. A. Ispol'zovanie videograficheskoy informacii dlya utochneniya dinamicheskoy stohasticheskoy modeli processa flotacii kaliynoy rudy // Naukovedenie. 2017. № 2 (39). S. 87.
5. Zatonskiy A. V., Varlamova S. A. Use of reflection flare spots for automatic recognition of froth parameters in potassium ores flotation // Obogashchenie Rud (Mineral processing). 2016. № 2. P. 49-56. DOI:https://doi.org/10.17580/or.2016.02.09.
6. Malysheva A. V. Vizual'nyy komp'yuternyy kontrol' processa proizvodstva bituma // Dostizheniya, problemy i perspektivy razvitiya neftegazovoy otrasli: materialy Mezhdunar. nauch.-prakt. konf., posvyasch. 60-letiyu vysshego neftegazovogo obrazovaniya v Respublike Tatarstan (Al'met'evsk, 28-29 oktyabrya 2016 g.). Al'met'evsk: Al'met'ev. gos. neftyanoy in-t, 2016. S. 171-172.
7. Outotec FrothSenseTM. URL: https://www.outotec.ru/globalassets/ru/brochures/outotec-sensornaya-sistema-frothsense-pdf (data obrascheniya: 21.12.2018).
8. Flotation optimization. URL: https://www.metso.com/services/optimization/flotation (data obrascheniya: 21.12.2018).
9. Romanenko S. A., Olennikov A. S. Opyt vnedreniya videokamer FROTHMASTER na obogatitel'noy fabrike «Zelenaya gora-2» // Obogaschenie rud. 2014. № 2. S. 23-28.
10. Xu D., Chen Y., Chen X., Xie Y., Yang C., Gui W. Multi-model soft measurement method of the froth layer thickness based on visual features // Chemometrics and Intelligent Laboratory Systems. 2016. V. 154. P. 112-121. DOI: https://doi.org/10.1016/j.chemolab.2016.03.029.
11. Sadr-Kazemi N., Cilliers J. J. An image processing algorithm for measurement of flotation froth bubble size and shape distributions // Minerals Engineering. 1997. V. 10. Iss. 10. P. 1075-1083.
12. Horn Z. C., Auret L., McCoy J. T., Aldrich C., Herbst B. M. Performance of Convolutional Neural Networks for Feature Extraction in Froth Flotation Sensing // IFAC-PapersOnLine. 2017. V. 50. Iss. 2. P. 13-18. DOI: https://doi.org/10.1016/j.ifacol.2017.12.003.
13. Tan J., Liang L., Peng Y., Xie G. The concentrate ash content analysis of coal flotation based on froth images // Minerals Engineering. 2016. V. 92. P. 9-20. DOI: https://doi.org/10.1016/j.mineng.2016.02.006.
14. Zhao L., Peng T., Zhao Y., Xia P., Xie Y. Features Extraction of Flotation Froth Based on Equivalent Binocular Stereo vision // IFAC-PapersOnLine. 2016. V. 49. Iss. 20. P. 90-95. DOI: http://dx.doi.org/https://doi.org/10.1016/j.ifacol.2016.10.102.
15. Jahedsaravani A., Massinaei M., Marhaban M. H. An Image Segmentation Algorithm for Measurement of Flotation Froth Bubble Size Distributions // Measurement. 2017. V. 111. P. 29-37. DOI: http://dx.doi.org/https://doi.org/10.1016/j.measurement.2017.07.023.
16. Zatonskiy A. V., Varlamova S. A. Primenenie tehnicheskogo zreniya dlya upravleniya processom flotacii // Virtual'noe modelirovanie, prototipirovanie i promyshlennyy dizayn: materialy II Mezhdunar. nauch.-prakt. konf. (Tambov, 17-19 noyabrya 2015 g.). Tambov: Tamb. gos. tehn. un-t, 2016. S. 313-318.
17. Aldrich C., Smith L. K., Verrelli D. I., Bruckard W. J., Kistner M. Multivariate image analysis of realgar - orpiment flotation froths // Mineral Processing and Extractive Metallurgy. 2017. V. 127 (3). P. 146-156. DOI:https://doi.org/10.1080/03719553.2017.1318570.
18. Tan J., Liang L., Xia W., Xie G. The effect of different flotation operating parameters on the froth properties and their relation to clean coal ash content // Separation Science and Technology. 2018. V. 53 (9). P. 1434-1444. DOI:https://doi.org/10.1080/01496395.2018.1444051.
19. Fu Y., Aldrich C. Froth image analysis by use of transfer learning and convolutional neural networks // Minerals Engineering. 2018. V. 115. P. 68-78. DOI:https://doi.org/10.1016/j.mineng.2017.10.005.
20. Myasnikov A. A., Gaynutdinova K. V. Ocenka processa flotacii s ispol'zovaniem komp'yuternogo zreniya // Pervyy shag v nauku. 2016. № 11. S. 32-36.
21. Zatonskiy A. V., Malysheva A. V. Modernization of algorithms for flare detection of froth layer parameters during flotation of potash ores // Obogashchenie Rud (Mineral processing). 2017. N. 2. P. 35-41. DOI:https://doi.org/10.17580/or.2018.02.07.
22. Malysheva A. V. Usovershenstvovannyy algoritm raspoznavaniya blikov binarizovannogo izobrazheniya peny // Sovremennaya nauka: aktual'nye problemy teorii i praktiki. 2018. № 5. S. 110-116.
23. Bekker V. F. Reshenie tehnologicheskih problem deystvuyuschego proizvodstva sredstvami avtomatizacii. T. 1: Proizvodstvo kaliynyh udobreniy. Perm': BF PNIPU, 2012. 312 s.
24. Open Source Computer Vision Library. URL: https://opencv.org (data obrascheniya: 30.08.2018).
