基于改进Canny算子的水域边界自动提取方法研究.pdf
《基于改进Canny算子的水域边界自动提取方法研究.pdf》由会员分享,可在线阅读,更多相关《基于改进Canny算子的水域边界自动提取方法研究.pdf(7页珍藏版)》请在咨信网上搜索。
1、第 22卷 第 9期2023年 9月Vol.22 No.9Sept.2023软 件 导 刊Software Guide基于改进Canny算子的水域边界自动提取方法研究郑云柯,王俊霖,张世涛(昆明理工大学 国土资源工程学院,云南 昆明 650093)摘要:选用云南省大理市洱海的Landsat-8 OIL影像作为基础数据,将Canny算子与遥感常用的水体指数法相结合,以自动提取出湖泊边界的矢量数据,再经高分影像目视判别和空间缓冲区分析对比各水体指数法的提取精度。针对传统Canny边缘检测算法需要人为设定阈值的情况,在分析水体影像的灰度直方图呈明显双峰分布后,利用类间方差最大阈值分割法优化Canny
2、算子中的高阈值设定原则,以保证高低阈值的自适应性。结果表明,将类间方差最大阈值分割法应用于Canny算子的门限选择能有效提高信噪比。基于遥感影像,利用归一化差分水体指数、改进的归一化差分水体指数、增强型水体指数与Canny算子相结合所提取的湖泊边界均达到了较理想的效果,其中改进的归一化差分水体指数的提取精度相较于其它两种水体指数更优。关键词:多源遥感影像;Canny算子;大津法;水边界提取DOI:10.11907/rjdk.231267开 放 科 学(资 源 服 务)标 识 码(OSID):中图分类号:TP751 文献标识码:A文章编号:1672-7800(2023)009-0167-07Re
3、search on Automatic Extraction Method of Water Boundary Based on Improved Canny OperatorZHENG Yunke,WANG Junlin,ZHANG Shitao(Faculty of Land Resources Engineering,Kunming University of Science and Technology,Kunming 650093,China)Abstract:Landsat-8 OIL images of the Erhai Sea in Dali,Yunnan Province,
4、are used as the base data.The Canny operator is combined with the water body index method commonly used in remote sensing to automatically extract the vector data of the lake boundary,and then the extraction accuracy of each water body index method is compared by visual discrimination of the high-re
5、solution images and spatial buffer analysis.To address the situation that the traditional Canny edge detection algorithm requires artificially set thresholds,the maximum threshold segmentation method with interclass variance is used to optimise the high threshold setting principle in the Canny opera
6、tor to ensure the adaptiveness of high and low thresholds after analysing the grey-scale histogram of the water body image with an obvious bimodal distribution.The results show that the application of the interclass variance maximum threshold segmentation method to the threshold selection of the Can
7、ny operator effectively improves the signal-to-noise ratio,and the lake boundaries extracted based on remote sensing imagery using the normalised difference water index(NDWI),the improved NDWI,and the enhanced water index(EWI)combined with the Canny operator all achieved satisfactory results,with th
8、e improved NDWI achieving better extraction accuracy compared to the other two indices.Key Words:multi-source remote sensing image;Canny operator;Otsu method;extraction method of water boundary0 引言水是人类生存必不可少的资源,是大自然赋予人类最宝贵的财富1。随着全球气候变化和人类活动的不断增加,地球上水体的分布和特征发生了显著变化,给人类社会和自然生态带来了严重挑战,水资源的可持续性问题引起人们越来越
9、多的关注。因此,及时、准确、全面地监测和收稿日期:2023-03-15基金项目:国家自然科学基金项目(U1502231)作者简介:郑云柯(1998-),男,昆明理工大学国土资源工程学院硕士研究生,研究方向为地层学、遥感生态环境;张世涛(1964-),男,博士,昆明理工大学国土资源工程学院教授、博士生导师,研究方向为地层学、地质调查、地质灾害、生态环境。本文通讯作者:张世涛。2023 年软 件 导 刊分析水体信息对于理解水循环过程、评估水资源状况、保障水资源的可持续开发利用至关重要。卫星遥感技术因其能够高效、准确地获取全球范围内的水体信息而备受关注。利用卫星遥感技术,可以获取大量水体信息,从多源
10、、多时相、多尺度的遥感影像中提取和反演水体信息,包括水体范围、深度、透明度、色度、温度等参数,并结合地面观测数据和数值模型进行验证与分析,可以实现对大范围、复杂区域内水体变化情况的动态监测和定量评估,为制定合理、有效的水资源开发利用和管理措施提供科学依据2。如今,快速、准确地从卫星遥感影像上提取水体信息已成为水资源调查、水资源宏观监测及湿地保护的重要手段。1 相关研究目前,利用遥感技术监测湖泊主要集中在湖泊的水质3、水位4、面积5等方面。常见的水体信息提取方法有单波段阈值法6、多波段谱间关系法7和水体指数法等。其中,水体指数法是一种简单、快速、高效的水体信息识别方法,如 Mcfeeters8在
11、 1996 年提出归一化差异水体指数(Normalized Difference Water Index,NDWI),利用绿光波段和近红外波段增强了影像中的水体特征;徐涵秋9在NDWI的基础上,提出改进的归一化差异水体指数(Modified Normalized Difference Water Index,MNDWI),采用中红外波段代替近红外波段,取得了比 NDWI 更好的提取效果;闫霈等10提出增强型水体指数(Enhanced Water Index,EWI),相较于NDWI和MNDWI,能更有效地区分河道与背景噪音。针对不同研究区的水体特性,以上水体指数法都取得了较好效果。湖泊作为具有
12、一定面积的相对封闭水域,为了获取湖泊的矢量边界,常用的方法包括手工数字化方法、区域分割法11与边缘检测法12。当前国内外学者对水边界提取的研究主要集中在海岸线与湖岸线,如Alesheik等13利用阈值分割法对乌尔米耶湖岸线进行提取;李秀梅等14利用Canny算子对渤海湾海岸带进行提取以监测其时空变化;Karantzalos 等15基于遥感影像,分别用 Laplacian 算子和 Canny算子提取海岸线;申家双等16在分析了边缘检测算法用于影像水边界提取的优缺点后,提出将Canny算子与GAC模型结合提取影像水边界的方法;Abolhassani等17利用边缘检测算法提取了美国东海岸线;魏东岚等
13、18在MATLAB平台上利用边缘检测算法、阈值分割法及小波变换法提取海岸线,结果表明边缘检测法中的Canny算子能准确提取出围填区域的海岸线。已有研究主要使用Canny边缘检测算子提取海岸线,且 Canny算子的阈值多是通过大量实验获取,这是 Canny算子自身的不足:高低阈值需人为设定,容易导致检测中出现大量虚假边缘19。因此,本文在已有研究基础上,以湖泊为研究对象,将Canny边缘检测算子与遥感水体指数法相结合以提取湖岸边界。针对Canny算子的不足,利用类间方差最大阈值分割法来计算其阈值,使Canny算子具有较理想的阈值,从而实现自动化提取湖泊水边界。2 研究区概况洱海位于云南省大理市郊
14、区,地跨大理市和洱源县。洱海是云贵高原九大湖泊之一,也是云南省第二大淡水湖,海拔1 980 m,经纬度为10005-10018E,2536-2558N。洱海北起洱源县,南至大理市下关镇,湖泊呈东西窄、南北长的条带状,因湖泊形状像耳朵而得名。洱海水域面积250 km2,平均水深10.5 m,平均水位1 974 m,湖水主要靠河流补给。3 数据来源及预处理本文采用的遥感数据是从USGS网站上下载的大理市洱海的Landsat-8 OIL影像,轨道号/行号为131/42,选取影像的云量少于 3%,遥感影像获取时间为 2019 年 11 月 29号,处于湖泊的枯水期。Landsat-8 OIL数据共有9
15、个波段,除全色波段空间分辨率为15 m外,其余波段都为30 m。另外,为验证水体边界提取精度,从中国资源卫星应用中心获取了一幅时期相近的大理市洱海的2 m全色/8 m多光谱GF-6号影像,遥感影像获取时间为2019年12月7号。对Landsat-8 OIL影像进行辐射定标、大气校正、图像裁剪等一系列预处理,对 GF-6号影像进行 2 m全色波段和8 m多光谱融合,形成空间分辨率为2 m的融合影像。4 研究方法4.1水体信息增强国内外学者根据水体在蓝绿波段吸收率较低、近红外波段吸收率较高这一光谱特征,开展了大量基于遥感影像的表面水体自动提取算法研究。本文将利用不同的影像波段,分别采用归一化差分水
16、体指数法、改进的归一化差分水体指数法和增强型水体指数法来增强洱海表面水体信息,如表1所示。Table 1Model of the water body index表1水体指数模型指数名称NDWIMNDWIEWI模型公式Green-NIRGreen+NIRGreen-MIRGreen+MIRGreen-NIR-MIRGreen+NIR+MIR 公式意义Green代表绿波段NIR代表近红外波段MIR代表中红外波段 168第 9 期郑云柯,王俊霖,张世涛:基于改进Canny算子的水域边界自动提取方法研究4.2水体边缘检测算法传统的边缘检测算法有Roberts算子、Prewitt算子、Sobel算子等
17、,这些算法简单,虽容易实现,但处理噪声的能力较差,裂纹边缘识别不完整,还容易出现伪边缘现象。能否选择合适的边缘检测算子,将直接影响到结果的精度。笔者对前人的研究结果进行比较,发现使用Canny算子处理边缘和噪声的效果最佳20。因此,本文选用Canny算子对水体增强的遥感影像进行边缘检测,主要分为以下4个步骤:(1)利用高斯函数对图像进行平滑处理。高斯函数表达式为:G(x,y)=122exp(-x2+y222)(1)平 滑 图 像 后 的 结 果I(x,y)表 达 式 为:I(x,y)=G(x,y)*f(x,y),其中G(x,y)为高斯函数,f(x,y)为原图像。(2)对平滑后的图像进行梯度幅值
18、和方向的计算。基于22模板,通过X、Y方向像素的一阶导数来确定梯度幅值。设Fx(x,y)为X方向的偏导数,Fy(x,y)为Y方向的偏导数,则梯度幅值M(x,y)和方向的表达式为:M(x,y)=Fx()x,y2+Fy()x,y2=arctanFy()x,yFx()x,y(2)(3)对梯度幅值的非极大值进行抑制处理,判断梯度幅值在其八邻域内是否为最大值,若为最大值则为边缘,否则置零。(4)对梯度选取两个门限阈值,即高门限阈值TH和低门限阈值TL,两者通常的关系为21:TL=0.5*TH(3)取出经过非极大值抑制后图像中的最大梯度幅值,标记大于高门限阈值的点(边缘点),将小于低门限阈值的点置0,即可
19、提取出完整边缘。4.3Canny算子阈值设定由Canny算法步骤可知,门限阈值的选取是进行图像边缘提取的关键,而Canny算子需要人为预先设定高低阈值,因此需要很多次反复实验才能找到合适的阈值。为解决阈值需要人为预先设定的问题,本文利用类间方差最大阈值分割算法来确定Canny算子的门限阈值。类间方差最大阈值分割算法别称大津(Otsu)法22,是由日本学者大津在 1979 年提出的一种自适应阈值确定法。其数学描述为对于给定的图像L(x,y),设图像的灰度级为M,灰度值h像素点的像素数为n,则其出现的概率为23:ph=1MNf(x,y)=hf(x,y)(4)采用阈值t为界限把图像区分为两个区域,假
20、设背景像素区域为A1,前景像素区域为A2,则两个区域的总概率为:A1=h=0tPhA2=h=t+1M-1Ph(5)设图像整体的灰度均值为,则前景像素区域A1与背景像素区域A2的灰度均值为:A 1=h=0thPhA 1A2=h=t+1M-1hPhA2(6)即当2为最大值时,对应的t为选择的最佳阈值。根据上述Otsu算法原理确定一个分割阈值t,将选取的分割阈值t作为Canny算子的高阈值TH,再用此高阈值乘以一个比例因子 0.5作为其低阈值TL24-25。得到阈值的表达式如下:TH=tTL=0.5t(7)由此 Canny 算子的高低阈值选择只与分割阈值t有关,不需要再人为多次反复实验设定,这就为
21、Canny算子中阈值的确定问题提供了一个较好的解决方法,也增强了Canny算子的分割能力和自适应性。总体技术路线如图1所示。5 结果与分析5.1水体边界自动提取5.1.1水体灰度差增强如图 2 所示,NDWI、MNDWI 和 EWI 水体指数法都能有效拉伸水体和非水体之间的灰度差异。从目视解译的效果来看,3种方法拉伸灰度差异的效果各不相同。相比于 MNDWI 法和 EWI法,NDWI 法对洱海南边的出水口即西洱河的提取有所欠缺。5.1.2边缘矢量化先对3幅水体增强影像进行归一化处理,由大津法的影像预处理应用3种水体指数法增强水体信息影像归一化利用Canny算子提取水体边界分析比较3种水体指数法
22、的边缘检测结果精度评价空间统计分析目视解译得到精度最高的水体提取方法Otsu算法得到分割阈值t结论代入Canny算子的高低阈值关系式0.5THtTLt=Fig.1Technical flow图1技术流程 1692023 年软 件 导 刊数学原理可知其特性26-27:当目标地物与背景地物的面积比例悬殊时,灰度直方图可能会无明显“双峰”或呈“多峰”形态,此时使用大津法的分割效果不佳。为得到最优阈值,需满足灰度频率直方图呈“双峰”的情况28。通过图2可看出目标地物水体与背景地物非水体的面积比例相差不大,研究发现 3 种水体指数法的灰度频率直方图均呈“双峰”分布,如图3所示。其中,横坐标Data Va
23、lue代表归一化后的灰度像素值。然后,使用 Ostu 算法得到水体增强影像的分割阈值(保留两位有效数字),最后利用关系式(7)得到 Canny算子的高阈值TH和低阈值TL,结果如表2所示。将提取好的栅格影像结果转换为矢量数据,再与原始遥感影像叠加显示。如图4所示,整体上,水体边界提取的边缘连续性及去噪效果都达到最佳。3种水体指数法均可提取出洱海的水体边界,对于洱海的北部、西部、东部的非目标地物类,也均提取了其边缘。就北部的目视效果而言,MNDWI+Canny算子提取的非目标地物略少一点。相比其它两种水体指数法,NDWI+Canny算子还提取出了位于洱海西南部的苍山积雪边缘及南部的部分建筑边缘。
24、因此,选取合适的边缘检测阈值,可减少目标地物的非必要干扰信息,为后续的处理节省时间。5.1.3边缘细部提取效果选取洱海湖岸的湿地、河滩及洱海出水口3种典型地区进行分析,湿地是水体、植被与裸地等按不同比例混合(a)NDWI+Canny operator(a)NDWI+Canny算子(b)MNDWI+Canny operator(b)MNDWI+Canny算子(c)EWI+Canny operator(c)EWI+Canny算子Fig.4Vector boundary overlay image of water body图4水体矢量边界叠加影像0.00.20.40.60.81.00.05.010
25、31.01041.51042.01042.51043.01043.5104FrequencyData Value(a)NDWI histogram(a)NDWI直方图0.00.20.40.60.81.00.01.01042.01043.01044.01045.01046.01047.0104FrequencyData Value(c)EWI histogram(c)EWI直方图0.00.20.40.60.81.00.05.01031.01041.51042.01042.51043.01043.51044.0104FrequencyData Value(b)MNDWI histogram(b)M
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 改进 Canny 算子 水域 边界 自动 提取 方法 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。