Изобретения в сфере сельского хозяйства, животноводства, рыболовства

 
Изобретения в сельском хозяйстве Обработка почвы в сельском и лесном хозяйствах Посадка, посев, удобрение Уборка урожая, жатва Обработка и хранение продуктов полеводства и садоводства Садоводство, разведение овощей, цветов, риса, фруктов, винограда, лесное хозяйство Новые виды растений или способы их выращивания Производство молочных продуктов Животноводство, разведение и содержание птицы, рыбы, насекомых, рыбоводство, рыболовство Поимка, отлов или отпугивание животных Консервирование туш животных, или растений или их частей Биоцидная, репеллентная, аттрактантная или регулирующая рост растений активность химических соединений или препаратов Хлебопекарные печи, машины и прочее оборудование для хлебопечения Машины или оборудование для приготовления или обработки теста Обработка муки или теста для выпечки, способы выпечки, мучные изделия

Способ вычисления запаса лесных массивов

 
Международная патентная классификация:       A01G G03B

Патент на изобретение №:      2242867

Автор:      Давыдов В.Ф. (RU), Корольков А.В. (RU), Гренц Н.В. (RU), Шалаев В.С. (RU)

Патентообладатель:      Московский государственный университет леса (RU)

Дата публикации:      27 Декабря, 2004

Начало действия патента:      10 Октября, 2003

Адрес для переписки:      141005, Московская обл., г. Мытищи-5, МГУЛ, патентный отдел


Изображения





Изобретение относится к лесному хозяйству, в частности к дистанционным методам таксации лесов на обширных площадях. Способ основан на получении изображения лесов в виде цифровой матрицы зависимости яркости I(х, у) от пространственных координат. Далее рассчитывают пространственный спектр Фурье матрицы. Находят среднюю частоту Fср, диаметр кроны среднего дерева Дср=I/F ср. По массовым таблицам находят высоту и ступень толщины среднего дерева. Матрицу последовательно разбивают на фрагменты /2×2/ смежных пикселей и делят четырехугольник по главной диагонали на два треугольника. Площадь рельефа лесного полога Sрп представляют в виде мозаики аппроксимирующих треугольников. По формуле Герона находят площадь каждого треугольника и их сумму. Вычисляют полноту насаждения Р по заданной формуле. Вычислением среднего количества деревьев на участке Nср по заданной формуле определяют прикрепляющую точку огивы анализируемого насаждения. Запас вычисляют по массовым таблицам как сумму 5 классов ступеней толщины Лорея по заданной формуле. Способ позволяет повысить точность и статистическую устойчивость результата вычисления запаса лесного массива. 7 ил.

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

За двести лет существования науки о лесах накоплены обширные статистические материалы о закономерностях соотношений между морфометрическими характеристиками как отдельных деревьев, так и их сообществ: высоты дерева, диаметре ствола, диаметре кроны, возрасте, полноте насаждения, бонитете, запасе и др. Выявленные статистические зависимости между характеристиками лесных ценозов зафиксированы в так называемых массовых таблицах (см., например, "Объемные таблицы Союзлеспрома" в учебнике Анучин Н.П. Лесная таксация, 5-е издание. - М.: Лесная промышленность, 1982 г., стр. 167, 207, 250, 273, 296).

Известен способ натурной таксации путем индивидуального пересчета деревьев на пробных площадках (см., например, Анучин Н.П. Лесная таксация, 5-е издание. - М.: Лесная промышленность, 1962 г., стр. 248-250 - аналог).

В способе-аналоге разбивают площадь лесного массива на мозаику участков, проводят измерения таксационных параметров на пробном (ключевом) участке: Нi - высоты дерева, di - толщины ствола на высоте 1,3 м, Ni - числа деревьев данной высоты и толщины, полноты насаждения Рo, fi - видовых чисел, а запас насаждения участка рассчитывают как сумму произведений

распространяют результаты расчета на всю площадь лесного массива.

Недостатком известного аналога является большая трудоемкость, связанная с необходимостью учета и обмера каждого дерева, неоперативность работ, недоступность горных и отдаленных районов, погрешность ( 20%) при распространении результатов измерений отдельных участков на весь таксируемый массив. Существующие космические средства позволяют оперативно получать изображения ранее не учтенных лесных массивов. Имеется ряд запатентованных способов (патенты РФ №2080051, №2130707, №2133565 - аналоги), позволяющие проводить дистанционную оценку запаса насаждений по их изображениям.

Ближайшим аналогом к заявляемому решению является "Способ таксации насаждений", патент РФ №2183847, 2002 г.) В способе ближайшего аналога получают изображение лесного массива в виде цифровой матрицы |m× n| зависимости яркости I(х, у) от пространственных координат, вычисляют пространственный спектр матрицы и его математическое ожидание (Fcp), осуществляют высокочастотную фильтрацию пространственного спектра с частотой среза, равной Fcp , проводят обратное Фурье-преобразование отфильтрованного спектра, рассчитывают таксационные параметры: полноту Р как отношение количества светлых пикселей nc к их общему количеству |m× n|, находят прикрепляющую точку огивы насаждения расчетом среднего количества деревьев в насаждении: где Дср=I/Fcp - средний диаметр крон древесного полога, So – площадь анализируемого участка, Н cp=7Д1,2cp - высота среднего дерева, а запас вычисляют по массовым таблицам, связывающим видовые числа fi, высоту Нi, площадь сечения ствола g i= d2icp/4 и количество деревьев в насаждении Ni данной ступени толщины из соотношения:

Недостатками ближайшего аналога являются:

- статистическая неустойчивость результата, основанная на использовании одного независимого признака изображения - пространственного спектра или производных от него; Fcp, nc;

- неточность расчета полноты Р насаждения через производные от пространственного спектра.

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

Решение поставленной задачи обеспечивается тем, что в способе вычисления запаса лесных мостов, содержащих пробные площадки, получает изображение лесов в виде цифровой матрицы зависимости яркости I(х, у) от пространственных координат, рассчитывают пространственный спектр Фурье матрицы, находят среднюю частоту Fср, диаметр кроны среднего дерева Дср=I/Fcp, а по массовым таблицам высоту и ступень толщины среднего дерева, последовательно разбивают матрицу на фрагменты |2× 2| смежных пикселей и делят четырехугольник по главной диагонали (слева сверху - направо вниз) на два треугольника, представляют площадь рельефа лесного полога Sрп в виде мозаики аппроксимирующих треугольников, находят по формуле Герона площади каждого треугольника и их сумму, вычисляют полноту насаждения Р через отношение площади рельефа пробной площадки Sпрп к удельной площади рельефа лесного полога Sпр/So домноженного на полноту насаждения пробной площадки Ро, вычислением среднего количества деревьев на участке Nсp=4· So· P/ Д2cр определяет прикрепляющую точку огивы анализируемого насаждения, вычисляют запас по массовым таблицам как сумму 5 классов ступеней толщины Лорея

где So - площадь изображения обрабатываемого лесного массива, га;

Ni - количестве деревьев i-го класса ступени толщины;

Нi - высота среднего дерева i-го класса ступени толщины, м;

gi - площадь сечения среднего дерева i-го класса ступени толщины, м2;

fi - видовое число среднего дерева i-го класса ступени толщины.

Изобретение поясняется чертежами, где

фиг.1 - исходное изображение анализируемого лесного массива;

фиг.2 - визуализированный пространственный спектр Фурье сигнала матрицы изображения;

фиг.3 - огибающая пространственного спектра Фурье и его средняя частота Fcp;

фиг.4 - последовательность разбиения матрицы на фрагменты |2× 2| из 4-х смежных пикселей и мозаику треугольников,

фиг.5 - зависимость между полнотой насаждения и рельефом (сечение древесного полога плоскостью X),

фиг.6 - огива насаждения для расчета запаса;

фиг.7 - функциональная схема устройства, реализующего способ.

Техническая сущность изобретения состоит в следующем. По определению, полнота насаждения характеризует совершенство использования растительным сообществом занимаемого пространства и оценивается сомкнутостью древесного полога, густотой. Даже в наземных условиях расчет полноты насаждения представляет определенные трудности (см., например, "Полнота насаждения" в учебнике Н.П.Анучин "Лесная таксация", 5-е издание. - М.: Лесная промышленность, 1982 г., стр. 234-241). Чаще всего расчет полноты проводят путем сравнения таксируемого участка с нормальным полным насаждением. Очевидно, что характеристикой полноты может быть объем пространства, занимаемый насаждением. На изображении вершина кроны дерева обладает большей яркостью, так как отражает падающий световой поток почти зеркально. Промежутки между деревьями из-за многократного диффузного отражения светового потока от соседних крон и лишь частичного отражения в сторону регистратора, имеют на изображении наименьшую яркость. Таким образом, на изображении древесный полог представляется в виде совокупности колоколообразных пиков разной высоты. Глубина перепада яркости отражает степень изрезанности рельефа древесного полога и связана с полнотой насаждения. Очевидно, что чем больше полнота, тем больше сомкнутость полога и тем меньше совокупная площадь "рельефа" функции яркости изображения I(х, у). Совокупную площадь рельефа вычисляют, аппроксимируя ее поверхность мозаикой треугольников. Чем выше разрешение снимка, тем точнее аппроксимирующая поверхность совпадает с истинной площадью рельефа древесного полога. На фиг.5 иллюстрируется зависимость полноты насаждения от площади рельефа его полога. Количественно полноту насаждения находят через отношение площади рельефа пробной площадки Sпpп к удельной площади рельефа лесного полога всего снимка, домноженного на полноту насаждения пробной площадки Ро

