【航空和卫星图像中检测建筑物】使用gabor特征和概率的城市区域和建筑物检测研究(Matlab代码实现)

简介: 【航空和卫星图像中检测建筑物】使用gabor特征和概率的城市区域和建筑物检测研究(Matlab代码实现)

💥💥💞💞欢迎来到本博客❤️❤️💥💥


🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。


⛳️座右铭:行百里者,半于九十。


📋📋📋本文目录如下:🎁🎁🎁


目录


💥1 概述


📚2 运行结果


🎉3 参考文献


🌈4 Matlab代码实现


💥1 概述

文献来源:


623fd52e272148debdaca3570df4e48b.png


从甚高分辨率 (VHR) 航空和卫星图像检测建筑物在地图制作、城市规划和土地利用分析中非常有用。虽然可以从这些VHR图像中手动定位建筑物,但此操作可能不可靠且快速。因此,需要自动化系统从VHR航空和卫星图像中检测建筑物。不幸的是,这种系统必须处理重大问题。首先,建筑物具有不同的特征,它们的外观(照明,视角等)在这些图像中不受控制。其次,城市地区的建筑物通常密集而复杂。很难从它们中检测到单独的建筑物。为了克服这些困难,我们提出了一种使用局部特征向量和概率框架的新型建筑检测方法。我们首先介绍四种不同的局部特征向量提取方法。提取的局部特征向量用作要估计的概率密度函数 (pdf) 的观测值。使用可变核密度估计方法,我们估计相应的pdf。换句话说,我们将图像中的建筑物位置(待检测)表示为联合随机变量并估计它们的pdf。使用估计密度的模式以及其他概率属性,我们检测图像中的建筑物位置。我们还引入了基于概率框架的数据和决策融合方法来检测建筑物位置。我们挑选了VHR全色航空和Ikonos卫星图像的某些作物来测试我们的方法。我们假设这些作物是使用我们以前的城市区域检测方法检测的。我们的测试图像由两个不同的传感器采集,它们具有不同的空间分辨率。此外,这些图像中的建筑物具有不同的特征。因此,我们可以在不同的数据集上测试我们的方法。广泛的测试表明,我们的方法可用于在Ikonos卫星和我们的航空图像中以稳健和快速的方式自动检测建筑物。


原文摘要:


Abstract:

Detecting buildings from very high resolution (VHR) aerial and satellite images is extremely useful in map making, urban planning, and land use analysis. Although it is possible to manually locate buildings from these VHR images, this operation may not be robust and fast. Therefore, automated systems to detect buildings from VHR aerial and satellite images are needed. Unfortunately, such systems must cope with major problems. First, buildings have diverse characteristics, and their appearance (illumination, viewing angle, etc.) is uncontrolled in these images. Second, buildings in urban areas are generally dense and complex. It is hard to detect separate buildings from them. To overcome these difficulties, we propose a novel building detection method using local feature vectors and a probabilistic framework. We first introduce four different local feature vector extraction methods. Extracted local feature vectors serve as observations of the probability density function (pdf) to be estimated. Using a variable-kernel density estimation method, we estimate the corresponding pdf. In other words, we represent building locations (to be detected) in the image as joint random variables and estimate their pdf. Using the modes of the estimated density, as well as other probabilistic properties, we detect building locations in the image. We also introduce data and decision fusion methods based on our probabilistic framework to detect building locations. We pick certain crops of VHR panchromatic aerial and Ikonos satellite images to test our method. We assume that these crops are detected using our previous urban region detection method. Our test images are acquired by two different sensors, and they have different spatial resolutions. Also, buildings in these images have diverse characteristics.Therefore, we can test our methods on a diverse data set. Extensive tests indicate that our method can be used to automatically detect buildings in a robust and fast manner in Ikonos satellite and our aerial images.


