Способ таксации насажденийПатент на изобретение №: 2183847 Автор: Давыдов В.Ф., Илларионов Г.П., Шалаев В.С., Комаров Е.Г., Мухин А.С. Патентообладатель: Московский государственный университет леса Дата публикации: 20 Июня, 2002 Начало действия патента: 28 Июня, 2000 Адрес для переписки: 141005, Московская обл., Мытищи-5, МГУЛ, патентный отдел Изображения    Изобретение относится к лесному хозяйству, методам дистанционного решения задач лесохозяйственного назначения. При реализации заявленного способа получают изображение лесного массива, преобразуют функцию яркости изображения I(x,y) в цифровую матрицу дискретных отсчетов размерностью |m  n| элементов, рассчитывают числовые характеристики матрицы и огибающую пространственного спектра, вычисляют таксационные параметры по статистическим зависимостям между числовыми характеристиками матрицы и элементами насаждений. Способ позволяет вычислить несколько независимых элементов насаждений по нескольким независимым признакам изображения. Технический результат - повышение точности и устойчивости результатов оценок расчетных параметров. 7 ил., 1 табл.
, , , , , , ,
ОПИСАНИЕ ИЗОБРЕТЕНИЯ К ПАТЕНТУИзобретение относится к лесному хозяйству, в частности к методам дистанционного решения лесохозяйственных задач на основе дешифрирования изображений лесов, получаемых в различных диапазонах электромагнитных волн. За период становления лесной науки разработано множество методов натурной таксации насаждений путем непосредственных измерений параметров деревьев на пробных площадках и последующего аналитического расчета таксационных характеристик по статистическим зависимостям (см. , например, Н.П. Анучин "Лесная таксация", учебник, 5-е изд., М., "Лесная промышленность", 1982 г., стр. 250, стр. 344-348 - аналоги). При всех методах натурной таксации вычисляют такие таксационные характеристики, как высоту насаждений Н, диаметр ствола деревьев d1,3 полноту насаждений Р, запас древесины на 1 га, м3, бонитет, возраст и другие характеристики. Недостатками методов натурной таксации являются: - неоперативность и большая трудоемкость; - ограниченность статистической выборки измерений и, как следствие, ошибки оценок при распространении результатов измерений пробных площадок на весь таксируемый массив. Дистанционные методы таксации на основе получения изображения лесов позволяют сформировать репрезентативную выборку измерений и рассчитать числовые характеристики матрицы изображения сразу всего лесного массива с точностью разрешения одного пикселя. Массив пикселей содержит информацию как о параметрах отдельных деревьев, так и о статистических характеристиках элементов насаждения в целом. Последнее обстоятельство обеспечивает повышение точности оценок. Известен способ оценки запаса насаждений путем получения изображения лесного массива, преобразования его в дискретную матрицу отсчетов яркости, нахождения всех локальных максимумов матрицы, расчета среднего расстояния между максимумами, выделения контуров вокруг каждого максимума методами пространственного дифференцирования и оценки запаса по аналитическим зависимостям между рассчитанными числовыми характеристиками матрицы и элементами дерева среднего класса Лорея (см., например, патент РФ 2133565, кл. А 01 G 23/00, 23/02, 1999 г. - аналог). Способ-аналог может быть использован только при наличии высокодетальных снимков с разрешением менее 1 м, когда крона одиночного дерева отображается на снимке несколькими десятками пикселов. При меньшем числе пикселов (отсчетов) методы пространственного дифференцирования не работоспособны. Получить снимки с разрешением менее 1 м затруднительно, как по режимным ограничениям, так и по техническим возможностям космических систем. Скрытые периодичности чередования крон древесного полога могут быть выявлены по снимкам среднего разрешения (более 1 м) путем вычисления спектра Фурье матрицы изображения. Ближайшим аналогом по технической сущности с заявляемым способом является "Способ определения запаса насаждений" (патент РФ 2080051, кл. А 01 G 23/00, 1997 г.). В способе ближайшего аналога получают изображение лесного массива, преобразуют функцию яркости изображения I(x, y) в матрицу |mn| дискретных отсчетов, вычисляют характеристики электрического сигнала матрицы: огибающую пространственного спектра, математическое ожидание, среднеквадратическое отклонение  , а запас насаждения рассчитывают по статистическим зависимостям между огибающей пространственного спектра, полнотой насаждения Р в функции от и параметрами отдельного дерева: высотой Н и сечением ствола g. Недостатками ближайшего аналога являются: - неточность расчета полноты насаждения Р через среднеквадратическое отклонение сигнала матрицы ; - статистическая неустойчивость результатов оценок, основанная на использовании фактически единственного независимого признака изображения - его текстуры и соответствующей ей огибающей пространственного спектра. Задача, решаемая данным изобретением, заключается в повышении точности и устойчивости результатов оценок расчетных параметров путем вычисления нескольких независимых элементов насаждения по нескольким независимым признакам изображения. Решение поставленной задачи обеспечивается тем, что в способе таксации насаждений, при котором получают изображение лесного массива, преобразуют функцию яркости изображения I(x,y) в цифровую матрицу дискретных отсчетов размерностью |mn| элементов, рассчитывают числовые характеристики матрицы и огибающую пространственного спектра, вычисляют таксационные параметры по статистическим зависимостям между числовыми характеристиками матрицы и элементами насаждений, рассчитывают математическое ожидание пространственного спектра (Fcp), осуществляют высокочастотную фильтрацию пространственного спектра с частотой среза, равной его математическому ожиданию, проводят обратное Фурье-преобразование отфильтрованного спектра, определяют полноту насаждения P как отношение количества выделенных в результате обратного Фурье-преобразования пикселей (nc) к их общему количеству в матрице, находят прикрепляющую точку огивы насаждения к его изображению расчетом среднего количества деревьев в насаждении  где Dcp=I/Fcp - средний размер крон древесного полога, S - площадь анализируемого участка насаждения, Н= 7 Dcp 1,2 - средняя высота деревьев, а запас насаждения вычисляют по массовым таблицам, связывающим видовые числа fi, высоту hi, ступени толщины (d1,3i), площади сечений  и количество деревьев в насаждении Ni данной ступени толщины из соотношения:  Изобретение поясняется чертежами, где: фиг.1 - исходное изображение анализируемого лесного массива; фиг. 2 - визуализированный пространственный спектр Фурье сигнала матрицы изображения на комплексной плоскости; фиг.3 - огибающая пространственного спектра Фурье: а) распечатка с дисплея в логарифмическом масштабе программного обеспечения ПЭВМ, б) реальная огибающая в масштабе исходного изображения; фиг. 4 - маска оптимального высокочастотного фильтра на фоне спектра Фурье; фиг. 5 - визуализированное изображение обратного Фурье-преобразования отфильтрованного спектра; фиг. 6 - огива закономерности распределения деревьев в насаждении по естественным ступеням толщины; фиг.7 - функциональная схема устройства, реализующего способ. Вновь введенные операции, образующие совокупность существенных признаков, обеспечивают достижение таких качественных свойств способа, как: - возможность обработки изображений с разрешением на уровне крон деревьев; - достоверность и статистическая устойчивость результатов за счет использования нескольких независимых признаков изображения; - адекватность оценок за счет точной привязки известных закономерностей в строении насаждений к числовым характеристикам матрицы изображения. Техническая сущность изобретения заключается в следующем. Многие таксационные показатели, определяемые при оценке насаждений, не получают своего непосредственного отображения на снимках. Скрытые закономерности в строении насаждения и отдельных деревьев могут быть выявлены на основе использования нескольких независимых признаков изображения, таких как цвет, тон, текстура, топология и др. Чем больше независимых признаков изображения использовано при расчете, тем выше точность метода. Известны статистические зависимости между параметрами деревьев и элементами насаждения. Так, например, полнота насаждений связана с сомкнутостью полога, а диаметр крон - с высотой древостоя. Также установлено (см. , например, Н.П. Анучин "Лесная таксация", учебник, 5-е изд., М., "Лесная промышленность", с. 271-273), что распределение количества деревьев в насаждении по естественным ступеням толщины не зависит ни от породы, ни от бонитета, ни от полноты насаждения. Общей для всех насаждений и статистически устойчивой характеристикой распределения количества деревьев по естественным ступеням толщины является так называемая огива насаждения. Таким образом, исходными независимыми элементами насаждения являются: количество деревьев в насаждении N и распределение размеров крон деревьев по диаметрам D, м. Остальные таксационные параметры могут быть вычислены по аналитическим зависимостям, которые отражают статистические закономерности в виде функций регрессии этих параметров от независимых N и D. Для расчета независимых элементов N и D в заявляемом способе используют два независимых признака изображения: тон и текстуру. Текстура изображения содержит информацию о степени изрезанности древесного полога, образуемой кронами деревьев. Крона отдельного дерева при съемке сверху представляется некоторой колоколообразной фигурой. Вершина кроны освещена ярче, чем промежутки между кронами деревьев. Поэтому, текстура изображения повторяет периодичности чередования крон деревьев. Скрытая информация о периодичности чередования и размерах крон деревьев может быть выявлена путем вычисления двумерного спектра Фурье от анализируемой матрицы изображения (см., например, "Анализ пространственных частот" в книге Р. Дуда, П. Харт "Распознавание образов и анализ сцен", перевод с англ., Мир, Москва, 1976 г. , стр. 319-321). Двумерный пространственный спектр L(Fx, Fy) от функции яркости изображения I(x, y) рассчитывают преобразованием Фурье:  Представленное аналитическое выражение является стандартной математической операцией, реализуемой алгоритмами быстрого Фурье-преобразования (БПФ), входящих в комплект специализированного программного обеспечения: "Пакет программ для обработки изображений в науках о Земле" (см., например, "ER MAPPER Reference", Earth Resource Mapping Pty Lid, Western Australia, 6005, 1998 г., р. 295). На фиг.1 представлено исходное изображение лесного массива, масштаб снимка М 1:10000. На фиг.2 представлен визуализированный двумерный пространственный спектр Фурье на комплексной плоскости после программной реализации операции двойного интегрирования I(x, y). На фиг.3 представлена огибающая пространственного спектра: а) распечатка с дисплея в логарифмическом масштабе программного обеспечения ER MAPPER, б) реальная функция в масштабе исходного изображения. Полученная функция по физическому смыслу представляет собой величину, обратную линейным размерам объектов подстилающей поверхности, отождествляемых с распределением диаметров крон деревьев D. Расчет огибающей реализуется стандартной программой, входящей в комплект программного обеспечения ER MAPPER 5,0 (см. там же, стр. 58, 286). Амплитуда А огибающей в каждой точке Fi определяет удельный вес соответствующей пространственной гармоники в общем спектре. Математическое ожидание (среднюю частоту) пространственного спектра находят из условия:  При этом расчетный средний диаметр крон анализируемого участка насаждения соответствует:  Другим независимым признаком изображения является тон. При преобразовании изображения в цифровую матрицу дискретных отсчетов функция яркости I(x, y) обычно квантуется в стандартной шкале от 0 до 256 уровней. Для расчета количества деревьев необходимо превратить исходное изображение в контурный рисунок, т.е. выделить контуры всех крон. Контур - это край участка изображения, где велик градиент изменения функции яркости I(x, y). Если разрешение снимка невелико и число пикселов недостаточно для расчета градиента, то эффективным методом выделения контурного рисунка является высокочастотная фильтрация (см., например, "Пространственная фильтрация" в книге Р. Дуда, П. Харт "Распознавание образов и анализ сцен", с. 329). Чтобы эффективно выделить контуры крон на анализируемом изображении, необходимо оптимальным образом рассчитать передаточную характеристику К(Fx,Fy) высокочастотного фильтра. Неидентифицируемые разрывы между кронами, протяженные прогалины, опушки, соответствуют низким частотам пространственного спектра изображения и расположены на комплексной плоскости в начале координат. Следовательно, для повышения точности оценок необходимо срезать низкие частоты (НЧ), и подчеркнуть высокие частоты (ВЧ). Границей раздела (частотой среза) между НЧ и ВЧ служит рассчитанная выше характеристика - математическое ожидание спектра Fcp. На фиг.4 представлена визуализированная маска высокочастотного фильтра с частотой среза Fcp. Высокочастотная фильтрация заключается в перемножении исходного пространственного спектра L(Fx,Fy) на передаточную функцию высокочастотного фильтра К(Fx, Fy), в результате чего получают отфильтрованный спектр L0(Fx,Fy): L0(Fx,Fy)=L(Fx,Fy)K(Fx,Fy). Высокочастотная цифровая фильтрация реализуется программным расчетом на ПЭВМ с использованием стандартной процедуры (см., там же ER MAPPER 5,0 стр. 299). После оптимальной фильтрации сигнала матрицы изображения осуществляют восстановление функции яркости изображения Iо(x, y) обратным Фурье-преобразованием: I0(x,y) = Ф-1[Lo(Fx,Fy)]. Данная операция также реализуется стандартной программной процедурой, см. ER MAPPER 5,0, стр. 299. Визуализированное изображение после операции обратного Фурье-преобразования представлено распечаткой с ПЭВМ на фиг.5. Как следует из фиг.5, визуально, операция ВЧ фильтрации состоит в ранжировании пикселей по яркости. Светлые участки изображения соответствуют площади, занимаемой кронами деревьев, темные участки - прогалинам. Отношение числа пикселей nс, принадлежащих к светлым участкам, к общему числу пикселей |mn| представляет собой характеристику насаждения - проекцию древесного полога на поверхность Земли или полноту (Р):  Располагая численным значением полноты, определяют среднее количество деревьев Ncp на анализируемом изображении через отношение площади участка S, домноженной на полноту P, к площади кроны среднего дерева:  Ncp связывает прикрепляющую точку огивы анализируемого насаждения с матрицей его изображения, для которой естественная ступень толщины, в долях среднего диаметра, равна единице (см. фиг.6). Количественное распределение деревьев (N) в насаждении по естественным ступеням толщины, задаваемое огивой, представлено таблицей . Зависимость высоты дерева Н от диаметра кроны D представляется степенной функцией вида (см. аналог, патент РФ 2133565, 1999 г.) H,м 7D1,2(м). Между высотой дерева Н и его толщиной d1,3 существует зависимость, представленная в массовых таблицах А.В. Тюрина, Д.И. Товстолеса, В.К. Захарова (см. , например, "Общесоюзные нормативы для таксации лесов", Справочник, М. , "Колос", 1992 г., стр. 67-73). Кроме того, существуют таблицы всеобщих видовых чисел (f) М.Е. Ткаченко (см., например, Н.П. Анучин "Лесная таксация", 5-е издание, М., "Лесная промышленность", 1982 г., стр. 148). Имея исходные данные массовых таблиц, запас насаждения рассчитывают по формуле  где fi - видовое число, gi - площади поперечных сечений, Ni, hi - количество и высота деревьев i-й ступени толщины на анализируемом изображении. Примеры реализации способа. Заявляемый способ может быть реализован на базе устройства по схеме фиг. 7. Функциональная схема устройства фиг.7 содержит орбитальный комплекс наблюдения 1 типа космического аппарата "Ресурс-Ф" с установленной на его борту фотокамерой 2 типа КФА-1000, осуществляющей съемку участков 3 подстилающей поверхности 4 по командам из Центра управления полетом 5, передаваемым по радиолинии управления 6. Заполненные кассеты с фотоматериалами съемки отстреливаются с борта космического аппарата по командам из ЦУПа и в виде капсул 7 доставляются в лабораторию обработки фотоматериалов 8. При обработке фотоматериалов осуществляется их предварительная геометрическая и радиометрическая коррекция. Материалы космической съемки в виде откорректированной серии снимков 9 поступают заказчику, в Центр тематической обработки 10. Посредством сканера 11 высокого разрешения типа "Panasonic" аналоговые снимки преобразуются в цифровые матрицы дискретных отсчетов и вводятся в ПЭВМ 12, включающей стандартный набор элементов в составе: процессора 13, оперативного запоминающего устройства 14, винчестера 15, клавиатуры 16, дисплея 17 и принтера 18. В качестве ПЭВМ 12 тематической обработки снимков 9 использовалась станция CVN со специальным программным обеспечением ER MAPPER 5.0, предварительно записанным на винчестер 15. Любое специализированное программное обеспечение в качестве единицы расстояния на плоскости изображения оперирует размерами одного элемента изображения - пикселя. Единицей пространственной частоты при этом считается величина, обратная элементу изображения  Следовательно, чтобы оперировать реальными линейными размерами объектов на изображении, необходимо связать масштаб снимка с пространственным разрешением одного пикселя. Масштаб исходного снимка соответствовал M 1:10000. Для ввода в ПЭВМ снимок был отсканирован с разрешением 400 точек на дюйм. Стандартная матрица изображения имела размерность 512512 отсчетов. Пространственное разрешение одного пикселя составило 0,64 м, а площадь анализируемого изображения S = |330330|, м  11 га. Программное обеспечение ER MAPPER 5,0 выдает функцию огибающей пространственного спектра в логарифмическом, по основанию 2, масштабе. На фиг. 3а иллюстрируется распечатка этой функции с экрана дисплея. На фиг.3б представлена огибающая пространственного спектра, пересчитанная из логарифмического масштаба в масштаб реального изображения. Математическое ожидание пространственного спектра исходного снимка составило Fср= 0,26, а средний диаметр крон древесного полога Dcp=I/Fcp=3,8 м. Средняя высота Н=7Dcp 1,2=34 м. После высокочастотной фильтрации с частотой Fcp=0,26 и обратного Фурье-преобразования, произошло ранжирование пикселей изображения по яркости. Восстановленное после ВЧ фильтрации изображение иллюстрировано на фиг.5. Численный расчет количества светлых пикселей (nc) к общему их числу в матрице дает значение полноты:  Количество деревьев в прикрепляющей точке огивы насаждения:  что составляет 66,3% от общего количества деревьев на 1 га. Расчетный диаметр ствола среднего дерева в прикрепляющей точке огивы насаждения равен d1,3=46 см. По массовым таблицам находят диаметры стволов и видовые числа деревьев насаждения всех ступеней толщины, а по огиве - их количество Ni в каждой ступени. Запас насаждения рассчитывают численными методами с использованием массовых таблиц. Расчетный запас насаждения для рассмотренного примера реализации составил V =  gihifiNi = 320 м3/га. Эффективность заявляемого способа оценивалась путем сравнения полученных таксационных характеристик насаждения с результатами наземной перечислительной таксации. Перечислительная таксация была выполнена на момент получения исходного изображения. Для анализируемого участка таксационные характеристики натурной таксации соответствовали: Р=0,4, d1,3=48 см, Н=33 м, V=340 м3/га, что не выходит за пределы допуска (см., например, "Точность определения таксационных показателей" в справочнике "Общесоюзные нормативы для таксации лесов", М., "Колос", 1992 г., с. 121). При внедрении в производство заявляемый способ измерений может рассматриваться как метрологический.
ФОРМУЛА ИЗОБРЕТЕНИЯСпособ таксации насаждений, при котором получают изображение лесного массива, преобразуют функцию яркости изображения 1(х,у) в цифровую матрицу дискретных отсчетов размерностью |m x n| элементов, рассчитывают числовые характеристики матрицы и огибающую пространственного спектра, вычисляют таксационные параметры по статистическим зависимостям между числовыми характеристиками матрицы и элементами насаждений, отличающийся тем, что рассчитывают математическое ожидание пространственного спектра (Fcp), осуществляют высокочастотную фильтрацию пространственного спектра с частотой среза, равной его математическому ожиданию, проводят обратное Фурье-преобразование отфильтрованного спектра, определяют полноту насаждения Р как отношение количества выделенных в результате обратного Фурье-преобразования пикселей (nс) к их общему количеству в матрице, находят прикрепляющую точку огивы насаждения к его изображению расчетом среднего количества деревьев в насаждении где Dср= I/Fср- средний размер крон древесного полога, S - площадь анализируемого участка насаждения, Н=7Dcр 1,2 - средняя высота деревьев, а запас насаждения вычисляют по массовым таблицам, связывающим видовые числа fi, высоту hi, ступени толщины (d1,3i), площади сечений  и количество деревьев в насаждении Ni данной ступени толщины из соотношения V = qihifiNi.Популярные патенты: 2027757 Способ получения растений - регенерантов in vitro ... - регенеранты получают in vitro, при этом для увеличения регенерации в питательные среды для инициации каллуса и формирования регенерантов дополнительно вводят ацетон в количестве 10000 - 20000 мг/л. Изобретение относится к биотехнологии, а именно к получению регенератов in vitro. Развитие методов биотехнологии, в основе которых лежит выращивание изолированных органов, тканей и клеток на искусственных питательных средах с последующей регенерацией целых растений, определяет все большее внедрение их в практику сельского хозяйства. Они способны значительно ускорить селекционный процесс и значительно расширить границы воздействия человека на живую природу. При разработке методов ... 2402211 Способ получения трансгенных кроликов, продуцирующих белки в молочную железу ... которые обладают следующими особенностями по сравнению с прототипом (К.Chang et al., 2002).Во-первых, эти антитела, как было известно, связываются со сперматозоидами человека. Во-вторых, нами экспериментально показано их способность связываться с ДНК (плазмидой гена человеческого гранулоцитарного колониестимулирующего фактора ч-ГКСФ) и ДНК (плазмидой гена человеческого гранулоцитарного макрофагового колониестимулирующего фактора ч-ГМКСФ).В-третьих, методом проточной цитометрии были выявлены человеческие моноклональные антитела (мат), специфически связывающиеся со сперматозоидами. Наиболее высокий процент связывания со сперматозоидами показали мат CD71 и CD25, матCD95, CD59, CD11b ... 2056755 Способ регулирования роста овощных культур ... нитратов в овощах при отсутствии токсичности для организма человека растения обрабатывают за 2 - 7 сут до уборки, до полного смачивания поверхности листьев 0,1 - 0,5%-ным водным раствором цитокининового препарата Сапфир, представляющего собой очищенную и сконцентрированную вытяжку из отходов производства пива, содержащих корешки проростков ячменя. 8 табл. , , ОПИСАНИЕ ИЗОБРЕТЕНИЯ К ПАТЕНТУ Изобретение относится к сельскому хозяйству, а именно к агротехнике выращивания овощных культур, в частности салата, томатов и др. Известен способ выращивания овощных культур, включающий обработку растений томатов в фазе формирования завязи плодов на первой и второй ... 2152151 Гербицидная водорастворимая гранулированная композиция ... либо Y и Z могут быть взяты вместе с образованием -СН=СН-СН=СН-. 7. Композиция по п.3, где второй гербицид представляет натриевую соль дикамба. 8. Композиция по п.6, где основание присутствует в соотношении около 1,1 - 1,8 мольных эквивалентов основания к мольным эквивалентам имидазолиноновых гербицидов. 9. Композиция по п.7, где имидазолиноновым гербицидом является 5-этил-2-(4-изопропил-4-метил-5-оксо-2-имидазолин-2-ил)никотиновая кислота. 10. Композиция по п.9, где основанием является карбонат ... 2407280 Устройство и способ для осушения воздуха в теплице и теплица ... по сравнению с другими известными системами и способами осушения и охлаждения относятся:Низкие затраты на оборудование, поскольку не требуются отдельные камеры конденсаторов или вентиляторы для проведения воздуха теплицы в конденсатор. Существенно более низкие эксплуатационные расходы, так как появилась возможность исключить потребляющие по существу наибольшую часть энергии устройства, используемые в прежних системах, то есть вентиляторы.Данный способ действует повсюду глобально, где имеется в наличии достаточно холодной воды, или вода может быть охлаждена посредством достаточно сухого наружного воздуха. Использование данного способа не создает тип проблем, связанных с ... |
Еще из этого раздела: 2171570 Устройство для группового учета надоев молока при доении 2066320 Производные тиазола, способ их получения и способ борьбы с грибками 2015654 Теплица для подземной выработки 2060624 Валкообразующий транспортер жатки-накопителя 2414114 Зерноуборочный комбайн 2423042 Электронно-оптический способ регулирования технологии производства агропродукции 2154938 Способ охлаждения молока на животноводческих фермах и устройство для его осуществления 2492650 Микроэмульсионная бактерицидная композиция 2005344 Способ облучения живых организмов или растений 2048055 Устройство для отрезания и погрузки сенажа и силоса |
Изобретения в сельском хозяйстве
Обработка почвы в сельском и лесном хозяйствах
Посадка, посев, удобрение
Уборка урожая, жатва
Обработка и хранение продуктов полеводства и садоводства
Садоводство, разведение овощей, цветов, риса, фруктов, винограда, лесное хозяйство
Новые виды растений или способы их выращивания
Производство молочных продуктов
Животноводство, разведение и содержание птицы, рыбы, насекомых, рыбоводство, рыболовство
Поимка, отлов или отпугивание животных
Консервирование туш животных, или растений или их частей
Биоцидная, репеллентная, аттрактантная или регулирующая рост растений активность химических соединений или препаратов
Хлебопекарные печи, машины и прочее оборудование для хлебопечения
Машины или оборудование для приготовления или обработки теста
Обработка муки или теста для выпечки, способы выпечки, мучные изделия
|
|
||