气体传感器鉴别花椒产地研究

2019-11-08 01:03陈晓燕陶怀亮李蒙良
农业工程学报 2019年18期
关键词:汉源产地花椒

庞 涛,杨 霄,陈晓燕,陶怀亮,李蒙良

气体传感器鉴别花椒产地研究

庞 涛1,杨 霄1,陈晓燕2,3※,陶怀亮1,李蒙良1

(1. 四川农业大学机电学院,雅安 625000; 2. 四川农业大学信息工程学院,雅安 625000;3. 四川农业大学农业信息工程四川省重点实验室,雅安 625000)

目前花椒产地鉴别基本以感官评定为主,缺乏客观性,在实施应用时难以做到量化和标准化,难以做出判断。因此设计研发一种快速鉴别花椒的智能装置。该装置以气体传感器阵列为核心,能够独立对花椒气味信息进行检测和鉴别,区分不同产地的同类花椒。利用主成分分析和Wilks Λ统计分析对检测数据进行处理。提取主成分5个,累积贡献率为94.41%,其对应Fisher判别模型训练集平均准确率达到88.6%,验证集90%,Wilks Λ统计分析最终选取8个变量,其对应判别Fisher模型训练集平均准确率91.82%,验证集95%。对Wilks Λ统计所选取变量建立细分类交叉验证的Fisher判别模型,平均正确率达到97.27%,将模型移植到采集装置,完成智能花椒品种鉴别装置。该方法是一种简便高效的花椒品种鉴别方法,可为今后进一步研究花椒产地、分级提供检测仪器和理论依据。

传感器;农作物;气体监测器;花椒产地鉴别;Fisher判别

0 引 言

花椒是中国的重要经济作物,中国西部地区如四川、陕西、甘肃、青海等地花椒种植的面积很大,成为重要的经济作物甚至是地方经济的支柱产业[1]。市面上花椒品种众多,外形、色泽、风味各不相同,这也造成了不同品种的花椒价格差异很大,中国约有39种花椒,14个变种[2],大部分花椒品种仍处于野生状态,人工栽培的花椒主要分为青花椒和红花椒,不同产地的花椒质量不同。但随着花椒深度开发和利用得到越来越多的重视,市场需求快速增长,随之而来的是市场上以劣充优、以陈冒新、掺假掺杂等现象层出不穷。一方面,消费者开始更多地关心他们所消费的农产品的来源与真实性,另一方面,企业与花椒原产地渴望寻求保护自身品牌的有效方法,因此研究开发出一种简单、快速、无损的花椒产地鉴别检测方法,具有重要的现实意义。

现有花椒品种鉴别的常用方法有感官分析技术、气相色谱法、与气质联用技术(GC-MS)等。感官辨识对于专业人员的要求较高,需要长期的经验,而这些专业人员很难做到将不同地区的花椒完全辨认。而如吴习宇等采用近红外光谱技术鉴别花椒产地[3],吴莉莉等采用机器视觉对花椒品种进行鉴别[4]。此外,也有众多研究发现不同种或不同品系的花椒间化学成分也有较大差异[5]。这一类的非感官的鉴别手段,测试周期较长,运行成本高,在操作上不利于推广。

气体传感器是一种模仿生物嗅觉的气体检测系统[6]。而气体传感器的核心部件为由多种气敏传感器组成的气敏传感器阵列,其原理是将多个传感器感知到的时间或空间上互补或冗余的信息,并在某种准则下进行综合与分析,以获得单个或单类传感器无法获得的有价值的综合信息,从而形成对观测对象客观的描述[7-9]。鉴于气体传感器对待测气体样品的信息综合分析能力,它已经在农业生产[10-12]、生物医学[13-15]、环境监测[16-18]、食品检测[19-21]等领域得到了广泛应用。

本文自行研发设计出一种基于气体传感器阵列的花椒信息采集装置,用于4种花椒气味信息检测,依据Wilks Λ统计量对气体传感器中的传感器阵列数据进行优化,剔除冗余信息。将优化数据利用Fisher判别建立分类器,结合花椒信息采集装置,建立一种快速、简单、准确且便于推广的花椒产地鉴别方法,为保证名品花椒、产地优选提出新的方法。