Р=Ро· Sпрп /Sрп/Sо.

Площадь рельефа древесного полога исходной матрицы изображения лесного массива вычисляют программным методом на ПЭВМ. Последовательность разбиения матрицы на фрагменты из 4-х смежных пикселей и мозаику треугольников по главной диагонали (1-4) иллюстрируется фиг.4.

Для вычисления площади поверхности обрабатываемого участка необходимо указать масштаб по координатам dх, dу и масштаб глубины dI. Алгоритм вычисления элементарной площади четырехточечного шаблона основан на триангуляции элемента поверхности, т.е. разбиение его по главной диагонали (слева-сверху-направо-вниз) на два треугольника. Длину сторон каждого из треугольников при известных dх, dу, dI вычисляют по теореме Пифагора, а площадь треугольников по формуле Герона

где a, b, c - длины сторон треугольника, р - полупериметр треугольника.

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

Пример реализации способа

Заявляемый способ может быть реализован на базе устройства по схеме фиг.7. Функциональная схема устройства фиг.7 содержит орбитальный комплекс наблюдения 1 типа космического аппарата "Ресурс-Ф" с установленной на его борту фотокамерой 2 типа КФА-1000, осуществляющей съемку участков подстилающей поверхности 3 по командам из Центра управления полетом 4, передаваемым по радиолинии управления 5. Заполненные кассеты с фотоматериалами съемки отстреливаются с борта космического аппарата по командам из ЦУПа 4 и в виде капсул 6 доставляются в лабораторию обработки фотоматериалов 7. При обработке фотоматериалов осуществляют их геометрическую фотометрическую коррекцию. Материалы космической съемки в виде откорректированной серии снимков 8 поступают заказчику в Центр тематической обработки 9. Посредством сканера 10 высокого разрешения типа "Panasonic" аналоговые снимки преобразуют в цифровые матрицы дискретных отсчетов и вводят в ПЭВМ 11, включающую стандартный набор элементов: процессор 12, оперативное запоминающее устройство 13, винчестер 14, дисплей 15, принтер 16, клавиатуру 17. В качестве ПЭВМ 11 использовалась станция типа IBM со специальным программным обеспечением MATH САД.

Масштаб снимка (Лосиный Остров, г.Москва) составлял 1:50000, после сканирования с разрешением 2000 точек на дюйм, разрешение одного пиксела матрицы изображения соответствовало 0,65 м. Исходное изображение для обработки иллюстрируется фиг.1.

Вычисление пространственного спектра матрицы является стандартной процедурой специализированного программного обеспечения MATH САД. Визуализированное изображение пространственного спектра (распечатка с ПЭВМ) иллюстрируется фиг.2.

На фиг.3 представлена нормированная А/Аmax функция огибающей пространственного спектра исходного изображения. Математическое ожидание пространственной частоты исходного спектра Fср=0,33; средний диаметр кроны Дcp=I/Fcp=3 м, а высота среднего дерева Нcp=7Д1,2ср=24 м. Расчетный диаметр ствола среднего дерева d1,3=42 см (см., например, Общесоюзные нормативы для таксации лесов. - М.: Колос, 1992 г., стр. 121). Таксационные параметры пробной площадки анализируемого лесного массива (по данным натурной таксации) составили d 1,3=44 см, Н=25 м, Ро=0,6, V=260 м3 /га, видовые числа по классам f1=0,575; f2 =0,57; f3=0,55; f4=0,53; f5=0,52. Вся площадь обрабатываемого лесного массива (стандартная матрица |512× 512| элементов, при разрешении на пиксел 0,65 м) составила S0=11 га. Площадь рельефа лесного полога, вычисленная программным методом Sрп=158 га, удельная площадь рельефа 158/11=14,4 га. Площадь рельефа пробной площадки, вычисленная программным методом 7,9 га. Расчетная полнота насаждения лесного массива

