现代水稻产业技术体系专题
  首页 > 新闻 > 科研新闻 > 精细农业中农情信息采集技术的研究进展  
精细农业中农情信息采集技术的研究进展  
作者:罗锡文 臧英 等
来源:见觚钿
发布时间:2006-2-15 11:51:29

0 引言

  精细农业是20世纪80年代初国际农业领域发展起来的一门新兴的跨学科综合技术。它是通过可变尺度的田块区别管理的方式,来确定经济、合理的投入,高效利用农业资源,从而获得经济、环境等方面最高回报的一种管理策略和技术体系。其本质是基于农田小区内作物生长环境和产量差异性的变量管理,因此,快速有效地采集和描述影响作物生长环境的空间变量信息,是开展精细农业实践的重要基础。

  农情信息主要包括地理环境、土壤环境、小气候、水环境、与作物生长状况相关的信息以及管理信息等六大要素,具有量大、多维(信息多种多样)、动态、不确定(系统的噪声或随机噪声)、不完整、稀疏性、时空变异性强等特点。在精细农业研究中,目前优先需要考虑的是土壤水分、土壤养分、土壤压实、耕作层深度和作物病虫草害及作物苗情分布信息等,要求能够定位、快速、精确、连续地测量。传统的实验室分析方法已很难满足这一要求,为此,需要开发适用于精细农业的农田信息快速采集技术。目前,农情信息快速采集技术的研究落后于支持精细农业的其他技术(如产量传感器技术等)的发展,已成为国际上众多单位研究的重要课题。

  本文主要对土壤环境、作物产量分布、病虫草害及作物长势等信息采集技术的研究情况进行分析,并提出今后进一步开展农情信息采集技术研究的方向。

1 土壤环境信息获取

  目前,国内外开展精细农业的实践研究,大多是从研究农田土壤特性的空间变异性开始的,研究内容主要集中在土壤养分( N、P、K、SOM 等)、土壤水分、电导率、pH 值、耕作阻力(圆锥指数)与耕作层深度等要素的快速采集方面。

1.1 土壤养分

  土壤养分的快速测量一直是精细农业信息采集技术的难题。目前采用的测量仪器主要有3类:一是基于光电分色等传统养分速测技术基础上的土壤养分速测仪,国内已有产品投入使用,其稳定性、操作性和测量精度虽然尚待改进,但对农田主要肥力因素的快速测量具有实用价值。河南农业大学开发的便携式 YN 型土壤养分速测仪,相对误差为5%~10%,尽管每个项目测试所需时间仍在40~50 min 之间,但较传统的实验室化学仪器分析在速度上提高了20倍;二是基于近红外技术( NIR )通过土壤或叶面反射光谱特性直接或间接进行农田肥力水平快速评估的仪器,已在试验中使用,如 Hummel 等通过 NIR 土壤传感器测量土壤在1603~2598 nm 波段的反射光谱,预测土壤的有机质和水分含量,测量相对误差分别为0.62%和5.31%;三是基于离子选择场效应晶体管( ISFET )集成元件的土壤主要矿物元素含量测量仪器,在国外已取得初步进展, S.J.Birrell 等研究的 ISFET/FIA 土壤分析系统,可以在1.25 s 的时间内完成土壤溶液硝酸钠浓度的分析,基本上可以满足实时采样分析的要求。

1.2 土壤水分

  土壤水分是土壤的重要组成部分,土壤水分的测量是精细农业中实施节水灌溉的基础。土壤水分传感器技术的研究与发展直接关系到精细农业变量灌溉技术的研究与发展。目前,土壤水分传感器的种类很多,根据其测量原理不同,土壤水分的测量方法可分为:

  1) 基于时域反射仪( TDR )原理的测量方法

  该方法目前在国际上比较流行,且在技术上较为成熟,已有商业化产品,仪器探头长度15~70 cm,精度可达98%。当前主要问题是降低生产成本。目前,中国土壤水分的测量研究主要是在 TDR 基础上开发经济实用的基于驻波比、频域法原理、近红外技术的快速测量仪。最常用的快速水分信息采集仪器有 SWR 型土壤水分传感器。裘正军等针对精细农业中节水灌溉的需要,采用 SWR-2 型土壤水分传感器,研制了基于 GPS 定位的土壤水分快速测量仪,并通过 GIS 软件生成农田的墒情图,但尚未实现实时在线测量。

  2) 基于中子法技术的测量方法

  中子土壤水分计能及时采集数据,不扰动被测土壤,测量深度不限,田间操作简单,携带方便。高国治等根据不同领域对实际测量的具体要求,开发了智能化中子土壤水分计。但该方法测定结果与传统方法测定结果之间存在一定差异,而且土壤表层因中子外溢问题使土壤水分含量的测定值偏差较大。

  3) 基于土壤水分张力的测量方法

  该方法的特点是产品价格低,但测量响应速度较慢。

  4) 基于电磁波原理的测量方法

  采用该方法的土壤水分传感器测量精度可达98%,且响应速度快,产品成本也较低,适用于便携式田间土壤水分采集,但该方法测量结果受土壤密度的影响比较大。