1 材料与方法

1.1 材料选取

试验花椒分为4个品种,分别是四川汉源红花椒、陕西韩城红花椒(大红袍)、四川汉源青花椒、云南鲁甸青花椒(后文简称云南青花椒)。4种花椒均取自原产地种植园区,均为干燥花椒果皮。其中汉源青、红花椒和云南青花椒各取60份,陕西红花椒取40份,每份均15 g,此220份花椒用作训练样本。各类花椒取样本各20份作为验证集。

1.2 试验设备

智能花椒品种鉴别装置。该装置系统结构图如图1a所示,作为整个设备的设计指导。图1b为实际设计时的气室结构设计剖图,从图中可以看出,气室,样品槽,传感器阵列,加热片,离心风机的组装结构,风机从上端吸收被加热片加热的花椒气体,由下侧向四周的传感器阵列散发,使得各个传感器均匀接触到花椒气体,气体再由外侧返回气室,达成气体循环。图1c为装置设计的外部结构,从外部只能看到气室和样品槽,无法看到气室内部的具体结构。图1d为拆下外壳的实物俯视图,俯视结构中可以看到本装置在设计时包含了气体流通通道,可以保证气体实现循环。图1e是在1d的基础上,拆除传感器阵列的保护外壳,传感器设置为倒置摆放,阵列通道内为线路连接,只将传感器探头露出,保证传感器传递线路不受到花椒气味的颗粒影响,确保实验安全进行。图1f为实物的侧放图,将侧面挡板取下,内部为搭建好的控制电路及外围电路等,挡板可拆卸保证了随时可以对装置程序,电路进行修改。

1.样品槽 2.气敏传感器阵列 3.离心风机 4.PTC加热片 5.隔离壳体 6.液晶显示屏 7.开关 8.风速调节旋钮 9.控制芯片(内部)

本装置所采用的传感器选型时参照花椒气体组成成分进行选择,并对传感器本身的灵敏性,恢复性,稳定性进行筛选,最终所选的7个传感器分别为MQ135、MQ5、MQ2、TGS2611、TGS2600、TGS2610和TGS2602。图2为本文所用气体传感器阵列实物图,表1为传感器所对应敏感响应特性。该装置有采集检测和鉴别种类2种工作模式,检测到的数据和鉴别结果会在显示屏上显示。

图2 试验所用传感器

表1 气体传感器响应特性

1.3 信息采集和选择

首先在采样前,将颗粒状活性炭置入样品槽中,进行清洗,直至传感器阵列的响应信号稳定,将此时传感器阵列的响应值作为基准值,清洗时间为10 min。将花椒样本送入样品槽,开启气味采集工作模式,先对气味采集装置进行约60 s的预热,使传感器阵列处于正常工作状态,气室内的温度稳定在26℃(±1℃)。进行数据采集,采样时间约为10 min,每个传感器分别采集50次数据,每次采样间隔10 s,取50次采样数值的平均值(aver),最大值(max),最小值(min)作为一个样本的记录值,每组花椒样本包含21个数据参数,构成完整的数据带。

1.4 数据分析

1.4.1 主成分分析

主成分分析(principal component analysis, PCA),是一种通过正交变换将一组可能存在相关性的变量转换为一组线性不相关的变量,从而代表所有变量的统计方法,转换后的这组变量叫主成分[22]。对测定花椒的21个参数进行主成分分析,选取方差累积贡献率达到90%以上的稳定水平的主成分,并确定其数量。

1.4.2 Wilks Λ统计分析

Wilks Λ统计量是常用的数据筛选方式,用于检验多个母体的判别效果和各个变量的判别能力[23]。其实质是,样本组内离差平方和与样本总体离差平方和之比。其比值越小说明样本贡献越大。

对2种分析方法选定的成分分别建立Fisher判别分析,选择最优方式并将最优判别分析移植到鉴别装置中。

2 结果与分析

2.1 主成分提取与Wilks Λ统计筛选

对220个样本所测试花椒21个特征参数进行主成分分析,方差累积结果如图3,主成分分析在第5个以后,累积方差逐增长渐趋于平稳,故提取主成分5个,方差累积贡献率达到94.41%。

