CN103745453A - 基于Google Earth遥感影像的城镇信息提取方法 - Google Patents
基于Google Earth遥感影像的城镇信息提取方法 Download PDFInfo
- Publication number
- CN103745453A CN103745453A CN201310676619.XA CN201310676619A CN103745453A CN 103745453 A CN103745453 A CN 103745453A CN 201310676619 A CN201310676619 A CN 201310676619A CN 103745453 A CN103745453 A CN 103745453A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- image
- pixel
- marginal density
- sigma
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
本发明公开一种基于Google Earth遥感影像的城镇信息提取方法,首先对遥感影像进行中值滤波等预处理,再根据sobel边缘检测算子生成影像的边缘图像,然后构造高斯空间域加权边缘密度生成算子,生成边缘密度图像,再对边缘密度图像进行OSTU最优阈值分割,最后对分割图进行后处理,最终获取遥感影像中的城镇信息。本发明的提取方法效果好、速度快、准确率高,适用于城镇尤其是城郊结合处的城镇信息的提取。
Description
技术领域
本发明涉及一种城镇信息提取方法,具体涉及一种基于Google Earth遥感影像的城镇信息提取方法。
背景技术
城镇作为独立于其周边农用地、植被、裸地或者水体的地物,对于土地规划、土地覆盖、利用制图都有着重要意义。二十一世纪以来,随着空间技术的发展和各类图像处理、模式识别算法的出现,基于遥感手段的城镇信息提取算法越来越丰富,在各项制图工作和研究中取得了广泛应用。目前,常用的城镇信息提取算法都是基于多光谱遥感影像如SPOT和TM遥感影像,利用特定波段的光谱值生成归一化建筑物指数(NDBI)或其改进形式,再基于监督分类或者阈值分割方法进行城镇信息提取,但是鲜见将Google Earth遥感影像作为数据源或使用空间结构信息的提取方法。
本发明针对Google Earth遥感影像中城镇区域和其他区域在光谱空间排列结构上的差异而设计了本算法,首先对遥感影像进行预处理,再通过sobel边缘检测算子生成边缘图像,然后根据空间域高斯边缘密度生成算子生成边缘密度图像,最后通过0STU阈值分割方法完成城镇信息提取。
关于城镇信息提取问题,鲜见基于空间结构化信息的方法,一般仅利用了光谱信息,未能考虑到空间地物丰富的结构信息,故未能够最大化利用遥感影像中丰富的信息;其次,Google Earth遥感影像获取容易,但针对这一数据源开展的研究并不多,因此有必要针对该数据源做相应的研究工作。
综上,提出基于Google Earth遥感影像结构信息的城镇信息提取算法非常有意义,也将为相近研究提供一种新思路。
发明内容
发明目的:本发明的目的是针对现有城镇信息提取技术中存在未利用空间信息和Google Earth遥感影像在该领域利用较少的情况,提供一种基于GoogleEarth遥感影像的城镇信息提取方法。
技术方案:本发明是一种基于Google Earth遥感影像的城镇信息提取方法, 包括以下步骤:
(1)将原始RGB图像转化为灰度图像,利用中值滤波算法进行图像去噪,中值滤波模板大小取为3×3,以每个像元为中心,取其3×3邻域内的9个像元灰度值的中值作为该像元的滤波结果;
(2)选择sobel边缘检测算子,使用横向和纵向两个模板对灰度图像进行卷积得到横向和纵向两个方向的梯度图像,累加两幅图像得到原始灰度图的近似梯度图像,通过下述公式梯度阈值g_thrd:g_thrd=k*mean(imgradient),其中imgradient为梯度图像,mean(*)为梯度均值运算操作,k为阈值乘系数,乘系数一般取值范围为[0.5,2],例如可以选取k=1,大于该阈值的认为是图像边缘,以此生成边缘特征图像;
(3)在上述步骤(2)中所得的边缘特征图像中,根据空间地物在空间范围内的相关性,越靠近某像元的边缘对该像元的边缘密度值贡献越大,由于二维高斯函数的旋转对称不变形和单值性,因此二维高斯函数可以反应此相关性,生成高斯空间域加权边缘密度生成算子Bm×m,其中Bm×m中元素Bij表示如下: 根据上述算子Bm×m生成边缘密度图像,Bm×m的行数以及列数均为m;
(4)基于步骤(3)中的边缘密度图像,通过OSTU方法选定阈值,对边缘密度图像进行阈值分割,选取大于给定面积阈值的图像斑块,最后进行闭运算以获取最终城镇信息。
进一步的,所述步骤(3)中边缘密度生成算子的生成步骤如下:
1)确定边缘密度生成算子B的大小m和参数σ2:m代表算子空间作用范围,σ2用于度量算子中各元素随着该元素至中心像元距离递增时的递减能力,σ2越小则算子的递减能力越强,越能体现中心像元的重要性,σ2越大则算子的递减能力越弱,越能体现邻域像元的重要性,σ2的选择能体现中心像元和其邻域像 元在空间范围内的重要程度,m和σ2可以通过多次试验根据最终分割结果与目视解译结果进行比较择优确定,例如,m=7,σ2=64其中
上述i和j分别代表算子B中各元素对应的行号和列号;
2)计算归一化参数c,
3)根据归一化参数计算各个元素b(i,j),组成边缘密度生成算子B,
4)生成大小为M×N的初始边缘密度图像DM×N,计算各像元的边缘密度值,
其中,im为原始遥感影像灰度图,col代表像元在原始遥感影像灰度图中的列数,row代表像元在原始遥感影像灰度图中的行数,所得边缘密度图像为:
因为边缘密度生成算子大小为m×m,因此对于图像靠近边界的行或列是未进行边缘密度计算的,对于每一个未参与计算的像元,在参与计算的像元中,寻找和该像元行列坐标距离最为接近的像元,并将此像元的边缘密度值赋给该像元,以此填补所有未计算像元的边缘密度值,最终得到的边缘密度图像为:
5)遍历边缘密度图像中各像元,使用冒泡法搜索边缘密度图像中各边缘密度值的最大值和最小值,首先初始化边缘密度最大值和最小值:dmax=0,dmin=100000;然后依次将i行j列像元的边缘密度值dij分别与dmax及dmin进行比较,若dmax<dij,则将dij赋给dmax,若dmin>dij,则将dij赋给dmin,以此获取边缘密度最大值和最小值,对各像元边缘密度值线性拉伸至0-255,并依据四舍五入法取整:
dij=round(255((dij-dmin)/(dmax-dmin)))
上述公式中,dij为第i行第j列像元的边缘密度值,dmax为边缘密度最大值,dmin为边缘密度最小值,round(*)为四舍五入取整算子。
进一步的,所述步骤(4)具体包括以下步骤:
i)将边缘密度图像中的各个像元的边缘密度值按照一定间隔分为若干等级,即边缘密度图像的灰度级,然后进行排列,可以使用集合P={p1,p2,…,pL}来表示,其中L为边缘密度图像的灰度等级,由于边缘密度图像已经进行线性拉伸至[0,255],则本发明中,L=256,边缘密度值为pi的像元数量为ri,则总的像元数量为其中若L<256,则相当于对边缘密度图像进行灰度级压缩;
ii)以第k个边缘密度值pk为分界处,将集合p分为两类,分别为P0={p1,p2,…,pk}和P1={pk+1,pk+2,…,pL},(k=1,2,…,L),依次计算两个类别之间的类间方差,类间方差计算公式如下所示:
其中, Pi=ri/R;
有益效果:本发明提出一种基于Google Earth遥感影像的城镇信息提取算方法,能够方便快速地提取出城镇信息,用于大尺度城市制图和城市扩张研究等。本发明可以充分利用遥感影像的空间结构信息,本发明引入的空间域高斯加权边缘密度加权算子充分考虑了地物之间的空间相关性,使得各地物的边缘密度生成更为合理,最后通过OSTU阈值分割方法方便快速地提取城镇信息。
附图说明
图1是本发明实施例中采集的城镇信息原始遥感影像的灰度图像;
图2是图1转化为灰度图的梯度图像;
图3是由梯度图像转化得到的边缘图像;
图4是图3的边缘密度图像;
图5是图4的OSTU分割结果示意图;
图6是图5中的OSTU分割结果叠加于边缘图像获得的结果的示意图;
图7是本发明的城镇信息提取流程图。
具体实施方式
下面对本发明技术方案结合附图进行详细说明。
本发明的一种基于Google Earth遥感影像的城镇信息提取方法,包括以下步骤:
(1)将原图像转化为灰度图像,利用中值滤波算法去除原始遥感影像中存在的噪声;
(2)使用sobel边缘检测算子生成边缘特征图像;
(3)根据空间域高斯函数,生成一定大小的边缘密度生成算子,然后生成边缘密度图像;
(4)基于步骤(3)中的边缘密度图像,通过OSTU方法选定阈值,对边缘密度图像进行阈值分割,选取大于给定面积阈值的图像斑块,最后进行闭运算以获取最终城镇信息。
进一步的,所述步骤(3)中边缘密度生成算子的生成步骤如下:
1)确定边缘密度生成算子B的大小m和参数σ2,其中
2)计算归一化参数c,
3)根据归一化参数计算各个元素b(i,j),组成边缘密度生成算子B,
4)生成大小为M×N的初始边缘密度图像DM×N,计算各像元的边缘密度值,
其中,im为原始遥感影像灰度图,col代表像元在原始遥感影像灰度图中的列数,row代表像元在原始遥感影像灰度图中的行数,所得边缘密度图像为:
因为边缘密度生成算子大小为m×m,因此对于图像靠近边界的行或列是未进行边缘密度计算的,对于每一个未参与计算的像元,在参与计算的像元中,寻找和该像元行列坐标距离最为接近的像元,并将此像元的边缘密度值赋给该像元,以此填补所有未计算像元的边缘密度值,最终得到的边缘密度图像为:
5)遍历所有像素点,搜索边缘密度图像的最大值和最小值,对各个像元进行线性拉伸至0-255,并依据四舍五入法取整:
dij=round(255((dij-dmin)/(dmax-dmin)))
进一步的,所述步骤(4)具体包括以下步骤:
i)将边缘密度图像中的各个像元的边缘密度值按照一定间隔(例如:间隔为1)分为若干等级,即边缘密度图像的灰度级,然后进行排列,可以使用集合P={p1,p2,…,pL}来表示,其中L为边缘密度图像的灰度等级,由于边缘密度图像已经进行线性拉伸至[0,255],本实施例中,L=256,边缘密度值为pi的像元数量为ri,则总的像元数量为
ii)以第k个边缘密度值pk为分界处,将集合p分为两类,分别为P0={p1,p2,…,pk}和P1={pk+1,pk+2,…,pL},(k=l,2,…,L),依次计算两个类别之间的类间方差类间方差计算公式如下所示:
其中, Pi=ri/R。
如上所述,尽管参照特定的优选实施例已经表示和表述了本发明,但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下,可对其在形式上和细节上作出各种变化。
Claims (3)
1.一种基于Google Earth遥感影像的城镇信息提取方法,其特征在于包括以下步骤:
(1)将原始RGB图像转化为灰度图像,利用中值滤波算法进行图像去噪,中值滤波模板大小取为3×3,以每个像元为中心,取其3×3邻域内的9个像元灰度值的中值作为该像元的滤波结果;
(2)选择sobel边缘检测算子,使用横向和纵向两个模板对灰度图像进行卷积得到横向和纵向两个方向的梯度图像,累加两幅图像得到原始灰度图的近似梯度图像,通过下述公式选取梯度阈值g_thrd:g_thrd=k*mean(imgradient),其中imgradient为梯度图像,mean(*)为梯度均值运算操作,k为阈值乘系数,乘系数的取值范围为[0.5,2];
(3)在上述步骤(2)中所得的边缘特征图像中,根据空间地物在空间范围内的相关性,越靠近某像元的边缘对该像元的边缘密度值贡献越大,由于二维高斯函数的旋转对称不变形和单值性,因此二维高斯函数可以反应此相关性,生成高斯空间域加权边缘密度生成算子Bm×m,其中Bm×m中元素Bij表示如下: 根据上述算子Bm×m生成边缘密度图像,Bm×m的行数以及列数均为m;
(4)基于步骤(3)中生成的边缘密度图像,通过OSTU方法确定最优阈值,对边缘密度图像进行阈值分割,选取大于给定面积阈值的图像斑块,最后进行闭运算以获取最终城镇信息。
2.根据权利要求1所述的基于Google Earth遥感影像的城镇信息提取方法,其特征在于:所述步骤(3)中边缘密度生成算子的生成步骤如下:
1)确定边缘密度生成算子B的大小m和参数σ2:m代表算子空间作用范围,σ2用于度量算子中各元素随着该元素至中心像元距离递增时的递减能力,m和σ2由最终分割结果与目视解译结果确定,其中
上述i和j分别代表算子B中各元素对应的行号和列号;
2)根据下式计算归一化参数c:
3)根据归一化参数计算各个元素b(i,j),组成边缘密度生成算子B:
4)生成大小为M×N的初始边缘密度图像DM×N,计算各像元的边缘密度值:
其中,im为原始遥感影像灰度图,row代表像元在原始遥感影像灰度图中的行数,col代表元在原始遥感影像灰度图中的列数,所得边缘密度图像为:
因为边缘密度生成算子大小为m×m,其中m为奇数,因此对于图像靠近边界的行或列是未进行边缘密度计算的,对于每一个未参与计算的像元,在参与计算的像元中,寻找和该像元行列坐标的欧氏距离最为接近的像元,并将此像元的边缘密度值赋给该像元,以此填补所有未计算像元的边缘密度值,最终得到的边缘密度图像为:
5)遍历边缘密度图像中各像元,使用冒泡法搜索边缘密度图像中各边缘密度值的最大值和最小值,首先初始化边缘密度最大值和最小值:dmax=0,dmin=100000;然后依次将i行j列像元的边缘密度值dij分别与dmax及dmin进行比较,若dmax<dij,则将dij赋给dmax,若dmin>dij,则将dij赋给dmin,以此获取边缘密度最大值和最小值,对各像元边缘密度值线性拉伸至0-255,并依据四舍五入法取整:
dij=round(255((dij-dmin)/(dmax-dmin)))
上述公式中,dij为第i行第j列像元的边缘密度值,dmax为边缘密度最大值,dmin为边缘密度最小值,round(*)为四舍五入取整算子。
3.根据权利要求1所述的基于Google Earth遥感影像的城镇信息提取方法,其特征在于所述步骤(4)具体包括以下步骤:
i)按照灰度间隔为1的标准,将边缘密度图像中的各像元边缘密度值分为若干等级,即边缘密度图像的灰度级,然后对灰度级进行排列,使用集合P={p1,p2,…,pL}来表示,其中L为边缘密度图像的灰度等级,由于边缘密度图像已经进行线性拉伸至[0,255],则L=256,记边缘密度值为pi的像元数量为ri,则总的像元数量为
ii)以第k个边缘密度值pk为分界处,将集合p分为两类,分别为P0={p1,p2,…,pk}和P1={pk+1,pk+2,…,pL},(k=1,2,…,L),依次计算两个类别之间的类间方差,类间方差计算公式如下所示:
其中, Pi=ri/R;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310676619.XA CN103745453B (zh) | 2013-12-11 | 2013-12-11 | 基于Google Earth遥感影像的城镇信息提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310676619.XA CN103745453B (zh) | 2013-12-11 | 2013-12-11 | 基于Google Earth遥感影像的城镇信息提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103745453A true CN103745453A (zh) | 2014-04-23 |
CN103745453B CN103745453B (zh) | 2016-08-17 |
Family
ID=50502468
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310676619.XA Expired - Fee Related CN103745453B (zh) | 2013-12-11 | 2013-12-11 | 基于Google Earth遥感影像的城镇信息提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103745453B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794478A (zh) * | 2015-05-04 | 2015-07-22 | 福建师范大学 | 一种用于遥感影像中具有均匀光谱特性的建筑物提取方法 |
CN108596103A (zh) * | 2018-04-26 | 2018-09-28 | 吉林大学 | 基于最佳光谱指数选择的高分辨率卫星遥感影像建筑物提取方法 |
CN109740489A (zh) * | 2018-12-27 | 2019-05-10 | 核工业北京地质研究院 | 一种利用casi图像识别赤铁矿的方法 |
CN110070545A (zh) * | 2019-03-20 | 2019-07-30 | 重庆邮电大学 | 一种城镇纹理特征密度自动提取城镇建成区的方法 |
CN110570427A (zh) * | 2019-07-19 | 2019-12-13 | 武汉珈和科技有限公司 | 一种融合边缘检测的遥感影像语义分割方法及装置 |
WO2020001631A1 (zh) * | 2018-06-29 | 2020-01-02 | 长城汽车股份有限公司 | 基于视觉摄像机的自阴影物体边缘识别方法、装置及车辆 |
CN114266138A (zh) * | 2021-11-29 | 2022-04-01 | 西南大学 | 一种利用云端数据进行城市边缘区识别与验证方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6738512B1 (en) * | 2000-06-19 | 2004-05-18 | Microsoft Corporation | Using shape suppression to identify areas of images that include particular shapes |
CN102945374A (zh) * | 2012-10-24 | 2013-02-27 | 北京航空航天大学 | 一种高分辨率遥感图像中民航飞机自动检测方法 |
-
2013
- 2013-12-11 CN CN201310676619.XA patent/CN103745453B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6738512B1 (en) * | 2000-06-19 | 2004-05-18 | Microsoft Corporation | Using shape suppression to identify areas of images that include particular shapes |
CN102945374A (zh) * | 2012-10-24 | 2013-02-27 | 北京航空航天大学 | 一种高分辨率遥感图像中民航飞机自动检测方法 |
Non-Patent Citations (5)
Title |
---|
BERIL SIRMACEK 等: "Urban Area Detection Using Local Feature Points and Spatial Voting", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
JUNPING ZHANG 等: "Multi-scale Segmentation in Change Detection for Urban High Resolution Images", 《GEOSCIENCE AND REMOTE SENSING SYMPOSIUM (IGARSS), 2011 IEEE INTERNATIONAL》 * |
张军 等: "基于边缘增强的遥感图像变化检测技术", 《计算机工程与应用》 * |
李启青 等: "基于遥感数据光谱和空间特征的边缘提取方法", 《计算机应用》 * |
李文庆 等: "基于Google Earth的ETM+遥感图像自动分类方法", 《江西农业学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794478A (zh) * | 2015-05-04 | 2015-07-22 | 福建师范大学 | 一种用于遥感影像中具有均匀光谱特性的建筑物提取方法 |
CN104794478B (zh) * | 2015-05-04 | 2017-12-19 | 福建师范大学 | 一种用于遥感影像中具有均匀光谱特性的建筑物提取方法 |
CN108596103A (zh) * | 2018-04-26 | 2018-09-28 | 吉林大学 | 基于最佳光谱指数选择的高分辨率卫星遥感影像建筑物提取方法 |
CN108596103B (zh) * | 2018-04-26 | 2021-03-19 | 吉林大学 | 基于最佳光谱指数选择的高分辨遥感影像建筑物提取方法 |
WO2020001631A1 (zh) * | 2018-06-29 | 2020-01-02 | 长城汽车股份有限公司 | 基于视觉摄像机的自阴影物体边缘识别方法、装置及车辆 |
CN109740489A (zh) * | 2018-12-27 | 2019-05-10 | 核工业北京地质研究院 | 一种利用casi图像识别赤铁矿的方法 |
CN110070545A (zh) * | 2019-03-20 | 2019-07-30 | 重庆邮电大学 | 一种城镇纹理特征密度自动提取城镇建成区的方法 |
CN110070545B (zh) * | 2019-03-20 | 2023-05-26 | 重庆邮电大学 | 一种城镇纹理特征密度自动提取城镇建成区的方法 |
CN110570427A (zh) * | 2019-07-19 | 2019-12-13 | 武汉珈和科技有限公司 | 一种融合边缘检测的遥感影像语义分割方法及装置 |
CN114266138A (zh) * | 2021-11-29 | 2022-04-01 | 西南大学 | 一种利用云端数据进行城市边缘区识别与验证方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103745453B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103745453B (zh) | 基于Google Earth遥感影像的城镇信息提取方法 | |
Li et al. | A review of remote sensing image classification techniques: The role of spatio-contextual information | |
Huang et al. | A new building extraction postprocessing framework for high-spatial-resolution remote-sensing imagery | |
Yin et al. | Hot region selection based on selective search and modified fuzzy C-means in remote sensing images | |
CN105528619B (zh) | 基于小波变换和svm的sar遥感影像变化检测方法 | |
CN103473764B (zh) | 一种遥感影像目标变化检测方法 | |
CN103871039B (zh) | 一种sar图像变化检测差异图生成方法 | |
CN102831598B (zh) | 多分辨率NMF和Treelet融合的遥感图像变化检测方法 | |
CN109583321A (zh) | 一种基于深度学习的结构化道路中小物体的检测方法 | |
CN102426700A (zh) | 基于局部和全局区域信息的水平集sar图像分割方法 | |
CN106156758B (zh) | 一种sar海岸图像中海岸线提取方法 | |
CN103824302B (zh) | 基于方向波域图像融合的sar图像变化检测方法 | |
Zhang et al. | Multi-focus image fusion based on non-negative matrix factorization and difference images | |
Kalkan et al. | Comparison of support vector machine and object based classification methods for coastline detection | |
Lu et al. | Cascaded multi-task road extraction network for road surface, centerline, and edge extraction | |
CN105138992A (zh) | 一种基于区域主动轮廓模型的海岸线检测方法 | |
CN105447452A (zh) | 一种基于地物空间分布特征的遥感亚像元制图方法 | |
CN103761717A (zh) | 一种基于全色遥感影像的城市水体提取方法 | |
CN105469428B (zh) | 一种基于形态学滤波和svd的弱小目标检测方法 | |
Zhang et al. | Impervious surface extraction from high-resolution satellite image using pixel-and object-based hybrid analysis | |
CN104933719A (zh) | 一种积分图块间距离检测图像边缘方法 | |
CN106504219A (zh) | 有约束的路径形态学高分辨率遥感影像道路增强方法 | |
CN107341449A (zh) | 一种基于云块特征变化的静止气象卫星降水估算方法 | |
CN104200472B (zh) | 基于非局部小波信息的遥感图像变化检测方法 | |
CN116883679B (zh) | 基于深度学习的地物目标提取方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 Termination date: 20181211 |
|
CF01 | Termination of patent right due to non-payment of annual fee |