甚高分辨率(VHR)航空和卫星图像提供了有价值的信息。特别是,从这些图像中检测建筑物需要特别考虑,因为这些信息可用于多种遥感应用,例如自动地图制作、城市规划和土地利用分析。不幸的是,出于多种原因,人类专家在给定的航空或卫星图像中手动标记建筑物是乏味的。首先,可以从不同的视角对建筑物进行成像。它们可能没有唯一的表示形式。其次,建筑物可能与环境有复杂的相互作用(例如树木的遮挡)。此外,它们可能会相互遮挡。第三,图像中的照明和对比度可能不足以可靠地检测建筑物。第四,这些图像可能覆盖了大片地理区域,其中有许多建筑物。分析图像可能需要一些时间。最后,建筑物没有标准的大小和形状。因此,在VHR航空和卫星图像上开发鲁棒且快速的建筑检测算法已成为必要。


在过去的二十年中,研究人员开发了使用航空和卫星图像的自动建筑检测方法。在航空和卫星图像中对建筑物检测有很好的评论[23],[38]。有兴趣的读者可以从这些研究中获得有关大多数自动化建筑检测方法的详细信息。接下来,我们总结了最近引入的建筑检测方法,重点介绍了与本文所提出的方法相关的方法。Kim和Muller[17]使用图论来检测航空图像中的建筑物。他们提取给定图像中的线性特征,并将它们用作图形的顶点。然后,他们通过将子图与模型建筑图进行匹配来提取建筑物。最后,他们使用强度和阴影信息来验证建筑物的外观。与我们不同的是,他们使用彩色航拍图像和线性特征。Krishnamachari和Chellappa[18]在航空图像中引入了一种基于马尔可夫随机场(MRF)的建筑物检测方法。它们受益于图像中的直线段,并在相互作用上形成基于MRF的检测方法。与我们的系统相比,这个系统更复杂。Segl和Kaufmann [29]以迭代方式将监督形状分类与无监督图像分割相结合。他们的方法允许在高分辨率卫星图像中搜索小物体(如建筑物)。Molinier等人[25]考虑通过训练自组织地图来检测卫星图像中人造结构的边界。Gamba等人[11]使用边界信息提取城市区域的地图。他们将边界和非边界数据提供给两个不同的分类器。然后,他们将结果结合起来,在VHR图像上检测城市地区的建筑物。在这些研究中,总是需要一个训练集。Benediktsson等人[6]使用数学形态学运算来提取结构信息,以检测卫星图像中的城市区域。此方法可用于检测影像中的建筑物。Ünsalan和Boyer[38]研究了多光谱卫星图像,以检测居民区的建筑物和街道网络。他们的方法使用植被指数,聚类,分解二进制图像和图论。虽然这种方法很有前途,但需要多光谱信息来检测建筑物。Akçay和Aksoy [1]还提出了一种在高分辨率卫星图像中进行无监督分割和目标检测的新方法。这种方法还需要多光谱信息。Idrissa等人[13]使用Gabor滤波器和SPOT5图像中的归一化差异植被指数提取人造结构(建筑物和道路)的边缘。通过比较从同一区域拍摄的两个图像序列的边缘,他们还检测到了变化。与我们不同,他们受益于多光谱信息。在最近的一项研究中,我们引入了一种使用尺度不变特征变换(SIFT)关键点和图论形式化来检测全色Ikonos卫星图像中建筑物的方法[32]。此方法给出了良好的结果,但它的计算负载很高。它还依赖于模板构建图像作为训练集。在类似的框架中,Xiong和Zhang[41]使用兴趣点进行卫星图像匹配。还有各种研究侧重于航空和卫星图像中的建筑物形状提取[5],[11],[15],[16],[40]。与建筑物检测相比,这是一个更复杂的问题。但是,检测建筑物位置可能有助于从图像中提取建筑物形状。


📚2 运行结果


2da742e59a4a4fac9a747f499c7fbdcf.png


放大:


93057ac9d1754f43bbddb0d2e3e7e485.png

ead14d27fe104d3988d43abacdd3fc4d.png

49653a89f2154a8a8579a7762176ec9c.png


放大:


efd34fd686c1450bb2be34b8d4ccbff9.png

7385e83d0b87400a8019ce478003d368.png


放大图:


9709876a189e4e1fa9a7771811cfa496.png


部分代码:

%
    % INPUT
    % =====
    % im     : the graylevel image
    %
    % OUTPUT
    % ======
    % points : the interest points extracted
    %
    % REFERENCES
    % ==========
    % C.G. Harris and M.J. Stephens. "A combined corner and edge detector",
    % Proceedings Fourth Alvey Vision Conference, Manchester.
    % pp 147-151, 1988.
    %
    % Alison Noble, "Descriptions of Image Surfaces", PhD thesis, Department
    % of Engineering Science, Oxford University 1989, p45.
    %
    % C. Schmid, R. Mohrand and C. Bauckhage, "Evaluation of Interest Point Detectors",
    % Int. Journal of Computer Vision, 37(2), 151-172, 2000.
    %
    % EXAMPLE
    % =======
    % points = kp_harris(im)
    % only luminance value
    im = double(im(:,:,1));
    sigma = 1.5;
    % derivative masks
    s_D = 0.7*sigma;
    x  = -round(3*s_D):round(3*s_D);
    dx = x .* exp(-x.*x/(2*s_D*s_D)) ./ (s_D*s_D*s_D*sqrt(2*pi));
    dy = dx';
    % image derivatives
    Ix = conv2(im, dx, 'same');
    Iy = conv2(im, dy, 'same');
    % sum of the Auto-correlation matrix
    s_I = sigma;
    g = fspecial('gaussian',max(1,fix(6*s_I+1)), s_I);
    Ix2 = conv2(Ix.^2, g, 'same'); % Smoothed squared image derivatives
    Iy2 = conv2(Iy.^2, g, 'same');
    Ixy = conv2(Ix.*Iy, g, 'same');
    % interest point response
%    cim = (Ix2.*Iy2 - Ixy.^2)./(Ix2 + Iy2 + eps);                % Alison Noble measure.
     %k=0.31;
     %k=0.26;
     %k=0.21;
     %k=0.16;
     %k=0.11; 
     k=0.06;
     %k=0.01;
     %k=0.001; 
     cim = (Ix2.*Iy2 - Ixy.^2) - k*(Ix2 + Iy2).^2;    % Original Harris measure.
    % find local maxima on 3x3 neighborhood
    [r,c,max_local] = findLocalMaximum(cim,3*s_I);
    % set threshold 1% of the maximum value
%    t = 0.01*max(max_local(:));
    t = 0.1*max(max_local(:));
    % find local maxima greater than threshold
    [r,c] = find(max_local>=t);
    % build interest points
    points = [r,c];


🎉3 参考文献

部分理论来源于网络,如有侵权请联系删除。


[1]B. Sirmacek and C. Unsalan, "A Probabilistic Framework to Detect Buildings in Aerial and Satellite Images," in IEEE Transactions on Geoscience and Remote Sensing, vol. 49, no. 1, pp. 211-221, Jan. 2011, doi: 10.1109/TGRS.2010.2053713.


🌈4 Matlab代码实现