图3 主成分方差累积贡献率

Wilks Λ统计分析筛选出8个变量,涉及试验传感器6个,如表2所示。

表2 Wilks Λ统计量选择

2.2 Fisher判别分析

根据筛选的变量PCA和Wilkss Λ统计筛选的变量,分别建立PCA-Fisher判别分析模型与Wilks-Fisher判别分析模型,2种模型均包含3个判别函数,PCA-Fisher判别函数如下:

PF1=1.655×PC1+0.101×PC2−0.247×PC3−

1.362×PC4+0.739×PC5(1)

PF2=0.224×PC1−0.376×PC2+0.073×PC3+

0.619×PC4−0.531×PC5(2)

PF3=0.072×PC1+0.498×PC2+0.308×PC3+

0.519×PC4−0.534×PC5(3)

式中PF1,PF2,PF3表示分类坐标值,PC1~PC5表示5个主成分。

Wilks-Fisher判别函数如下:

WF1=−0.001×1−0.031×2+0.028×3−0.045×4+

0.015×5+0.026×6+0.005×7−0.004×8−5.56(4)

WF2=−0.013×1−0.049×2+0.002×3+0.093×4−

0.048×5+0.066×6−0.001×7+0.036×8−2.417(5)

WF3=−0.002×1−0.008×2−0.024×3−0.074×4−

0.138×5−0.005×6+0.034×7+0.116×8+5.101 (6)

式中WF1,WF2,WF3表示分类坐标值,1~8表示表2中序号1~8的8个筛选变量。将输入导入2种模型,分别与其对应质心相比较,求得距离最短即为分析结果。分析结果如表3所示。

表3 Fisher判别分类结果

由表3知,Wilks-Fisher判别分析的整体正确率相对较高。其中在PCA-Fisher判别模型中,云南青花椒与汉源红花椒的识别率明显较低,这是因为,不同花椒的挥发性成分相对含量不同。文献[24]中指出,汉源青花椒和汉源红花椒的挥发性气体成分含量差异极大,汉源青花椒醇类化合物含量较高,而汉源红花椒酯类化合物较高,这是导致传感器在响应2种花椒时数值差异的主要原因,以表2中出现最多参数的传感器TGS2602为例,图4a是传感器TGS2602每组所采集样本的平均值(aver)折线图,从中可以明显得看出4种花椒的响应差异值较为明显;图4b是表2中未筛选传感器TGS2611每组采集样本的平均值折线图,可以发现,该传感器对4种花椒的响应值差异较低,印证了该传感器所提供信息存在较多冗余,而PCA会提取这一部分的冗余信息,这是导致PCA-Fisher模型识别精度较低的原因,也证明了Wilks统计分析对剔除冗余数据的有效性。图4 c是4类花椒在Fisher判别分析坐标下的坐标点,从图中可以看出,陕西红花椒相对其他3种花椒区分度明显,云南青花椒和汉源青花椒、汉源红花椒在一定程度上具有重叠部分,这可能是因为汉源鲁甸2地相隔较近,在地理位置,土壤等因素较为接近引起的,这也导致了云南青花椒和汉源红花椒的判别率相对较低,故仅仅使用Fisher判别分析,并不能完全将产地鉴别,还需要对函数或数据进行优化分类。

2.3 细分类交叉验证下的Wilks—Fisher判别

对Wilks—Fisher判别模型进行改进,将样本集进行区分,分别训练出红花椒的判别函数和青花椒的判别函数。为降低交互数据的干扰性,对分析案例进行交叉验证。交叉验证是在机器学习建立模型和验证模型参数时常用的办法。交叉验证,就是重复的使用数据,把得到的样本数据进行切分,组合为不同的训练集和测试集,用训练集来训练模型,用测试集来评估模型预测的好坏。在此基础上可以得到多组不同的训练集和测试集,某次训练集中的某样本在下次可能成为测试集中的样本,即所谓“交叉”。

图4 传感器部分响应及判别坐标

采用S折交叉验证,根据样本数量S,将对应训练集样本分成S-1份,剩下1份作为测试机,循环S次,直到每个样本均作为测试集对象进行测试。在判别时,将数据分为青红花椒进行分别验证,结果如表4所示。

