CN118038280A - 基于航拍影像的建筑施工进度监测及预警方法 - Google Patents
基于航拍影像的建筑施工进度监测及预警方法 Download PDFInfo
- Publication number
- CN118038280A CN118038280A CN202410444932.9A CN202410444932A CN118038280A CN 118038280 A CN118038280 A CN 118038280A CN 202410444932 A CN202410444932 A CN 202410444932A CN 118038280 A CN118038280 A CN 118038280A
- Authority
- CN
- China
- Prior art keywords
- degree
- window
- pixel point
- dust
- building construction
- 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
- 238000009435 building construction Methods 0.000 title claims abstract description 117
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000012544 monitoring process Methods 0.000 title claims abstract description 39
- 239000000428 dust Substances 0.000 claims abstract description 177
- 230000002146 bilateral effect Effects 0.000 claims abstract description 23
- 238000001914 filtration Methods 0.000 claims abstract description 23
- 238000010276 construction Methods 0.000 claims abstract description 18
- 238000010606 normalization Methods 0.000 claims description 14
- 238000004891 communication Methods 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000002708 enhancing effect Effects 0.000 description 3
- 239000000779 smoke Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/176—Urban or other man-made structures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/08—Construction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/36—Applying a local operator, i.e. means to operate on image points situated in the vicinity of a given point; Non-linear local filtering operations, e.g. median filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/17—Terrestrial scenes taken from planes or by drones
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Primary Health Care (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Economics (AREA)
- Nonlinear Science (AREA)
- Health & Medical Sciences (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
本发明涉及建筑施工航拍图像增强技术领域,具体涉及一种基于航拍影像的建筑施工进度监测及预警方法。该方法根据建筑施工航拍图像中的灰度分布情况,获得每个预设窗口的灰尘影响程度,进而对图像进行区域划分,获得所有灰尘区域;根据灰尘区域内像素点点的梯度分布获得灰尘区域中每个像素点的显著程度,利用灰尘区域中的整体灰尘影响程度与每个像素点的显著程度获得灰尘区域中像素点的调整像素范围域权重,并进行双边滤波,获得图像增强结果;利用图像增强结果对建筑施工进度进行监测预警。本发明通过调整像素范围域权重,提高了双边滤波的鲁棒性,能够获得理想的建筑施工图像的图像增强结果,使建筑施工进度监测预警的准确性。
Description
技术领域
本发明涉及建筑施工航拍图像增强技术领域,具体涉及一种基于航拍影像的建筑施工进度监测及预警方法。
背景技术
在进行建筑施工时,需要根据建筑的施工进度进行监测,避免建筑施工时由于施工问题导致其逾期,影响后续的各项工程和收益情况。在建筑施工过程中通过航拍获取到施工过程中的建筑影像信息,分析采集的影像信息中的施工进度和预计设定的施工进度,分析两者之间进度的差异性,对现在施工的情况进行判断对其进行预警,避免施工进度导致逾期问题。
在现有技术中,常使用双边滤波算法对采集到的图像进行图像增强以保留更多的图像信息。对采集图像进行增强的过程中,由于建筑会产生烟尘,进而导致获取的建筑施工航拍图像边缘过于平滑,导致双边滤波中的空间域权重过大,而像素范围域权重过小,进而导致增强后的建筑施工航拍图像会损失大量边缘信息,无法获得理想的图像增强结果,进而影响建筑施工进度监测的准确性。
发明内容
为了解决由于建筑产生烟尘导致获取的图像边缘过于平滑,导致双边滤波中的空间域权重过大,而像素范围域权重过小,进而导致增强后的图像会损失大量边缘信息,无法获得理想的图像增强结果,进而影响建筑施工进度监测的准确性的技术问题,本发明的目的在于提供一种基于航拍影像的建筑施工进度监测及预警方法,所采用的技术方案具体如下:
一种基于航拍影像的建筑施工进度监测及预警方法,所述方法包括:
获取建筑施工航拍图像;
根据所述建筑施工航拍图像中每个预设窗口及其预设第一邻域范围内的邻域窗口内的灰度分布情况,获得所述建筑施工航拍图像中每个窗口的灰尘影响程度;依据所有窗口之间的距离与每个窗口的所述灰尘影响程度对所述建筑施工航拍图像中的窗口进行区域划分,获得所有灰尘区域;根据每个所述灰尘区域内所有窗口的所述灰尘影响程度获取每个所述灰尘区域的整体灰尘影响程度;
根据每个所述灰尘区域内像素点的梯度分布特征,获得每个所述灰尘区域内每个像素点的显著程度;根据每个所述灰尘区域的整体灰尘影响程度与每个所述灰尘区域内每个像素点的所述显著程度,获得每个灰尘区域内每个像素点的调整像素范围域权重;
利用所述调整像素范围域权重与空间域权重对所述建筑施工航拍图像进行双边滤波,获得图像增强结果;
利用所述图像增强结果对建筑施工进度进行监测预警。
进一步地,所述灰尘影响程度的获取方法包括:
将每个窗口内像素点的灰度方差的倒数作为每个窗口内的灰度平缓程度;
根据每个窗口内像素点的灰度位置分布特征获得每个窗口内的灰度混乱程度;
根据每个窗口与所有邻域窗口内的所述灰度混乱程度,获得每个窗口与所有邻域窗口之间的灰度混乱差异程度;
根据所述灰度平缓程度、所述灰度混乱程度与所述灰度混乱差异程度获得每个窗口的灰尘影响程度。
进一步地,所述灰度混乱程度的获取方法包括:
对每个窗口进行阈值分割,获得每个窗口受到灰尘影响的连通域;计算获取连通域的周长与面积;
将每个窗口内的连通域周长与连通域面积之间的比值作为每个窗口内的灰度混乱程度。
进一步地,所述灰度混乱差异程度获取方法包括:
将每个窗口的连通域数量与所述灰度混乱程度之间的乘积作为窗口混乱程度,将每个窗口及其对应的所有邻域窗口的所述窗口混乱程度进行累加求和获得第一混乱程度;
计算每个窗口与所有邻域窗口组成的整体中连通域的数量与所述灰度混乱程度之间的乘积作为第二混乱程度;
将所述第一混乱程度与所述第二混乱程度之间的比值进行归一化,作为每个窗口与所有邻域窗口之间的灰度混乱差异程度。
进一步地,根据所述灰度平缓程度、所述灰度混乱程度与所述灰度混乱差异程度获得每个窗口的灰尘影响程度,包括:
将所述灰度平缓程度、所述灰度混乱程度与所述灰度混乱差异程度的乘积归一化处理,获得每个窗口的灰尘影响程度。
进一步地,所述显著程度的获取方法包括:
获取每个像素点预设第二邻域范围内的所有邻域像素点;
根据每个像素点所有方向上的梯度幅值变化程度获得每个像素点的梯度变化程度;
根据每个像素点的所有邻域像素点的所述梯度变化程度,及每个像素点与所有邻域像素点的梯度特征,获得每个像素点与所有邻域像素点之间的梯度关联程度;
根据每个像素点的所述梯度变化程度及每个像素点与所有邻域像素点之间的梯度关联程度,获得每个像素点的显著程度。
进一步地,所述梯度变化程度的获取方法包括:
将每个像素点在所有方向中每两个方向之间的梯度幅值差异最大值进行归一化处理,获得每个像素点的梯度变化程度。
进一步地,所述梯度关联程度的获取方法包括:
将每个像素点的梯度方向与梯度方向上的梯度幅值之间的乘积作为每个像素点的梯度特征;
将每个像素点的梯度特征与每个像素点的每个邻域像素点的梯度特征之间的比值进行归一化处理,并与每个邻域像素点的所述梯度变化程度相乘获得第一乘积,将每个邻域像素点的所述第一乘积累加求和,获得每个像素点与所有邻域像素点之间的梯度关联程度。
进一步地,根据每个像素点的所述梯度变化程度及每个像素点与所有邻域像素点之间的梯度关联程度,获得每个像素点的显著程度,包括:
将每个像素点的所述梯度变化程度,和每个像素点与所有邻域像素点之间的梯度关联程度的乘积进行归一化处理,作为每个像素点的显著程度。
进一步地,所述调整像素范围域权重的获取方法包括:
将每个像素点的所述显著程度与每个像素点所在灰尘区域的整体灰尘影响程度之间的和进行归一化处理,获得每个像素点的所述调整像素范围域权重。
本发明具有如下有益效果:
本发明获取建筑施工航拍图像,为了后续对建筑施工进度进行监测;根据建筑施工航拍图像中每个窗口与每个窗口的所有邻域窗口内的灰度分布情况,获得建筑施工航拍图像中每个窗口的灰尘影响程度,反映了不同窗口内细小边缘信息需要保留的程度;依据所有窗口之间的距离与每个窗口的灰尘影响程度对建筑施工航拍图像中的窗口进行区域划分,获得所有灰尘区域,灰尘区域内包含这灰尘影响程度相似的窗口,这些窗口内对细小边缘的保留程度相同;根据每个灰尘区域内所有窗口的灰尘影响程度获取每个灰尘区域的整体灰尘影响程度;根据每个灰尘区域内像素点的梯度分布特征,获得每个灰尘区域内每个像素点的显著程度,显著程度可以区分边缘信息是否为建筑施工航拍图像背景区域的有效边缘信息;由于相同显著程度的像素点由于所在灰尘区域的整体灰尘影响程度不同,对像素范围域权重的调整情况也不相同,根据每个灰尘区域的整体灰尘影响程度与每个灰尘区域内每个像素点的显著程度,获得每个灰尘区域内每个像素点的调整像素范围域权重;利用调整像素范围域权重与空间域权重对建筑施工航拍图像进行双边滤波,获得图像增强结果;利用图像增强结果对建筑施工进度进行监测预警。本发明通过调整像素范围域权重,提高了双边滤波的鲁棒性,能够获得理想的建筑施工图像的图像增强结果,使建筑施工进度监测预警的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案和优点,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它附图。
图1为本发明一个实施例所提供的一种基于航拍影像的建筑施工进度监测及预警方法流程图。
具体实施方式
为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明提出的一种基于航拍影像的建筑施工进度监测及预警方法,其具体实施方式、结构、特征及其功效,详细说明如下。在下述说明中,不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外,一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。
一种基于航拍影像的建筑施工进度监测及预警方法实施例:
下面结合附图具体的说明本发明所提供的一种基于航拍影像的建筑施工进度监测及预警方法的具体方案。
请参阅图1,其示出了本发明一个实施例提供的一种基于航拍影像的建筑施工进度监测及预警方法流程图,该方法包括:
步骤S1:获取建筑施工航拍图像。
本发明实施例提供了一种基于航拍影像的建筑施工进度监测及预警方法,针对于建筑施工进度进行监测,首先需要获取建筑施工场所的图像,再进行后续的图像分析。为了使图像能够更加完整的反映建筑施工场所的实际情况,本发明一个实施例中,利用无人机携带的可见光摄像头拍摄建筑施工场所的图像,并对采集到的图像进行预处理操作,将采集到的图像通过无线设备传输到数据处理终端中,获得建筑施工航拍图像。
本发明实施例中,对采集图像进行灰度化处理的预处理操作。
需要说明的是,图像采集过程可以由实施人员依据具体实施场景自行设置,在此不做限定。
步骤S2:根据建筑施工航拍图像中每个预设窗口及其预设第一邻域范围内的邻域窗口内的灰度分布情况,获得建筑施工航拍图像中每个窗口的灰尘影响程度;依据所有窗口之间的距离与每个窗口的灰尘影响程度对建筑施工航拍图像中的窗口进行区域划分,获得所有灰尘区域;根据每个灰尘区域内所有窗口的灰尘影响程度获取每个灰尘区域的整体灰尘影响程度。
由于在进行建筑施工的过程中会产生大量灰尘,会导致建筑施工航拍图像中很多边缘被灰尘影响,使得边缘信息不清晰。因为灰尘出现位置具有随机性,在建筑航拍图像中不同区域受到灰尘的影响可能不同,即灰度分布可能在不同的局部区域有较为明显的变化,为了研究不同区域受到灰尘的影响,需要将建筑施工航拍图像分割成许多窗口以保留图像中受到灰尘影响后的更多局部灰度特征。不同窗口受到灰尘影响的程度不同,此时不同窗口内细小边缘信息需要保留的程度不同,为了对灰尘影响的图像窗口内的边缘信息进行保留,需要获取窗口受到的灰尘影响程度,由于不同灰尘影响程度的窗口内灰度分布情况有很大不同,所以需要首先分析窗口内像素点灰度分布情况。故本发明实施例中,根据建筑施工航拍图像中每个窗口及其预设第一邻域范围内的邻域窗口内的灰度分布情况,获得建筑施工航拍图像中每个窗口的灰尘影响程度。
优选地,本发明一个实施例中,灰尘影响程度的获取方法包括:
将每个窗口内像素点的灰度方差的倒数作为每个窗口内的灰度平缓程度,由于受灰尘影响较大的窗口内灰度变化较为平缓,灰度方差越大,此时窗口内灰度变化越剧烈,说明窗口内灰尘影响程度越低,此时窗口内的灰度平缓程度越小;根据每个窗口内像素点的灰度位置分布特征获得每个窗口内的灰度混乱程度,若窗口内相似灰度值的分布越混乱,越有可能受到灰尘影响,此时灰度混乱程度越大;根据每个窗口与所有邻域窗口内的灰度混乱程度,获得每个窗口与所有邻域窗口之间的灰度混乱差异程度,灰度混乱差异程度越大,每个窗口与所有邻域窗口之间的灰度关联性越差,越有可能受到灰尘影响;根据灰度平缓程度、灰度混乱程度与灰度混乱差异程度获得每个窗口的灰尘影响程度。
本发明一个实施例中,预设窗口设置为边长为的正方形。需要说明的是,在本发明其他实施例中,窗口的参数可以由实施人员依据具体实施场景自行设置,在此不做限定。
本发明一个实施例中,预设第一邻域范围设置为以预设窗口中心点为圆心,半径为6的圆形,则邻域窗口为预设窗口在第一邻域范围内除自身外的所有窗口,且邻域窗口的参数与预设窗口相同。需要说明的是,在本发明其他实施例中,第一邻域范围与邻域窗口的设置可以由实施人员依据具体实施场景自行设置,在此不做限定。
优选地,本发明一个实施例中,灰度混乱程度的获取方法包括:
对每个窗口进行阈值分割,由于灰尘影响区域灰度值变化较为平缓,而建筑施工的背景区域灰度值会发生较大改变,阈值分割可以获得灰度值相似的区域,即灰尘影响的区域,获得每个窗口受到灰尘影响的连通域;计算获取连通域的周长与面积,连通域的周长和面积反映出窗口内受到灰尘影响的区域特征;将每个窗口内的连通域周长与连通域面积之间的比值作为每个窗口内的灰度混乱程度,比值越大,说明连通域组成的形状越复杂,灰度混乱程度越大。
优选地,本发明一个实施例中,灰度混乱差异程度获取方法包括:
将每个窗口的连通域数量与所述灰度混乱程度之间的乘积作为窗口混乱程度,将每个窗口及其对应的所有邻域窗口的所述窗口混乱程度进行累加求和获得第一混乱程度;计算每个窗口与所有邻域窗口组成的整体中连通域的数量与所述灰度混乱程度之间的乘积作为第二混乱程度,若每个窗口与所有邻域窗口组成的整体中连通域数量与分开时的连通域数量越接近,说明窗口与邻域窗口之间关联性越差,此时整体的连通域的复杂情况越高,此时受到的灰尘影响程度越大;将所述第一混乱程度与所述第二混乱程度之间的比值进行归一化,作为每个窗口与所有邻域窗口之间的灰度混乱差异程度,灰度混乱差异程度越大,说明窗口与邻域窗口之间的灰度分布越混乱,灰度分布关联性越差,此时窗口受到的灰尘影响程度越大。
优选地,本发明一个实施例中,将灰度平缓程度、灰度混乱程度与灰度混乱差异程度的乘积归一化处理,获得每个窗口的灰尘影响程度。本发明一个实施例中,灰尘影响程度计算公式如下所示:
式中,表示第/>个窗口内的灰尘影响程度;/>表示第/>个窗口内像素点的灰度方差;/>表示第/>个窗口内受到灰尘影响的连通域的周长;/>表示第/>个窗口内受到灰尘影响的连通域的面积;/>表示第/>个窗口与所有邻域窗口的数量;/>表示第/>个窗口与所有邻域窗口的序号;/>表示第/>个窗口与所有邻域窗口内第/>个窗口中连通域的数量;/>表示第/>个窗口与所有邻域窗口内第/>个窗口中连通域的周长;/>表示第/>个窗口与所有邻域窗口内第/>个窗口中连通域的面积;/>表示第/>个窗口与所有邻域窗口组成的整体内连通域的数量;表示第/>个窗口与所有邻域窗口组成的整体内连通域的周长;/>表示第/>个窗口与所有邻域窗口组成的整体内连通域的面积;/>表示第/>个窗口的灰度平缓程度;/>表示第/>个窗口的灰度混乱程度;/>表示第/>个窗口与所有邻域窗口之间的灰度混乱差异程度;/>表示以自然常数为底数的指数函数;/>表示归一化函数。
在灰尘影响程度计算公式中,越大,第/>个窗口的灰度平缓程度越大,说明此时第个窗口内建筑施工场所的背景区域越不明显,第/>个窗口受到的灰尘影响程度越大;由于灰尘出现位置是随机的,相较于建筑施工航拍图像的背景区域而言,灰尘影响下窗口的灰度分布要更加混乱,/>越大,第/>个窗口内的连通域周长相对于连通域面积越大,此时连通域的边缘越复杂,说明此时第/>个窗口的灰度混乱程度越大,第/>个窗口受到的灰尘影响程度越大;在实际情况中,第/>个窗口与所有邻域窗口组成的整体内的连通域数量、连通域周长和连通域面积的最大值等于第/>个窗口与所有邻域窗口每个窗口内连通域数量、连通域周长和连通域面积,/>越大,第/>个窗口与所有邻域窗口之间的灰度混乱差异程度越大,此时第/>个窗口与所有邻域窗口的第一混乱程度与第/>个窗口与所有邻域窗口整体的第二混乱程度越接近,说明此时第/>个窗口与所有邻域窗口分开时的连通域数量和连通域的区域分布,和第/>个窗口与所有邻域窗口组成的整体内的连通域数量和连通域的区域分布越接近,说明第/>个窗口与所有邻域窗口之间的灰度分布关联性越差,第个窗口受到的灰尘影响程度越大。
由于在受到不同灰尘影响程度的窗口内,对细小边缘信息的保留程度不同,所以需要对受到不同灰尘影响程度的窗口进行分类,所以本发明实施例中,依据灰尘影响程度与窗口之间的距离对建筑施工航拍图像进行区域划分,获得所有灰尘区域,并获得每个灰尘区域的整体灰尘影响程度。
本发明一个实施例中,利用窗口的灰尘影响程度与窗口间的距离对每个窗口进行聚类,聚类结果可以将灰尘影响程度相似且距离较为接近的窗口分为一类,根据聚类结果将建筑施工航拍图像进行区域划分,获得建筑施工航拍图像中的所有灰尘区域,将每个灰尘区域中所有窗口的灰尘影响程度求平均,获得每个灰尘区域的整体灰尘影响程度。
需要说明的是,聚类算法可以采用K-Means聚类算法、DBSCAN聚类算法等多种聚类算法,在此不做限定,由于各种聚类算法均为本领域技术人员所熟知的技术手段,在此不做赘述。
步骤S3:根据每个灰尘区域内像素点的梯度分布特征,获得每个灰尘区域内每个像素点的显著程度;根据每个灰尘区域的整体灰尘影响程度与每个灰尘区域内每个像素点的显著程度,获得每个灰尘区域内每个像素点的调整像素范围域权重。
由于后续利用双边滤波对建筑施工航拍图像进行图像增强时要保留边缘信息,然而在建筑施工航拍图像中的某些区域内细小边缘信息是由噪声或其他因素引起的,实际情况中并不需要由噪声或其他因素引起的细小边缘信息,由于在建筑施工航拍图像中建筑施工场所的背景包含的有效边缘信息,其梯度分布特征与噪声或其他因素引起的无效边缘信息的梯度分布特征不同,即有效边缘信息与无效边缘信息包含的像素点的显著程度不同。所以本发明实施例中,根据每个灰尘区域内像素点的梯度分布程度,获得每个灰尘区域内每个像素点的显著程度。
优选地,本发明一个实施例中,显著程度的获取方法包括:
获取每个像素点预设第二邻域范围内的所有邻域像素点,通过每个像素点与第二邻域范围内邻域像素点的梯度特征可以反映出每个边缘像素点是否属于有效的边缘信息;根据每个像素点所有方向上的梯度幅值获得每个像素点的梯度变化程度,每个方向的梯度幅值可以反映出每个像素点在所有方向的纹理特征;每个像素点的梯度方向与梯度幅值反映出每个像素点的梯度特征,根据每个像素点的所有邻域像素点的梯度变化程度,及每个像素点与所有邻域像素点的梯度特征,获得每个像素点与所有邻域像素点之间的梯度关联程度,梯度关联程度越大,说明每个像素点与邻域像素点的梯度特征越相似;根据每个像素点的梯度变化程度及每个像素点与所有邻域像素点之间的梯度关联程度,获得每个像素点的显著程度。
本发明一个实施例中,预设第二邻域范围设置为以像素点为中心,半径为5的圆形。需要说明的是,在本发明其他实施例中,预设第二邻域范围的参数可以由实施人员依据具体实施场景自行设置,在此不做限定。
优选地,本发明一个实施例中,梯度变化程度的获取方法包括:
将每个像素点在所有方向中每两个方向之间的梯度幅值差异最大值进行归一化处理,获得每个像素点的梯度变化程度,像素点的梯度幅值差异越大,说明该像素点越不可能是灰尘影响造成的梯度变化,此时的梯度变化程度越大,像素点的显著程度越高。
优选地,本发明一个实施例中,梯度关联程度的获取方法包括:
将每个像素点的梯度方向与梯度方向上的梯度幅值之间的乘积作为每个像素点的梯度特征;将每个像素点的梯度特征与每个像素点的每个邻域像素点的梯度特征之间的比值进行归一化处理,并与每个邻域像素点的梯度变化程度相乘获得第一乘积,将每个邻域像素点的第一乘积累加求和,获得每个像素点与所有邻域像素点之间的梯度关联程度。
本发明一个实施例中,对建筑施工航拍图像建立笛卡尔坐标系,梯度方向则为建筑施工航拍图像中每个像素点的梯度方向向量与坐标横轴之间夹角的余弦值。需要说明的是,在本发明其他实施例中,梯度方向可以由其他数学运算形式进行表达,在此不做限定。
优选地,本发明一个实施例中,将每个像素点的梯度变化程度,和每个像素点与所有邻域像素点之间的梯度关联程度的乘积进行归一化处理,作为每个像素点的显著程度。本发明一个实施例中,显著程度计算公式如下所示:
式中,表示第/>个像素点的显著程度;/>表示第/>个像素点的梯度变化程度;/>表示第/>个像素点的邻域像素点的数量;/>表示第/>个像素点的邻域像素点的序号;/>表示第/>个像素点的第/>个邻域像素点的梯度变化程度;/>表示第/>个像素点的梯度方向;/>表示第/>个像素点的梯度幅值;/>表示第/>个像素点的第/>个邻域像素点的梯度方向;表示第/>个像素点的第/>个邻域像素点的梯度幅值;表示第/>个像素点与所有邻域像素点之间的梯度关联程度;/>表示归一化函数;/>表示以自然常数为底数的指数函数。
在显著程度计算公式中,第个像素点的梯度变化程度越大,说明第/>个像素点在某两个方向上梯度幅值差异越明显,说明第/>个像素点的变化越剧烈,越不可能是受到灰尘影响,此时第/>个像素点的显著程度越高;/>表示第/>个像素点的第/>个邻域像素点的梯度变化程度,同理,梯度变化程度越大,第/>个像素点的第/>个邻域像素点的显著程度越高;表示第/>个像素点的梯度特征,/>表示第/>个像素点的第/>个邻域像素点的梯度特征,当第/>个像素点的梯度特征与第/>个邻域像素点的梯度特征越接近时,第/>个像素点与第/>个邻域像素点越有可能属于同一条纹理,此时/>越大,其中如果第/>个邻域像素点的梯度变化程度越大,说明第/>个邻域像素点变化越显著,此时与第/>个邻域像素点在同一条纹理的第/>个像素点变化也会越显著,即第/>个像素点与所有邻域像素点之间的梯度关联程度越大,说明第/>个像素点的延展程度越大,越有可能属于建筑施工航拍图像背景区域的边缘信息,此时第/>个像素点的显著程度越高。
在本发明实施例中,使用双边滤波算法对建筑施工航拍图像进行图像增强。由于灰尘无法影响建筑施工航拍图像中空间上的变化,所以在双边滤波中,无需改变空间域权重,但在灰尘影响下,建筑施工航拍图像中的灰度信息发生很大改变,所以需要重新调整像素范围域权重,由于相同显著程度的像素点由于所在灰尘区域的整体灰尘影响程度不同,对像素范围域权重的调整情况也不相同,所以本发明实施例中,根据每个灰尘区域的整体灰尘影响程度与每个灰尘区域内每个像素点的显著程度,获得每个灰尘区域内每个像素点的调整像素范围域权重。
优选地,本发明一个实施例中,调整像素范围域权重的获取方法包括:
将每个像素点的显著程度与每个像素点所在灰尘区域的整体灰尘影响程度之间的和进行归一化处理,获得每个像素点的调整像素范围域权重。本发明一个实施例中,调整像素范围域权重计算公式如下所示:
式中,表示第/>个像素点的调整像素范围域权重;/>表示第/>个像素点的显著程度;/>表示第/>个像素点所在灰尘区域的整体灰尘影响程度;/>表示归一化函数。
至此,获得双边滤波中的调整像素范围域权重。
步骤S4:利用调整像素范围域权重与空间域权重对建筑施工航拍图像进行双边滤波,获得图像增强结果。
利用调整像素范围域权重与空间域权重对建筑施工航拍图像进行双边滤波,通过调整像素范围域权重与空间域权重达到对建筑施工航拍图像的平坦区域进行平滑,对受到灰尘影响的边缘信息进行保留的目的,最终获得建筑施工航拍图像的图像增强结果。
需要说明的是,双边滤波算法为本领域技术人员所熟知的技术手段,在此不做赘述。
步骤S5:利用图像增强结果对建筑施工进度进行监测预警。
本发明一个实施例中,通过建筑施工航拍图像的图像增强结果进行分析,获得多个视角下的建筑施工航拍图像的图像增强结果后,利用三维重构算法进行建筑施工模型的重构,将重构后的模型与建筑施工过程中的BIM建筑施工模型进行比较,计算在相同时刻下建筑施工进度的差异性,根据差异性对建筑施工进度进行监测预警。
至此,完成对建筑施工进度进行监测预警。
综上所述,本发明获取建筑施工航拍图像;根据建筑施工航拍图像中每个窗口与每个窗口的所有邻域窗口内的灰度分布情况,获得建筑施工航拍图像中每个窗口的灰尘影响程度;依据窗口之间的距离与灰尘影响程度对建筑施工航拍图像中的窗口进行区域划分,获得所有灰尘区域;根据每个灰尘区域内像素点的梯度分布程度,获得每个灰尘区域内每个像素点的显著程度;根据每个灰尘区域的整体灰尘影响程度与每个灰尘区域内每个像素点的显著程度,获得每个灰尘区域内每个像素点的调整像素范围域权重;利用调整像素范围域权重与空间域权重对建筑施工航拍图像进行双边滤波,获得图像增强结果;利用图像增强结果对建筑施工进度进行监测。本发明通过调整像素范围域权重,提高了双边滤波的鲁棒性,能够获得理想的建筑施工航拍图像的图像增强结果,使建筑施工进度监测预警的准确性。
一种建筑施工航拍图像的图像增强方法实施例:
在现有技术中,利用双边滤波算法对采集到的建筑施工航拍图像进行增强的过程中,会出现由于建筑会产生烟尘,进而导致获取的建筑施工航拍图像边缘过于平滑,导致双边滤波中的空间域权重过大,而像素范围域权重过小,进而导致增强后的建筑施工航拍图像会损失大量边缘信息,无法获得理想的图像增强结果的技术问题。为了解决该技术问题,提出一种建筑施工航拍图像的图像增强方法实施例,该实施例包括:
步骤S1:获取建筑施工航拍图像。
步骤S2:根据建筑施工航拍图像中每个预设窗口及其预设第一邻域范围内的邻域窗口内的灰度分布情况,获得建筑施工航拍图像中每个窗口的灰尘影响程度;依据所有窗口之间的距离与每个窗口的灰尘影响程度对建筑施工航拍图像中的窗口进行区域划分,获得所有灰尘区域;根据每个灰尘区域内所有窗口的灰尘影响程度获取每个灰尘区域的整体灰尘影响程度。
步骤S3:根据每个灰尘区域内像素点的梯度分布特征,获得每个灰尘区域内每个像素点的显著程度;根据每个灰尘区域的整体灰尘影响程度与每个灰尘区域内每个像素点的显著程度,获得每个灰尘区域内每个像素点的调整像素范围域权重。
步骤S4:利用调整像素范围域权重与空间域权重对建筑施工航拍图像进行双边滤波,获得图像增强结果。
由于步骤S1-步骤S4在上述一种基于航拍影像的建筑施工进度监测及预警方法中已经提及,在此不再赘述。
本实施例的技术效果:本实施例获取建筑施工航拍图像,为了后续对建筑施工进度进行监测;根据建筑施工航拍图像中每个窗口与每个窗口的所有邻域窗口内的灰度分布情况,获得建筑施工航拍图像中每个窗口的灰尘影响程度,反映了不同窗口内细小边缘信息需要保留的程度;依据所有窗口之间的距离与每个窗口的灰尘影响程度对建筑施工航拍图像中的窗口进行区域划分,获得所有灰尘区域,灰尘区域内包含这灰尘影响程度相似的窗口,这些窗口内对细小边缘的保留程度相同;根据每个灰尘区域内所有窗口的灰尘影响程度获取每个灰尘区域的整体灰尘影响程度;根据每个灰尘区域内像素点的梯度分布特征,获得每个灰尘区域内每个像素点的显著程度,显著程度可以区分边缘信息是否为建筑施工航拍图像背景区域的有效边缘信息;由于相同显著程度的像素点由于所在灰尘区域的整体灰尘影响程度不同,对像素范围域权重的调整情况也不相同,根据每个灰尘区域的整体灰尘影响程度与每个灰尘区域内每个像素点的显著程度,获得每个灰尘区域内每个像素点的调整像素范围域权重;利用调整像素范围域权重与空间域权重对建筑施工航拍图像进行双边滤波,获得图像增强结果;利用图像增强结果对建筑施工进度进行监测预警。本发明通过调整像素范围域权重,提高了双边滤波的鲁棒性,能够获得理想的建筑施工图像的图像增强结果。
需要说明的是:上述本发明实施例先后顺序仅仅为了描述,不代表实施例的优劣。在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。
Claims (10)
1.一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,所述方法包括:
获取建筑施工航拍图像;
根据所述建筑施工航拍图像中每个预设窗口及其预设第一邻域范围内的邻域窗口内的灰度分布情况,获得所述建筑施工航拍图像中每个窗口的灰尘影响程度;依据所有窗口之间的距离与每个窗口的所述灰尘影响程度对所述建筑施工航拍图像中的窗口进行区域划分,获得所有灰尘区域;根据每个所述灰尘区域内所有窗口的所述灰尘影响程度获取每个所述灰尘区域的整体灰尘影响程度;
根据每个所述灰尘区域内像素点的梯度分布特征,获得每个所述灰尘区域内每个像素点的显著程度;根据每个所述灰尘区域的整体灰尘影响程度与每个所述灰尘区域内每个像素点的所述显著程度,获得每个灰尘区域内每个像素点的调整像素范围域权重;
利用所述调整像素范围域权重与空间域权重对所述建筑施工航拍图像进行双边滤波,获得图像增强结果;
利用所述图像增强结果对建筑施工进度进行监测预警。
2.根据权利要求1所述的一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,所述灰尘影响程度的获取方法包括:
将每个窗口内像素点的灰度方差的倒数作为每个窗口内的灰度平缓程度;
根据每个窗口内像素点的灰度位置分布特征获得每个窗口内的灰度混乱程度;
根据每个窗口与所有邻域窗口内的所述灰度混乱程度,获得每个窗口与所有邻域窗口之间的灰度混乱差异程度;
根据所述灰度平缓程度、所述灰度混乱程度与所述灰度混乱差异程度获得每个窗口的灰尘影响程度。
3.根据权利要求2所述的一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,所述灰度混乱程度的获取方法包括:
对每个窗口进行阈值分割,获得每个窗口受到灰尘影响的连通域;计算获取连通域的周长与面积;
将每个窗口内的连通域周长与连通域面积之间的比值作为每个窗口内的灰度混乱程度。
4.根据权利要求3所述的一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,所述灰度混乱差异程度获取方法包括:
将每个窗口的连通域数量与所述灰度混乱程度之间的乘积作为窗口混乱程度,将每个窗口及其对应的所有邻域窗口的所述窗口混乱程度进行累加求和获得第一混乱程度;
计算每个窗口与所有邻域窗口组成的整体中连通域的数量与所述灰度混乱程度之间的乘积作为第二混乱程度;
将所述第一混乱程度与所述第二混乱程度之间的比值进行归一化,作为每个窗口与所有邻域窗口之间的灰度混乱差异程度。
5.根据权利要求2所述的一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,根据所述灰度平缓程度、所述灰度混乱程度与所述灰度混乱差异程度获得每个窗口的灰尘影响程度,包括:
将所述灰度平缓程度、所述灰度混乱程度与所述灰度混乱差异程度的乘积归一化处理,获得每个窗口的灰尘影响程度。
6.根据权利要求1所述的一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,所述显著程度的获取方法包括:
获取每个像素点预设第二邻域范围内的所有邻域像素点;
根据每个像素点所有方向上的梯度幅值变化程度获得每个像素点的梯度变化程度;
根据每个像素点的所有邻域像素点的所述梯度变化程度,及每个像素点与所有邻域像素点的梯度特征,获得每个像素点与所有邻域像素点之间的梯度关联程度;
根据每个像素点的所述梯度变化程度及每个像素点与所有邻域像素点之间的梯度关联程度,获得每个像素点的显著程度。
7.根据权利要求6所述的一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,所述梯度变化程度的获取方法包括:
将每个像素点在所有方向中每两个方向之间的梯度幅值差异最大值进行归一化处理,获得每个像素点的梯度变化程度。
8.根据权利要求6所述的一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,所述梯度关联程度的获取方法包括:
将每个像素点的梯度方向与梯度方向上的梯度幅值之间的乘积作为每个像素点的梯度特征;
将每个像素点的梯度特征与每个像素点的每个邻域像素点的梯度特征之间的比值进行归一化处理,并与每个邻域像素点的所述梯度变化程度相乘获得第一乘积,将每个邻域像素点的所述第一乘积累加求和,获得每个像素点与所有邻域像素点之间的梯度关联程度。
9.根据权利要求6所述的一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,根据每个像素点的所述梯度变化程度及每个像素点与所有邻域像素点之间的梯度关联程度,获得每个像素点的显著程度,包括:
将每个像素点的所述梯度变化程度,和每个像素点与所有邻域像素点之间的梯度关联程度的乘积进行归一化处理,作为每个像素点的显著程度。
10.根据权利要求1所述的一种基于航拍影像的建筑施工进度监测及预警方法,其特征在于,所述调整像素范围域权重的获取方法包括:
将每个像素点的所述显著程度与每个像素点所在灰尘区域的整体灰尘影响程度之间的和进行归一化处理,获得每个像素点的所述调整像素范围域权重。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410444932.9A CN118038280B (zh) | 2024-04-15 | 2024-04-15 | 基于航拍影像的建筑施工进度监测及预警方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410444932.9A CN118038280B (zh) | 2024-04-15 | 2024-04-15 | 基于航拍影像的建筑施工进度监测及预警方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN118038280A true CN118038280A (zh) | 2024-05-14 |
CN118038280B CN118038280B (zh) | 2024-06-14 |
Family
ID=90991837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410444932.9A Active CN118038280B (zh) | 2024-04-15 | 2024-04-15 | 基于航拍影像的建筑施工进度监测及预警方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118038280B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106651783A (zh) * | 2016-09-29 | 2017-05-10 | 厦门美图之家科技有限公司 | 图像滤波方法及装置 |
WO2019144581A1 (zh) * | 2018-01-29 | 2019-08-01 | 江苏宇特光电科技股份有限公司 | 一种智能化红外图像场景增强方法 |
WO2020057062A1 (zh) * | 2018-09-19 | 2020-03-26 | 北京图森未来科技有限公司 | 在fpga中实现图像双边滤波的方法和装置、fpga |
CN115272315A (zh) * | 2022-09-27 | 2022-11-01 | 南通致和祥智能装备有限公司 | 基于无人机航拍的市政园林面积测量方法 |
CN115953421A (zh) * | 2022-12-22 | 2023-04-11 | 郑州大学 | 一种用于蜂窝结构规整度检测的Harris蜂窝顶点提取方法 |
US20230305572A1 (en) * | 2021-12-29 | 2023-09-28 | China University Of Mining And Technology | Method for drivable area detection and autonomous obstacle avoidance of unmanned haulage equipment in deep confined spaces |
CN117036353A (zh) * | 2023-10-09 | 2023-11-10 | 东莞市盈圆合创新材料有限公司 | 一种新能源电池用的耐温泡棉涂层检测方法 |
CN117197005A (zh) * | 2023-09-26 | 2023-12-08 | 中国十七冶集团有限公司 | 一种炼钢工程区域图像的图像增强系统及其安全监测系统 |
CN117422712A (zh) * | 2023-12-15 | 2024-01-19 | 青岛合丰新材料有限公司 | 基于图像滤波处理的塑料母粒视觉检测方法及系统 |
CN117830139A (zh) * | 2024-03-01 | 2024-04-05 | 天津世亚模具股份有限公司 | 一种高强度冲压件缺陷检测方法 |
CN117853510A (zh) * | 2023-12-14 | 2024-04-09 | 北京航天控制仪器研究所 | 基于双边滤波和自适应阈值的Canny边缘检测方法 |
CN117876971A (zh) * | 2024-03-12 | 2024-04-12 | 武汉同创万智数字科技有限公司 | 基于机器视觉的建筑施工安全监测及预警方法 |
-
2024
- 2024-04-15 CN CN202410444932.9A patent/CN118038280B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106651783A (zh) * | 2016-09-29 | 2017-05-10 | 厦门美图之家科技有限公司 | 图像滤波方法及装置 |
WO2019144581A1 (zh) * | 2018-01-29 | 2019-08-01 | 江苏宇特光电科技股份有限公司 | 一种智能化红外图像场景增强方法 |
WO2020057062A1 (zh) * | 2018-09-19 | 2020-03-26 | 北京图森未来科技有限公司 | 在fpga中实现图像双边滤波的方法和装置、fpga |
US20230305572A1 (en) * | 2021-12-29 | 2023-09-28 | China University Of Mining And Technology | Method for drivable area detection and autonomous obstacle avoidance of unmanned haulage equipment in deep confined spaces |
CN115272315A (zh) * | 2022-09-27 | 2022-11-01 | 南通致和祥智能装备有限公司 | 基于无人机航拍的市政园林面积测量方法 |
CN115953421A (zh) * | 2022-12-22 | 2023-04-11 | 郑州大学 | 一种用于蜂窝结构规整度检测的Harris蜂窝顶点提取方法 |
CN117197005A (zh) * | 2023-09-26 | 2023-12-08 | 中国十七冶集团有限公司 | 一种炼钢工程区域图像的图像增强系统及其安全监测系统 |
CN117036353A (zh) * | 2023-10-09 | 2023-11-10 | 东莞市盈圆合创新材料有限公司 | 一种新能源电池用的耐温泡棉涂层检测方法 |
CN117853510A (zh) * | 2023-12-14 | 2024-04-09 | 北京航天控制仪器研究所 | 基于双边滤波和自适应阈值的Canny边缘检测方法 |
CN117422712A (zh) * | 2023-12-15 | 2024-01-19 | 青岛合丰新材料有限公司 | 基于图像滤波处理的塑料母粒视觉检测方法及系统 |
CN117830139A (zh) * | 2024-03-01 | 2024-04-05 | 天津世亚模具股份有限公司 | 一种高强度冲压件缺陷检测方法 |
CN117876971A (zh) * | 2024-03-12 | 2024-04-12 | 武汉同创万智数字科技有限公司 | 基于机器视觉的建筑施工安全监测及预警方法 |
Non-Patent Citations (3)
Title |
---|
ANNA GABIGER-ROSE等: "Image denoising using bilateral filter with noise-adaptive parameter tuning", IECON 2011 - 37TH ANNUAL CONFERENCE OF THE IEEE INDUSTRIAL ELECTRONICS SOCIETY, 2 January 2012 (2012-01-02) * |
常戬;白佳弘;: "基于回转对称双边滤波的Retinex图像增强算法", 计算机工程, no. 06, 15 June 2016 (2016-06-15) * |
韩震;王红斌;余正涛;朱映柔;: "双边非局部均值滤波图像去噪算法", 传感器与微系统, no. 06, 20 June 2016 (2016-06-20) * |
Also Published As
Publication number | Publication date |
---|---|
CN118038280B (zh) | 2024-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114418957B (zh) | 基于机器人视觉的全局和局部二值模式图像裂缝分割方法 | |
US7359555B2 (en) | Detecting roads in aerial images using feature-based classifiers | |
US8749658B2 (en) | Data processing device, image matching method, program, and image matching system | |
US8861853B2 (en) | Feature-amount calculation apparatus, feature-amount calculation method, and program | |
CN109741356B (zh) | 一种亚像素边缘检测方法及系统 | |
CN106485724B (zh) | 一种基于组合感受野和朝向特征调制的轮廓检测方法 | |
CN117593592B (zh) | 车底异物智能扫描识别系统及方法 | |
CN112819094A (zh) | 一种基于结构相似度度量的目标检测与识别方法 | |
CN116137036A (zh) | 基于机器学习的基因检测数据智能处理系统 | |
CN108108740B (zh) | 一种主动毫米波人体图像性别识别方法 | |
CN116129195A (zh) | 图像质量评价装置、方法、电子设备和存储介质 | |
CN118155149B (zh) | 一种智慧城市道路智能监控系统 | |
CN116152115A (zh) | 基于计算机视觉的垃圾图像去噪处理方法 | |
CN118038280B (zh) | 基于航拍影像的建筑施工进度监测及预警方法 | |
CN107977608B (zh) | 一种应用于公路视频图像道路区域提取的方法 | |
CN115797374B (zh) | 基于图像处理的机场跑道提取方法 | |
CN110222749B (zh) | 一种可见光图像与红外图像匹配方法 | |
CN115661464B (zh) | 图像分割方法、装置、设备及计算机存储介质 | |
CN116486480A (zh) | 一种基于点云的人体跌倒检测方法及装置 | |
CN116188990A (zh) | 一种基于无人机遥感的地表植被识别方法及系统 | |
CN115937311A (zh) | 摄像头偏移检测方法及其相关装置 | |
CN115100510A (zh) | 一种轮胎磨损度识别方法 | |
CN112734666B (zh) | 基于相似值的sar图像相干斑非局部均值抑制方法 | |
CN114049337A (zh) | 一种基于人工智能的隧道变形检测方法及系统 | |
CN115512310A (zh) | 一种基于视频监控下车脸特征的车型识别方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |