整合与物种分布模型遥感;映射红柳入侵使用软件进行辅助建模居（SAHM）

摘要

We demonstrate the utility of remotely sensed data and the newly developed Software for Assisted Habitat Modeling (SAHM) in predicting invasive species occurrence on the landscape. An ensemble of predictive models produced highly accurate maps of tamarisk (Tamarix spp.) invasion in Southeastern Colorado, USA when assessed with subsequent field validations.

摘要

入侵植物物种的早期检测是对自然资源和生态系统过程的保护管理是至关重要的。使用用于映射侵入植物分布卫星遥感的正变得越来越普遍，但是传统的成像软件和分类方法已被证明是不可靠的。在这项研究中，我们测试和评估采用沿阿肯色河东南部科罗拉多五个品种分布模型技术契合卫星遥感数据映射侵入柽柳（ 柽柳 ）。测试的模型包括升压回归树（BRT），随机森林（RF），多元自适应回归样条（MARS），广义线性模型（GLM）和Maxent模型。这些分析采用了新开发的软件包叫做软件辅助建模居（SAHM）进行。所有型号都带有499的存在点10000伪缺席点，并预测变量acqu受训来自八个月期间陆地卫星5专题制图仪（TM）传感器，红外发光二极管，从使用的物候差异检测原生河岸植被区分红柳。从陆地卫星的场景中，我们使用各个波段和计算的归一化植被指数（NDVI），土壤，调节植被指数（SAVI），并加盖流苏转换。确定对景观当前红柳分配所有五款车型成功的基础上独立的门槛和阈值相关的评价指标与位置无关的数据。为了说明模型的具体差异，我们生产的所有五款车型的合奏地图输出高亮和协议方面不确定性的领域。我们的研究结果表明物种分布模型的有效性在分析遥感数据和测绘合奏的效用，并展示在预处理SAHM的能力，执行多个复杂的模型。

引言

在整个美国西南部滨河湿地生态系统正由红柳的入侵的威胁（ 柽柳 ），欧亚大陆在1800年¹引入了非本土木本灌木。红柳具有许多生理机制，允许类利用水资源，外竞争本地物种，并改变生态系统过程^1-2。评估环境影响，并制定有效的控制策略映射红柳分布是为资源管理者高度优先事项。虽然地面调查仍经常使用，它们是不切实际的非常大的地区，由于人力，时间，和物流相关成本。

卫星遥感起到了红柳侵扰检测和映射一个重要的，但有限，角色。传统的分类分析和遥感软件有边际的成功^3-5。最近的几项研究具备检测利用遥感数据^1,6入侵植物探索非传统的方法。红柳，像许多外来入侵植物，展品在整个生长季节，从本地河岸物种的物候的不同物候的变化。在一些地区，例如，红柳叶出一些原生河岸植物面前，红柳保留其枝叶比其他本地物种更长的时间。通过使用光谱波段，并从整个生长季节时间序列卫星数据得出光谱指数，我们可以区分基于这些物候差异^1,6原生植物红柳。建立在Evangelista 等人的工作。 2009年¹月，在这项研究中，我们从时间序列的Landsat 5专题制图仪（TM）的卫星图像纳入各个波段1-7和衍生归一化植被指数（NDVI），土壤调节植被指数（SAVI）和缨帽从这些频带转换。归differen策植被指数（NDVI）是用于估计植被的生物量，篷盖，和叶面积指数^8-9中最常用的光谱指数之一，并且是可见的（红色）和近之间的比率的一个非线性变换红外波段^10。土壤调节植被指数（SAVI）是用来减少土壤背景对植被指数¹¹影响的修改NDVI。缨帽变换进行加权六个陆地卫星频带的复合材料成测量土壤亮度（缨帽，频带1），植被绿色（缨帽，频带2），和土壤/植被湿润（缨帽，带3）和三个正交频段经常被用来区分植被组合物，龄级，和结构^12-14。我们使用了克里斯特^{（1985）15日}报道^，为所有缨帽变换系数。