表4 青红花椒分类交叉验证结果

由表4可知,在分类判别时,汉源青花椒的正确率为98.3%,误判1组,云南青花椒的正确率为93.3%,误判4组,陕西红花椒的正确率为100%,汉源红花椒的正确率为98.3%,误判1组,平均正确率达到97.27%,所有验证结果较原模型检测均有所提高。该验证结果表明分类交叉验证的Wilks—Fisher判别模型可以较好地鉴别花椒的产地,将判别函数写入采集装置并调试程序,即构成完整的智能花椒品种鉴别系统。

3 结 论

本文采用自制的智能花椒品种鉴别装置,利用气体传感器阵列对花椒数据进行采集,通过主成分分析、WilksΛ统计分析提取有效变量,建立对应的Fisher判别模型。试验表明,Wilks—Fisher具有更好的结果,训练集的平均判别正确率达到91.82%,验证集平均判别正确率达到95%。在细分类交叉验证的下的该模型判别率达到97.27%。该模型能够正确地识别4个产地的花椒,其检测结果较为理想。而将模型函数写入装置后,构成智能花椒品种鉴别系统装置,实现花椒产地的无损智能鉴别。该装置无需样品预处理,不需要使用检测理化值所需的精密仪器,检测成本较高光谱等精细设备极具优势,为花椒产地的识别及追溯提供了技术支持,可以有效控制以次充优、假冒伪劣等掺假手段,为名优名品花椒,具有地理标志保护花椒的鉴别提供了简便快捷的新思路,具有广阔的应用前景。

[1] 毕君,赵京献,王春荣,等. 国内外花椒研究概况[J]. 经济林研究,2002,20(1):46-49.

[2] Huang Chengjiu, Rutaceae, Dicotyledoneae, Angiospermae. Flora Reipublicae Popularis Sinicae[M]. Beijing: Science Press, 1997.13.

[3] 吴习宇,祝诗平,黄华,等. 近红外光谱技术鉴别花椒产地[J]. 光谱学与光谱分析,2018(1):68-72.

Wu Xiyu, Zhu Shiping, Huang Hua, et al. Near infrared spectroscopy for determination of the geographical origin of Huajiao[J].Spectroscopy and Spectral Analysis, 2018(1): 68-72. (in Chinese with English abstract)

[4] 吴莉莉,邢玉清,郑宝周,等. 机器视觉在花椒品种鉴别中的初步研究[J]. 传感技术学报,2016,29(1):136-140.

Wi Lili, Xing Yuqing, Zheng Baozhou, et al. A preliminary study of zanthoxylurn bungeanum maxim varieties discriminating by computer vision[J]. Chinese Journal of Sensors and Actuators, 2016, 29(1): 136-140. (in Chinese with English abstract)

[5] Ito C, Katagiri H, Sato A, et al. Pharm acogno stical Studies on the Sino-Japanes crude drugs “Huajiao” and “Sansho”: Part5.On essential oils and pungent priciples of pericarps of subgen Zanthoxylum plants and commercial samples[J].Natural Medicines, 1997, 51(3): 249.

[6] Chen F Q. Means for rapid detection of toxic gases summary[J]. Autom Instrum, 2014(10): 79.

[7] 韩崇昭,朱洪艳,段战胜,等. 多源信息融合[M]. 北京:清华大学出版社,2006:15-16.

[8] 王欣. 多传感器数据融合问题的研究[D]. 吉林:吉林大学,2006:13-15.

[9] Hall D L, Llinas J. Handbook of Multisensor Data Fusion[M]. New York: CRC Press, 2001: 10-13.

[10] 周亦斌,王俊. 基于电子鼻的番茄成熟度及贮藏时间评价的研究[J]. 农业工程学报,2005,21(4):113-117.

