CN104794724A - 一种基于非线性尺度滤波的遥感影像建筑物提取方法 - Google Patents
一种基于非线性尺度滤波的遥感影像建筑物提取方法 Download PDFInfo
- Publication number
- CN104794724A CN104794724A CN201510220753.8A CN201510220753A CN104794724A CN 104794724 A CN104794724 A CN 104794724A CN 201510220753 A CN201510220753 A CN 201510220753A CN 104794724 A CN104794724 A CN 104794724A
- Authority
- CN
- China
- Prior art keywords
- band
- sensing image
- remote
- remote sensing
- scale filtering
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Landscapes
- Image Processing (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明涉及一种基于非线性尺度滤波的遥感影像建筑物提取方法。包括如下步骤:步骤1,分离影像波段;步骤2,非线性尺度滤波;步骤3,判断迭代停止条件;步骤4,二值化处理;步骤5,波段叠加。可克服目前高空间分辨率遥感影像中建筑物提取困难的问题,充分利用遥感影像的光谱特征,可以检测遥感影像中的建筑物目标,且准确地保留其边缘,无需人工干预,自动化程度高。
Description
技术领域
本发明涉及一种遥感影像处理领域,具体说是一种基于非线性尺度滤波的遥感影像建筑物提取方法。
背景技术
建筑物是城市主要地理要素之一,是各种城市专题图的重要内容,研究建筑物的提取对综合考察城市地理信息环境具有重要意义。随着高分辨率遥感影像获取技术的快速发展,遥感影像的处理、分析和应用有了更好的数据源,其数字产品则有了更广泛、更深入的应用。计算机图像处理技术、模式识别、人工智能等方面的都取得不同程度的进展,为高效地提取海量影像中的有效信息提供了可能。但建筑物信息的提取要比其他信息如道路、水体的获取难得多,主要原因如下:
(1)数据源主要是二维的遥感影像,大多数情况下缺少直接的三维数据;
(2)不同的遥感影像常因为光谱范围、分辨率、传感器的几何图像以及成像条件等因素的不同而有较大的差异;
(3)不同种类的建筑物其所表现出来的外观和纹理细节等千变万化,表现在遥感图像上差异很大,统一的建筑物模型库难以建立,这使得信息的自动提取变得相当困难;
(4)建筑物所处场景的复杂性,如对比度较低时、房屋相互遮挡、建筑物自身的阴影以及处在其它地物的阴影等,所以想自动地从背景中提取出边界清晰的建筑物较为困难。
发明内容
本发明提供了一种基于非线性尺度滤波的遥感影像建筑物提取方法,可克服目前高空间分辨率遥感影像中建筑物提取困难的问题,充分利用遥感影像的光谱特征,可以检测遥感影像中的建筑物目标,且准确地保留其边缘,无需人工干预,自动化程度高。
为实现本发明的目标所采用的技术方案是:方法包括以下步骤:
步骤1:分离输入遥感影像image1的各个波段,分别记为:band1,band2,…,bandn,其中n为image1中的波段总数;
步骤2:利用各向异性的热传递方程迭代计算每个波段band1,band2,…,bandn中每个像元的取值,得到每个波段的滤波结果,分别记为:F_band1,F_band2,…,F_bandn;
步骤3:对F_band1,F_band2,…,F_bandn进行二值化处理,分别得到:BW_F_band1,BW_F_band2,…,BW_F_bandn;
步骤4:叠加BW_F_band1,BW_F_band2,…,BW_F_bandn。
所述的各向异性的热传递方程为:
It+1=It+λ(cNx,y▽N(It)+cSx,y▽S(It)+cEx,y▽E(It)+cWx,y▽W(It))
k和λ均为平滑控制系数,控制平滑程度越大,越平滑,越不易保留边缘。
▽k(Ix,y)表示在四个方向上对当前像素求偏导:
▽N(Ix,y)=Ix,y-1-Ix,y
▽S(Ix,y)=Ix,y+1-Ix,y
▽E(Ix,y)=Ix-1,y-Ix,y
▽W(Ix,y)=Ix+1,y-Ix,y
若▽k(Ix,y)越大,则变化较大,就表示边界。
cKx,y表示四个方向上的导热系数:
cNx,y=e-(||▽N(I)||2/k2)
cSx,y=e-(||▽S(I)||2/k2)
cEx,y=e-(||▽E(I)||2/k2)
cWx,y=e-(||▽W(I)||2/k2)
方向变化越大,求得的cKx,y的取值越小,从而达到保留边界的目的。
所述的各向异性的热传递方程迭代的停止条件为对应图像的标准差曲线出现第一个谷点。
所述的二值化处理使用的阈值由OSTU方法自动确定。
本发明的有益效果是:解决了高空间分辨率遥感影像中建筑物边缘检测困难的问题,达到完全自动化的效果。可以用于遥感影像制图、地理信息系统的数据获取和自动更新。
附图说明
图1是本发明的总体处理流程图。
具体实施方式
下面结合附图详细描述本发明的具体实施方式。
在步骤101,输入的待处理遥感影像image1为包含3个波段的高空间分辨率遥感影像,传感器为Quick bird,且已经进行辐射校正和几何校正等预处理。
在步骤102,分离输入遥感影像image1的各个波段,分别得到:band1,band2和band3。
在步骤103,利用各向异性的热传递方程分别对band1,band2和band3进行非线性尺度滤波,迭代计算每个波段band1,band2和band3中每个像元的取值,得到每个波段的滤波结果,分别记为:F_band1,F_band2和F_band3。
在步骤104,分别计算F_band1,F_band2和F_band3中全部像素取值的标准差,如果某一波段的标准差取值从小变大,则其满足迭代停止条件,进入步骤105;否则进入步骤103。
在步骤105,通过OSTU方法自动确定阈值,对F_band1,F_band2和F_band3进行二值化处理,得到BW_F_band1,BW_F_band2和BW_F_band3。
在步骤106,叠加BW_F_band1,BW_F_band2和BW_F_band3。
在步骤107,得到边界准确的建筑物对象。
Claims (4)
1.一种基于非线性尺度滤波的遥感影像建筑物提取方法,其特征在于包括以下步骤:
步骤1:分离输入遥感影像image1的各个波段,分别记为:band1,band2,…,bandn,其中n为image1中的波段总数;
步骤2:利用各向异性的热传递方程迭代计算每个波段band1,band2,…,bandn中每个像元的取值,得到每个波段的滤波结果,分别记为:F_band1,F_band2,…,F_bandn;
步骤3:对F_band1,F_band2,…,F_bandn进行二值化处理,分别得到:BW_F_band1,BW_F_band2,…,BW_F_bandn;
步骤4:叠加BW_F_band1,BW_F_band2,…,BW_F_bandn。
2.根据权利要求1所述的一种基于非线性尺度滤波的遥感影像建筑物提取方法,其特征在于各向异性的热传递方程为:
It+1=It+λ(cNx,y▽N(It)+cSx,y▽S(It)+cEx,y▽E(It)+cWx,y▽W(It))
k和λ均为平滑控制系数,控制平滑程度越大,越平滑,越不易保留边缘。
▽k(Ix,y)表示在四个方向上对当前像素求偏导:
▽N(Ix,y)=Ix,y-1-Ix,y
▽S(Ix,y)=Ix,y+1-Ix,y
▽E(Ix,y)=Ix-1,y-Ix,y
▽W(Ix,y)=Ix+1,y-Ix,y
若▽k(Ix,y)越大,则变化较大,就表示边界。
cKx,y表示四个方向上的导热系数:
方向变化越大,求得的cKx,y的取值越小,从而达到保留边界的目的。
3.根据权利要求1所述的一种基于非线性尺度滤波的遥感影像建筑物提取方法,其特征在于各向异性的热传递方程迭代的停止条件为对应图像的标准差曲线出现第一个谷点。
4.根据权利要求1所述的一种基于非线性尺度滤波的遥感影像建筑物提取方法,其特征在于二值化处理使用的阈值由OSTU方法自动确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510220753.8A CN104794724A (zh) | 2015-05-04 | 2015-05-04 | 一种基于非线性尺度滤波的遥感影像建筑物提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510220753.8A CN104794724A (zh) | 2015-05-04 | 2015-05-04 | 一种基于非线性尺度滤波的遥感影像建筑物提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104794724A true CN104794724A (zh) | 2015-07-22 |
Family
ID=53559504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510220753.8A Pending CN104794724A (zh) | 2015-05-04 | 2015-05-04 | 一种基于非线性尺度滤波的遥感影像建筑物提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104794724A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102032903A (zh) * | 2010-12-01 | 2011-04-27 | 南京大学 | 基于Landsat数据源的珊瑚岛礁遥感信息自动提取方法 |
CN102509262A (zh) * | 2011-10-17 | 2012-06-20 | 中煤地航测遥感局有限公司 | 一种遥感图像薄云去除方法 |
CN103955896A (zh) * | 2014-04-22 | 2014-07-30 | 中国资源卫星应用中心 | 一种提升卫星图像视觉效果的真彩色增强方法 |
-
2015
- 2015-05-04 CN CN201510220753.8A patent/CN104794724A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102032903A (zh) * | 2010-12-01 | 2011-04-27 | 南京大学 | 基于Landsat数据源的珊瑚岛礁遥感信息自动提取方法 |
CN102509262A (zh) * | 2011-10-17 | 2012-06-20 | 中煤地航测遥感局有限公司 | 一种遥感图像薄云去除方法 |
CN103955896A (zh) * | 2014-04-22 | 2014-07-30 | 中国资源卫星应用中心 | 一种提升卫星图像视觉效果的真彩色增强方法 |
Non-Patent Citations (3)
Title |
---|
PIETRO PERONA等: "Scale-Space and Edge Detection Using Anisotropic Diffusion", 《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELIGENCE》 * |
包文林等: "《模型黄土高原建设方略纵论》", 1 December 2005 * |
杨先武等: "基于高通滤波的高分辨率遥感影像城市道路提取", 《测绘科学》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Isikdogan et al. | RivaMap: An automated river analysis and mapping engine | |
Boonpook et al. | Deep learning-based multi-feature semantic segmentation in building extraction from images of UAV photogrammetry | |
Du et al. | Analysis of Landsat-8 OLI imagery for land surface water mapping | |
Ramiya et al. | Segmentation based building detection approach from LiDAR point cloud | |
CN108369651B (zh) | 一种提取天空区域的方法、系统及非暂时性计算机可读存储介质 | |
Yue et al. | A new accuracy evaluation method for water body extraction | |
GB2513702A (en) | Methods and systems for automatic and semi-automatic geometric and geographic feature extraction | |
Liu et al. | Interactive geospatial object extraction in high resolution remote sensing images using shape-based global minimization active contour model | |
Zani et al. | Multitemporal Landsat data applied for deciphering a megafan in northern Amazonia | |
Liu et al. | Using SPOT 5 fusion-ready imagery to detect Chinese tamarisk (saltcedar) with mathematical morphological method | |
Li et al. | Automatic cloud-shadow removal from flood/standing water maps using MSG/SEVIRI imagery | |
CN103218811A (zh) | 一种基于统计分布的卫星多光谱图像波段配准方法 | |
ZhiYong et al. | A modified mean filter for improving the classification performance of very high-resolution remote-sensing imagery | |
CN109165553B (zh) | 基于深度学习的石油设施遥感自动识别方法及装置 | |
CN104050674A (zh) | 一种显著性区域检测方法及装置 | |
Yang et al. | Improved multi-scale image matching approach for monitoring Amery ice shelf velocity using Landsat 8 | |
Bian et al. | Individual tree delineation in deciduous forest areas with LiDAR point clouds | |
CN110738134B (zh) | 面向无人机可见光影像的土壤信息提取方法及装置 | |
Tanathong et al. | Using GPS/INS data to enhance image matching for real-time aerial triangulation | |
Sharma et al. | A new water-resistant snow index for the detection and mapping of snow cover on a global scale | |
CN104778714A (zh) | 一种基于阴影和分割的高空间分辨率遥感影像建筑物检测方法 | |
CN105761250A (zh) | 一种基于模糊场景分割的建筑物提取方法 | |
CN104794724A (zh) | 一种基于非线性尺度滤波的遥感影像建筑物提取方法 | |
WO2014165787A1 (en) | System and method for detection of structures | |
CN103955936A (zh) | 一种基于堆栈式去噪自编码机的显著性物体检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150722 |
|
WD01 | Invention patent application deemed withdrawn after publication |