在这项研究中，我们测试五个品种分布模型与时间序列光谱波段和蔬菜从陆地卫星5 TM衍生etation指数映射沿下阿肯色河红柳在科罗拉多州东南部，美国。阿肯色河，跨越2364公里（1,469英里），是密苏里州的密西西比系的第二大支流。它涵盖分水岭在科罗拉多落基山脉的源头435123公里^2（168002英里^2）。从其在2965米起源，阿肯色州大幅下降的高程，拉平普韦布洛，CO附近，并通过农田和短草草原蜿蜒。这条河是受季节性洪水，并依靠在洛矶福特，拉詹塔，和拉马尔市政和农业用水，持续到堪萨斯州，俄克拉何马州，阿肯色州和它流入密西西比河之前。红柳首次在阿肯色河由R. Niedrach观察到1913年近¹⁶拉马尔的当今镇。今天，据估计，柽柳覆盖普韦布洛和堪萨斯州林间100多^平方公里即，与另外60 ^平方公里 ，沿阿肯色河¹⁷条支流。研究区域包括灌溉沟渠，湿地，农田等几个支流汇流;所有伴有不同程度的红柳为患。牧场和农业是主要的土地利用毗邻主要包括苜蓿，干草，玉米和小麦的河岸走廊。

物种分布模型依赖于地理参考事件（ 例如，纬度，经度）来确定一个物种的发生与环境¹⁸之间的关系。环境数据可以包括多个遥感和其他空间层。我们测试了五种分布模型包括升压回归树^{（BRT）19，}随机森林^{（RF）20，}多元自适应回归样条^{（MARS）21，}广义线性模型^{（GLM）22，}和Maxent模型^23。这五个MOD埃尔算法是最常用的物种分布建模之间，并且许多研究已证明其功效^24-25。我们使用的软件辅助建模居（SAHM）诉2.0模块来执行五款机型，其中包含在VisTrails v.2.2.2 ²⁶可视化和处理软件。有几个优点使用SAHM为比较建模。除了形式化和建模过程的记录温顺，SAHM允许用户与个别，有不同的接口，软件和文件格式²⁷多个物种分布模型的算法工作。 SAHM产生一致的门槛独立和阈值相关的评价指标，评价模型的性能。其中之一是地区根据ROC曲线（AUC），计算结果为背景，从²⁸歧视存在的模型的能力的门槛独立指标。一个AUC VAL0.5或更低UE表示模型的预测并不比随机更好或更坏; 0.5和0.70之间的值表示性能差;和值从0.70到1.0增加表明逐步更高的性能。另一个指标是百分之正确分类（PCC），即重灵敏度和特异性基于度量的用户定义的阈值的阈值依赖度量;灵敏度测量观察派驻人员划分为适宜和特异性的百分比衡量的背景地点列为不适宜的百分比。然而，另一个指标是真正的技术统计（TSS =灵敏度+特异- 1），这使更多的重量比特异性模型的灵敏度，值的范围-1到1，其中值> 0表示比²⁹的机会更好的模型性能之间。

使用模式输出映射红柳，我们构建二进制分类使用均衡的灵敏度和特异性来定义在p阈resence与否红柳。然后，这些个体模型推导出的地图进行相加，形成一个集成的地图^30。乐团地图结合个别物种分布模型的预测产生分类图，排列测试的模型的集体协议。例如，一个合奏细胞值表示只有一个型号划分该小区为合适的栖息地，而五个值指示所有五款车型分类的小区作为合适的栖息地。一个优点这种方法是合奏的地图产生较低的平均误差比任何单独的模型。它也允许用户直观地比较测试每个模型的性能。我们的总体目标是提供的这些方法可适合于物种对景观的电流分布进行建模的详细描述。

研究方案

1.野外数据采集

1. 导出字段的数据来自于2005年和2006年¹⁷全州调查由红柳联盟收集了矢量多边形数据红柳。
   注:通过密集的地面调查，在现场技术人员映射的所有红柳使用全球定位系统（GPS）和航拍照片沿着阿肯色河站获得的数据。
2. 这些多边形内，生成随机的499点（ 即派驻）来训练模型。放置两个约束集的随机点:（1）各点需要为"30米任何多边形边，以确保它是红柳内站在一个陆地卫星5 TM分辨率;和（2）所需的每个点可以从任何相邻点≥60米，以确保在整个研究站点健壮分布和最小化空间自相关。
3. 在三列的MS Excel文件编译现场数据与标题标有"反应"中，"X"，和"Y"，这里的响应值是:（1）对存在，UTM东向为X，以及为Y.保存该文件.csv格式用于SAHM（FieldData模块）UTM北向。
4. 生成红柳联盟的红柳多边形内的多边形作为独立的测试数据外额外的100随机点和另一个100个随机点来评估模型结果。以.csv格式保存该文件用于SAHM（FieldData模块）。