Zhou Yibin, Wang Jun. Evaluation of maturity and shelf life of tomato using an electro nic nose[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2005, 21(4): 113-117. (in Chinese w ith Eng lish abstr act)

[11] 殷勇,郝银凤,于慧春. 基于多特征融合的电子鼻鉴别玉米霉变程度[J]. 农业工程学报,2016,32(12):254-260.

Yin Yong, Hao Yinfeng, Yu Huichun. Identification method for different moldy degrees of maize using electronic nose coupled with multi-features fusion[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(12): 254-260. (in Chinese with English abstract)

[12] 张婷婷,孙群,杨磊,等. 基于电子鼻传感器阵列优化的甜玉米种子活力检测[J]. 农业工程学报,2017,33(21):275-281.

Zhang Tingting, Sun Qun, Yang Lei, et al. Vigor detection of sweet corn seeds by optimal sensor arraybased on electronic nose[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(21): 275-281. (in Chinese with English abstract)

[13] 郝燕,王永清,王悦,等. 检测细胞代谢气体成分诊断肺癌的方法[J]. 浙江大学学报:工学版,2008,42(2):294.

[14] Brekelmans M P, Fens N, Brinkman P, et al. Smelling the diagnosis: The electronic nose as diagnostic tool in inflammatory arthri-tis[J]. PLOS One, 2016, 11(3): e0151715.

[15] 张斌,吴哲,何宏建,等. 结合磁共振成像和脑机接口的新型在体生物电子鼻的研究[J]. 中国生物医学工程学报,2018,37(1):57-63.

Zhang Bin, Wu Zhe, He Hongjian, et al. Characterization of in vivo bioelectronic nose with combined manganese- enhanced mri and brain-computer interface[J]. Chinese Journal of Biomedical Engineering. 2018, 37(1): 57-63.

[16] 方向生,施汉昌,何苗,等. 气体传感器在环境监测中的应用与进展[J]. 环境科学与技术,2011,34(10):112.

[17] Tian F C, Zhang J, Yang S X, et al. Suppression of strong back-ground interference on E-nose sensors in an open country environ-ment[J]. Sensors, 2016, 16(2): 233.

[18] 俞守华,董绍娴,区晶莹. 猪舍有害气体NH3、H2S的电子鼻定量识别[J]. 农业工程学报,2009,25(7):153-157.

Yu Shouhua, Dong Shaoxian, Ou Jingying. Quantitative identification of pernicious gases NH3and H2S in piggery using electronic nose[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(7): 153-157. (in Chinese with English abstract)

[19] 殷勇,赵玉珍,于慧春. 基于多种变量分析方法鉴别食醋种类电子鼻信号特征筛选[J]. 农业工程学报,2018,34(15):290-297.

Yin Yong, Zhao Yuzhen, Yu Huichun. Feature selection of electronic nose signal for vinegar discrimination based on multivariable analysis [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(15): 290-297. (in Chinese with English abstract)

[20] 沈飞,刘鹏,蒋雪松,等. 基于电子鼻的花生有害霉菌种类识别及侵染程度定量检测[J]. 农业工程学报,2016,32(24):297-302.

Shen Fei, Liu Peng, Jiang Xuesong, et al. Recognition of harmful fungal species and quantitative detection of fungal contamination in peanuts based on electronic nose technology[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(24): 297-302. (in Chinese with English abstract)

[21] 裴高璞,史波林,赵镭,等.典型掺假蜂蜜的电子鼻信息变化特征及判别能力[J]. 农业工程学报,2015,31(增刊1):325-331.

Pei Gaopu, Shi Bolin, Zhao Lei, et al. Information variation feature and discriminant capabilities of electronic nose for typical adulteration honey identification[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(Supp.1): 325-331. (in Chinese with English abstract)

[22] 金慧芳,史东梅,陈正发,等. 基于聚类及 PCA 分析的红壤坡耕地耕层土壤质量评价指标[J]. 农业工程学报,2018,34(7):155-164.

Jin Huifang, Shi Dongmei, Chen Zhengfa, et al. Evaluation indicators of cultivated layer soil quality for red soil slope farmland based on cluster and PCA analysis[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(7): 155-164. (in Chinese with English abstract)

[23] 高惠璇. 应用多元统计分析[M]. 北京:北京大学出版社,2005:199-205.

[24] 蒲凤琳,孙伟峰,车振明,等. 水蒸气蒸馏结合GC-MS法分析比较四川汉源青、红花椒挥发性香气成分[J]. 中国调味品,2017,42(1):23-27.

Pan Fenglin, Sun Weifeng, Che Zhenming, et al. Analysis of volatile aroma components in hanyuan green and redby steam distillation and GC-MS[J]. China Condiment, 2017, 42(1): 23-27. (in Chinese with English abstract)

Identification oforigin based on gas sensor

Pang Tao1, Yang Xiao1, Chen Xiaoyan2,3※, Tao Huailiang1, Li Mengliang1

(1.,,625000,; 2.,625000,;3.,625000,)

At present, the identification of the origin ofis basically based on sensory evaluation, lack of objectivity, and it is difficult to quantify standardize when applying, and is difficult for non-professionals to make judgments. Therefore, in this paper, a smart device to quickly identifywas designed and developed. The device was based on the gas sensor array, including a control module, a temperature module, a data storage module, a fan module, and a display module, it could not only independently detect and identify the odor information of the, but also distinguish the same kind offrom different places. The sensor array contained seven gas sensors, which could respond to irritating gases emitted bysuch as benzene, alkanes, alcohols, and aldehydes. When the temperature was stable at about 26 degrees Celsius, it could effectively collect information on the odor emitted by. Each group ofsamples was collected 50 times, and the average value, the maximum value, and the minimum value were taken as sample recording parameters. In this paper, four kinds ofwere selected as experimental subjects. Two kinds of greenwere from Ludian in Yunnan and Hanyuan in Sichuan. At the same time, the two kinds of redwere from Hancheng in Shaanxi and Hanyuan in Sichuan. A total of 220 samples were collected as training sets, including 40 redin Shaanxi and 60 samples in the remaining three samples. Another 80 samples were taken as the verification set, the number of samples for eachwas 20 in the verification set as well. The detection data were processed using principal component analysis (PCA) and Wilks statistical analysis. Five principal components were extracted, and the cumulative contribution rate was 94.41%. The average accuracy rate of the training model corresponding to the Fisher discriminant model was only 88.6%, and the verification set was 90%. As a comparison, the Wilks statistical analysis finally eliminated 13 variables as well as selected 8 variables, and only TGS2611 sensor acquisition was not used. The average accuracy of the Fisher model training set was 91.82%, and the validation set was 95%. The results of the comparison of the two models indicate that the variables screened by Wilks are more effective in discriminating thefield. Among the four kinds of, the recognition rate of Yunnan greenand Hanyuan redwas relatively lower than the others, and there was a phenomenon that the boundary data overlaps in the discrimination result graph. Then, to solve the problem, a Fisher discriminant model with cross-validation was established for the variables selected by Wilks statistic. In addition, the average accuracy rate reached 97.27%. Finally, the model was transplanted to the collection device to complete the identification device of intelligentvariety. It was a simple and efficient method for identifyingvarieties and could provide a testing instrument and theoretical basis for further research on the origin and classification of.

sensors; crops; gas detectors; identification offield; fisher discriminant

庞 涛,杨 霄,陈晓燕,陶怀亮,李蒙良. 气体传感器鉴别花椒产地研究[J]. 农业工程学报,2019,35(18):267-272.doi:10.11975/j.issn.1002-6819.2019.18.032 http://www.tcsae.org

Pang Tao, Yang Xiao, Chen Xiaoyan, Tao Huailiang, Li Mengliang.Identification oforigin based on gas sensor[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(18): 267-272. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.18.032 http://www.tcsae.org

2019-03-18

2019-04-26

四川省教育厅一般项目(自然科学)立项编号:17ZB0333 基于高光谱图像技术的苹果品质无损检测方法研究

庞 涛,讲师,主要从事农业信息检测。Email:349380993@qq.com

陈晓燕,教授。主要从事农产品无损检测。Email:chenxy@sicau.edu.cn

10.11975/j.issn.1002-6819.2019.18.032

S-3

A

1002-6819(2019)-18-0267-06

猜你喜欢
汉源产地花椒
乌斯河汉源
推动产地农产品追溯进入下游流通消费环节
当前花椒价格走低 椒农如何积极应对
国内花椒产业进入低谷期
印尼燕窝产地探秘
警惕“洗产地”暗礁
食物离产地越远越好
乡村振兴的汉源实践
花椒泡脚好处多
第六届阳光汉源赏花旅游文化月开幕