《一种基于全色遥感影像的城市水体提取方法.pdf》由会员分享，可在线阅读，更多相关《一种基于全色遥感影像的城市水体提取方法.pdf（15页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 103761717 A (43)申请公布日 2014.04.30 CN 103761717 A (21)申请号 201410037069.1 (22)申请日 2014.01.26 G06T 5/00(2006.01) G06T 7/00(2006.01) (71)申请人 河海大学 地址 210098 江苏省南京市鼓楼区西康路 1 号 (72)发明人 徐南 马符讯 马成 严朝霞 徐大钊 贾涛 李代云 艾斯卡尔阿不力米提 黄卓 (74)专利代理机构 南京苏高专利商标事务所 ( 普通合伙 ) 32204 代理人 柏尚春 (54) 发明名称 一种基于全色遥感影像的城市水体提取。
2、方法 (57) 摘要 本发明公开一种基于全色遥感影像的城市水 体提取方法， 通过遥感影像预处理、 最优尺度选 取、 均变纹理生成和最优阈值分割四个步骤最终 获得所需提取的城市地区的水体信息。本发明计 算方法简单， 运算量较小， 并且最终提取结果较为 精准， 能够应用于城市规划、 环境科学、 地理信息 制图等多个领域。 (51)Int.Cl. 权利要求书 2 页 说明书 5 页 附图 7 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书5页 附图7页 (10)申请公布号 CN 103761717 A CN 103761717 A 1/2 页 2 1. 一种基。
3、于全色遥感影像的城市水体提取方法， 其特征在于包括以下步骤 ： （1） 遥感影像预处理 ： 利用中值滤波抑制遥感影像中的噪声 ； （2） 最优尺度选取 ： 在经步骤 （1） 预处理的图像中， 选定最小尺度、 最大尺度以及尺度 变化步长， 依次计算全局平均方差， 求得最优空间尺度 ； （3） 均变纹理生成 ： 依据步骤 （2） 中所得的最优空间尺度， 遍历遥感影像中的各个像 素， 分别计算均变方差， 进而生成遥感影像的均变方差纹理图 ； （4） 最优阈值分割 ： 基于步骤 （3） 中所得遥感影像的均变方差纹理图， 人工选定分割阈 值， 对方差纹理图进行分割， 选取大于给定面积阈值的斑块， 进行闭。
4、运算， 获得最终水体信 息。 2. 根据权利要求 1 所述的基于全色遥感影像的城市水体提取方法， 其特征在于 ： 所述 步骤 （2） 中确定最优尺度的详细步骤如下 ： （1） 确定最小尺度 scale_min， 最大尺度 scale_max 以及尺度变化步长 scale_step ； （2）对 于 第 k 个 空 间 尺 度 scalek=scale_min+scale_step*(k-1)， 其 中， k=1,2, ,(scale_max-scale_min)/scale_step+1，在 原 始 遥 感 影 像 上 设 置 一 个 大 小 为 (2*scalek+1)(2*scalek+1。
5、) 的窗口， 不断移动该窗口， 对原始图像 im 中第 i 行第 j 列的 像元 im(i,j)， 窗口中的所有像元光谱值记为 im(i-scalek:i-scalek,j-scalek:j-scalek)， 然后计算窗口中所有像元光谱值的方差 imvar(i,j,k) ： 上述公式中， ii 代表某个像元在窗口中的行号， jj 代表某个像元在窗口中的列号 ； 将 imvar(i,j,k) 作为该窗口中心像元 im(i,j) 在空间尺度 scalek下的方差， 并按照上述方法 计算处全图中所有像元的方差， 并求其平均值 imvar0； （3） 比较 k 个空间尺度下的 imvar0， 选择 i。
6、mvar0最大时所对应的空间尺度 scalebest作为 遥感影像的最优空间尺度， scalebest=1， 则移动窗口大小为 33。 3. 根据权利要求 1 所述的基于全色遥感影像的城市水体提取方法， 其特征在于 ： 所述 步骤 （3） 中遥感影像均变方差纹理图的具体生成方法如下 ： （1）依据所得的最优空间尺度 scalebest， 在原始遥感影像上设置一个大小为 (2*scalebest+1)(2*scalebest+1) 的窗口 ； （2） 不断移动该窗口， 对于像元 im(i,j)， 窗口中的所有像元光谱值记为 im(i-scalebes t:i-scalebest,j-scaleb。
7、est:j-scalebest)， 每个窗口中的像元总个数记作 N=(2*scalebest+1)*( 2*scalebest+1)， 将其按照列方向排列记作 imN， 将窗口中所有像元的行坐标和列坐标依次记 作 rowN和 colN； （3）依此对每个窗口中的所有像元的光谱值进行最小二乘平面拟合， 拟合方程为 ： Ax+By+Cz+1=0， 将待拟合的 N 个像元表示成以下矩阵形式 ： 权 利 要 求 书 CN 103761717 A 2 2/2 页 3 依据最小二乘法则得到平面拟合系数为 ： 上述公式中， h 代表 N 个像元中的第 h 个像元 ； （4） 计 算 窗 口 中 所 有 像 。
8、元 光 谱 值 至 拟 合 平 面 的 距 离 ： 其中 abs(*) 是绝对值算子 ； 计算该窗口中的各个像元所对应的距离， 求得方差 dvar 将 dvar作为该窗口中心像元的均变方差纹理， 遍历所有点， 即生成均变方差纹理图 imvar_texture； 遍历所有像素点， 搜索均变方差纹理图 imvar_texture的最大值和最小值， 分别记作 imvar_ texture_max 和 imvar_texture_min， 对每个像素点作如下操作用于将纹理值拉伸至 0,255 ： 4. 根据权利要求 1 所述的基于全色遥感影像的城市水体提取方法， 其特征在于 ： 所述 步骤 （4） 中。
9、最优阈值分割的具体步骤如下 ： （1） 将方差纹理图中， 由均变方差纹理图的灰度直方图确定影像中的水体和非水体的 纹理值分布， 从而获得分割阈值 ； （2） 在步骤 （1） 中确定的分割阈值中， 根据实际情况确定一个面积阈值用于抑制微小斑 块， 依次遍历每一个斑块， 计算其面积， 如果该斑块的面积小于该阈值， 则去除该斑块， 否则 保留该斑块 ； （3） 基于上述结果， 进行闭运算， 用于填充水体中的微小空洞， 至此， 已经得到最终的水 体信息 ； 通过另设一个面积阈值， 保留面积大于给定阈值的水体区域， 进而显示面积较大的 水体。 权 利 要 求 书 CN 103761717 A 3 1/5。
10、 页 4 一种基于全色遥感影像的城市水体提取方法 技术领域 0001 本发明涉及一种图像信息提取方法， 具体涉及一种基于全色遥感影像的城市水体 提取方法。 背景技术 0002 城市水体作为一种重要的空间地物和自然资源， 对于城市规划、 环境科学和城市 遥感制图都发挥着不可或缺的作用。近年来， 随着卫星遥感技术的发展和各种空间地物提 取算法的出现， 水体提取方法不断进步， 在各项研究中得到了广泛应用。目前， 常用的水体 提取算法都是基于多光谱遥感影像如 SPOT 和 TM 遥感影像， 使用光谱变换生成各种水体指 数， 然后通过监督分类或者阈值分割进行水体提取， 但是少有基于全色遥感影像和针对城 。
11、市地区水体提取的算法。 0003 针对水体提取问题， 目前少有针对全色遥感影像和城市区域的提取算法， 并且一 般仅利用了光谱信息， 没有考虑到纹理特征， 没能够充分利用遥感影像的结构化特征， 未能 够最大化利用遥感中丰富的信息。 0004 因此， 提出基于全色遥感影像并且针对城市区域的水体提取算法非常有意义， 也 会对同类研究起到积极作用。 发明内容 0005 发明目的 ： 为解决现有技术中存在的不足， 本发明提供了一种基于全色遥感影像 的城市水体提取方法。 0006 技术方案 ： 本发明的一种基于全色遥感影像的城市水体提取方法， 包括以下步 骤 ： 0007 （1） 遥感影像预处理 ： 利用。
12、中值滤波抑制遥感影像中的噪声 ； 0008 （2） 最优尺度选取 ： 在经步骤 （1） 预处理的图像中， 选定最小尺度、 最大尺度以及 尺度变化步长， 依次计算全局平均方差， 求得最优空间尺度 ； 0009 （3） 均变纹理生成 ： 依据步骤 （2） 中所得的最优空间尺度， 遍历遥感影像中的各个 像素， 分别计算均变方差， 进而生成遥感影像的均变方差纹理图 ； 0010 （4） 最优阈值分割 ： 基于步骤 （3） 中所得遥感影像的均变方差纹理图， 人工选定分 割阈值， 对方差纹理图进行分割， 选取大于给定面积阈值的斑块， 进行闭运算， 获得最终水 体信息。 0011 进一步的， 所述步骤 （2。
13、） 中确定最优尺度的详细步骤如下 ： 0012 （1）确定最小尺度 scale_min， 最大尺度 scale_max 以及尺度变化步长 scale_ step， 考虑到尺度过大之后会造成图像中出现明显的边缘效应， 因此过大的尺度不予考虑， 例如本发明中各项参数选取如下 ： scale_min=1， scale_max=11， scale_step=2 ； 0013 （2） 对 于 第 k 个 空 间 尺 度 scalek=scale_min+scale_step*(k-1)，其 中， k=1,2,(scale_max-scale_min)/scale_step+1， 在原始遥感影像上设置一个。
14、大小为 说 明 书 CN 103761717 A 4 2/5 页 5 (2*scalek+1)(2*scalek+1) 的窗口， 不断移动该窗口， 对原始图像 im 中第 i 行第 j 列的 像元 im(i,j)， 窗口中的所有像元光谱值记为 im(i-scalek:i-scalek,j-scalek:j-scalek)， 然后计算窗口中所有像元光谱值的方差 imvar(i,j,k) 0014 0015 0016 上述公式中， ii 代表某个像元在窗口中的行号， jj 代表某个像元在窗口中的列 号， 将imvar(i,j,k)作为该窗口中心像元im(i,j)在空间尺度scalek下的方差， 并。
15、按照上述 方法计算处全图中所有像元的方差， 并求其平均值 imvar0 0017 （3） 比较 k 个空间尺度下的 imvar0， 选择 imvar0最大时所对应的空间尺度 scalebest 作为遥感影像的最优空间尺度 ； 0018 进一步的， 上述步骤 （3） 中遥感影像均变方差纹理图的具体生成方法如下 ： 0019 （1） 依据上述步骤得到的最优空间尺度 scalebest， 在原始遥感影像上设置一个大 小为 (2*scalebest+1)(2*scalebest+1) 的窗口 ； 0020 （2） 不断移动该窗口， 对于像元 im(i,j)， 窗口中的所有像元光谱值记为 im(i-sc。
16、a lebest:i-scalebest,j-scalebest:j-scalebest)， 每个窗口中的像元总个数记作 N=(2*scalebest+ 1)*(2*scalebest+1)， 将其按照列方向排列记作 imN， 将窗口中所有像元的行坐标和列坐标依 次记作 rowN和 colN； 0021 （3） 依此对每个窗口中的所有像元的光谱值进行最小二乘平面拟合， 拟合方程为 ： Ax+By+Cz+1=0， 将待拟合的 N 个点表示成以下矩阵形式 ： 0022 0023 依据最小二乘法则得到平面拟合系数为 ： 0024 0025 其中， h 代表 N 个像元中的第 h 个像元 ； 0026。
17、 （4） 计 算 窗 口 中 所 有 像 元 光 谱 值 至 拟 合 平 面 的 距 离 ： 其中abs(*)是绝对值算子， 计算 该窗口中的各个像元所对应的距离， 求得方差 dvar 说 明 书 CN 103761717 A 5 3/5 页 6 0027 0028 0029 将dvar作为该窗口中心像元的均变方差纹理， 遍历所有点， 即生成均变方差纹理图 imvar_texture； 0030 遍历所有像素点， 搜索均变方差纹理图 imvar_texture的最大值和最小值， 分别记作 imvar_texture_max 和 imvar_texture_min， 对每个像元进行如下操作用于将。
18、纹理值拉伸至 0,255 ： 0031 0032 进一步的， 上述步骤 （4） 中最优阈值分割的具体步骤如下 ： 0033 （1） 在所得的方差纹理图中， 确定一个分割阈值 ； 0034 （2） 在步骤 （1） 中确定的分割阈值中， 根据实际情况确定一个面积阈值用于抑制微 小斑块， 依次遍历每一个斑块， 计算其面积， 如果该斑块的面积小于该阈值， 则去除该斑块， 否则保留该斑块 ； 0035 （3） 基于上述结果， 进行闭运算， 用于填充水体中的微小空洞， 至此， 已经得到最终 的水体信息， 如果有必要需要显示面积较大的水体， 可以再设置一个面积阈值， 保留面积大 于给定阈值的水体区域。 00。
19、36 有益效果 ： 本发明的一种基于全色遥感影像的水体提取方法， 基于水体因其流动 性而造成的水体区域光谱响应局部变化率较小的原理， 并提出均变方差纹理的概念用以度 量光谱值局部变化率的均一性， 充分利用纹理信息用以消除和水体相似光谱地物诸如阴影 和居民地等的影响， 使得水体部分的纹理特征和其他区域的纹理值之间差异更大， 最后利 用最优阈值分割方法高效快速地提取水体信息， 为城市水体信息的提取提出了一种高效、 快速的新方法。本发明计算方法简单， 运算复杂度较低， 并且最终提取结果较为精准， 能够 应用于城市规划、 环境科学、 地理信息制图等多个领域。 附图说明 0037 图 1 为本发明的流程。
20、图 ； 0038 图 2 为本发明中提取的水体的原始遥感影像 ； 0039 图 3 为图 1 的均变方差纹理图 ； 0040 图 4 为图 3 的分割结果示意图 ； 0041 图 5 为图 4 去除微小噪声后的结果示意图 ； 0042 图 6 为对图 5 结果进行闭运算得到的结果示意图 ； 0043 图 7 为对图 6 中较大面积水体的显示结果示意图。 具体实施方式 说 明 书 CN 103761717 A 6 4/5 页 7 0044 下面对本发明技术方案进行详细说明， 但是本发明的保护范围不局限于所述实施 例。 0045 下面对本发明技术方案结合附图进行详细说明。 0046 本发明的一种基。
21、于全色遥感影像的城市水体提取方法， 包括以下步骤 ： 0047 （1） 遥感影像预处理 ： 利用中值滤波抑制遥感影像中的噪声 ； 0048 （2） 最优尺度选取 ： 在经步骤 （1） 预处理的图像中， 选定最小尺度、 最大尺度以及 尺度变化步长， 依次计算全局平均方差， 求得最优空间尺度 ； 0049 （3） 均变纹理生成 ： 依据步骤 （2） 中所得的最优空间尺度， 遍历遥感影像中的各个 像素， 分别计算均变方差， 进而生成遥感影像的均变方差纹理图 ； 0050 （4） 最优阈值分割 ： 基于步骤 （3） 中所得遥感影像的均变方差纹理图， 人工选定分 割阈值， 对方差纹理图进行分割， 选取大。
22、于给定面积阈值的斑块， 进行闭运算， 获得最终水 体信息。 0051 本发明中步骤 （2） 中确定最优尺度的详细算法如下 ： 0052 （1）确定最小尺度 scale_min， 最大尺度 scale_max 以及尺度变化步长 scale_ step ； 0053 （2） 对 于 第 k 个 空 间 尺 度 scalek=scale_min+scale_step*(k-1)，其 中， k=1,2,(scale_max-scale_min)/scale_step+1， 在原始遥感影像上设置一个大小为 (2*scalek+1)(2*scalek+1) 的窗口， 不断移动该窗口， 对于像元 im(i,。
23、j)， 窗口中的所有 像元光谱值记为 im(i-scalek:i-scalek,j-scalek:j-scalek)， 然后计算窗口中所有像元光 谱值的方差 imvar(i,j,k) 0054 0055 0056 将 imvar(i,j,k) 作为该窗口中心像元 im(i,j) 在空间尺度 scalek下的方差， 并按 照上述方法计算处全图中所有像元的方差， 并求其平均值 imvar0 0057 （3） 比较 k 个空间尺度下的 imvar0， 选择 imvar0最大时所对应的空间尺度 scalebest 作为遥感影像的最优空间尺度 ； 0058 本发明中步骤 （3） 遥感影像均变方差纹理图的。
24、具体生成方法如下 ： 0059 （1） 依据上述步骤得到的最优空间尺度 scalebest， 在原始遥感影像上设置一个大 小为 (2*scalebest+1)(2*scalebest+1) 的窗口 ； 0060 （2） 不断移动该窗口， 对于像元 im(i,j)， 窗口中的所有像元光谱值记为 im(i-sca lebest:i-scalebest,j-scalebest:j-scalebest)， 每个窗口中的像元总个数记作 N=(2*scalebest+ 1)*(2*scalebest+1)， 将其按照列方向排列记作 imN， 将窗口中所有像元的行坐标和列坐标依 次记作 rowN和 colN。
25、； 0061 （3） 依此对每个窗口中的所有像元的光谱值进行最小二乘平面拟合， 拟合方程为 ： 说 明 书 CN 103761717 A 7 5/5 页 8 Ax+By+Cz+1=0， 将待拟合的 N 个点表示成以下矩阵形式 ： 0062 0063 依据最小二乘法则得到平面拟合系数为 ： 0064 0065 计 算 窗 口 中 所 有 像 元 光 谱 值 至 拟 合 平 面 的 距 离 ： 0066 对该窗口中的所有距离求得方差 dvar 0067 0068 0069 将dvar作为该窗口中心像元的均变方差纹理， 遍历所有点， 即生成均变方差纹理图 imvar_texture； 0070 遍历。
26、所有像素点， 搜索均变方差纹理图 imvar_texture的最大值和最小值， 分别记作 imvar_texture_max 和 imvar_texture_min， 对每个像素点作如下操作 ： 0071 0072 本发明中步骤 （4） 中最优阈值分割的具体步骤如下 ： 0073 （1） 在所得的方差纹理图中， 确定一个分割阈值 ； 0074 （2） 在步骤 （1） 中确定的分割阈值中， 根据实际情况确定一个面积阈值用于抑制微 小斑块， 依次遍历每一个斑块， 计算其面积， 如果该斑块的面积小于该阈值， 则去除该斑块， 否则保留该斑块 ； 0075 （3） 基于上述结果， 进行闭运算， 用于填充。
27、水体中的微小空洞， 至此， 已经得到最终 的水体信息， 如果有必要需要显示面积较大的水体， 可以再设置一个面积阈值， 保留面积大 于给定阈值的水体区域。 说 明 书 CN 103761717 A 8 1/7 页 9 图 1 说 明 书 附 图 CN 103761717 A 9 2/7 页 10 图 2 说 明 书 附 图 CN 103761717 A 10 3/7 页 11 图 3 说 明 书 附 图 CN 103761717 A 11 4/7 页 12 图 4 说 明 书 附 图 CN 103761717 A 12 5/7 页 13 图 5 说 明 书 附 图 CN 103761717 A 13 6/7 页 14 图 6 说 明 书 附 图 CN 103761717 A 14 7/7 页 15 图 7 说 明 书 附 图 CN 103761717 A 15 。