Р=0,6× 7,9/14,4=0,32.

Ошибка расчета полноты насаждения аналога при распространении полноты пробной площадки на весь таксируемый массив достигает 40%. Лесной наукой выявлена закономерность, что в любом насаждении, независимо от возраста, состава пород, вида насаждения - однородное или смешанное, распределение количества деревьев по ступеням толщины подчиняется одному и тому же закону, носящему название огивы насаждения. График огивы насаждения представлен на фиг.6. При полученных расчетных параметрах определяют среднее число деревьев в анализируемой массиве

где Р - расчетная полнота насаждения лесного массива (Р=0,32);

Д2ср/4 - проекция кроны среднего дерева на земную поверхности (Дср=3 м).

Разбивая все количество деревьев на 5 классов ступеней толщины Лорея, вычислено, что расчетный запас всего лесного массива равен 218 м3 /га.

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

Формула изобретения

Способ вычисления запаса лесных массивов, содержащих пробные площадки, при котором получают изображения лесов в виде цифровой матрицы зависимости яркости I(х, у) от пространственных координат, рассчитывают пространственный спектр Фурье матрицы, находят среднюю частоту Fср, диаметр кроны среднего дерева Дcр =I/Fcр, а по массовым таблицам - высоту и ступень толщины среднего дерева, отличающийся тем, что последовательно разбивают матрицу на фрагменты /2× 2/ смежных пикселей и делят четырехугольник по главной диагонали на два треугольника, представляют площадь рельефа лесного полога Sрп в виде мозаики аппроксимирующих треугольников, находят по формуле Герона площадь каждого треугольника и их сумму, вычисляют полноту насаждения Р через отношение площади рельефа пробной площадки Sпрп к удельной площади рельефа лесного полога S рп/Sо, домноженного на полноту насаждения пробной площадки Ро, вычислением среднего количества деревьев на участке Ncр=4· Sо· Р/ Д2ср определяют прикрепляющую точку огивы анализируемого насаждения, вычисляют запас по массовым таблицам как сумму 5 классов ступеней толщины Лорея

где Sо - площадь изображения обрабатываемого лесного массива, га;

Ni - количество деревьев i-го класса ступени толщины;

Нi - высота среднего дерева i-го класса ступени толщины, м;

gi - площадь сечения среднего дерева i-го класса ступени толщины, м;

fi - видовое число среднего дерева i-го класса ступени толщины.

MM4A - Досрочное прекращение действия патента СССР или патента Российской Федерации на изобретение из-за неуплаты в установленный срок пошлины за поддержание патента в силе

Дата прекращения действия патента: 11.10.2005

Извещение опубликовано: 20.06.2007        БИ: 17/2007





Популярные патенты:

2161391 Комбинированная почвообрабатывающая посевная машина

... , , , , , ОПИСАНИЕ ИЗОБРЕТЕНИЯ К ПАТЕНТУ Изобретение относится к сельскохозяйственной технике, в частности к машинам для обработки почвы одновременно с посевом за один проход агрегата. Известна комбинированная машина - устройство для обработки почвы и внесения удобрений (SU а.с. N 1512499, МКИ A 01 B 13/04 от 07.07.87), включающее рамы рыхлящей части и прикатывающей, соединенные серьгой, дисковые батареи, рабочие органы для внесения удобрений в почву, трубопроводы, опорные колеса, прицепное устройство. К недостаткам известной машины следует отнести то, что она неудовлетворительно адаптируется на смену почвенных условий и не производит одновременно с обработкой почвы ...


2073513 Способ профилактики технологических стрессов молодняка крупного рогатого скота

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


2012206 Инсектицидная композиция для борьбы с тараканами

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


2233582 Устройство для охлаждения молока

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


2168887 Машина для добычи корней

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


Еще из этого раздела:

2402189 Роликовая сортировальная машина

2121787 Устройство для регулирования температуры воздуха в теплице

2488263 Система механической подачи недомолота для вторичного обмолота на возвратную доску

2090040 Машина для возделывания корнеклубневых культур

2121263 Способ лесоводственной оценки технологического комплекса машин

2010519 Способ биологической борьбы с вредителями растений

2267924 Способ стимулирования роста растений

2189718 Пневматический высевающий аппарат

2307495 Пневматический высевающий аппарат

2470922 Сокристаллы