1.3 土壤电导率

  土壤学的研究结果表明,土壤电导率能不同程度地反映土壤中的盐分、水分、有机质含量、土壤质地结构和孔隙率等参数的大小。有效地获取土壤电导率值,对于确定各种田问参数时空分布的差异具有一定意义。快速测量土壤电导率的方法有:

  1) 电流—电压四端法

  该方法属于接触式测量方法,虽为接触测量但却不需要取样,基本不用扰动土体,而且在作物生长前和生长期间都可以实现实时测量,可测量不同深度的土壤电导率。但要求土壤和电极之间接触良好,在测量干或多石土壤电导率时可靠性差。李民赞等人开发了一种基于电流—电压四端法的适合中国较小地块应用的便携式土壤电导率实时分析仪。试验结果表明,不仅可以满足大田裸地实时测量要求,也可以适应作物生长期和设施栽培中的实时测量;在 Schlumberger 组态( a=15,b=30 )结构下,传感器可获得最高的测量精度;非线性回归建模结果还显示幂函数模型具有最优的预测结果,测定精度达到0.994,模型的常数项与传感器结构有关,乘幂项则不随传感器结构参数变化,而受土壤质地、土壤密度等因素的影响需要通过田间试验进一步验证。Carter,Rhoades 等人开发了车载式电流—电压四端法测量土壤电导率的测量设备,但它只适于大块田地使用,而且只能在作物种植前测量。

  2) 基于电磁感应原理的测量方法

  该方法属于非接触式测量方法,测量精度高,可测量不同深度的土壤电导率,但测量深度不如电流—电压四端法容易确定。Myers 等利用电磁感应的方法实现了电导率的非接触式检测,该传感器与土壤生产力指数联系起来,能综合反映土壤容重、持水量、盐度、pH 值等参数。Carter,Rhoades 等人开发了车载式基于电磁感应原理测量土壤电导率的测量设备,可在作物生长期进行测量。

  3) 基于时域反射仪原理的测量方法

  该方法测量精度高,可在同一点连续自动测量。但由于是点对点测量,绘制的土壤电导率图空间分辨率比上述两种方法低。

1.4 土壤 pH 值

  土壤 pH 值检测大多采用 pH 试纸或 pH 玻璃电极,前者只能定性检测,后者不能在线测量,因此,这两种方法均无法满足精细农业对农情信息获取技术的要求。目前适合精细农业要求的方法有光纤 pH值传感器和 pH-ISFET 电极。其中光纤 pH 值传感器虽然易受环境光干扰,但在精度和响应时间上基本能满足田间实时快速采集的需要。现已研制出适用于 pH1~14 范围内不同区间 pH 测量的光极,测量精度可达±0.001个 pH 单位,响应时间短,可进行连续、自动测量。如:张靖等研究的以多孔聚合物为载体的光纤 pH 值传感器,可将 pH 值测量的平衡时间缩短到1.3 s;据荆淼等报道,Kopelman 等研制的亚微米光纤 pH 传感器可以在很少量的样品内进行 pH 测量,并且响应时间达到毫秒级。基于 pH-ISFET 电极测量的方法具有良好的精度和较短的响应时间,美国利诺公司( Leeds & Northrup )的 Durafet pH 电极,在2.8个 pH 单位的范围内变化时,其稳定响应的时间在1 s 左右。但基于 pH-ISFET 电极测量的方法易受温度影响,需要温度补偿,且电极的寿命较短,在田间实时快速采集时要加以注意。

