一种基于RS/GIS的冬虫夏草产量预测模型制造技术

一种基于RS(Remote Sensing，遥感)/GIS(Geographic Information System，地理信息系统)技术的冬虫夏草产量预测模型，解决了如何科学预测冬虫夏草这种野生中药材的年产量的问题。该模型基于RS(Remote Sensing，遥感)/GIS(Geographic Information System，地理信息系统)技术，利用雪线高程，采收期(4月和5月)的平均气温，降水量和日照时长等环境因子和虫草产量之间的关系，采用加权几何平均方法，在每年的6月上旬预测当年的虫草产量，经过检验，本模型的精度可达82.16％以上。此研究可为虫草产业的健康、可持续发展提供基础数据和信息。

A prediction model of Cordyceps sinensis yield based on RS/GIS

A method based on RS (Remote Sensing, remote sensing (Geographic) /GIS Information System, geographic information system) technology forecasting model of Cordyceps yield, how to solve the scientific prediction of wild Chinese herbal medicine Cordyceps annual production problems. This model is based on RS (Remote Sensing, remote sensing (Geographic) /GIS Information System, geographic information system) technology, the snowline elevation, harvest period (April and May) the relationship between the average temperature, precipitation and sunshine duration and other environmental factors and yield of Cordyceps, the weighted geometric mean method, in a year the early June forecast year Cordyceps production, after inspection, the model accuracy can reach more than 82.16%. This study can provide basic data and information for the health and sustainable development of the caterpillar fungus industry.

