基于MCR和HY模型的莆田市生态廊道构建与优化.pdf

引文格式：金博闻，谢祯，柯I.基于MCR和HY模型的莆田市生态廊道构建与优化 .西南林业大学学报（自然科Jan.2024JOURNAL OF SOUTHWEST FORESTRY UNIVERSITY2024年1月Vol.44No.1大南第44卷西第1期报学学业林DOI:10.11929/j.swfu.202210015学）,2 0 2 4,44(1):7 9-8 7.基于MCR和HY模型的莆田市生态廊道构建与优化金博闻1，2谢祯？柯山1耿建伟2 潘辉！（1.闽江学院，福建福州35 0 10 8；2.福建农林大学风景园林与艺术学院，福建福州35 0 10 6）摘要：以莆田市主城区为例，采用形态学空间格局方法（MSPA）结合景观连通性分析筛选重要生态源地，运用最小累积阻力模型（MCR）与重力模型方法，提取研究区潜在生态廊道；同时结合水文分析模型（HY）识别研究区辐射廊道和生态节点，最终构建并优化莆田市生态廊道。结果表明：基于MSPA与景观连通性指数，共识别出14个核心源地；采用MCR与重力模型，共识别出2 9 条廊道，一级廊道7 条，二级廊道10 条，三级廊道12 条；结合水文分析法共识别出2 0 5 条辐射廊道和5 3个生态功能节点，包括一级节点5 个，二级节点17 个，三级节点31个，确定各类型生态修复的关键区域，提取生态节点建设策略。本研究可为莆田市生态廊道构建及生态节点保护时序性研究提供科学依据。关键词：生态廊道；最小累积阻力模型；水文分析；辐射廊道中图分类号：X171.1文献标志码：A文章编号：2 0 9 5-19 14(2 0 2 4)0 1-0 0 7 9-0 9Construction and Optimization of Ecological Corridor Based onMCR and HY Analysis in PutianJin Bowen 2,Xie Zhen2,Ke Shan,Geng Jianwei?,Pan Hui(1.Minjiang University,Fuzhou Fujian 350108,China;2.Fujian Agriculture and Forestry University,Fuzhou Fujian 350106,China)Abstract:In this study,taking the main urban area of Putian City as an example,the morphological spatialpattern approach(MSPA)combined with landscape connectivity analysis is used to screen important ecologicalsource sites,and the minimum cumulative resistance model(MCR)and gravity model approach are applied to ex-tract potential ecological corridors in the study area,while the hydrological analysis model(HY)is combined toidentify radiation corridors and ecological function nodes in the study area,and finally construct and optimize eco-logical corridors in Putian City.The results show that based on MSPA and landscape connectivity index,a total of14 core source sites were identified;using MCR and gravity model,a total of 29 corridors were identified,includ-ing 7 primary corridors,10 secondary corridors and 12 tertiary corridors;combined with the hydrological analysismethod,a total of 205 radiating corridors and 53 ecological function nodes were identified,including 5 primarynodes,17 secondary nodes and 31 tertiary nodes.key areas for various types of ecological restoration were identi-fied.Ecological node construction strategies were extracted.The results of this study provide a scientific basis forthe study of ecological corridor construction and ecological node protection temporality in Putian City.收稿日期：2 0 2 2-10-0 9；修回日期：2 0 2 2-12-0 6基金项目：福州市社科规划重大项目（2 0 17 FZA02）资助。第1作者：金博闻（19 9 7 一），男，硕士研究生。研究方向：风景园林规划与设计、旅游地理。Email:。通信作者：潘辉（19 6 8 一），男，博士，教授，享受国务院政府特殊津贴专家。研究方向：环境资源管理、生态旅游等。Email:。80西南林业大学学报第44卷Key words:ecological corridor;minimum cumulative resistance model;hydrological analysis;radiation cor-ridor随着城镇化进程的快速推进，高密度的土地开发不断侵占城市中的生态用地，导致了生态用地的斑块侵蚀、连通性降低和景观破碎化，对城市生态安全与整体景观格局造成严重威胁与破坏 ，产生了如城市热岛效应 2 、洪涝 3 和生物多样性丧失 4 等系列问题。研究表明，生态斑块破碎化是造成区域生物多样性减少的症结之一。作为连接城市主要生态斑块的线性要素，生态廊道通过在生态规划中的布局、整合和衔接，利用其线性特征连接和组织生态源地，不仅有利于促进城市生物流的良性循环，也对城市土地价值提升和城市可持续发展具有重要意义。然而，主城区具有人为干扰大，土地覆被复杂等特点，导致生态源地较为破碎，对生态廊道的构建造成了诸多困难。因此，如何在高人口密度的主城区，构建并优化生态廊道显得尤为重要。目前，生态廊道已形成了“源地-阻力面-廊道”的研究范式 5 ，其中，如何准确地识别生态廊道，成为当今生态学及风景园林学等关联学科讨论研究的重点。当前对于生态廊道的研究主要集中在连通性阈值设置、阻力面构建、廊道的模拟等方面，涉及渗透理论、图论、阻力模型理论、电流理论等。杨娅楠等 6 基于渗透理论并结合模型模拟等方法，讨论了当前在研究景观连通性过程中存在的问题，侧重于合理构建廊道对于野生动物生境保护的重要性。Urban等 7 首次将图论理论应用于解决有关景观连通性的实际问题，重点为保护墨西哥斑点猫头鹰种群及其栖息地提供了方法。阻力模型理论认为生物体在穿过不同景观空间的过程中需要耗费一定的成本或克服相应阻力，物种最佳的生境通道则为累积成本最小或阻力最小的生境通道 8 ，其操作过程复杂且易生成杂廊道。近年来又发展出许多以阻力模型理论为基础，计算最佳迁移路径的方法和廊道构建模型，例如最小成本路径模型（LCP）、最小累积阻力模型（MCR）9-1 等。最小累积阻力模型（MCR）运用ArcGIS的成本路径工具识别源地间最小耗费路径，确定物种迁徙与扩散的所有可能路径，忽略了源地与周围基质的连接和连接潜在源地等属性，以及斑块内部的物种运动路径，对生态廊道与节点的识别具有局限性12 。而水文分析模型（HY）能有效的识别辐射廊道，即在低阻力区域形成辐射型谷线 3-15 ，能够有效补充区域内潜在生态廊道；MCR与水文分析模型（H Y）结合所产生的潜在廊道与辐射廊道的交点，能有效补充斑块内部的生态节点。因此，水文分析模型在一定程度上弥补了MCR模型对于潜在廊道提取与生态节点选择中的不足，对于加强源地与周围基质的连接以及连接潜在源地具有重要作用 12 基于此，本研究以莆田市主城区为例，采用形态学空间格局（MSPA）结合景观连通性识别、分析、筛选生态源地，运用最小累积阻力模型（MCR）与重力模型，构建研究区各级潜在生态廊道；同时结合水文分析模型（HY），提取研究区辐射廊道，并识别出潜在廊道与辐射廊道的交点为生态节点，最终构建并优化莆田市生态廊道。本研究可为莆田市生态廊道构建提供科学依据。1研究区概况莆田市是福建省辖地级市，位于我国华东地区（2 45 9 2 5 46 N，118 2 7 119 6 E），东临台湾海峡，背倚戴云山脉。境内地势呈西北高、东南低，西北部中低山众多，海拔约5 0 0 18 0 0 m。中部为低山与丘陵，平原与盆谷错综其间，地势较平坦，是研究区中人口主要分布的建成区。本研究划定莆田市主城区范围：荔城区、城厢区、涵江区、秀屿区。由于莆田市主城区内人口密集，城市化程度较高，造成内部生态源地数量多，面积小，与城郊的大型生态源地联系性差。因此，合理地构建莆田市生态廊道，对提升城市生态环境，促进城市可持续发展具有重要的意义。2材料与方法2.1数据来源与预处理本研究所采用的数据为：1）莆田市10 m81第1期金博闻等：基于MCR和HY模型的莆田市生态廊道构建与优化10m精度的土地利用栅格数据【欧洲航天局（ESA）的全球陆地覆盖数据集（https:/www.esa-landcover-cci.org），根据ESA数据分类结果，在研究区内提取8 类土地利用类型：林地、灌木地、草地、水体、湿地、耕地、裸地和建设用地；2）30 m30m的DEM据【地理空间数据云（h t t p:/w w w.g s c l o u d.c n）l；3）2 0 2 1年莆田市行政边界数据（莆田市自然资源和规划局）；4）道路矢量数据（https:/www.openstreetmap.org/），提取莆田市的铁路、国道和高速公路3类道路信息。2.2重要生态源地识别与评价形态学空间格局分析（MSPA）是基于数学形态学原理来识别、测量和分析栅格影像空间的图像处理方法 16 ，能够高效识别研究区范围内的景观类型。本研究通过将土地利用数据中具有生态价值的要素（前景）和其他要素（背景）进行二值化处理后，以林地、灌木地、草地、水体、湿地要素作为前景数据，借鉴相关学者的参数设定 17 ，景观类型划分参考GuidosToolbox手册 18 ,将前景数据分为7 种景观类型。最后对分析结果进行面积、占比等数据统计，并提取核心区斑块作为连通性分析的生态源地 19-2 0 。其次，利用连通性指数（PC）与整体连通性指数（IC），定量评价斑块之间的连通性。数值越大说明斑块的连通性越高，对物种的交流和迁徙越有利，说明斑块具有更高的生态功能价值。相关计算公式如下：a:a,Pij4j1PC,(0PC 1)(1)A2a;ajIC(0IIC 1)(2)式中：PC为可能连通性指数；IIC为整体连通性指数；AL为景观总面积；n为斑块的总量；a;与a,为斑块i与j的面积；P,为斑块在i与j之间连通的最大概率，nl，指斑块在i与j之间的连接数量。PC-PC,removedPC100%(3)PCIC-IIC,removedIIC100%(4)IC式中：dPC与dIIC值越高，代表斑块的连通性越好，作用越大；PCremove与IICremove为剔除某一斑块后连接指数值。本研究综合相关研究 2 1-2 2 并结合相关领域专家建议，分别赋予dPC、d I C为0.5 的权重值，结果大于0.0 2 的生态斑块划定为核心源地。2.3基于MCR模型的生态网络构建2.3.1综合阻力面构建物种的空间传播和迁移可以促进生物多样性，有利于生态系统的稳定；然而，迁移的难易程度受到空间阻力的影响，因此合理地构建综合阻力面是提取生态廊道的关键。其中，土地利用类型是廊道构建的基础，坡度和高程影响植被生长和城镇开发强度 13，而与铁路、公路 16 之间的距离可以反映人类活动的影响程度。综上所述，本研究选取土地利用类型、MSPA景观类型、高程、坡度、距离国道16 、距离水域、距离铁路、距离高速公路作为生态廊道提取的主要阻力因子。结合研究区实际情况 2 3 和相关研究成果 19.2 4,并通过咨询相关专家和AHP层次分析法赋予莆田市的阻力因子相应权重和阻力值，借助ArcGIS10.8栅格加权叠加等方法构建综合阻力面（表1）。表1阻力层权重及因子阻力值Table1Resistance layer weights and factor resistance values阻力类型阻力因子阻力值权重阻力类型阻力因子阻力值权重土地利用类型林地10.350MSPA景观类型核心区10.221灌木1桥接1草地1环道3耕地5支线5建设用地9孤岛5裸地7边缘7水体3孔隙7湿地3背景9高程/m(80,200)10.112距离水域/m0,10010.070200,4003100,2003400,6005200,500582西南林业大学学报第44卷续表1阻力类型阻力因子阻力值权重阻力类型阻力因子阻力值权重高程/m600,80070.112距离水域/m500,100070.070800,+o01000,+09坡度/)0,810.112距离铁路/m2000,+o010.0450,831500,2000315,2551000,1500525,357500,1000735,900,500距离国道/m2000,+8010.045距离高速公路/m2000,+0010.0451500,200031500,200031000,150051000,15005500,10007500,100070,50090,50092.3.2生态廊道的提取物种的迁徙受多种阻力的影响，结合研究区生态环境特征，以Yu等 2 5-2 6 最小累积阻力模型（M CR）相关研究为依据，基于ArcGIS10.8中的成本路径工具，以构建的综合阻力面为基础计算生态斑块之间最小成本路径，得到潜在生态廊道。公式如下：MCR=fmZ(D,R.)(5)式中：f是模型与变量间的正函数关系；D；是生态流从i到i所消耗的空间距离；R，是加权阻力值。2.3.3生态廊道的分级重力模型又称为引力模型 2 7 ，该方法的关键在于能够将不同斑块之间的相互影响程度计量化，常被应用于空间相互作用的研究中。斑块之间的相互影响程度越高，表明斑块之间联系越紧密，以此作为筛选相对重要廊道 2 8 的关键依据。重力模型的计算公式如下：N.N,LmaxlnV,lnV,F=GabD2LaP.P,ab(6)式中：F为重力模型式；Gab为生态斑块a与b间的相互作用力；Na、N,为斑块a和b的权重；Dab为与b之间潜在廊道阻力的标准化值；Pa、P,是斑块a、b 的阻力值；Va、Vb 是斑块a、b 的面积；Lab是斑块与b之间廊道的累加阻力值；Lmax是研究区内所有廊道累积的最大阻力值。2.3.4基于水文分析的生态廊道优化基于MCR与重力模型所提取的仅为斑块间的廊道，只能定量描述物种在斑块间的迁徙最大可能路径，而忽视了物种在斑块内部的运动路径，而通过水文分析模型所提取辐射廊道 13 能够有效弥补前者的不足。水文分析模型（HY）常用于提取河流网络、划分流域等水文研究中 2 8-2 9 。本研究采用ArcGIS10.8中水文分析工具，提取综合阻力面的谷线，进而构成低阻力值网络，再通过河网分级工具中STRAHLER方法 15 识别潜在的辐射廊道。利用ArcGIS10.8的距离测算工具分析得到研究区累积耗费距离表面，以此为基础，结合研究区实际情况，利用填洼、流向、流量等工具，通过反复测试阈值，确定汇流累积的值，补充识别研究区的辐射廊道。结合研究区现状测试与相关研究 13,15 ，设定大于15 0 0 的阈值范围，得到物种迁移的低阻力线路，最终确定研究区辐射廊道的分布位置。借鉴生态节点的相关研究 15 ,在最小耗费路径即生态廊道与辐射廊道的交叉处提取生态节点 2 9 3结果与分析3.1生态源地识别将二值化后的前景与背景数据导入Guid-osToolbox工具，分析得到MSPA分类（图1），计算后得出各景观类型的面积与占比数据（表2）。结果表明，前景要素中面积最大的为核心区（9 6 3.32 k m），景观类型占比为7 8.7 3%，景观类型主要以林地为主，其次为水体、草地等，最小的为环道（2 6.5 3km），仅占2.17%。从空间分布来看，核心区主要分布于莆田市北部戴云山余脉和西部林地。其中，北部核心区呈块状组团，整体连通性良好，包括瑞云山森林公园、壶公山、囊山等；西部呈条块状，连通性较好，包括九龙谷国家森林公园、东圳水库和木兰溪上游；而中心城区内的核心区数量上较多且面积小，破碎化程度高，连通性差，不利于物种的迁徒扩散和物质循环。83第1期金博闻等：基于MCR和HY模型的莆田市生态廊道构建与优化N核心区孤岛区孔隙边缘环道景051020km图1MSPA景观类型Fig.1MSPAlandscapetype基于MSPA分析结果，提取面积前15 的核心区进行连通性分析，采用Conefor2.6设置斑块间的连通距离阈值为12 0 0 m，连通率为0.5 30 。通过各斑块的整体连通性指数（dIC）、可能连通性指数（dPC）计算结果，选择dIIC、d PC均值大于0.0 2 的生态斑块作为核心源地，总源地数量为14个。由于14号斑块面积过大会导致数据偏差，故计算其余的13个源地的连通性指数（dIIC、dPC）。由图2、表3分析可知，城郊西北部戴云山的大型山体与西部壶公山构成研究区生态源地的主体，面积较大。其余城区内的源地面积小且分散，主要包括城厢区的凤凰山公园、绶溪公园、秀屿区的土海湿地公园以及木兰溪流域等。主要由于城区往沿海地区的土地资源稀缺、人类活动密集以及沿海地区工业用地规模的扩张和侵蚀。表2 MSPA景观类型统计Table2Statistics on MSPA landscape type林地灌木地草地水体湿地景观类型面积/km占比/%面积/km?占比/%面积/km?占比/%面积/km占比/%面积/km2占比/%面积/km2占比/%核心区963.3278.73905.1293.966.750.0116.651.7341.374.290.120.01孤岛39.103.2026.9668.9512.870.3310.8527.741.012.590.160.40孔隙56.944.6550.9989.560.450.014.167.311.763.090.020.04边缘67.715.5350.0873.976.570.106.729.9310.6715.760.160.24环道26.532.1721.5681.281.440.053.6713.851.164.360.120.46桥接33.552.7426.2778.303.690.114.8414.442.276.770.130.38支线36.462.9826.2972.135.040.147.3720.232.526.920.210.58总和1223.59100.001107.2790.4936.810.0354.274.4460.764.970.920.07N研究区边界核心源地2051020km图2核心源地分布图Fig.2Core source distribution map表3核心生态源地连通性指数Table3Connectivity index of core ecological source斑块编号斑块重要性指数1098.74611.9292.0641.4420.8530.62130.21120.0670.05110.0510.0580.0350.0284西南林业大学学报第44卷3.2生态网络构建分析3.2.1综合阻力面分析基于ArcGIS10.8，对各个阻力因子进行加权叠加计算，生成综合阻力面栅格图（图3）。综合阻力面高：7.47 2：1051020km图3综合阻力面Fig.3Comprehensiveresistance surface由图3分析可知，阻力值较大的区域主要位于莆田市中心区域与东部沿海工业区，这些地区建筑、人口和产业密集，干扰较大，不利于生物的迁徙和交流。其中木兰溪、雁阵山、汀港山等中小型生态源地，维持着一定的生态服务功能和区域内生态系统平衡，促进了莆田市南北物种的交流，但明显数量不足且较为分散。3.2.2关键生态廊道提取分析重力模型能高效筛除MCR模型提取出的余廊道，保留下的2 9 条廊道在各生态源地间均有直接或间接联系。计算各核心源地之间的相互作用力，参考相关研究的阅值设置 31，提取作用力值大于5 0 0 0 的廊道定义为一级廊道，共计7 条；值大于10 0 0 而小于5 0 0 0 的定义为二级廊道，共计10 条；值大于5 0 0 小于10 0 0 的为三级廊道，共计12 条，构建斑块相互作用矩阵（表4）和莆田市生态廊道（图4）。如图4所示，一级廊道和二级廊道构成大型生态源地与中小型源地物质循环、能量流动的重要载体，主要分布在外围城郊。三级廊道主要为促进城市中小型生态源地的内部联系。中心城区阻力值较高，一级和二级廊道数目较少，三级廊道则多沿着水系分布。一级和二级廊道主要依托城郊山地的自然资源，大型面状水域或绿地，而三级廊道大多经过木兰溪、延寿溪、芦溪等带状水域。表4斑块相互作用矩阵Table 4Patch interaction matrix斑块序号1234567891011121314105012493467862409231472323951352559236420278210199219181732487172780878285410259830113511411112819265804288562502844864018248112587814608453747409955850150261627418624924427512414660202709633157675929739570585279366335588501156807210028535131906490652968360112523510094172714344811012385611205102221306614085第1期金博闻等：基于MCR和HY模型的莆田市生态廊道构建与优化N研究区边界核心源地级廊道级廊道三级廊道12051020km图4基于MCR模型提取结果Fig.4ExtractionresultsbasedonMCRmodel3.2.3廊道优化与节点的识别辐射廊道作为低阻力区域形成的辐射型谷线 13，对于加强生态源地与周围基质的物质能量交流具有重要作用。本研究基于水文分析法共提取2 0 5 条辐射廊道（图5），辐射廊道分布多呈西北至东南走向，分布较为均匀。由图4可知，基于MCR与重力模型提取的生态廊道分布受核心生态源地空间分布的影响主要集中于研究区西侧，分布不均。而辐射廊道的构建，一方面加强核心区于周围基质的物质能量交流，另一方面可缓解核心区空间分布不均的缺点，在一定程度上增强了研究区西部与中心城区的生态稳定，为潜在生态源地的连接提供了可能。基于MCR与重力模型所提取的各级潜在生态廊道，借鉴相关研究 13,2 9，识别辐射廊道与各级廊道的交点作为生态节点，共提取5 3个生态节点，其中一级节点5 个，二级节点17 个，三级节点31个，得到优化后的莆田市生态廊道（图6）。生态节点保护时序上，应将一、二级生态节点划分为优先保护区域，严格控制对节点区域的生态侵扰，在保护基本农田的基础上，可采取适当的退耕还林、退耕还草等方式，合理规划建设用地数量，保证公园绿地建设；建立完善的监测体系，有效了解修复过程存在的问题，及时调整修复策略。空间分布上，生态节点在生态廊道中的分布较为均匀，对生态廊道的稳定连接起到一定程度的促进作用。N研究区边界核心源地辐射廊道051020km图5基于HY模型提取结果Fig.5Extractionresults basedonHYmodelN研究区边界核心源地级廊道二级廊道三级廊道辐射廊道一级节点二级节点三级节点051020km图6莆田市生态廊道Fig.6Ecological corridor in Putian City4结论与讨论4.1结论基于MSPA、C o n f e r 工具评价识别核心源地，利用MCR及重力模型提取潜在生态廊道，并运用水文分析法优化研究区生态廊道。通过上述研究得出以下结论：1）本研究基于MSPA识别出莆田市7 种景观类型，核心区面积为96 3.32 km，景观类型占比为7 8.7 3%；综合选取dPC与dIIC连通性指数和斑块面积定量筛选出14个生态源地。主要分布在86西南林业大学学报第44卷研究区西部、北部中低山丘陵，中部平原与东南沿海区域分布稀疏，生态源地连通性较差。2）在阻力因子选取方面，依据莆田市土地利用类型，综合8 类要素作为主要阻力因子，进行栅格叠加分析后得到综合阻力面。空间布局上看，阻力值分布呈西北低东南高，分布阻力值较大的区域主要位于莆田市主区与东部沿海工业区，这些地区建筑、人口和产业密集，不利于物种的迁徙。3）基于MCR与重力模型提取2 9 条生态廊道，其中一级7 条，二级10 条，三级12 条。一、二级廊道主要集中在莆田市西部较高阻力值汇集的区域，需在这些区域加强生态建设，降低阻力值以促进生态廊道的连通性。而三级廊道大多途经水系，因此区域内水源涵养是未来工作的重点。4）基于水文分析模型，共提取辐射廊道205条，生态节点5 3个，一级节点5 个，二级节点17 个，三级节点31个。该方法为物种迁移及保护提供了一定补充，辐射廊道的轨迹能够提升生态廊道的合理延展性，促进区域内生态系统的优化组合和生态资源的保护。而生态节点有利于加强源地与周围基质的连接和潜在源地的连接，节点等级的划分对于莆田市主城区未来规划建设时序优先级判定具有重要意义。4.2讨论本研究旨在利用水文分析法，对传统MCR与重力模型的廊道识别与筛选方法进行优化，在前人研究基础上进一步提取辐射廊道与生态节点，补充对生态廊道以及辐射廊道研究的不足，提升研究区生态廊道的稳定性和整体性。在生态源地的识别中，MSPA对于研究区域像元大小十分敏感，不同粒度、边缘宽度对生态廊道构建的影响尚未纳人考虑范围之内；在Con-efor2.6连通性计算分析的过程中，参考操作手册设定与前人的研究结果 30-31，设置宽度阅值为12 0 0,连通概率为0.5，用以评价源地的重要程度，宽度阈值的大小对连通性评价结果影响的大小还可以进一步研究。阻力面构建是生态廊道的关键步骤，本次研究借鉴已有的研究，选取的阻力因子侧重陆生物种生态廊道构建，未考虑到水生物种对生态网络的敏感性。当前，阻力因子赋值主要来源于专家经验数据，本研究采用AHP层次分析法与德尔菲法主观性较强，阻力值的合理设置则要考虑多方面的因素，包括研究目的、景观格局特征等。已有研究表明，结合研究目的，建立多套阻力赋值，对研究区做阻力值影响分析 32 ，筛选出合适的综合阻力。由于莆田市主城区综合阻力值较大，生态廊道及源地难以通过与提取。本研究在高密度的主城区补充构建了辐射廊道及生态节点，优化了生态廊道的组合形式，但对高密度人口特征区域的筛选仍存在不足。因此，在今后研究中，如何准确、合理地提取人类活动影响较大的城区生态源地及如何构建高阻力城市中心区域的生态廊道，是未来应深入探讨的关键问题。参考文献1邓金杰，陈柳新，杨成，等.高度城市化地区生态廊道重要性评价探索：以深圳为例 .地理研究，2 0 17,36(3):573582.2裔传祥，胡继超，李小军.土地覆盖类型对城市热岛效应的响应分析 J.测绘通报,2 0 18(1):7 2-7 6,116.3吴豪，许刚，虞孝感.关于建立长江流域生态安全体系的初步探讨 .地域研究与开发，2 0 0 1,2 0(2):34-37.4武晶，刘志民.生境破碎化对生物多样性的影响研究综述 .生态学杂志,2 0 14,33(7):1946-195 2.5刘壮壮，吴未，刘文锋，等.基于“源地-廊道”生态安全格局构建逻辑范式的建设用地减量化研究 .生态学报,2 0 2 0,40(2 2):8 2 30-8 2 38.6杨娅楠，王金亮，习晓环，等.野生动物生境景观连通性综述 .环境科学导刊,2 0 15,34(3):10-15.7Urban D,Keitt T.Landscape connectivity:a graph-the-oreticperspective.Ecology,2001,82(5):1205-1218.8LaRue M A,Nielsen C K.Modelling potential dispers-al corridors for cougars in Midwestern North Americausing least-cost path methods J.Ecological Model-ling,2008,212(3/4):372-381.9单楠，周可新，潘扬，等.生物多样性保护廊道构建方法研究进展 .生态学报，2 0 19,39(2)：411-42 0.10沈泽昊.景观生态学的实验研究方法综述 .生态学报,2 0 0 4,2 4(4):7 6 9-7 7 4.11张美丽，齐跃普，张利，等.基于LinkageMapper与粒度反推法的太行山中北段生态节点识别与分析：以河北省阜平县为例 .生态与农村环境学报，2 0 2 0,36(12):1569-1578.12王越，赵雯琳，刘纯青.基于MSPA-Conefor-MCR路陆驰）（责任编辑87第1期金博闻等：基于MCR和HY模型的莆田市生态廊道构建与优化径的生态网络优化及其构建：以彭泽县为例 .江西农业大学学报，2 0 2 2,44(2)：5 0 4-5 18.13韦宝婧，苏杰，胡希军，等.基于“HY-LM”的生态廊道与生态节点综合识别研究 .生态学报，2 0 2 2,42(7):2995-3009.14王原，何成，刘荣国，等.宁夏沙坡头国家自然保护区鸟类景观生态安全格局构建 .生态学报，2 0 17,37(16):5531-5541.15朱军，李益敏，余艳红.基于GIS的高原湖泊流域生态安全格局构建及优化研究：以星云湖流域为例 .长江流域资源与环境，2 0 17,2 6 8)：12 37-12 5 0.16Soille P.Morphological segmentation of binary pat-terns.Pattern Recognition Letters,2009,30(4):456-459.17崔凤奎，王晓强，张丰收，等.二值图像细化算法的比较与改进 .洛阳工学院学报，1997，18(4):48-5 2.18Vogt P,Ritters K.GuidosToolbox:universal digitalimage object analysis J.European Journal of RemoteSensing,2017,50(1):352-361.19吴健生，张理卿，彭建，等.深圳市景观生态安全格局源地综合识别 .生态学报，2 0 13,33(13):4125-4133.20 Mitchell M G E,Bennett E M,Gonzalez A.Linkinglandscape connectivity and ecosystem service provi-sion:current knowledge and research gaps J.Ecosys-tems,2013,16(5):894-908.21郭家新，胡振琪，李海霞，等.基于MCR模型的市域生态空间网络构建 .农业机械学报，2 0 2 1,5 2(3)：275-284.22陈南南，康帅直，赵永华，等.基于MSPA和MCR模型的秦岭（陕西段）山地生态网络构建 .应用生态学报,2 0 2 1,32(5):15 45-15 5 3.23陈菁，廖克.基于空间数据库的生态环境综合信息图谱的指标体系研究：以莆田市为例 .云南师范大学学报（自然科学版），2 0 0 5,2 5(1)：6 0-6 4.24李怡欣，李菁，陈辉，等.基于MSPA和MCR模型的贵阳市2 0 0 8 2 0 17 年景观连通性评价与时空特征 J.生态学杂志，2 0 2 2,41(6)：12 40-12 48.25Yu K J.Ecological security patterns of landscapes:concept,method and a caseC/Proceedings for Inter-nationalSymposiumofGeoinformatics.HongKong,1995:396-405.26Knaapen J P.Estimating habitat isolation in landscapeplanningJ.Landscape and Urban Planning,1992,23(1):116.27艾婧文，余坤勇，黄茹鲜，等.风电项目对潜在生态廊道的影响一一基于MSPA-MCR模型 .生态学报，2023,43(9):3665-3676.28张豆，渠丽萍，张桀.基于生态供需视角的生态安全格局构建与优化：以长三角地区为例 .生态学报,2 0 19,39(2 0):7 5 2 5-7 5 37.29杨亮洁，王晶，魏伟，等.干旱内陆河流域生态安全格局的构建及优化：以石羊河流域为例 .生态学报，2020,40(17):5915-5927.30耿建伟，余坤勇，谢祯，等.基于生境质量模型的福州市绿色基础设施网络优化研究 .西南林业大学学报（自然科学），2 0 2 3,43(2)：118-12 5.31黄河，余坤勇，高雅玲，等.基于MSPA的福州绿色基础设施网络构建 .中国园林,2 0 19,35(11):7 0-7 5.32陈春娣，吴胜军，MeurkColinDouglas，等.阻力赋值对景观连接模拟的影响 .生态学报，2 0 15,35(2 2)7367-7376.