1.5 土壤耕作层深度和耕作阻力

  圆锥指数 CI( Cone Index )可以综合反映土壤机械物理性质,研究表明,它可用以表征土壤耕作层深度和耕作阻力。

  圆锥指数 CI 是用圆锥贯入仪(简称圆锥仪)来测定的。圆锥仪的研制工作不断发展,从手动贯入到机动贯入,从目测读数到电测记录,出现了多种多样的圆锥仪。美国研制成功了空投圆锥仪,德国正准备利用卫星测定土壤的 CI 值。中国最初都是采用先从现场采样再在室内进行试验的方法。20世纪70年代初,中国农机部门研制了 SR-2 型土壤贯入仪。南京工程兵学院于1981年研制成功工兵圆锥仪 Ⅰ 型。1984年,该学院又研制出工兵圆锥仪 Ⅱ 型。这两种仪器质量轻、结构可靠、使用方便,但仍采用手动贯入,目测读数的方式。华中农业大学和武汉工学院研制的圆锥指数仪采用转动手柄推动圆锥贯入和滚筒自记的方式。1985年,华南农业大学研制成功了电测圆锥、剪切仪。1987年,又对圆锥仪部分进一步进行了改进,这种电测圆锥仪采用单板计算机进行采样、分析、计算和打印,在圆锥指数的自动记录方面前进了一大步。但贯入方式仍采用手动的方法,影响圆锥指数的测量精度。滚筒记录方式虽然解决了自记的问题,但整理记录曲线的工作量仍很大。1996年,罗锡文等研制了一种以电动机为贯入动力,用微型计算机( IBMPC/XT )控制、记录和处理数据的机械式电测圆锥仪,可保证贯入速度均匀,减少对测量精度的影响,该仪器可用于室内土槽、野外土壤的测试,但仍是逐点测量,尚未实现连续测量。

2 作物产量分布信息采集技术

  作物产量是众多影响作物生产因素综合作用的结果,是变量作业管理的重要依据。获取作物产量信息,建立产量分布图,是实现作物生产过程中科学调控投入和制定管理决策措施的重要基础,也是实施精细农业的先决条件。

  在精细农业研究中,国外已商品化的产量监视系统产品集中于谷类作物收获机械方面,主要有美国 CASE IH 公司的 AFS ( Advanced Farming System )系统,英国 AGCO 公司的 FieldStar 系统,美国 John Deree 公司的 Greenstar 系统,美国 AgLeader 公司 PF ( Precision Farming )系统及英国 RDS 公司的产量监测系统等。这些系统都具有功能较强的 GIS 综合功能,能自动完成产量监测和生成产量分布图。中国对谷物产量测产系统的研究才刚刚起步,目前还没有商品化谷物产量监测系统产品。

  谷物流量传感器是整个测产系统的核心,目前应用的谷类作物产量传感器主要有冲量式流量传感器、γ 射线式流量传感器及光电式容积流量传感器。

  1) 冲量式流量传感器

  冲量式流量传感器基于冲击原理,当谷物流冲击感力板时会改变运动方向,造成冲量的变化,在感力板上反映为力的变化,检测该变化即可得到谷物流量值。该测量方法的精度取决于集谷搅笼的速度、谷物类型和粮食的湿度。动态试验表明,该方法测量精度达96%。由于冲量式流量传感器结构简单,使用安全,测量精度较高,被认为是较实用的传感器类型。

  2) γ 射线式流量传感器

  γ 射线式流量传感器虽然可以得到较高的测量精度,但是由于它造价高,且对操作者的人身健康可能造成潜在伤害,限制了其推广。

  3) 光电式容积流量传感器

  光电式容积流量传感器直接测量谷物的容积流量,但测量结果受谷物密度、含水率、机器倾斜度、探头污染等的影响,需经常清洗和标定,性能不稳。测量误差为±3%。

  影响产量的另一主要因素是谷物水分,因此谷物水分传感器也是测产系统的重要组成部分。常见的谷物水分传感器有电阻式、电容式、红外式、微波式和中子式传感器。

  1) 电阻式水分传感器

  电阻式水分传感器已利用于谷物等固体物料的水分连续测量,其特点是结构简单,价格便宜,但对低水分或高水分测量不准,对被测物料的接触状态要求较高。

  2) 电容式水分传感器

  电容式水分传感器是使用较广泛的一种测水装置,其结构简单,易于进行连续测水,但测量精度不高,稳定性差,不同品种谷物对测量的影响较明显。

  3) 红外式水分传感器

  红外式水分传感器为非接触测量形式,易于快速连续测量水分,但因受被测物料的形状、大小、密度的影响,不能测量物料内部水分。

  4) 微波式水分传感器

  微波式水分传感器优点是灵敏度高,非接触测量,可实现快速连续测定,其缺点是受被测物料的形状、密度的影响,不能测量物料内部水分。

  5) 中子式水分传感器

  中子式水分传感器是非接触式测量,对动态物料可进行快速连续测量,并能测量物料内部水分,不受被测物料的形状、密度、水的形态影响。

