EA035343B1 - Способ определения объема зерна в зернохранилищах на основе бесконтактных 3d-измерений - Google Patents

Способ определения объема зерна в зернохранилищах на основе бесконтактных 3d-измерений Download PDF

Info

Publication number
EA035343B1
EA035343B1 EA201700618A EA201700618A EA035343B1 EA 035343 B1 EA035343 B1 EA 035343B1 EA 201700618 A EA201700618 A EA 201700618A EA 201700618 A EA201700618 A EA 201700618A EA 035343 B1 EA035343 B1 EA 035343B1
Authority
EA
Eurasian Patent Office
Prior art keywords
grain
granary
volume
points
camera
Prior art date
Application number
EA201700618A
Other languages
English (en)
Other versions
EA201700618A1 (ru
Inventor
Кайрат Темирбаевич Кошеков
Наталья Владимировна Астапенко
Александр Валентинович Демьяненко
Original Assignee
Кайрат Темирбаевич Кошеков
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Кайрат Темирбаевич Кошеков filed Critical Кайрат Темирбаевич Кошеков
Publication of EA201700618A1 publication Critical patent/EA201700618A1/ru
Publication of EA035343B1 publication Critical patent/EA035343B1/ru

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F22/00Methods or apparatus for measuring volume of fluids or fluent solid material, not otherwise provided for
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

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

Description

Изобретение относится к сельскому хозяйству, а именно к обеспечению количественнокачественной сохранности зерна. Для обеспечения количественно-качественной сохранности зерна традиционно проводится только визуальный осмотр мест хранения зерна. Однако практика показала, что традиционные формы организации работы зернохранилищ выявляют такие глобальные проблемы как хищение зерна и сохранение качества зерна.
При определении наполненности зернохранилища необходимо учитывать некоторые особенности работы с сыпучими продуктами, такие как наличие сложной формы поверхности продукта (наличие горки, воронки, кусковые слипнувшиеся комки).
Известны два способа измерения объемов неровных тел: вычисление объема тела при помощи емкости с водой и вычисление объема тела через разбиение его на правильные геометрические фигуры (http//u ztes.t ru abstracts/idabst ract=216487 (просмотр 22.06.2017). Недостаток этих способов расчета в том, что первый способ очень сложно организовать с габаритными объектами, реализация же второго представляет собой сложную задачу, т.к. нельзя явно выделить правильные геометрические фигуры.
Известен способ хранения зерна (патент РК №31368. Опубл. 2016, №8). В известном способе для контроля процесса хранения зерна в зернохранилище используют автоматизированную систему управления для определения температуры зерновой массы и при наличии местного перегрева производят изъятие его, не определяя при этом объема оставшегося качественного зерна. Недостаток известного способа в трудоемкости процессса и отсутствии показателя объема зерна после изъятия перегревшего.
Задача, на решение которой направлено заявляемое изобретение, заключается в том, что оно позволяет осуществлять ежедневный контроль заполнения зернохранилища, проводить анализ полученных данных, выявляя возможную недостачу зерна, осуществлять контроль качества зерна.
Технический результат достигается за счет использования видеокамер и набора лазерных диодов, позволяющих сформировать массив 3D точек поверхности зерна на регулярной основе.
Сущность изобретения заключается в осуществлении способа измерения объема зерна в зернохранилище на основе бесконтактных 3D измерений и получения координат точек поверхности, образующейся в зернохранилище в результате его заполнения.
Способ осуществляется в виде последовательности следующих этапов:
Этап 1 Этап 2 Этап 3
Вычисление объема поверхности
В 1 этап входит аналого-цифровое преобразование данных: сканирование аналоговых оригиналов с последующей векторизацией, а также импорт готовых цифровых данных. Требуется получить данные о дальности для некоторого набора различных фрагментов, которые в совокупности содержат всю поверхность объекта. Далее происходит регистрация или совмещение карт глубин каждого фрагмента. Облако точек представляет собой набор точек поверхности, заданных тремя значениями координат х, у, z. Чем плотнее облако точек, тем точнее будет результат вычислений.
На этапе получения данных предполагается использование видеокамеры и набора лазерных диодов, позволяющих в дальнейшем получить облако точек поверхности зерна в зернохранилище. Важным условием на данном этапе является принятие решения о требованиях к видеокамере, количестве и расположении лазерных диодов.
1. Расположение видеокамеры и ее параметры (фиг. 1).
В качестве основного параметра при выборе видеокамеры учитывается угол обзора. Используя геометрические правила для прямоугольных треугольников можно определить отношение между углом обзора камеры (УО), расстоянием до поверхности (D) и шириной области (W), попадающей в обзор камеры при заданном расстоянии.
( W^l ( W
УО = 2 · arctg --- , в радианах, или УО = 2 · arctg --\2D) \2D
Полученное отношение позволяет принять решение о выборе и расположении камеры/камер для зернохранилища. Рассмотрим данный процесс на примере прототипа зернохранилища инновационного типа, вместимостью 200 т (фиг. 2,3).
Случай 1. Качество распознавания поверхности зерна требует попадания всей ширины зернохранилища на расстоянии минимум метра от камеры. Архитектура зернохранилища не позволяет поместить камеру в центре секции, так как там расположена загрузочная щель. Таким образом, на расстоянии D=1 м в обзор камеры должно попадать W=4,5 м, то есть УО должен быть более 132°.
Использование одной широкоугольной камеры в секции наиболее простой и экономичный вариант решения, однако при полном заполнении зернохранилища образуется слепая зона, которая в случае многократного перемешивания зерна может быть плохо прогнозируемой (фиг. 4).
Случай 2. Наилучший случай - это расположение двух камер по обе стороны от загрузочной щели (фиг. 5). Такое расположение предполагает, что на расстоянии D=1 м в обзор камеры должно попадать
180°
-----, в градусах, π
- 1 035343
W=2,25 м, то есть УО должен быть более 97°. В формате 4:3 вертикальный обзор на том же расстоянии составит немногим более метра.
В этом случае образуется слепая зона в границах ширины загрузочной щели в дальности менее метра. Данные в этой области можно восстановить, используя методы интерполяции данных.
Рассмотрим реализацию способа на примере сборно-разборного зернохранилища инновационного типа. Ранее мы определили, что необходимо зеркальное расположение двух камер относительно линии верхнего транспортера. Для полного сканирования зернохранилища требуется размещение камер с углом обзора 92,8° на каждом метре по длине зернохранилища или на каждые два метра камеры с углом обзора 120°.
Луч лазерного диода в виде точки проектируется на поверхность зерна, расположенную в области видения камеры (фиг. 6). Луч лазерного диода параллелен оптической оси камеры. Простой алгоритм ищет на изображении яркие пиксели, так как освещенная точка является яркой на фоне более темной обстановки в силосе. Расчет расстояния до поверхности основывается на том, что чем ближе лазерная точка к центру изображения, тем дальше находится поверхность.
Как мы видим из фиг. 6, расстояние (D) может быть рассчитано по формуле:
D = —.
tgX
То есть необходимо точно знать h - фиксированное расстояние между лазерным диодом и камерой. Знаменатель рассчитывается так:
λ = pfc · rpc + го где pfc - номер светлого пикселя от центра изображения, rpc - радиан на пиксель, ro - компенсация. h
Таким образом: D = tg(pfc · rpc + го)
Количество пикселей от центра изображения до центра яркой точки может быть просто рассчитано с картинки. Для получения других параметров необходимо выполнить калибровку. Для калибровки системы собираем серию измерений, где нам известна дальность до цели, а также количество пикселей центра изображения до точки, освещенной лучем лазерного диода.
Видеокамера крепится к потолку зернохранилища таким образом, чтобы линия обзора была направлена перпендикулярно поверхности крепления. Направление лазерных диодов - параллельно линии обзора видеокамеры. Расстояния между диодами устанавливаются так, что точки их лучей будут представлять собой узлы регулярной сетки. Условная схема размещения блока лазерных диодов относительно одной видеокамеры представлена на фиг. 7.
Рекомендуемое нами расстояние между лазерными диодами - 0,2 м. Количество лазеров определяется в соответствии с шириной обзора камеры на расстоянии метра. Итак, получаем М точек по вертикали и N точек по горизонтали. Расстояния от каждого лазерного диода до камеры статичны, их храним в двумерном массиве Н.
Алгоритм получения координат точек поверхности для одной камеры.
Шаг 1. Определяем количество точек на изображении M'xN'.
Шаг 2. Если количество точек на изображении равно MxN, то устанавливаем координаты для каждой точки, начиная с левого верхнего угла:
Хи (ОД), Хи (ОД,2), Χί3 (ОДД), Хш (0; N*0,2-1)
Х21 (0,2;0), Х22 (ОДД,2), Х23 (0,2Д,4), ..., X2n (0,2; N*0,2-1)
Х31 (ОДД), Хз2 (ОДД,2), Хзз (ОДДД), ..., X3n (ОД; N*0,2-1)
Χμι (М*0,2-1Д), Хм2 (М*0,2-1 Д,2), Хмз (ОДД), ..., ΧΜΝ (М*0,2-1; N*0,2-1)
Иначе, если точек меньше MxN, то устанавливаем координаты для каждой точки, начиная с левого верхнего угла:
Хп ([М-М’] *0,2-1; [N-Ν’] *0,2-1), Χί2 ([М-М’] *0,2-1; [N-Ν’] *0,2-1+0,2), Х13 ([М-М’] *0,2-1; [N-Ν’] *0,2-1+0Д), Хш’ ([М-М’] *0,2-1; [N-Ν’] *0,2-1+Ν’*0,2-1)
Χ2ι ([М-М’] *0,2-1+0,2; [N-Ν’] *0,2-1), Х22 ([М-М’] *0,2-1+0,2; [N-Ν’] *0,2-1+0,2), Х2з ([М-М’] х0,2-1+0,2; [N-Ν’] *0,2-1+0,4), Х21+ ([М-М’] *0,2-1+0,2; [N-Ν’] *0,2-1+Ν’*0,2-1)
Хз1 ([М-М’] *0,2-1+0Д; [N-Ν’] *0,2-1), Х32 ([М-М’] *0,2-1+0Д; [N-Ν’] *0,2-1+0,2), Хзз ([М-М’] *0,2-1+0Д; [N-Ν’] *0,2-1+0,4), ..., Χ’ ([М-М’] *0,2-1+0Д; [N-Ν’] *0,2-1+Ν’*0,2-1)
Χμί ([М-М’] *0,2-1+ М’*0,2-1; [N-Ν’] *0,2-1), ХМ2 ([М-М’] *0,2-1+М’*0,2-1; [N-Ν’] *0,2-1+0,2), Хм’з ([М-М’] *0,2-1+ М’*0,2-1; [N-Ν’] *0,2-1+0Д), Хм+г ([М-М’] *0,2-1+ М’х0,2-1; [N-Ν’] *0,2-1+N’*0,2-1)
Шаг 3. Для каждой точки определяем расстояние
- 2 035343 tg(pfc·rpc + ro) по методу, описанному выше, hy - расстояние между лазерным диодом и камерой из массива Н.
Предложенный алгоритм позволяет сформировать массив 3D точек поверхности зерна на регулярной сетке.
2. Способ получения 3D измерений полной поверхности зерна, образованной при заполнении зернохранилища.
Для получения 3D измерений полной поверхности зерна в зернохранилище размещают набор видеокамер с комплектом лазерных диодов для каждой на одном реле. Количество диодов для камеры и общее количество таких комплектов с камерой зависит от технических характеристик конкретной выбранной модели.
Например, для зернохранилища инновационного типа (фиг. 2, 3) в соответствии с расчетами, предложенными выше, требуется 18 комплектов камера+набор лазерных диодов, где камера с горизонтальным углом обзора 120° и вертикальным 90°, количество лазерных диодов в комплекте составляет 10 по горизонтали и 10 по вертикали.
Расположение и индексация комплектов камера+набор лазерных диодов, прикрепленных к потолку опытного зернохранилища, представлены на фиг. 8. Ориентация всех камер одинакова и выставлена таким образом, что левый нижний угол фиг. 8 является верхним левым углом для камеры 11.
Получение изображений выполняется последовательно. Активируется комплект лазерных диодов и фиксируется изображение с камеры комплекта, далее активируется комплект и фиксируется изображение с его камеры и так далее. При формировании 3d измерений необходимо только будет учитывать фиксированный сдвиг по вертикали и горизонтали для каждого комплекта. Для точек, полученных от лучей лазерных диодов, закрепленных на потолочных скатах, рассчитывая расстояние до поверхности необходимо вычесть компенсацию по сдвигу лазерного диода относительно поверхности закрепления камеры.
В результате выполнения описанной процедуры получаем упорядоченное облако точек с тремя координатами (X, Y, Z), где X - координата по горизонтали относительно верхнего левого угла комплекта^ Y - координата по вертикали относительно верхнего левого угла комплекта^ Z - расстояние до поверхности зерна.
Далее происходит регистрация или совмещение карт глубин каждого фрагмента. Если облако точек не упорядочено, необходимо привести его к упорядоченности. В зависимости от метода сбора исходные данные представляются либо в виде модели DEN, либо в виде модели TIN.
Модель DEM (Digital Elevation Model) (фиг. 9) более популярна и удобна для практического использования, так как представлена на регулярной сетке со сторонами, параллельными координатным осям X и Y системы местности. Это наиболее простой случай, значения точек поверхности могут распологаться в узлах регулярной сетки. Пример показан на фиг. 9 - представлен фрагмент объекта размером метр на метр, растояние между узлами регулярной сетки 1дм.
Модель TIN (Triangulated Irregular Network), или модель по триангуляционной нерегулярной сетке, или свободная модель (фиг. 10). Построение цифровой модели рельефа с использованием модели данных TIN сводится к созданию оптимальной сети треугольников, элементы которой стремятся быть как можно ближе к равносторонним. При этом любая точка двумерного пространства обладает только одной высотной координатой.
Предлагаемый в данной работе способ получения облака точек позволяет построить модель, близкую к DEM, но тем не менее предлагаем алгоритм расчета объемов, независимый от вида полученной модели.
Традиционно при вычислении объема поверхности суммируют объемы прямых призм с треугольным основанием и высотой, соответствующей среднему арифметическому высот вершин призмы. В разрабатываемой подсистеме мониторинга важна точность вычислений, поэтому мы предлагаем рассчитать объем зерновой массы как сумму объемов прямых косоусеченных призм с треугольным основанием, вершины которых принадлежат исходному облаку точек, образующих поверхность зерна. Данное решение позволит повысить точность вычислений.
Общий процесс функционирования подсистемы мониторинга объема зерна в силосе реализован в соответствии с фиг. 11.
Исходные данные хранятся в текстовых файлах в виде значений точек в трехмерной системе координат. В том случае, когда исходное облако точек является нерегулярным, осуществляется их упорядочивание в соответствии с заданной координатной сеткой с помощью существующих методов триангуляции. В поставленной задаче точность расчетов играет немаловажную роль, поэтому вычисление объема поверхности осуществляется самым ресурсозатратным способом, путем суммирования объемов косоусеченных призм. Для уменьшения погрешности вычисления в случае разреженности облака исходных точек могут использоваться методы интерполяции.
- 3 035343
Метод расчета объема прямой косоусеченной призмы - элементарного элемента, составляющего поверхность
Рассмотрим механизм нахождения точного объема прямой косоусеченной призмы с треугольным основанием, заданной тремя точками в трехмерном пространстве. Пусть косоусеченная призма задана точками A(ai, а2, а3), B(b1, b2, b3) и C(ci, c2, с3) (см. фиг. 12а). Точки A'(a1, а2, 0), B'(b1, b2, 0) и C'(c1, c2, 0) являются проекциями точек А, В и С соответственно на плоскость XOY.
Рассматриваемая фигура может быть представлена как совокупность трех фигур: прямая призма с треугольным основанием и две пирамиды.
Описание метода расчета объема прямой косоусеченной призмы
Шаг 1. Сначала упорядочиваем точки таким образом, чтобы точкой А обозначалась точка с наименьшим значением координаты по оси OZ, а С обозначало точку с наибольшей координатой по оси OZ.
Шаг 2. Разбиваем исходную фигуру на прямую призму (АВСА'В'С') и две пирамиды (ВАВС и САВС) в соответствии с фиг. 12b.
Шаг 3. Находим объем прямой призмы АВСА'В'С'. Объем представленной прямой призмы равен площади основания, умноженной на высоту. Высота рассматриваемой прямой призмы равна а3. Вычисление площади основания для общего случая наиболее удобно представить с помощью векторов: 3α.β.£.=1|ΑΈ·χΑΌΊ (1)
Определяем координаты векторов:
Ж -*pb2 - а2,0)
А С (cj aL2 а2;0}
Векторное произведение двух векторов:
a = {ax;ay;az} b = {bx;by;bz} в декартовой системе координат - это вектор, значение которого можно вычислить, используя следующие формулу (2):
i j k
ax ay az
bx by bz
- i(aybz - azby) - j(axbz - azbx) + k(axby - aybx) a x b — {aybz azby; azbx - axbz; axby - aybx}
Таким образом, по формуле (2) получаем:
(2)
ΑΈ'χΑ'σ= (0,0; (Ц - a,)(c2 - a2) - (b2 - a2)(C1 - аД)
Соответственно площадь треугольника и объем прямой призмы составят:
°А’В’С’
-аДс2 -a2)-(b2 -a2)(Cj -аД)2 I
-•а.
Шаг 4. Нахождение объема пирамиды ВАВС реализуем также с помощью векторов: ВЁА (0,0, а5-ЬД ВА{ ai-bi; a2-b2; a3-b5}? ВС{ Ci-bi; c2-b2; а33} Объем пирамиды по трем векторам определяется формулой (5):
νΕΛΕ.,.. Дввхвахвс|
О
Таким образом:
Vbab-c = 710 + 0 + (а33)(а! -^)(ο2 -b2)-(a3 -ЬДа2 -ЬД^ -bJ-O-OE = 7 I (а 3 -b3)[(a! - b!)(c2 -b2)-(a2 -b2)(Cj -ЬД I
-аД^ -a2)-(b2 -a2)(Cj -aj)2 I (5) (3) (4)
Шаг 5. Аналогично по формуле (5) вычислим объем пирамиды CABC с помощью векторов:
- 4 035343
СС”{0;0; а33), СА{ ai-сь а22; а33), СВ{ bi-ci; Ъ22; (Ъ33)-(с33))
V('aik'’ - — 10 + 0 + (а3С3)(аi ’ ci)(h2 ~ с2) — (аз - с3 )(bj - сj )(а2 - с2) — 0 — 01— = ||(а33)[(а! -ЬДЬ, -с2) - (Ь1 -сДа, -с2)]1
Шаг 6. Объединяем результаты, полученные на шагах 3-4. Суммируем объемы прямой призмы и двух пирамид и получаем общий объем заданной прямой косоусеченной призмы:
v = l, + -1(а3 о
+ -1(а3 о
-аДр -a2)-(b22)(с1 -аД)2 1 +
3)[(а1 -ЬДр -Ь2)-(а22)(с1 -ЬД]1 +
3)[(а1 -bj(b2 -cJ-Cy -сД(а22)]1 (6)
Полученная формула несколько тяжеловесна, однако дает точное значение объема рассматриваемой фигуры.
Алгоритм вычисления объема по результатам 3D измерения поверхности
Объем всей зерновой массы в процессе мониторинга можно получить, обходя последовательно все исходные точки поверхности и суммируя объемы прямых косоусеченных призм, прилегающих к точке. В результате мы получим объем пустот зернохранилища, ограниченных поверхностью зерна. Алгоритм вычисления объема сырья в зернохранилище по точкам поверхности реализован следующим образом.
Шаг 1. Получение облака точек с поверхности размером NxM. Обозначим количество точек как Num.
Шаг 2. Поиск min_length - минимального расстояния между точками в заданном облаке точек поверхности [14].
Шаг 3. Создание динамического двумерного массива, который является виртуальной проекцией изучаемой поверхности. Количество элементов в столбце таблицы определяется как целая часть от деления длины поверхности на минимальное расстояние между точками - int(N/min_length). Количество элементов в строке таблицы определяется как целая часть от деления ширины поверхности на минимальное расстояние между точками - int(M/min_length). Таким образом, мы эмулируем регулярную сетку с размером ячейки, соответствующей минимальному расстоянию между заданными точками.
Шаг 4. Помещаем исходные точки в ячейки, полученные на шаге 3 регулярной сетки. Благодаря тому, что размер ячеек соответствует минимальному расстоянию между заданными точками, каждой точке единственным образом соответствует элемент массива.
Шаг 5. Незаполненные на шаге 4 элементы массива заполняем с помощью линейной интерполяции.
Шаг 6. Обходим последовательно все элементы массива, за исключением элементов последнего столбца и последней строки. Для каждого элемента массива A[i][j] (где i=[0; int(N/min_length)) и j=[0; int(M/min_length))) по формуле 6 рассчитываем объем прямой косоусеченной призмы, определенный по точкам
а также объем прямой косоусеченной призмы, определений по точкам
В результате получаем объем пустот, ограниченных поверхностью.
Шаг 7. Объем зерновой массы определяем как разницу объемов зернохранилища и объема пустот, рассчитанного на шаге 6.
Таким, образом, в результате работы алгоритма получаем объем сырья в зернохранилище. Описанный алгоритм реализован в экспериментальном приложении, разработанном в среде программирования Visual Studio.NET 2010.
В результате проведенных исследований был получен алгоритм расчета объема зерновой массы в зернохранилище по точкам поверхности, который не зависит от регулярности полученных точек. В процессе реализации алгоритм был оптимизирован следующим образом.
Для исключения избыточных вычислений в случае min_length->0, ее значение устанавливается в 1. В этом случае в ячейке виртуальной регулярной сетки сохраняется только одно из значений облака, соответствующих координатам ячейки. Точность алгоритма при этом не страдает, так как разница между значениями в пределах 1см2 незначительна.
В алгоритме предусмотрена возможность задать размер ячейки сетки, соответствующий необходимой точности. Таким образом, исключается шаг 2 рассмотренного алгоритма, и сложность в целом уменьшается на O(Num).
Получение промежуточных значений можно выполнить с помощью линейной интерполяции, ин- 5 035343 терполяции полиномами Лагранжа или интерполяцией сплайнами на выбор.
После разработки приложения на базе описанного алгоритма были проведены эксперименты для расчета объема по данным пятнадцати поверхностей. Анализ результатов экспериментов показал, что отклонение значения, полученного в результате работы алгоритма, от фактического значения составляет в среднем 0,07%.
Вычисление объема и массы зерна в зернохранилище выполняется достаточно просто по остаточному принципу. На 2 этапе мы получили объем пустот в зернохранилище, таким образом, объем зерна составит разница между общим объемом зернохранилища и объемом образовавшихся пустот. Масса зерна определяется как отношение общего объема зерна к насыпной плотности соответствующего вида зерна.
По полученным в результате проведенного эксперимента данным можно сделать вывод о том, что точность вычислений тем выше, чем больше точек задействовано в вычислениях.
В целом проведенный эксперимент показал, что с помощью данного алгоритма можно рассчитать объем зерна, хранимого на зернохранилище.
Эффективное функционирование современных хранилищ различного типа предполагает решение ряда задач для обеспечения надежной сохранности сырья и оптимального использования складских площадей. Таким образом, актуальной является задача мониторинга хранящегося сырья. Некоторые трудности мониторинга возникают при хранении сыпучих габаритных объектов с неровной поверхностью.
Внедрение данного способа позволит осуществить непрерывное наблюдение за процессом работы автоматизированного зернохранилища и осуществлять контроль объема зерна, таким образом можно будет обеспечивать сохранность зерна от хищения и неправильного хранения, а результаты измерений будут способствовать формированию рекомендаций по управлению технологическими процессами.
Результаты исследования найдут применение при создании беспроводных программно-технических комплексов для других приложений, где необходимо обеспечение непрерывного контроля сложных поверхностей.

Claims (3)

1. Способ определения объема зерна в зернохранилище на основе бесконтактных 3П-измерений, включающий получение облака точек, обработку полученных данных, вычисление объема зерновой массы, отличающийся тем, что для получения облака точек используют видеокамеры, прикрепленные к потолку и направленные точно вниз, и комплекты лазерных диодов, закрепленных на фиксированном расстоянии от соответствующей камеры и направленных параллельно оптической оси камеры, причем предварительно определяют количество и расположение камер и лазерных диодов для получения облака точек необходимой плотности в зернохранилище так, чтобы точки их лучей формировали на поверхности зерна узлы регулярной сетки.
2. Способ по п.1, отличающийся тем, что включение лазерных диодов и получение изображения с каждой камеры выполняют автоматически и последовательно, исключая необходимость дополнительного упорядочивания значений и экономя вычислительные ресурсы.
3. Способ по п.1, отличающийся тем, что для расчета объема 3D точек используется суммирование косоусеченных призм с треугольным основанием.
EA201700618A 2017-07-12 2017-11-13 Способ определения объема зерна в зернохранилищах на основе бесконтактных 3d-измерений EA035343B1 (ru)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KZ20170587 2017-07-12

Publications (2)

Publication Number Publication Date
EA201700618A1 EA201700618A1 (ru) 2019-01-31
EA035343B1 true EA035343B1 (ru) 2020-05-29

Family

ID=65137784

Family Applications (1)

Application Number Title Priority Date Filing Date
EA201700618A EA035343B1 (ru) 2017-07-12 2017-11-13 Способ определения объема зерна в зернохранилищах на основе бесконтактных 3d-измерений

Country Status (1)

Country Link
EA (1) EA035343B1 (ru)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111815704B (zh) * 2020-07-16 2024-05-03 中国农业机械化科学研究院 一种基于双目相机的粮仓谷堆体积测算方法
CN115136799A (zh) * 2022-08-15 2022-10-04 安徽荣夏智能科技有限责任公司 一种粮食存储智能管理***
CN116358500B (zh) * 2023-06-01 2023-08-18 中粮资本科技有限责任公司 一种基于图像分析检测粮堆形变的方法及***
CN118196336B (zh) * 2024-05-17 2024-07-26 中山大学 随机颗粒堆积体的构建方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130144565A1 (en) * 2011-12-01 2013-06-06 Harris Corporation Accuracy-based significant point derivation from dense 3d point clouds for terrain modeling
CN103307980A (zh) * 2013-06-04 2013-09-18 中国科学院遥感与数字地球研究所 粮堆体积自动测量装置及其测量方法
CN105931238A (zh) * 2016-04-19 2016-09-07 北京理工大学 一种粮仓储粮体积测量的方法和***
CN106097318A (zh) * 2016-06-06 2016-11-09 北京理工大学 一种粮食体积测量***和方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130144565A1 (en) * 2011-12-01 2013-06-06 Harris Corporation Accuracy-based significant point derivation from dense 3d point clouds for terrain modeling
CN103307980A (zh) * 2013-06-04 2013-09-18 中国科学院遥感与数字地球研究所 粮堆体积自动测量装置及其测量方法
CN105931238A (zh) * 2016-04-19 2016-09-07 北京理工大学 一种粮仓储粮体积测量的方法和***
CN106097318A (zh) * 2016-06-06 2016-11-09 北京理工大学 一种粮食体积测量***和方法

Also Published As

Publication number Publication date
EA201700618A1 (ru) 2019-01-31

Similar Documents

Publication Publication Date Title
EA035343B1 (ru) Способ определения объема зерна в зернохранилищах на основе бесконтактных 3d-измерений
Fawzy 3D laser scanning and close-range photogrammetry for buildings documentation: A hybrid technique towards a better accuracy
CN102042814B (zh) 大型堆场三维形貌的投影辅助摄像测量方法
KR102559661B1 (ko) 부피 측량 방법, 시스템, 설비 및 컴퓨터 판독이 가능한 저장매체
US10582188B2 (en) System and method for adjusting a baseline of an imaging system with microlens array
Rupnik et al. 3D reconstruction from multi-view VHR-satellite images in MicMac
US6922234B2 (en) Method and apparatus for generating structural data from laser reflectance images
CN110230998B (zh) 基于线激光和双目相机的快速精密三维测量方法和装置
US20120114181A1 (en) Vehicle pose estimation and load profiling
ES2885230T3 (es) Procedimiento e instalación de control dimensional en línea de objetos fabricados
US10930062B2 (en) 3D view model generation of an object utilizing geometrically diverse image clusters
US20190065824A1 (en) Spatial data analysis
US10002437B2 (en) Method and electronic device of identifying redundant data
Charbonnier et al. Accuracy assessment of a canal-tunnel 3D model by comparing photogrammetry and laserscanning recording techniques
Dolloff et al. An assessment of WorldView-1 positional accuracy based on fifty contiguous stereo pairs of imagery
CN113989199A (zh) 一种基于深度学习的双目窄对接焊缝检测方法
CN116993942A (zh) 锥形物体三维重建及体积测量方法、***及无人机
CN101196396A (zh) 基于物方投影几何约束的线阵推扫式影像最佳扫描线搜索方法
Hosseininaveh et al. A low-cost and portable system for 3D reconstruction of texture-less objects
Peng et al. View selection with geometric uncertainty modeling
CN113160416A (zh) 一种用于煤流检测的散斑成像装置及方法
GB2569609A (en) Method and device for digital 3D reconstruction
Muratov et al. Estimation of distance to objects by stereovision
Putra et al. Approximate volume of sand materials stockpile based on structure from motion (sfm)
Singh et al. Digital photogrammetry for automatic close range measurement of textureless and featureless objects

Legal Events

Date Code Title Description
MM4A Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s)

Designated state(s): AM AZ KG TJ TM