【技术实现步骤摘要】
一种基于RS/GIS的冬虫夏草产量预测模型
本专利技术属于野生药用植物资源量预测领域
技术介绍
目前世界上已知的广义虫草属有500多种，其中冬虫夏草是我国特有的名贵中药资源，是由冬虫夏草菌(Ophiocordycepssinensis)侵染鳞翅目蝙蝠蛾幼虫后发育而成的由子座和菌核组成的菌体。冬草夏草一般生长在海拔3000m以上的山地阴坡、半阴坡的灌丛和草甸中，在我国主要分布于青藏高原及边缘地区，零星分布于境外的尼泊尔、不丹的喜马拉雅山脉高寒草甸，我国冬虫夏草的采挖量占世界总采挖量的98％以上。由于天然冬虫夏草的生长有严格的寄生性及要求特殊的生态地理环境，野生资源十分有限。另一方面，冬虫夏草还不能进行人工或半人工的培育，只能自然生长，于是其蕴藏量在日益趋高的价格逼迫下迅速下降。据统计，20世纪50年代以前全国产量曾达到100吨以上，60年代初为50～80吨，90年代就仅有5～15吨。自70年代以来，因全球气候变化、生态环境破坏，以及掠夺式的采挖，加之其自身生长缓慢，自然资源更新能力低下，更使天然蕴藏量迅速减少。据80年代中期第三次全国中草药资源普查结果显示其野生蕴藏量仅约400吨，总体还呈逐年下降趋势。目前冬虫夏草天然资源已濒于灭绝，被列为国家二级重点保护野生植物。由于各种原因从第三次中药资源普查至今，没有对虫草进行过全面的调查，特别是虫草在西藏的具体分布和蕴藏量，有关部门只能从各乡群众采挖范围上获取一些基本情况，基本上是经验数据。因此，准确掌握虫草的产量信息是可持续开发冬虫夏草资源的基础，也是一个亟待解决的问题。
技术实现思路
本专利技术基于冬虫夏草的分布和气候环境变化之间的定量关系，选择雪线高程，采收期(4月和5月)的平均气温，降水量和日照时长等环境因子，采用加权几何平均方法，在每年的6月上旬预测当年的虫草产量。该模型包括如下步骤：步骤一：基于RS/GIS的气候环境因子构建。气候环境因子包括：虫草所在地区的5月份雪线高程，4月、5月、6月的降水量，4月、5月、6月的平均气温，4月、5月、6月的日照时长共10个因子。步骤二：建立气候环境因子与虫草产量的关系。利用相关分析工具，分别确定各因子与虫草产量之间的相关系数。根据经验，按照相关系数大小排序来确定入选模型的因子，各因子与虫草产量的相关系数，就是该因子的权重。步骤三：利用加权几何平均模型来构建虫草产量的预测模型。并利用历史数据进行模型精度检验。步骤四：预测当年虫草产量。在当年6月上旬，把当年的气候环境因子数据输入模型，得到当年的虫草产量预测值。上述的步骤一中，具体过程包括如下步骤：步骤一一：虫草所在地区的雪线高程的构建方法如下：利用NDSI指数(Normalizeddifferencesnowindex，归一化雪被指数)提取雪线[15]。NDSI是观测冰雪定量指标，积雪有很强的可见光反射和很强的短波红外吸收特性。积雪在0.5m附近有高反射率，在1.6m和2.1m附近反射率较低。计算公式是：NDSI＝(Ref0.555um-Ref1.640um)/(Ref0.555um+Ref1.640um)(1)公式(1)中Ref0.555um、Ref1.640um分别为0.555um和1.64um处的反射率，NDSI＞＝0.4的区域即为冰雪覆盖区。把提取出来的雪线高程的“平均值-2*标准差”作为雪线的最低高程。步骤一二：降水量、平均气温、日照时长相关因子的构建方法如下：根据虫草所在地区的经纬度，找到该区域内的气象台站，利用这些气象台站的平均值作为环境因子的值。所述的步骤一一中，计算雪线高程的具体过程包括如下步骤：步骤a：把DN值图像转化成反射率图像根据遥感影像的来源，有两种情况，一种是针对LandsatTM5/LandsatTM7影像，另一种是针对LandsatTM8影像。①对于LandsatTM5和LandsatTM7影像：0.555um和1.64um波段分别对应着影像的第2和第5波段。计算过程如下：第一步：分别计算各个波段每个像元的辐射亮度L值：公式(2)中，QCAL为某一像元的DN值(DigitalNumber，像元亮度值)，即QCAL＝DN。QCALmax为像元可以取的最大值，QCALmin为像元可以取的最小值。Lmax为该波段的最大辐射亮度，Lmin为该波段的最小辐射亮度，可以从遥感影像的头文件中获取以上参数的信息。第二步：计算各波段反射率(反照率、反射率)p公式(3)中，D为日地距离，可从影像的头文件中获取该参数信息。ESUN(SolarExoatmosphericSpectralIrradiances，或ESUNI)为大气顶层太阳辐照度可从遥感权威单位定期测定并公布的信息中获取，一般如表1所示：表1大气顶层太阳辐射照度(ESUNI)Table1SolarExoatmosphericSpectralIrradiances(ESUNI)(w/cm2·ster·μm)波段TM1TM2TM3TM4TM5TM7L5TM1957182915571047219.374.52L7ETM+1969184015511044225.782.07θ为太阳的天顶角，θ＝90°-β，β为太阳高度角，从遥感影像的头文件可获取太阳高度角信息。②对于LandsatTM8影像：0.555um和1.64um波段分别对应着影像的第3和第6波段。计算过程如下：第一步：分别计算各个波段的反射率ρλ′ρλ′＝Mρ*Qcal+Aρ(4)公式(4)中，ρλ′是未经太阳高度角校正的反射率，Mρ各波段的乘子，Aρ是各波段的修正数，ρ是波段数，Qcal是该波段的像元值。第二步：用太阳高度角进行修正ρλ′公式(5)中，ρλ是修正后的反射率，θSZ是太阳天顶角。θSE是太阳高度角，可从从遥感影像的头文件可获取。步骤b：在生成的反射率图像上计算NDSI对于LandsatTM5/7，NDSI＝(B2-B5)/(B2+B5)；对于LandsatTM8，NDSI＝(B3-B6)/(B3+B6)。式中B2，B5，B3，B7分别代表各波段反射率图像。步骤c：提取NDSI上大于等于0.4的区域根据MODIS产品说明，认为NDSI大于等于0.4的区域被积雪覆盖。利用ENVI软件里的bandmath工具来提取积雪覆盖的区域。步骤d：提取积雪区域的边界由于提取的积雪覆盖的区域是栅格图形，利用ENVI软件转化为矢量后是多边形，接着利用Arctools工具把多边形转化为线类型的矢量图形，至此雪线已提取完成。步骤e：利用空间分析提取雪线处最低高程利用Arcmap10软件，加载DEM(DigitalElevationModel，数字高程模型)数据，与雪线分布数据进行空间分析，提取雪线处的高程信息，最后去除图像边缘的异常值。上述的步骤二中，具体过程包括如下步骤：步骤二一：利用线性相关分析工具，分别计算各因子与虫草产量之间的相关系数R。步骤二二：按照相关系数从大到小的顺序排序，选择前4个因子作为模型的输入因子。各因子与虫草产量的相关系数，就是该因子的权重。上述的步骤三中，具体过程包括如下步骤：步骤三一：基于时间序列的虫草产量数据和各因子数据，建立虫草产量与各因子的单因子线性回归模型。步骤三二：利用加权几何平均模型来构建虫草产量的预测模型，本文档来自技高网...

