卫星遥感辅助的大型海藻养殖...态对比监测——以威海市为例_侯英卓.pdf

第 卷，第 期自 然 资 源 遥 感，年 月 ，：引用格式：侯英卓，纪灵，邢前国，等 卫星遥感辅助的大型海藻养殖动态对比监测 以威海市为例 自然资源遥感，（）：（，（）：）卫星遥感辅助的大型海藻养殖动态对比监测 以威海市为例侯英卓，纪 灵，邢前国，盛德志，（中国科学院烟台海岸带研究所，中国科学院海岸带环境过程与生态修复重点实验室，烟台；山东省海岸带环境过程重点实验室，烟台；中国科学院大学，北京；国家海洋局烟台海洋环境监测中心站，烟台）摘要：大型海藻养殖的时空动态变化监测对其环境管理至关重要，目前关于不同品种的大型海藻养殖的对比监测鲜有研究报道。文章基于 卫星影像，利用归一化植被指数（，）与支持向量机（，），对山东省威海市文登区南部海域紫菜养殖区与荣成市南部海域海带养殖区的动态特征进行了监测。研究结果表明：威海市文登区的紫菜养殖在 年遥感影像上首次出现，与该市历史上首次出现紫菜养殖的年份相符；基于文章方法提取的紫菜养殖区与海带养殖区的整体提取效果较好，总体精度可达 以上；年度紫菜养殖区的遥感监测面积整体呈逐年增加趋势，空间上养殖区呈现远离岸边的分布趋势；紫菜与海带养殖区监测面积总体呈冬高、夏低的冷水型藻类养殖季节变化特征，但紫菜养殖区监测面积最小值与最大值出现时间较海带养殖区早 个月。卫星遥感较统计年鉴能提供更精确的时间与空间信息，研究可为中国北方海岸带大型海藻养殖管理提供监测技术与数据上的参考与借鉴。关键词：；大型海藻；紫菜；海带；山东半岛；黄海中图法分类号：文献标志码：文章编号：（）收稿日期：；修订日期：基金项目：中国科学院 类战略性先导科技专项“地球大数据科学工程”（编号：）、国家自然科学基金项目“海表漂浮大型藻空间分布特征及其在覆盖面积卫星遥感估算中的应用研究”（编号：）、“面向冬季大型藻类高分遥感的海表耀光消减与利用研究”（编号：）和中国科学院仪器设备研制重点项目“沿海水色环境污染和资源机载高光谱成像探测仪”（编号：）共同资助。第一作者：侯英卓（），男，硕士研究生，研究方向为海洋遥感。：。通信作者：邢前国（），男，博士，研究员，研究方向为海岸带 海洋环境遥感、数值模拟与评估。：。引言大型海藻养殖是海水养殖的重要组成部分，主要分布在我国东部沿海地区，包括江苏、福建、浙江、辽宁和山东等地。我国的大型海藻养殖种类众多，主要有紫菜、海带和裙带菜等。海藻不仅具有巨大经济价值，亦具有调节大气、净化水质和固碳等生态功能。然而，海藻养殖不当也会产生环境危害，已有研究显示，我国近海的海洋大型藻等灾害与海藻养殖有着密切的联系。开展藻类养殖区的动态监测，不仅可以掌握养殖区的动态变化特征，还可以辅助分析周围水环境的变化，为相关海洋灾害的产生机理研究提供基础数据。山东半岛，特别是威海市沿岸，是我国北方最重要的海藻养殖基地之一。近年，受市场对养殖品种的需求、海水环境与灾害等因素影响，威海市近海的海藻养殖品种与空间布局正发生深刻变化。对本地区海藻养殖动态进行监测，对其管理及灾害应对有重要的现实意义。卫星遥感具有周期短、范围大、观测连续等优势，可以提供传统调查不具备的精细化空间信息。国内外众多学者利用遥感影像对大型海藻养殖展开了广泛研究。在海藻养殖区的提取方法上，目前主要有人工目视解译、归一化植被指数（，）阈值法、最大似然分类、支持向量机（，）与面向对象分类等，。武易天等提出了一种耦合光谱与纹理信息的近海藻类养殖区自动提取方法；等提出了一种基于极化熵信息的水下紫菜养殖区提取方法。在提取水色背景相对稳定的大型海藻养殖区时，阈值法的监测效果较好。但当海水中混有大量悬浮第 期侯英卓，等：卫星遥感辅助的大型海藻养殖动态对比监测 以威海市为例泥沙、水色背景较复杂时，相比 阈值法，有更高的监测精度，。在大型海藻养殖的遥感动态监测上，国内学者研究较多，而国外目前相关研究较少。等对苏北浅滩的紫菜养殖动态进行了长时序精细化的遥感监测，获取了紫菜养殖年际变化以及每年的养殖设施回收进程的时空动态特征，并阐明了其对黄海浒苔绿潮的影响；郭瑞宏基于 卫星的时序数据，分析了江苏省紫菜养殖区的年际变化特征；等利用 期 影像，对连云港紫菜养殖区的空间变化特征展开研究；曹万云利用遥感影像与海带产量数据，分析了荣成市海带养殖的时空变化特征等。但是上述研究仅针对一种藻类养殖进行动态监测，缺乏不同类型大型海藻养殖的对比监测，且山东沿海的紫菜养殖动态研究鲜有报道。因此，本文基于 卫星影像，结合目前较常用的 与 方法提取山东省威海市文登区南部海域紫菜养殖区与荣成市典型海域海带养殖区信息，分析紫菜养殖的年际变化，并对比紫菜与海带养殖的季节变化特征。研究区概况研究区位于山东省威海市文登区与荣成市的南部海域。文登区（，）北连威海市环翠区，东邻荣成市，南濒黄海。荣成市（，）东、南、北三面濒临黄海。文登区与荣成市的近海是山东省海藻养殖的重要区域，如图 所示。据资料显示，威海市文登区于 年开始试养条斑紫菜；荣成市海域是山东省海带养殖的主要区域。（）真彩色影像（）区域 方形紫菜养殖区（）区域 条带状海带养殖区 （）区域 条带状海带养殖区图 研究区地理位置及大型藻类养殖区分布 数据源本文以 卫星影像为基础数据，数据在欧洲航空局官网（：）下载。所用卫星的传感器为多光谱成像仪（，），波段为 空间分辨率的蓝光、绿光、红光和近红外波段。经过筛选，共选取 年的 景影像用于监测文登区南部海域的紫菜养殖区，年间的 景影像用于监测荣成市南部海域的海带养殖区。所用影像时间如表 所示。统计年鉴数据来自中国渔业统计年鉴山东渔业统计年鉴和威海统计年鉴。表 所用影像介绍 年份影像日期 年，年，年，年，年，年，自 然 资 源 遥 感 年 研究方法本文技术路线如图 所示，数据预处理主要有辐射定标、大气校正和影像裁剪等，提取后处理包括面积筛选与人工修正。其中，大气校正使用欧洲航空局提供的 插件进行。图 技术路线 紫菜富含叶绿素、胡萝卜素和叶黄素等色素，具有与陆地绿色植物相似的光谱特征，即紫菜光谱存在红边效应。本文随机在 影像的图（）区域 范围内选取典型紫菜养殖区样本与海水样本，各 个像元，样本分布如图（）所示。分别统计紫菜养殖区与海水样本的反射率均值，如图（）所示。紫菜养殖区的光谱曲线在近红外波段处的值高于海水，而在红光波段处的值低于海水。所以，采用 可以区分紫菜养殖区与海水。的表达式为：，（）式中 和分别 为近红外和红光波段反射率，分别对应 影像第 与第 波段的反射率。针对紫菜养殖区提取的 阈值确定，在 影像上的紫菜养殖区域选取典型紫菜养殖区样本与海水样本，各 个像元，样本分布如图（）所示，底图为 值。统计样本的 分布频率，如图（）所示。根据频率分布情况，可初步确定区分紫菜养殖区与海水的阈值为 ，针对不同影像进行阈值调整。此外，经统计发现每块紫菜养殖方的面积均大于 ，所以将筛选出面积大于 的斑块作为紫菜养殖区。故本文通过设置 阈值初步提取紫菜养殖区，再对阈值提取结果进行面积筛选及人工修正，最终实现紫菜养殖区的监测。（）反射光谱样本分布（）样本的光谱曲线（型误差线表示标准差）（）及样本分布（）频率分布图 样本分布及相应的光谱曲线与 频率分布 第 期侯英卓，等：卫星遥感辅助的大型海藻养殖动态对比监测 以威海市为例 经多次实验发现，在利用 提取荣成市南部海域的海带养殖区时，阈值法提取的效果较差。海带养殖区与周围海水均混有悬浮泥沙，导致海带养殖区与海水难以通过 阈值法进行区分。故本文采用 机器学习分类方法，结合目视解译、面积筛选（将面积大于 的斑块作为海带养殖区）与人工修正实现海带养殖区的监测。其中，的核函数为线性核函数，且分类样本中典型海带养殖区样本与海水样本的分离度在 以上（研究表明样本分离度在，之间，尤其大于 时最佳）。在监测结果的精度评价中，分别目视选取养殖区与非养殖区验证样本，各 个像元，计算混淆矩阵，利用总体分类精度（，）与 系数进行评价。每个季节选一景进行评价，选取 年 月、月、月和 月的紫菜与海带养殖区的遥感监测结果进行精度评价。结果与分析 监测结果与精度评价图（）中区域 范围中紫菜养殖的年度遥感监测结果如图 所示。紫菜的年度养殖方与年度监测面积，由前一年下半年下苗和第二年上半年的收割构成的一个养殖周期内的各期结果合并而成，如 年度的紫菜养殖周期为 年 月 年 月。文登区南部海域的紫菜养殖块呈方形，与江苏连云港的紫菜养殖方形状类似；整体上，紫菜养殖区在五垒岛湾近海呈南北向分布，南北向跨度大，东西向跨度小。海带养殖区的监测结果如图 所示（仅展示 年 月 日的监测结果，此时海带养殖的监测面积达 年最大），荣成市南部海域的海带养殖呈条带状，紧沿岸边分布。通过目视观测，整体上遥感提取的紫菜与海带养殖区可以较好地覆盖底图显示的养殖方与养殖条带。（）年度（）年度（）年度（）年度（）年度图 年度紫菜养殖区的遥感监测结果 图 年 月 日海带养殖区的遥感监测结果 ，基于 阈值法，利用 遥感影像追溯监测文登区南部海域大型藻类养殖历程发现，年 月 日及以前影像上未监测到养殖区，年 月 日影像上开始有紫菜养殖区的分布；年 月 日，监测到的紫菜养殖区总面积达到极大值（图（），约 ；卫星监测可清晰识别养殖方，其最大数量为 个。据相关资料显示，威海市文登区于 年秋季开始在南部海域开展插杆式紫菜养殖，与卫星遥感结果一致；个养殖方（每方约 ），与遥感监测结果间存在差异，这可能由于实际生产与原计划有一定出入。基于验证样本计算 与 系数，计算结果如表 所示。针对 期紫菜养殖区的遥感监测结果，本文紫菜养殖区的 均大于，系数均大于 。针对 期海带养殖区的监测结果，均大于，除 年 月 日外，其余 期海带养殖区的 系数均大于 。自 然 资 源 遥 感 年表 精度评价结果 指标紫菜养殖区海带养殖区 紫菜养殖年际变化遥感监测的紫菜养殖区面积与养殖方个数如图 所示，年，文登区紫菜养殖监测面积的逐月变化具有明显的周期性，每个周期内的遥感监测面积先增后减，故将每一养殖周期定义为一个养殖年度。文登区自 年度开展养殖紫菜，年度便大幅增加，其后呈逐年增加趋势。年度紫菜养殖的监测面积比 年度约多 ，养殖方个数多 个。年度紫菜养殖的监测面积与养殖方个数变化最大；在 年以后，面积与养殖方数目的变化速率均减小，年 度 变 化 最 小。在 空 间 上，如 图 所 示，年度紫菜养殖的空间分布整体上逐渐变大，且变化明显，其中 年度变化最大，年度与 年度变化相对较小。文登区自开始养殖紫菜，养殖区域呈先向岸扩展后离岸扩展。（）年养殖区监测面积变化（）年度养殖区监测面积与养殖方个数变化图 遥感监测的紫菜养殖区面积与养殖方个数变化 紫菜与海带养殖季节变化文登区的紫菜养殖在每一年度均具有明显的季节性，如图（）所示，每一年度的逐月监测面积整体呈先增后减的趋势；年度，紫菜养殖监测面积在 月、月份监测面积较大，在 月份监测面积较小。紫菜养殖的月均值曲线亦具有明显的季节性，如图（）所示，在 月份监测面积最大，在 月份监测面积最小。整体上，文登区紫菜养殖区的监测面积冬季最大，夏季最小。（）逐月变化（）月均值变化（型误差线表示标准差）图 紫菜养殖遥感监测面积的季节变化 第 期侯英卓，等：卫星遥感辅助的大型海藻养殖动态对比监测 以威海市为例 选取 年与 年的 月的遥感影像，每月选用一景，开展紫菜与海带养殖区的对比监测，养殖区监测面积的逐月变化如图（）所示。在 年，紫菜养殖的监测面积分别于 年 月 日、年 月 日出现极小值。在 年 月 日和 年 月 日，海带养殖的监测面积处于极小值，面积分别约为 和。可以发现，紫菜养殖监测面积的极小值出现的时间比海带的大致早 个月。在 年，紫菜养殖的监测面积于 年 月 日和 年 月 日出现极大值，海带养殖的监测面积于 年 月 和 年 月 日出现极大值。可见紫菜养殖监测面积的极大值出现的时间也比海带的大致早 个月。（）逐月变化（）月均值变化（型误差线表示标准差）图 年紫菜与海带养殖的监测面积变化 将 年与 年相应月份的紫菜与海带养殖的监测面积分别取均值，如图（）所示，因夏季的影像云量较多，所以共得到 个月份的紫菜监测面积均值和 个月份的海带监测面积均值。可以看到，在 月份紫菜养殖监测面积的月均值整体先减后增，而海带养殖监测面积的月均值整体先增后减再增。紫菜养殖监测面积的月均值的最小值出现在 月，而海带的出现在 月，时间间隔为 个月。此外，图 中的海带养殖监测面积的曲线整体较紫菜的曲线均存在右移现象。上述表明，在同一养殖周期中，紫菜养殖监测面积的最大值出现在冬季，而海带养殖出现在初春；紫菜养殖监测面积的最小值出现在夏季，而海带养殖出现在秋季。整体而言，紫菜养殖监测面积最大值与最小值出现的时间比海带的早 个月。与统计年鉴数据的对比分析统计年鉴中紫菜与海带养殖面积的统计结果如图 所示。（）山东省、威海市与文登区的紫菜养殖面积变化 （）山东省、威海市与荣成市的海带养殖面积变化图 年统计年鉴中紫菜与海带的养殖面积变化 对比图（）与图（）发现，整体上，紫菜年度养殖面积在遥感监测与统计年鉴中均呈现增加趋势。紫菜养殖面积的大幅度变化在本文遥感监测的结果中出现在 年度，而在统计年鉴中出现在 年；这是因为统计年鉴中的紫菜年度养殖是由每年 月份的紫菜养殖组成，其自 然 资 源 遥 感 年与本文所定义的紫菜年度养殖不同。如图（）与图（）所示，本文遥感监测的荣成市部分海域的海带养殖面积整体变化不明显，同样统计年鉴中 年荣成市的海带养殖面积整体变化也不明显。此外，图 显示，威海市的紫菜养殖面积约占山东省总养殖面积的，海带养殖面积约占山东省总养殖面积的；威海市的紫菜养殖均在文登区，而海带养殖均在荣成市，紫菜与海带养殖面积变化与山东省的总养殖面积变化基本一致。针对养殖周期跨年的大型海藻养殖，统计年鉴按日历周期进行统计会导致时间信息不明确。由图 可知，遥感能提供精细的季节变化特征，可实现按养殖周期进行结果统计与分析；并且，遥感能提供精准的空间信息，可对养殖区大型海藻生长的空间差异进行即时评估。如前所述，考虑到研究区的紫菜与海带养殖在山东省占主要比重，基于遥感的监测结果对威海市及山东省的紫菜与海带养殖状态有更好的指示作用。结论与展望本研究基于 卫星数据，采用 阈值法和 方法对威海市文登区南部海域的紫菜养殖区以及荣成市南部海域的海带养殖区的动态进行了对比监测。研究发现并证实了，威海市及山东省紫菜养殖的大幅增加主要源于 年秋季开始的威海市文登区南海海域的紫菜养殖扩张；研究区的紫菜与海带养殖区遥感监测面积均呈冬高夏低的周期性变化特征，紫菜养殖遥感监测结果的变化周期较海带早 个月。卫星遥感数据能提供紫菜与海带养殖精准的生长期时空动态数据，可以弥补统计年鉴数据时间、空间信息不全的不足，这对威海市乃至山东省的大型海藻养殖管理有重要参考价值。在大型海藻养殖区遥感提取的基础上，可进一步利用遥感监测估算其生长状态及生物量，结合海温、水质等遥感监测对海藻产量及可能的灾害进行预警评估。参考文献（）：刘福利，牛建峰，隋正红，等 中国经济海藻养殖技术概况与展望 中国农业科技导报，（）：，（）：卢 霞，顾 杨，王晓静，等 连云港近海紫菜养殖区遥感识别、空间变异和驱动分析 海洋科学，（）：，（）：李 涵，青 平，张晓恒，等 我国藻类产业发展现状、问题及对策建议 中国农业科技导报，（）：，（）：岳冬冬，王鲁民，耿 瑞，等 中国近海藻类养殖生态价值评估初探 中国农业科技导报，（）：，（）：魏振宁，邢前国，郭瑞宏，等 基于遥感的 年南黄海紫菜养殖空间分布变化研究 海洋技术学报，（）：，（）：，：刘材材，徐 韧，何培民，等 南黄海绿潮暴发与紫菜养殖的关系 海洋科学，（）：，（）：许海蓬，张彦彦，郑 剑 基于 和中值滤波的连云港海域紫菜筏架自动提取研究 海洋技术学报，（）：，（）：郭瑞宏 江苏沿海紫菜养殖时空动态遥感（）北京：中国科学院大学，（）：，初佳兰，邵光辉，赵建华，等 高分一号的浮筏养殖信息提取方法 测绘科学，（）：，（）：武易天，陈 甫，马 勇，等 基于 数据的近海养殖区自动提取方法研究 国土资源遥感，（）：，（）：，第 期侯英卓，等：卫星遥感辅助的大型海藻养殖动态对比监测 以威海市为例 ，（）：胡姣婵，黄梦迪，于浩洋，等 基于哨兵二号遥感影像的近海养殖区提取方法研究 海洋环境科学，（）：，（）：，（）：曹万云 荣成海带养殖时空变化与适宜性评价研究 烟台：鲁东大学，：，海洋网 威海海洋食品又添新品 文登千亩紫菜收割上市（）：？（）：？有贺祐胜，于廷林 条斑紫菜（）的颜色和色素 国外水产，（）：，（）：吴健生，潘况一，彭 建，等 基于 决策树的遥感影像土地利用分类 以云南省丽江市为例 地理研究，（）：，：，（）：，（，；，；，；，）：，（）（），：，；，；，；（），：；（责任编辑：陈 理）