相关文章
|
2月前
|
算法 数据安全/隐私保护
织物图像的配准和拼接算法的MATLAB仿真,对比SIFT,SURF以及KAZE
本项目展示了织物瑕疵检测中的图像拼接技术,使用SIFT、SURF和KAZE三种算法。通过MATLAB2022a实现图像匹配、配准和拼接,最终检测并分类织物瑕疵。SIFT算法在不同尺度和旋转下保持不变性;SURF算法提高速度并保持鲁棒性;KAZE算法使用非线性扩散滤波器构建尺度空间,提供更先进的特征描述。展示视频无水印,代码含注释及操作步骤。
|
3月前
|
算法 数据可视化 数据安全/隐私保护
基于LK光流提取算法的图像序列晃动程度计算matlab仿真
该算法基于Lucas-Kanade光流方法,用于计算图像序列的晃动程度。通过计算相邻帧间的光流场并定义晃动程度指标(如RMS),可量化图像晃动。此版本适用于Matlab 2022a,提供详细中文注释与操作视频。完整代码无水印。
|
5月前
|
算法
基于粒子群优化的图像融合算法matlab仿真
这是一个基于粒子群优化(PSO)的图像融合算法,旨在将彩色模糊图像与清晰灰度图像融合成彩色清晰图像。在MATLAB2022a中测试,算法通过PSO求解最优融合权值参数,经过多次迭代更新粒子速度和位置,以优化融合效果。核心代码展示了PSO的迭代过程及融合策略。最终,使用加权平均法融合图像,其中权重由PSO计算得出。该算法体现了PSO在图像融合领域的高效性和融合质量。
|
6月前
|
机器学习/深度学习 算法 固态存储
m基于深度学习的卫星遥感图像轮船检测系统matlab仿真,带GUI操作界面
在MATLAB 2022a中,使用GoogLeNet对卫星遥感图像进行轮船检测,展示了高效的目标识别。GoogLeNet的Inception架构结合全局平均池化增强模型泛化性。核心代码将图像切块并分类,预测为轮船的部分被突出显示,体现了深度学习在复杂场景检测中的应用。
409 8
|
6月前
|
算法 计算机视觉 异构计算
基于FPGA的图像一维FFT变换IFFT逆变换verilog实现,包含tb测试文件和MATLAB辅助验证
```markdown ## FPGA 仿真与 MATLAB 显示 - 图像处理的 FFT/IFFT FPGA 实现在 Vivado 2019.2 中仿真,结果通过 MATLAB 2022a 展示 - 核心代码片段:`Ddddddddddddddd` - 理论:FPGA 实现的一维 FFT/IFFT,加速数字信号处理,适用于高计算需求的图像应用,如压缩、滤波和识别 ```
|
6月前
|
算法 计算机视觉
基于Chan-Vese算法的图像边缘提取matlab仿真
**算法预览展示了4幅图像,从边缘检测到最终分割,体现了在matlab2022a中应用的Chan-Vese水平集迭代过程。核心代码段用于更新水平集并显示迭代效果,最后生成分割结果及误差曲线。Chan-Vese模型(2001)是图像分割的经典方法,通过最小化能量函数自动检测平滑区域和清晰边界的图像分割,适用于复杂环境,广泛应用于医学影像和机器视觉。**
|
7月前
|
算法 数据安全/隐私保护 C++
基于二维CS-SCHT变换和扩频方法的彩色图像水印嵌入和提取算法matlab仿真
该内容是关于一个图像水印算法的描述。在MATLAB2022a中运行,算法包括水印的嵌入和提取。首先,RGB图像转换为YUV格式,然后水印通过特定规则嵌入到Y分量中,并经过Arnold置乱增强安全性。水印提取时,经过逆过程恢复,使用了二维CS-SCHT变换和噪声对比度(NC)计算来评估水印的鲁棒性。代码中展示了从RGB到YUV的转换、水印嵌入、JPEG压缩攻击模拟以及水印提取的步骤。
|
6月前
|
算法 计算机视觉 异构计算
基于FPGA的图像直方图均衡化处理verilog实现,包含tb测试文件和MATLAB辅助验证
摘要: 在FPGA上实现了图像直方图均衡化算法,通过MATLAB2022a与Vivado2019.2进行仿真和验证。核心程序涉及灰度直方图计算、累积分布及映射变换。算法旨在提升图像全局对比度,尤其适合低对比度图像。FPGA利用可编程增益器和查表技术加速硬件处理,实现像素灰度的均匀重分布,提升视觉效果。![image preview](https://ucc.alicdn.com/pic/developer-ecology/3tnl7rfrqv6tw_a075525027db4afbb9c0529921fd0152.png)
|
7月前
|
算法 异构计算
基于直方图的图像曝光量分析FPGA实现,包含tb测试文件和MATLAB辅助验证
该内容包括了算法的运行效果展示、软件版本信息、理论概述和核心程序代码。在正常图像中,`checkb`位于`f192b`和`f250b`之间,而多度曝光图像中`checkb`超出此范围,判断为曝光过度。使用的软件为Vivado 2019.2和MATLAB 2022a。理论依据未详细给出,但提及主要方法。提供的Verilog代码段用于处理图像数据,包括读取文件、时钟控制及图像histogram计算等,其中模块`im_hist`似乎是关键部分。
|
7月前
|
存储 算法 数据可视化
基于harris角点和RANSAC算法的图像拼接matlab仿真
本文介绍了使用MATLAB2022a进行图像拼接的流程,涉及Harris角点检测和RANSAC算法。Harris角点检测寻找图像中局部曲率变化显著的点,RANSAC则用于排除噪声和异常点,找到最佳匹配。核心程序包括自定义的Harris角点计算函数,RANSAC参数设置,以及匹配点的可视化和仿射变换矩阵计算,最终生成全景图像。

相关课程

更多