3 病虫害、杂草识别

  目前,病虫草害信息的自动快速采集主要是基于计算机图像处理和模式识别技术,以研究作物植株的根、茎、冠层(叶、花、果实)等的形态特征作为诊断判读的目标。主要分析方法有:光谱特征分析法、纹理特征分析法、形状特征分析法等。用于识别的图像有普通的灰度图像、彩色图像、光谱图像和遥感图像等。在宏观上大尺度病虫草害板块信息的快速采集多采用遥感图像的分析处理,而小区内病虫草害的精确定位识别主要借助于机器视觉识别方法。

3.1 杂草识别

  Yonekawa 等人和 Franz 等人采用形状特征分析法,利用作物冠层的形状特征因子(紧密度、圆度、伸长度、叶状度、粗糙度、叶片边缘的弯曲率等)进行杂草的识别研究,该方法分类识别效果较好,且识别效率高,基本能满足田间在线实时识别时对响应速度的要求,但未能很好地解决叶片相互遮挡造成识别率下降的问题。Burks 、Tang 及 Wang 等人采用光谱特征分析法,通过植物的光谱反射特性进行杂草的识别研究,该方法效率高,可以满足田间在线实时识别时对响应速度的要求,但干扰因素多,对设备要求高,且设备昂贵。纪寿文等根据叶片的投影面积、叶长、叶宽等特征量,研究了玉米苗期田间杂草识别方法,该方法能滤除土壤背景,识别出玉米苗期田间细长的单子叶杂草,测量精度达到80%以上,但易受光照等外因的影响,离实际应用尚有距离。王月青等研究了麦田杂草的实时识别系统,根据分割的图像统计杂草的密度,基本达到了对杂草进行实时探察的目的,但对杂草识别精度没有进行深入的试验研究。Goel 等利用航空遥感技术,对玉米和大豆地里的杂草进行了光谱分析,发现位于红外区(675.98~685.17 nm )和近红外区(743.93~830.43 nm )的反射率对识别玉米地里的杂草最有效,而只有位于近红外区(811.40 nm )的反射率对识别大豆地里的杂草有效。

3.2 病虫害识别

  利用基于 SOM 和多层感知的神经网络技术进行小麦叶面的黄萎病害识别,分类精度达到99%,但该识别算法的实时性还需要进行进一步的试验研究验证。陈佳娟等根据棉花叶片的孔洞及叶片边缘的残缺,进行了棉花受虫害程度的自动测定研究,测定精度可达95%。田有文等人根据植物病害的颜色特征,采用统计模式识别的方法,对植物病害斑块进行了识别研究,根据植物病害的纹理特征,提出将支持向量机和色度矩分析方法应用于病害斑块的识别,该方法处理速度较快,准确度高,但未考虑田间复杂的环境干扰因素,其实用性尚需进一步验证。马骏等应用数学形态学算法进行植物病虫识别,采用膨胀和腐蚀快速算法用于提取病虫骨架特征。根据不同种类病虫骨架的几何矩特征,通过神经网络进行病虫分类识别。试验结果表明,这种基于数学形态学和神经网络的植物病虫识别方法正确识别率高,识别速度快,适合于植物病虫种类的识别,但对病虫害程度则无法进行分类识别。M.Zhang 等利用多光谱遥感图像进行西红柿作物的病虫害检测。试验结果表明,近红外波段区(0.7~1.3 μm )比可见光波段在监测作物病虫害方面更有意义,在0.7~0.9 μm 波段对监测西红柿晚期枯萎病效果显著,但目前研究方法仅局限棚内试验阶段,有待通过田间试验进一步完善。Malthus 等用地物光谱仪研究大豆受蚕豆斑点葡萄孢感染后的反射光谱,发现其一阶导数反射率比原始的反射率要高,可用来监测病虫害的感染发生情况,但无法进行病虫害程度的识别。黄木易等对冬小麦条锈病光谱特征进行了研究,结果表明在条锈病害田间病叶率为5%时的最佳防治时期内,高光谱可以对冬小麦条锈病害做出早期诊断。利用多时相的高光谱航空图像可以快速、大面积的监测冬小麦条锈病病害程度,但对病虫害的诊断性探测、评估还需要进行进一步的研究。