【技术保护点】
一种基于RS(Remote Sensing，遥感)/GIS(Geographic Information System，地理信息系统)技术的冬虫夏草产量预测模型，其特征是基于RS(Remote Sensing，遥感)/GIS(Geographic Information System，地理信息系统)技术，利用雪线高程，采收期(4月和5月)的平均气温，降水量和日照时长等环境因子和虫草产量之间的关系，采用加权几何平均方法，在每年的6月上旬预测当年的虫草产量。该模型包括如下步骤：步骤一：基于RS/GIS的气候环境因子构建。气候环境因子包括：虫草所在地区的5月份雪线高程，4月、5月、6月的降水量，4月、5月、6月的平均气温，4月、5月、6月的日照时长共10个因子。步骤二：建立气候环境因子与虫草产量的关系。利用相关分析工具，分别确定各因子与虫草产量之间的相关系数，按照相关系数大小排序来确定入选模型的因子，各因子与虫草产量的相关系数，就是该因子的权重。步骤三：利用加权几何平均模型来构建虫草产量的预测模型。并利用历史数据进行模型精度检验。步骤四：预测当年虫草产量。在当年6月上旬，把当年的生态因子数据输入模型，得到当年的虫草产量预测值。...

【技术特征摘要】
1.一种基于RS(RemoteSensing，遥感)/GIS(GeographicInformationSystem，地理信息系统)技术的冬虫夏草产量预测模型，其特征是基于RS(RemoteSensing，遥感)/GIS(GeographicInformationSystem，地理信息系统)技术，利用雪线高程，采收期(4月和5月)的平均气温，降水量和日照时长等环境因子和虫草产量之间的关系，采用加权几何平均方法，在每年的6月上旬预测当年的虫草产量。该模型包括如下步骤：步骤一：基于RS/GIS的气候环境因子构建。气候环境因子包括：虫草所在地区的5月份雪线高程，4月、5月、6月的降水量，4月、5月、6月的平均气温，4月、5月、6月的日照时长共10个因子。步骤二：建立气候环境因子与虫草产量的关系。利用相关分析工具，分别确定各因子与虫草产量之间的相关系数，按照相关系数大小排序来确定入选模型的因子，各因子与虫草产量的相关系数，就是该因子的权重。步骤三：利用加权几何平均模型来构建虫草产量的预测模型。并利用历史数据进行模型精度检验。步骤四：预测当年虫草产量。在当年6月上旬，把当年的生态因子数据输入模型，得到当年的虫草产量预测值。2.根据权利要求1所述的一种基于RS/GIS技术的冬虫夏草产量预测模型，其特征在于，所述的步骤一中，具体过程包括如下步骤：步骤一一：虫草所在地区的雪线高程的构建方法如下：利用NDSI指数(Normalizeddifferencesnowindex，归一化雪被指数)提取雪线[15]。ND...

【专利技术属性】
技术研发人员：黄璐琦，朱寿东，郭兰萍，
申请(专利权)人：中国中医科学院中药研究所，
类型：发明
国别省市：北京,11