2.预测变量

1. 下载陆地卫星5表面反射L4-5 TM图片来自美国地质调查局的全球可视化浏览器/ E​​arthExplorer（路径32，排34）（http://earthexplorer.usgs.gov/）。这些场景包括由科罗拉多红柳联盟在2005年和2006年¹⁷取样的程度。为了确定用于模型的几个月里，选择场景是基本上无云（ 即 <10％的云层），并代表Source认为有区别的柽柳候个电子;这些都是2004年10月，2005年4月，2005年5月，2005年6月，2005年7月，2005年9月，2006年4月，2006年5月，2006年6月，2006年7月，2006年8月，2006年9月和2006年11月。
2. 请从遥感指数衍生工具
   https://github.com/rander38/Remote-Sensing-Indices-Derivation-Tool。
3. 运行该工具的Python脚本，无论是GDAL或ArcPy中的版本;建议GDAL。
4. 选择适当的卫星传感器，期望的指数，并设定输入的图像文件并在文件将被存储在输出文件夹（ 图1）。我们出口的各个频段和使用的NDVI，SAVI，和缨帽亮度，绿色，以及从每一个Landsat TM的场景衍生湿度指数。请注意，任何指标可以被修改或改变Sensors_Formulas_RSIDT.ini文件加入。

图1.遥感指数衍生工具的GUI。

5. 运行该工具，并在ArcMap v直观地验证输出文件。10.0（ESRI，雷德兰兹，CA）或其他地理信息系统软件。

3.软件辅助建模居（SAHM）（图2）

图2.整个SAHM工作流程包括输入数据，预处理，初步的模型分析和决策，相关模型，和输出过程。

1. 在https://my.usgs.gov/catalog/RAM/SAHM运行SAHM，来自美国地质调查局网站先下载文件（包括VisTrails）。请参阅用户手册在同一网站上下载和安装SAHM详细说明。请注意，该网站也有一个SAHM tutoriaL和附加说明的数据。
2. 为了开发红柳物种分布模型，使用随包下载SAHM_tutorial_2.0.vt文件（在SAHM例子文件夹）。在历史视图中，选择独立位置的工作流程。其它工作流的例子，可以选择并取决于研究的目标;描述提供与每个。选择管道 。
3. 设置输出文件夹，通过将包，然后SAHM> 更改会话文件夹 。整个工作流开发过程中，每一个步骤和选项的详细描述可通过选择位于SAHM查看器画面的右侧的文件标签被发现。在下列方法中列出的所有模块可在SAHM窗口的SAHM标签下的左侧找到。
4. 接着，直接SAHM到将成为字段数据用来训练的物种分布模型。
   1. 点击TemplateLayer模块。浏览到将被用作掩模，并以限定分析的投影，细胞的大小和程度的栅格。
   2. 点击在工作流的左侧FieldData模块。浏览到.csv（ 即 training.csv）字段数据模块中的现场数据（存在点或有无分）的文件。
   3. 点击PredictorListFile模块上浏览到.csv文件列表（ 例如 ，包含完整路径预测的所有文件，在使用建模参考用户手册）。
5. 接下来，进行预处理步骤。
   1. 单击FieldDataQuery模块上，并且响应柱与列标题填补反应（ 即 ，在FieldData.csv的列名）中，X和Y列。
   2. 点击MDSBuilder模块 。设置backgroundPointField 10000。
      注:如果您使用的品种存在和不存在数据被模仿，你并不需要更改backgroundPointField;你将包括在该领域Data.csv的响应（0）这些位置。它是可选的设置backgroundProbSurf如果希望通过指向值的范围从0到100的光栅表面，以限制的区域内的背景点选择（这些值代表了一个随机生成的将被保留下来应该它落入一个概率特定的细胞）。在这项研究中，使用具有100值的backgroundProbSurf一个5000米缓冲器阿肯色河和0为地区，这种缓冲液（基于由红柳联合取样的总面积）的外侧范围内。
6. 要使用，然后指定物种分布建模算法。
   1. 需要注意的是BoostedRegressionTree，GLM， 火星和镭ndomForest模块在独立的地点工作流程已经建立。在MAXENT模块添加到工作流程，测试所有五种型号。它连接到CovariateCorrelationAndSelection模块。
      注:为所有型号的默认设置开始;这些可以根据研究目的（详情请参阅模型文档）进行修改。
   2. 添加ModelOutputViewer模块，并将其连接到MAXENT模块;改变柱5和行为1。ModelOutputViewer产生一个可用于模型的结果进行比较的电子表格。
   3. 点击OutputName模块上，然后键入子文件夹名称。