4 作物长势信息采集

  作物长势信息是调控作物生长、进行作物营养缺素诊断、分析和预测作物产量的重要基础和根据。目前作物长势信息采集技术方面的研究主要集中在以下两方面:

  一是从宏观角度利用 RS 遥感的多时相影像信息研究植被生长发育的节律特征。目前,国内外许多学者已涉足高光谱遥感在植被生物物理和生物化学信息提取方面的研究,但大部分还停留在从光谱信息反演作物冠层生化参数这一步。杨敏华等以国产机载成像光谱数据为基础,依据田间同步采样数据的光谱响应,建立了小麦冠层理化参数估计模型,实现了航空高光谱遥感田间小麦冠层理化参量的整体反演。国家农业信息化工程技术研究中心联合中国科学院遥感应用研究所利用 OMIS ( Operative Modular Imaging Spectrometer,128个波段)和 pHI ( Push-broom Hyper spectral Imager,244个波段)在北京小汤山国家精准农业示范基地进行了农业遥感试验,并在农作物冠层生化参数反演研究方面取得一定的进展,进一步证明了高光谱遥感在农田信息提取中的巨大前景。

  二是在区域或田块的尺度上,近距离直接观测分析作物的长势信息。目前比较有代表性的采集方法有:

  1) 根茎的测量:传统的方法有钉板法、容器法、同位素测定法、地下根室法、基于计算机视觉技术的离土取样扫描分析法等。罗锡文、周学成等提出了基于多层螺旋 CT 技术的根系原位形态的可视化研究方法,并进行了试验研究,为准确、快速、无损地实现植物根系原位形态的定性观察和定量测量提供了新思路。

  2) 植株高度测量:Searcy 和 Beck 利用安装在机械手上的红外光束传感器实现了实时估测棉花活体植株的高度。高巍等采用40 kHz 超声波测量作物高度,可用于多种农业机械的自动调节系统、自动采摘机器手爪、农业机器人的控制系统、作物生长情况监测等方面,但该装置存在测量误差较大、采集数据点少、不适合室外田间实时批量采集等问题。

  3) 植物生理信息:白广存等对植物生理信息采集处理系统进行了研究,选用微电极传感器,实现了植物电位信号的多点采集;应用微电极传感器和恒流源技术,实现了植物电阻信号的多点采集;采用微重量传感器,实现了植物蒸腾量和蒸腾速率信号的采集;应用红外线气体分析传感器,实现了光合速率信号的采集。

  4) 作物水分胁迫指数( CWSI ):王卫星等以柳叶菜心为研究对象,利用红外测温仪、土壤含水量传感器等采集不同条件下的土壤含水量、空气温湿度和冠层温度等参数,以确定作物缺水指数经验模型的参数,并针对太阳辐射强度对模型做了改进,研制了蔬菜旱情检测系统,可实时测得田间作物的水分亏缺状态。

  5) 营养缺素胁迫:Anom 等开发的实时检测分光光度计,用于作物水分、病害、盐度、营养缺素等胁迫指数分布的在线测量;Ahmad 等用基于 DGPS 的叶绿素检测仪对玉米的氮素胁迫分布进行了研究。李云梅利用美国 ASD ( Analytical Spectral Device )公司的 Fieldspec 光谱仪,对水稻冠层二向反射特性进行了研究,总结了水稻冠层二向反射率随不同氮素水平而变化的规律,并集成了植被光谱模拟系统,通过输入冠层二向反射率等参数,可反演出作物的冠层结构参数和叶片生物化学含量,对监测水稻长势及预测产量很有意义。王人潮等对水稻氮素营养水平与光谱特征的关系进行系统深入的研究,通过光谱测定及其变换运算,可以区分水稻的不同氮素营养水平。赵春江等人开发的归一化差异植被指数( NDVI )仪,在近红外和红光两个特征波段,通过对太阳入射光和植被的反射光进行探测,得到 NDVI 值,反演出叶面积指数、植被覆盖度、发育程度、生物量等指标,可以对作物长势、营养诊断等方面作出评估。

5 农情信息采集系统

  目前,国外已经出现了一批商品化的农情信息采集软硬件产品,软件设计功能强大,尤其是在图形处理、可视化分析方面技术比较成熟,硬件设计上考虑到系统应用的实际环境,可满足机载或徒步采集的不同要求。如 Trimble 公司的 Ag160 Portable Computer 可用于田间徒步的定位数据采集,它具有数据采集、田间绘图、预定义模式导航及指导土壤采样的功能,能实现田间成图,并现场记录各种作物及其生长环境属性信息和各种田间环境传感器的信息,所有数据能通过 ESRI 的 Shape 文件格式与外界的 GIS 分析软件进行信息交换。 Field Worker 公司的基于掌上电脑的 Field Worker 信息采集软件能实时采集带有空间位置属性的田间作物生长状态信息,并能进行相应的计算处理,为作物管理提供科学依据。

  国内基于 GPS 和 GIS 的农情信息定位采集系统的研究较国外起步晚,多数还处于手工采集阶段。目前,一些科研院所正积极着手开发适于中国国情的农情信息采集系统,主要有:

  1) 以单片机为核心进行开发,结合 GPS 技术实现对某几个农情信息要素(如土壤水分、养分、电导率等)的定位测量,然后将采集的信息在 GIS 平台上进行事后管理和分析,这类设备功能简单、实用、开发成本较低,但功能过于专一。

  2) 在便携式计算机上进行开发,可实现大量农情信息现场快速采集及可视化分析处理。目前已开发了能够实现农情信息采集的通用型软件,如中国农业大学开发的基于便携式计算机的农用 GPS 数据采集系统,具有 GPS 数据接收、存储、坐标转换、定位采样和平行导航等基本功能,初步实现了对农田信息的半自动化采集。但便携式计算机存在体积和质量较大、强光下操作不便等问题。

  3) 基于掌上电脑的嵌入式农情信息采集系统的开发。该系统体积小、功耗低、功能强、可靠性高。如饶晓燕等利用桑夏 HC200O ( 863HC )开发了基于掌上电脑的 GPS 数据采集器,把 PDA ( Personal Digital Assistant )技术与农用 GPS 结合起来,实现了数据采集、存储等基本应用功能,体积小,携带方便,但由于 CPU 速度、存取容量,操作系统等方面的限制,应用程序开发仍然受到一定程度的影响。随着掌上电脑的快速发展,其相关的技术也在不断成熟,开发基于掌上电脑的农情信息采集系统必将成为今后的主要发展方向。

6 结语

  综上所述,关于农情信息采集技术,在研究内容方面,最为热门的是土壤环境信息的采集与处理,其次是作物的长势及病虫草害等相关的信息;在时效性方面,大多基于传统的物理原理,现场手工采集与事后实验室分析相结合,能够实现现场快速实时自动采集的报道较少;在采集设备的综合性能方面,一个要素对应一种采集设备的单一功能模式是主流,基于多源信息融合、能消除数据冗余、增强数据互补、集多种要素于一体的多功能采集设备,尚处于酝酿研究当中。

  纵观国内外农情信息采集技术的研究状况,近阶段重点在以下几个方面开展研究:

  1) 结合新的物理原理及学科移植的方法,把相关领域的新理论和新技术融合到农情信息的采集技术研究中,提高设备的采集精度和采集速度。

  2) 重点开展联合收获机械测产系统以及土壤环境信息快速测量仪器的研究与产品开发。

  3) 开发集多种测量要素于一体的多功能采集设备,以提高数据采集效率,降低数据采集的成本。

    注:
    (1)文章来源:农业工程学报,2006年 第22卷 第1期;
    (2)作者单位:华南农业大学工程学院。

(浏览次数:)
 
上篇文章:作物稳定性研究进展
下篇文章:中国杂交粳稻研究现状与对策
打印 | 关闭
相关文章