7. 接下来，添加创建模型输出的合奏的模块。这个模块产生两个输出映射;愉快，包括输出的平均连续概率和与模型的数量的计数的第二带正仓元概率。
   1. 一个EnsembleBuilder模块添加到工作流程。这是可选的设置阈值度量;本次研究中，选择AUC和0.75设为默认的阈值。这将确保只有机型的AUC值大于或等于0.75包括在合奏地图输出。连接BoostedRegressionTree，GLM，MARS， 随机森林 和MAXENT模块的EnsembleBuilder。
8. 接着，直接模式，以独立的测试数据。
   1. 点击其他FieldData模块（工作流的右侧）并浏览到包含模型验证数据的.csv文件。这些都是在协议的步骤1.4中产生的200有无分。
   2. 点击FieldDataQuery模块上，确保响应，x和y列相匹配的列我n个字段数据.csv格式。
   3. 添加ApplyModel模块，并将其连接到MAXENT模块。添加ModelOutputViewer模块，并将其连接到这个ApplyModel模块;更改列5和行1.在菜单中，选择套餐 - > SAHM - > 更改处理模式 。既然你要运行一个以上的型号，选择单一机型顺序（N - 1每个内核）。这将通过以多个计算机的核心优势加速模型的执行时间。
9. 接下来，执行物种分布模型。
   1. 保存文件.vt，然后点击执行 。
   2. 当CovariateCorrelationAndSelection插件（ 图3）显示，去选择每个相关对变量，其中相关系数为一的| r | ≥0.7（基于％ 越轨从单变量干事解释源化此小部件（ 图3）和知情生态决定的左侧看到的添加剂模型;在这项研究中的优先级的至少一个协变量的选择每个月捕获红柳候变化）。在观察窗图解的数量可以通过在数字输入（默认为8）和点击更新来改变。
   3. 敲定协后，选择在CovariateCorrelationAndSelection小部件（ 图3）的底部确定。在这项研究中，以下9个变量被保留:July_30_2006_Brightness，June_09_2005_SAVI，Sept_16_2006_SAVI，May_24_2005_B4，Oct_28_2004_NDVI，April_22_2005_Brightness，April_09_2006_SAVI，Aug_31_2006_B4和Nov_19_2006_SAVI。该物种分布模型将选择确定后执行。

JPG"/>
图3.协变量的相关性和选择SAHM接口。

10. 结果输出。
    注:该机型完成后，VisTrails电子表格出现模型比较（ 图4）
    1. 比较AUC地块，文本输出，响应曲线，校准，混淆矩阵，残差和响应曲线模型之间。
       注:从SAHM模型输出包括CovariateCorrelationOutputMDS，混淆矩阵，残差图，校正曲线，模型评估的情节，变量重要性的情节，响应曲线的文件夹，扩大输出文件夹bin地图，乱七八糟的地图，MOD地图，概率地图，残图，和一个output.txt的;看到SAHM用户指南，详细的描述。

图4. VisTrails电子表格可以用来评估模型输出。这是训练数据的AUC模型比较;从左至右分别为车型是BRT，GLM，MARS，RF和Maxent模型。

结果

基于独立测试数据集BRT，RF，MARS，GLM和Maxent模型的统计评价显示所有五个模型检测红柳表现相对较好;有模型中的阈值的独立和阈值取决于评价指标相差不大。 AUC值均> 0.88，％的正确分类值均> 77％，敏感性和特异性均> 0.77和TSS均> 0.54（表1）。二元模型输出的合奏透露沿阿肯色河（图5）地区多模协议。残局（多元环境的相似性面）每个模型映射输出指示研究区的现有环境下以及取样（图6），进一步提高我们的集成方法的信心。

方法"> 模型 AUC PCC 灵敏度特异性 TSS BRT 0.91 85 0.85 0.85 0.70 RF 0.92 85 0.85 0.85 0.70 火星 0.90 82 0.82 0.82 0.64 GLM 0.88 77 0.77 0.77 0.54 MAXENT 0.92 84 0.83 0.84 0.67

表1。独立的阈值（AUC）和阈值依赖（PCC，敏感性，特异性和TSS）评价指标BRT，RF，MARS，GLM和Maxent模型模型适合独立的测试红柳存在和不存在的数据集。

图5.合奏结果在ArcGIS结合BRT，GLM，MARS，RF和Maxent模型二进制输出地图。区域通过的模型数目在协议着色，从0（无颜色）至5（红色）。注意在预测西北角的彩色区域;这条线是陆地卫星图像的假象;因此，模型结果应该在这个区域谨慎服用。 请点击此处查看该图的放大版本。

图6.多元环境相似表面（MESS）输出。目标="_空白">点击此处查看该图的放大版本。

使用九个预测的，2006年6月30日的亮度是所有五个模型（表2）的最重要的变量。这是基于逐步赤池信息标准由GLM保留的唯一变量（AIC;这是在SAHM GLM模型选择默认值），不过需要注意的这款机型还包括这个变量的平方项是很重要的。 RF和Maxent模型默认保留所有变量。

预报器	BRT	RF	火星	GLM	MAXENT
2006年7月30日亮度	41.60	34.11	76.78	100	67.27
2006年8月31日四级	6.35	5.87	5.16	0	2.82
2005年6月9日SAVI	13.67	14.09	9.14	0	9.75
2005年4月22日亮度	6.29	6.30	0	0	0.43
2004年10月28日NDVI	5.66	8.25	0	0	2.94

表2.每个模型预测的相对重要性。

讨论

我们的结果表明拟合BRT，RF，MARS，GLM和Maxent模型与红柳存在点和时间序列遥感的Landsat卫星影像数据可在景观区分红柳，是一种有效的替代传统的单一场景分类方法。这是从我们的结果六月是我们研究范围内检测红柳一个特别重要的时间明确;此同意与2009年¹斯塔等。这表明六月旱涝是红柳发生在这方面基于一个Maxent模型模型拟合以时间序列陆地卫星图像的最重要的预测。

被纳入BRT，RF，火星和Maxent模型的模型，其他光谱指数和乐队可以从土壤基质，其他落叶乔木，包括三叶（ 杨树 ）和柳树进一步区分红柳（ 柳树 ），或灌溉农业的是在较低的共同阿肯色河流域。其他GIS图层，如地形，土壤类型和气候数据也被视为协变量和包含在这些模型中，但我们建议保持这些到最低限度，如果目的是检测当前的物种分布的景观，而不是预测的潜力发生或适宜的栖息地。

为我们的研究测试的模型提供了强有力的分析能力以及评估结果的多个选项。让所有在单一框架内，这些相关模型，如SAHM的，允许正规化和建模过程的记录温顺。前和响应和预测变量的后处理被标准化在SAHM，允许更好和更有效的模型进行比较，而工作流记录的分析提供便利的修改，迭代和复制的每一个步骤。

合奏映射的目的是几个相关模型的优势结合起来，同时尽量减少在W任何一个模型³⁰ eakness。我们认为，这是在我们的研究的情况;然而，我们告诫说是表现不佳（ 即下预测或过度预测）模型可以削弱整体效果。在文献中的有限使用合奏映射已有利的结果，但大多数这些方法试图"预测"物种发生，而不是"检测"。此外，乐团映射允许不确定性的不同的建模方法之间的视觉评估，确定示范协议的水平。最常见的是建模方法的选择（ 例如 ，GLM对BRT）是对模型的结果，而不是在建模过程中其他决定，如位置数据的不确定性最大的³¹量化的影响。虽然我们相信，我们最好的红柳地图是所有五种模式是一致的，进一步的测试和使用Ensemble映射的各种方法建议（ 例如 ，通过AUC加权^{）32 <}/ SUP>，并通过独立的野外观测的最佳验证。总之，这些方法可以很容易地适合于使用衍生于SAHM一个给定的研究区域的环境变量其他物种的分布进行建模。

材料

Name	Company	Catalog Number	Comments
Earth Explorer	USGS	http://earthexplorer.usgs.gov	Open Access: Yes
Remote Sensing Indices Derivation Tool	github	https://github.com/rander38/Remote-Sensing-Indices-Derivation-Tool	Open Access: Yes
Software for Assisted Habitat Modeling	USGS	https://my.usgs.gov/catalog/RAM/SAHM	Open Access: Yes
ArcGIS v.10.3	Esri	https://www.arcgis.com/features/	Open Access: No