Nothing Special   »   [go: up one dir, main page]

CN107578419A - A Stereo Image Segmentation Method Based on Consistent Contour Extraction - Google Patents

A Stereo Image Segmentation Method Based on Consistent Contour Extraction Download PDF

Info

Publication number
CN107578419A
CN107578419A CN201710820101.7A CN201710820101A CN107578419A CN 107578419 A CN107578419 A CN 107578419A CN 201710820101 A CN201710820101 A CN 201710820101A CN 107578419 A CN107578419 A CN 107578419A
Authority
CN
China
Prior art keywords
image
pixel
contour
stereo
segmentation
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
Application number
CN201710820101.7A
Other languages
Chinese (zh)
Other versions
CN107578419B (en
Inventor
厉旭杰
邵靖凯
沈顺辉
黄辉
胡明晓
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wenzhou University
Original Assignee
Wenzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wenzhou University filed Critical Wenzhou University
Priority to CN201710820101.7A priority Critical patent/CN107578419B/en
Publication of CN107578419A publication Critical patent/CN107578419A/en
Application granted granted Critical
Publication of CN107578419B publication Critical patent/CN107578419B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

The embodiment of the invention discloses a kind of stereo-picture dividing method based on uniformity contours extract, including input stereo-picture left view, left view segmentation figure picture and stereo-picture right view;According to stereo-picture or so view, Stereo matching item is calculated;Stereo-picture right view segmentation probability graph is estimated according to stereo-picture left view, left view segmentation figure picture and stereo-picture right image, and extracts and correct the lines of outline of estimating of segmentation probability graph, and further calculates the profile item based on image patch;According to Stereo matching item and the profile item based on image patch, profile probability graph is synthesized;The optimal profile line of profile probability graph is calculated using geodesic curve distance;Optimal profile line is scanned, obtains the right view segmentation figure picture of stereo-picture.Implement the present invention, the stereo-picture segmentation effect of high quality can be obtained, there is faster operational efficiency again.

Description

一种基于一致性轮廓提取的立体图像分割方法A Stereo Image Segmentation Method Based on Consistent Contour Extraction

技术领域technical field

本发明涉及计算机图形图像处理技术领域,尤其涉及一种基于一 致性轮廓提取的立体图像分割方法。The invention relates to the technical field of computer graphics and image processing, in particular to a stereoscopic image segmentation method based on consistent contour extraction.

背景技术Background technique

最近几年,立体媒体越来越受到人们的关注,使得立体图像的编 辑也变得越来越重要。然而,立体图像分割是立体图像编辑的基础, 因此进行精确的立体图像分割就显得尤为迫切。In recent years, stereo media has attracted more and more attention, making the editing of stereo images more and more important. However, stereoscopic image segmentation is the basis of stereoscopic image editing, so it is particularly urgent to perform accurate stereoscopic image segmentation.

到目前位置,单个图像的分割方法已经非常成熟,例如Carsten Rother,VladimirKolmogorov,and Andrew Blake(2004)在ACM Trans. Graph.23(3):309-314上发表的论文“grabcut”:Interactive fore-ground extraction using iterated graph cuts,其提出的GrabCut方法,用户仅仅 需要在物体附近画一个矩形框或者使用着色线条表明要求的前景和 背景区域,就能获得高质量的图像分割结果,并已在OpenCV上采用 GPU加速获得了实时的性能。So far, the segmentation method of a single image has been very mature, such as the paper "grabcut" published by Carsten Rother, Vladimir Kolmogorov, and Andrew Blake (2004) on ACM Trans. Graph.23(3):309-314: Interactive fore- ground extraction using iterated graph cuts, the GrabCut method proposed by it, the user only needs to draw a rectangular frame near the object or use colored lines to indicate the required foreground and background areas, and can obtain high-quality image segmentation results, and has been used on OpenCV Real-time performance is achieved with GPU acceleration.

但是,立体图像分割不仅仅局限于单个图像分割,而是建立在对 立体图像左右视图分别采用单个图像分割操作基础上,可是这种最直 接的立体图像分割方法不仅需要大量的用户交互,而且很难保证获得 左右一致性的立体分割结果。However, stereoscopic image segmentation is not limited to single image segmentation, but is based on a single image segmentation operation for the left and right views of the stereoscopic image, but this most direct stereoscopic image segmentation method not only requires a lot of user interaction, but also very It is difficult to guarantee the stereo segmentation results with left and right consistency.

为了克服上述方法所存在的问题,众多学者已开始研究采用立体 视差图同时对立体图像左右视图进行图像分割的方法。目前,众多学 者所研究出的立体图像分割方法大致可以分为三类:基于区域的立体 图像分割方法、基于着色线条的立体图像分割方法和基于轮廓提取的 立体图像分割方法。In order to overcome the problems of the above methods, many scholars have begun to study the method of image segmentation for the left and right views of the stereo image using the stereo disparity map. At present, the stereoscopic image segmentation methods developed by many scholars can be roughly divided into three categories: stereoscopic image segmentation methods based on regions, stereoscopic image segmentation methods based on colored lines, and stereoscopic image segmentation methods based on contour extraction.

基于区域的立体图像分割方法采用左右视图像中每个像素之间 的立体相关性,根据立体图像左视图的分割结果来推导出对应的立体 图像右视图分割结果,但是缺点在于:不仅分割效果严重依赖于立体 视差图的质量,而且计算立体视差图也将消耗大量的时间。例如, Wan-Yen Lo,Jeroen van Baar,Claude Knaus,Matthias Zwicker,and MarkusGross(2010)在ACM Trans.Graph.29(6):147:1-147:10上发 表的论文Stereoscopic 3dcopy&paste,通过立体视差图建立了立体图 像左右视图的一致性关系,并根据左视图分割结果来获得一致性的立 体图像右视图分割结果。The region-based stereo image segmentation method uses the stereo correlation between each pixel in the left and right view images, and derives the corresponding right view segmentation result of the stereo image according to the segmentation result of the left view of the stereo image, but the disadvantage is that not only the segmentation effect is serious It depends on the quality of the stereo disparity map, and calculating the stereo disparity map will also consume a lot of time. For example, Wan-Yen Lo, Jeroen van Baar, Claude Knaus, Matthias Zwicker, and Markus Gross (2010) published the paper Stereoscopic 3dcopy&paste on ACM Trans.Graph.29(6):147:1-147:10, through stereo disparity The figure establishes the consistency relationship between the left and right views of the stereo image, and obtains the consistent right view segmentation result of the stereo image according to the segmentation results of the left view.

基于着色线条的立体图像分割方法不同于基于区域的立体图像 分割方法,该方法通过立体一致性,仅仅迁移用户的着色线条,然后 分别对立体图像左右视图单独进行立体图像分割,但是依然需要计算 立体相关性,以及分别对立体图像左右视图进行图像分割操作,需要 消耗大量的时间。例如,Wei Ma,Luwei Yang,Yu Zhang,and Lijuan Duan(2016)在Multimedia Tools and Applications,75(18):10935-10948 上发表的论文Fastinteractive stereo image segmentation,提出采用立 体图像的稀疏立体一致性来代替原先的稠密立体视差图来提高立体 图像分割的效率;又如,W H.Emrah Tasli andA.Aydin Alatan(2013) 在Image Commun,28(10):1374-1389上发表的论文User assisteddisparity remapping for stereo images,提出利用稀疏立体一致性,把立 体图像左视图中的着色线条自动迁移到立体图像右视图中,随后利用 基于图割的能量最优化模型对立体图像左右视图分别进行图像分割 操作来提高立体图像分割的效果。The stereoscopic image segmentation method based on colored lines is different from the region-based stereoscopic image segmentation method. This method only transfers the user's colored lines through stereo consistency, and then performs stereoscopic image segmentation on the left and right views of the stereoscopic image separately, but still needs to calculate the stereoscopic image. Correlation, as well as performing image segmentation operations on the left and right views of the stereo image respectively, consume a lot of time. For example, the paper Fastinteractive stereo image segmentation published by Wei Ma, Luwei Yang, Yu Zhang, and Lijuan Duan (2016) on Multimedia Tools and Applications, 75(18): 10935-10948 proposed to use the sparse stereo consistency of stereo images to Replace the original dense stereo disparity map to improve the efficiency of stereo image segmentation; another example, W H.Emrah Tasli and A.Aydin Alatan (2013) published a paper on Image Commun,28(10):1374-1389 User assisted disparity remapping for Stereo images, proposed to use sparse stereo consistency to automatically migrate the colored lines in the left view of the stereo image to the right view of the stereo image, and then use the energy optimization model based on graph cuts to perform image segmentation operations on the left and right views of the stereo image to improve The effect of stereoscopic image segmentation.

基于轮廓提取的立体图像分割方法还是相对比较少。例如,Ran Ju,Tongwei Ren,and GangshanWu(2015)The IEEE International Conference on Computer Vision(ICCV)上发表的论文Stereosnakes: Contour based consistent object extractionfor stereo images,提出根据立 体图像左视图分割物体的轮廓,结合立体相关性和图像色彩直方图来 推导立体图像右视图分割物体轮廓,该方法虽然提高了图像分割的效 率,但是由于采用图像色彩直方图来预测图像遮挡区域,在遮挡区域 的图像分割结果很难获得满意的质量,而且该方法假设立体图像右视 图轮廓图的像素和左视图中的轮廓像素数量是一致的。由于存在遮挡 区域,这个假设是不成立的,因此,该方法在提取遮挡区域的轮廓效 果很难令人满意。Stereo image segmentation methods based on contour extraction are still relatively few. For example, Ran Ju, Tongwei Ren, and GangshanWu (2015) published the paper Stereosnakes: Contour based consistent object extraction for stereo images on The IEEE International Conference on Computer Vision (ICCV), which proposes to segment the outline of the object according to the left view of the stereo image, combined with the stereo Correlation and image color histogram are used to derive the object contour of the right view of the stereo image. Although this method improves the efficiency of image segmentation, it is difficult to obtain the image segmentation results in the occluded area because the image color histogram is used to predict the image occlusion area. Satisfactory quality, and the method assumes that the number of pixels in the contour map of the right view of the stereo image is the same as the number of contour pixels in the left view. Due to the existence of occluded areas, this assumption is not valid, so the method is not satisfactory in extracting the contours of occluded areas.

综上所述,众多学者所研究的立体图像分割方法都无法获得高质 量的立体图像分割效果,也无法具有较快的运行效率。因此,有必要 提供一种既能获得高质量的立体图像分割效果,又有较快的运行效率 的立体图像分割方法。To sum up, the stereoscopic image segmentation methods studied by many scholars cannot obtain high-quality stereoscopic image segmentation effects, nor can they have fast operating efficiency. Therefore, it is necessary to provide a stereoscopic image segmentation method that can not only obtain high-quality stereoscopic image segmentation effect, but also has faster operating efficiency.

发明内容Contents of the invention

本发明实施例的目的在于提供一种基于一致性轮廓提取的立体 图像分割方法,既能获得高质量的立体图像分割效果,又有较快的运 行效率。The purpose of the embodiment of the present invention is to provide a stereoscopic image segmentation method based on consistent contour extraction, which can not only obtain high-quality stereoscopic image segmentation effect, but also have faster operating efficiency.

为了解决上述技术问题,本发明实施例提供了一种基于一致性轮 廓提取的立体图像分割方法,包括步骤:In order to solve the above-mentioned technical problem, the embodiment of the present invention provides a kind of stereoscopic image segmentation method based on consistent contour extraction, comprising steps:

S1、输入立体图像左视图、左视图分割图像及立体图像右视图;S1. Input the left view of the stereo image, the segmented image of the left view and the right view of the stereo image;

S2、根据所述立体图像左右视图,计算立体匹配项;S2. Calculate a stereo matching item according to the left and right views of the stereo image;

S3、根据所述立体图像左视图、左视图分割图像及立体图像右视 图,预估立体图像右视图分割概率图,并根据所述预估的立体图像右 视图分割概率图,提取及修正预估轮廓线条,且进一步根据所述修正 后的预估轮廓线条,计算基于图像补丁的轮廓项;S3. Estimating the right view segmentation probability map of the stereo image according to the left view of the stereo image, the segmented image of the left view, and the right view of the stereo image, and extracting and correcting the estimate according to the estimated right view segmentation probability map of the stereo image Contour lines, and further calculate contour items based on image patches according to the revised estimated contour lines;

S4、根据所述计算出的立体匹配项及基于图像补丁的轮廓项,合 成轮廓概率图;S4. Synthesize a contour probability map according to the calculated stereo matching item and the contour item based on the image patch;

S5、采用测地线距离计算出所述合成的轮廓概率图的最优轮廓线;S5. Calculate the optimal contour line of the synthesized contour probability map by using the geodesic distance;

S6、扫描所述计算出的最优轮廓线,获得立体图像的右视图分割 图像。S6. Scan the calculated optimal contour line to obtain the right view segmented image of the stereo image.

其中,所述步骤S2中立体匹配项是通过公式 来实现;Wherein, in the step S2, the stereo matching item is obtained by the formula to fulfill;

其中,pi,pj为像素位置索引,pi为左视图分割图像轮廓线上的 像素,为以pi为中心的m×m的局部区域,m=5,d(pi)为像素pi的视差值,CAD为绝对误差和算法,且c表示彩色图像的三个通道。Among them, p i , p j are the pixel position index, p i is the pixel on the contour line of the left view segmented image, is the m×m local area centered on p i , m=5, d(p i ) is the parallax value of pixel p i , C AD is the absolute error sum algorithm, and c represents the three channels of the color image.

其中,所述步骤S3中“预估立体图像右视图分割概率图”的具体 步骤包括:Wherein, the specific steps of "estimated stereoscopic image right view segmentation probability map" in the step S3 include:

分别构建立体图像左视图和右视图立体图像的特征向量,特征向 量的构建如下:式中,为以pi为中心的m×m的 局部区域,为以pi为中心的m×m图像补丁,m=5;Construct the eigenvectors of the stereo image left view and the right view stereo image respectively, and the eigenvectors are constructed as follows: In the formula, is a m×m local area centered on p i , is an m×m image patch centered on p i , m=5;

采用预设的搜索策略在立体图像左视图中搜索与立体图像右视 图中每个像素点最相似的图像补丁;Using a preset search strategy to search for the image patch most similar to each pixel in the stereo image right view in the left view of the stereo image;

通过公式IRmask(pi)=ILmask(pj),预估立体图像右视图分割概率图; 其中,ILmask为左视图分割图像,pi,pj为像素位置索引,pj为右视图 中像素pi在立体图像左视图中找到的最相似图像补丁的中心位置。By the formula I Rmask (p i )=I Lmask (p j ), estimate the right view segmentation probability map of the stereo image; where, I Lmask is the left view segmentation image, p i , p j are the pixel position index, p j is the right In-view pixel p i is the center position of the most similar image patch found in the left view of the stereo image.

其中,所述采用预设的搜索策略在立体图像左视图中搜索与立体 图像右视图中每个像素点最相似的图像补丁的步骤具体包括:Wherein, the step of searching for the image patch most similar to each pixel in the right view of the stereo image in the left view of the stereo image using a preset search strategy specifically includes:

在立体图像右视图中的被遮挡像素采用全局范围搜索策略搜索 最相似补丁;其中,所述全局范围搜索策略是将搜索范围确定为立体 图像左视图整个图像;The occluded pixel in the right view of the stereo image adopts a global range search strategy to search for the most similar patch; wherein, the global range search strategy is to determine the search range as the entire image of the left view of the stereo image;

在立体图像右视图中的非遮挡像素采用扫描线搜索策略搜索最 相似补丁;其中,所述扫描线搜索策略是将搜索范围确定为立体图像 左视图沿着相同水平线在[0,dmax]视差范围内;dmax为预估立体图像最 大视差值;The non-occluded pixels in the right view of the stereo image use a scan line search strategy to search for the most similar patch; wherein, the scan line search strategy is to determine the search range as the left view of the stereo image along the same horizontal line at [0,d max ] parallax Within the range; d max is the maximum parallax value of the estimated stereo image;

其中,立体图像右视图中被遮挡像素和非遮挡像素的区分是通过 以下步骤实现,具体包括:采用扫描线搜索策略在立体图像左视图中 搜索和立体图像右视图中某一像素Pr最相似补丁的像素,搜索到的 像素设为Pl,然后同样采用扫描线搜索策略在立体图像右视图中搜索Among them, the distinction between the occluded pixel and the non-occluded pixel in the right view of the stereo image is realized through the following steps, which specifically include: using the scan line search strategy to search in the left view of the stereo image and the most similar pixel P r in the right view of the stereo image The pixel of the patch, the searched pixel is set to P l , and then the scan line search strategy is also used to search in the right view of the stereo image

其中,所述步骤S3中“根据所述预估的立体图像右视图分割概率 图,提取及修正预估轮廓线条”的具体步骤包括:Wherein, in the step S3, the specific steps of "according to the estimated right view segmentation probability map of the stereoscopic image, extracting and correcting the estimated contour lines" include:

采用逐行扫描立体图像右视图分割概率图,并待扫描到立体图像 右视图分割概率图中某一像素为非零且以该像素为中心的3×3区域 内至少包含一个非零像素后,则保留该像素点为轮廓点,且进一步将 所保留的轮廓点绘制成预估轮廓线条;Use progressive scanning of the right-view segmentation probability map of the stereoscopic image, and after scanning to a certain pixel in the right-view segmentation probability map of the stereoscopic image is non-zero and contains at least one non-zero pixel in the 3×3 area centered on this pixel, The pixel point is reserved as a contour point, and the retained contour point is further drawn into an estimated contour line;

统计预估轮廓线条中封闭曲线的个数,并按照封闭曲线围成的面 积从小到大排列;Count and estimate the number of closed curves in the contour lines, and arrange them according to the area enclosed by the closed curves from small to large;

设定面积阈值Ah=βAk,移除预估轮廓线条中面积小于阈值Ah的 封闭曲线,得到修正后的预估轮廓线条;其中,β为权重参数,设定 为0.8,Ak为已排完序的第k个封闭曲线的面积,k为左视图分割图 像中封闭曲线的个数。Set the area threshold A h = βA k , remove the closed curves whose area is smaller than the threshold A h in the estimated contour lines, and obtain the revised estimated contour lines; where β is the weight parameter, set to 0.8, and A k is The area of the kth closed curve that has been sorted, k is the number of closed curves in the left view segmented image.

其中,所述步骤S3中“根据所述修正后的预估轮廓线条,计算基于图像补 丁的轮廓项”是通过公式来 实现;其中,dc(qi)为像素点qi到预估轮廓线条的最短距离,D为权重参数,设 定为D=10,Td为阈值,设定为20;Cr为修正后的预估轮 廓线条,d(qi,qj)为两个像素之间的欧拉距离。Wherein, in the step S3 "according to the estimated contour line after the correction, calculate the contour item based on the image patch" is through the formula To achieve; wherein, d c (q i ) is the shortest distance from the pixel point q i to the estimated contour line, D is the weight parameter, set to D=10, T d is the threshold, set to 20; C r is the corrected estimated contour line, d(q i , q j ) is the Euler distance between two pixels.

其中,所述步骤S4中“轮廓概率图”是通过公式 CC(pi,d(pi))=CS(pi,d(pi))×CO(pi,d(pi))来实现;其中,CS为立体匹配项,CO为基 于图像补丁的轮廓项。Wherein, the "contour probability map" in the step S4 is calculated by the formula C C (p i ,d(p i ))=C S (p i ,d(p i ))×C O (p i ,d(p i ) i )) to realize; among them, C S is a stereo matching item, and C O is a contour item based on an image patch.

其中,所述步骤S5中的“最优轮廓线”是通过求解公式最优解来实现;其中,Cm为 左视图分割图像轮廓线,N(pi,pi-1)为轮廓连续性项,当像素pi-1为像 素pi的8邻域时,轮廓连续性项的值为0,否则轮廓连续性项的值为 ∞。Wherein, the "optimum contour line" in the step S5 is obtained by solving the formula The optimal solution is realized; where, C m is the contour line of the left view segmentation image, N(p i , p i-1 ) is the contour continuity item, when the pixel p i-1 is the 8 neighbors of the pixel p i , The value of the contour continuity item is 0, otherwise the value of the contour continuity item is ∞.

其中,所述步骤S6具体包括:Wherein, the step S6 specifically includes:

依次逐行从左到右扫描,如果像素值为0,表示该点为背景,直 到扫描到像素值为1的点,表示该点为左轮廓点,继续扫描直到扫描 到像素值为0的点,表示左轮廓线结束,该点为前景,接着继续扫描 直到再次扫描到像素值为1的点且其对应下一个像素点的像素值为0 时,表示该点为右轮廓点,Scan from left to right line by line. If the pixel value is 0, it means that the point is the background. Until the point with a pixel value of 1 is scanned, it means that the point is a left contour point. Continue to scan until the point with a pixel value of 0 is scanned. , indicating the end of the left contour line, the point is the foreground, and then continue to scan until a point with a pixel value of 1 is scanned again and the pixel value corresponding to the next pixel point is 0, indicating that the point is a right contour point,

将左轮廓点和右轮廓点之间的像素均填充为1,重复此步骤,直 到整个最优轮廓线填充完整,获得立体图像的右视图分割图像。Fill the pixels between the left contour point and the right contour point as 1, repeat this step until the entire optimal contour line is filled completely, and obtain the right view segmentation image of the stereo image.

实施本发明实施例,具有如下有益效果:Implementing the embodiment of the present invention has the following beneficial effects:

本发明利用立体图像的极线约束,在能量最优化框架中除了加入 了立体匹配项,同时加入了基于图像补丁的轮廓项,根据这两项合成 了轮廓概率图,随后采用测地线距离计算最优轮廓线。提出的基于一 致性轮廓提取的立体图像分割方法能对立体图像进行一致性的图像 分割,特别是在图像遮挡区域,也能取得很好的立体分割效果,从而 达到既能获得高质量的立体图像分割效果,又有较快运行效率的目的。The present invention utilizes the epipolar constraints of the stereo image, and in the energy optimization framework, in addition to the stereo matching item, the contour item based on the image patch is added at the same time, and the contour probability map is synthesized according to these two items, and then the geodesic distance is used to calculate optimal contour. The proposed stereoscopic image segmentation method based on consistent contour extraction can perform consistent image segmentation on stereoscopic images, especially in image occlusion areas, and can also achieve good stereoscopic segmentation results, so as to achieve both high-quality stereoscopic images Segmentation effect, and the purpose of faster operation efficiency.

附图说明Description of drawings

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面 将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而 易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域 普通技术人员来讲,在不付出创造性劳动性的前提下,根据这些附图 获得其他的附图仍属于本发明的范畴。In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only These are some embodiments of the present invention. For those of ordinary skill in the art, obtaining other drawings based on these drawings still belongs to the scope of the present invention without any creative effort.

图1为本发明实施例提供的一种基于一致性轮廓提取的立体图 像分割方法的流程图;Fig. 1 is a flow chart of a stereoscopic image segmentation method based on consistent contour extraction provided by an embodiment of the present invention;

图2为本发明立体图像分割技术流程框架图;Fig. 2 is a frame diagram of the process flow of the stereoscopic image segmentation technology of the present invention;

图3为本发明提出的两种搜索最相似的图像补丁的策略图,全局 范围搜索策略和扫描线搜索策略;Fig. 3 is the strategy diagram of two kinds of search most similar image patches that the present invention proposes, global range search strategy and scan line search strategy;

图4为本发明采用不同代价匹配项的最优化结果图;Fig. 4 is the optimization result figure that the present invention adopts different cost matching items;

图5为本发明与目前主流的立体图像分割方法效果对比图;Fig. 5 is a comparison diagram between the present invention and the current mainstream stereoscopic image segmentation method;

图6为本发明对多组立体图像进行图像分割的结果图,其中第一 至三列依次为立体图像左视图、右视图和左视图分割图像,第四列和 第五列为本发明方法获得右视图最优轮廓图和最优右视图分割图像, 第六列为标准的右视图分割图像。Fig. 6 is the result figure that the present invention carries out image segmentation to multiple groups of stereoscopic images, wherein the first to third columns are the stereoscopic image left view, right view and left view segmentation images successively, the fourth column and the fifth column are obtained by the method of the present invention Right view optimal contour map and optimal right view segmented image, the sixth column is the standard right view segmented image.

具体实施方式detailed description

为了使本发明的目的、技术方案及优点更加清楚明白,以下结合 附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描 述的具体实施例仅仅用以解释本发明,并不用于限定本发明。In order to make the object, technical solution and advantages of the present invention clearer, the present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.

如图1所示,为本发明实施例中,提出的一种基于一致性轮廓提 取的立体图像分割方法,包括步骤:As shown in Figure 1, in the embodiment of the present invention, a kind of stereoscopic image segmentation method that proposes based on consistent outline extraction, comprises steps:

步骤S1、输入立体图像左视图、左视图分割图像及立体图像右 视图;Step S1, input stereoscopic image left view, left view segmentation image and stereoscopic image right view;

具体过程为,输入立体图像左视图IL,左视图分割图像ILmask, 立体图像右视图IRThe specific process is as follows: input stereo image left view I L , left view segmented image I Lmask , and stereo image right view I R .

步骤S2、根据所述立体图像左右视图,计算立体匹配项;Step S2, calculating stereo matching items according to the left and right views of the stereo image;

具体过程为,通过公式计算立体 匹配项CSThe specific process is, through the formula Calculate the stereo matching term C S ;

其中,pi,pj为像素位置索引,pi为左视图分割图像轮廓线上的 像素,为以pi为中心的m×m的局部区域,m=5,d(pi)为像素pi的视差值,CAD为绝对误差和算法,且c表示彩色图像的三个通道。Among them, p i , p j are the pixel position index, p i is the pixel on the contour line of the left view segmented image, is the m×m local area centered on p i , m=5, d(p i ) is the parallax value of pixel p i , C AD is the absolute error sum algorithm, and c represents the three channels of the color image.

步骤S3、根据所述立体图像左视图、左视图分割图像及立体图 像右视图,预估立体图像右视图分割概率图,并根据所述预估的立体 图像右视图分割概率图,提取及修正预估轮廓线条,且进一步根据所 述修正后的预估轮廓线条,计算基于图像补丁的轮廓项;Step S3: Estimating the right view segmentation probability map of the stereo image according to the left view of the stereo image, the segmented image of the left view and the right view of the stereo image, and extracting and correcting the prediction probability map according to the estimated right view segmentation probability map of the stereo image. Estimating contour lines, and further calculating contour items based on image patches according to the revised estimated contour lines;

具体过程为,(a)预估立体图像右视图分割概率图;The specific process is, (a) predicting the right view segmentation probability map of the stereo image;

(i)分别构建立体图像左视图和右视图立体图像的特征向量, 特征向量的构建如下:式中,为以pi为中心的m×m 的局部区域,为以pi为中心的m×m图像补丁,m=5;(i) Construct the feature vectors of the left view of the stereo image and the stereo image of the right view respectively, and the construction of the feature vector is as follows: In the formula, is a local area of m×m centered on p i , is an m×m image patch centered on p i , m=5;

(ii)采用预设的搜索策略(两种搜索策略:全局范围搜索策略 和扫描线搜索策略)在立体图像左视图中搜索与立体图像右视图中每 个像素点最相似的图像补丁;(ii) using a preset search strategy (two search strategies: a global range search strategy and a scan line search strategy) to search for the image patch most similar to each pixel in the stereo image right view in the left view of the stereo image;

右视图中的被遮挡的像素采用全局范围搜索策略搜索最相似补 丁,右视图中的非遮挡像素采用扫描线搜索策略搜索最相似补丁;即 全局范围搜索策略是将搜索范围确定为立体图像左视图整个图像;扫 描线搜索策略是将搜索范围确定为立体图像左视图沿着相同水平线 在[0,dmax]视差范围内;dmax为预估立体图像最大视差值;The occluded pixels in the right view use the global range search strategy to search for the most similar patch, and the non-occluded pixels in the right view use the scan line search strategy to search for the most similar patch; that is, the global range search strategy is to determine the search range as the left view of the stereo image The whole image; the scan line search strategy is to determine the search range as the left view of the stereo image along the same horizontal line within the parallax range of [0, d max ]; d max is the maximum parallax value of the estimated stereo image;

立体图像右视图中被遮挡像素和非遮挡像素的区分是通过以下 步骤实现,具体包括:采用扫描线搜索策略在立体图像左视图中搜索 和立体图像右视图中某一像素Pr最相似补丁的像素,搜索到的像素 设为Pl,然后同样采用扫描线搜索策略在立体图像右视图中搜索和立 体图像左视图中像素Pl最相似补丁的像素,搜索到的像素设为P′r;如 果像素Pr和像素P′r为同一点,则像素Pr为非遮挡像素;反之,如果 像素Pr和像素P′r为不同点,则像素Pr为被遮挡像素。The distinction between occluded pixels and non-occluded pixels in the right view of the stereo image is achieved through the following steps, which specifically include: using the scan line search strategy to search for the most similar patch to a certain pixel P r in the left view of the stereo image Pixel, the searched pixel is set as P l , and then the same scan line search strategy is used to search for the pixel of the most similar patch to the pixel P l in the left view of the stereo image in the right view of the stereo image, and the searched pixel is set as P′ r ; If the pixel P r and the pixel P' r are the same point, then the pixel P r is a non-occluded pixel; otherwise, if the pixel P r and the pixel P' r are different points, then the pixel P r is an occluded pixel.

(iii)通过公式IRmask(pi)=ILmask(pj),预估立体图像右视图分割概 率图;其中,ILmask为左视图分割图像,pi,pj为像素位置索引,pj为 右视图中像素pi在立体图像左视图中找到的最相似图像补丁的中心 位置。(iii) By the formula I Rmask (p i )=I Lmask (p j ), estimate the right view segmentation probability map of the stereo image; where, I Lmask is the left view segmentation image, p i , p j are the pixel position index, p j is the center position of the most similar image patch found by pixel p i in the right view in the left view of the stereo image.

(b)根据预估立体图像右视图分割概率图,提取预估轮廓线条 以及修正预估轮廓线条;(b) According to the estimated stereo image right view segmentation probability map, extract the estimated contour lines and correct the estimated contour lines;

(i)采用逐行扫描立体图像右视图分割概率图,并待扫描到立 体图像右视图分割概率图中某一像素为非零且以该像素为中心的 3×3区域内至少包含一个非零像素后,则保留该像素点为轮廓点,且 进一步将所保留的轮廓点绘制成预估轮廓线条;(i) The right-view segmentation probability map of the stereoscopic image is scanned progressively, and a certain pixel in the right-view segmentation probability map of the stereoscopic image to be scanned is non-zero and contains at least one non-zero pixel in the 3×3 area centered on the pixel After the pixel, the pixel point is reserved as a contour point, and the reserved contour point is further drawn into an estimated contour line;

(ii)统计预估轮廓线条中封闭曲线的个数,并按照封闭曲线围 成的面积从小到大排列;(ii) Statistically estimate the number of closed curves in the contour lines, and arrange them according to the area enclosed by the closed curves from small to large;

(iii)设定面积阈值Ah=βAk,移除预估轮廓线条中面积小于阈值 Ah的封闭曲线,得到修正后的预估轮廓线条;其中,β为权重参数, 设定为0.8,Ak为已排完序的第k个封闭曲线的面积,k为左视图分 割图像中封闭曲线的个数。(iii) Set the area threshold A h = βA k , remove the closed curves whose area is smaller than the threshold A h in the estimated contour lines, and obtain the revised estimated contour lines; where β is a weight parameter, set to 0.8, A k is the area of the kth closed curve that has been sorted, and k is the number of closed curves in the left view segmented image.

(c)根据预估轮廓线条,计算基于图像补丁的轮廓项;(c) Calculate the contour item based on the image patch according to the estimated contour line;

根据预估轮廓线条,计算基于图像补丁的轮廓项COAccording to the estimated contour line, calculate the contour item C O based on the image patch:

式中,dc(qi)为像素点qi到预估轮廓线条的最短距离,D为权重 参数,设定为D=10,Td为阈值,设定为20;如果像素点qi到预估轮 廓线条的最短距离超过了设定的阈值Td,意味着该像素点不可能成为 最终的轮廓点,dc(qi)的值为无穷大;In the formula, d c (q i ) is the shortest distance from the pixel point q i to the estimated contour line, D is the weight parameter, set to D=10, T d is the threshold, set to 20; if the pixel point q i The shortest distance to the estimated contour line exceeds the set threshold T d , which means that the pixel point cannot become the final contour point, and the value of d c (q i ) is infinite;

像素点qi到预估轮廓线条的最短距离dc(qi)的计算方式如下:The calculation method of the shortest distance d c (q i ) from the pixel point q i to the estimated contour line is as follows:

式中,Cr为立体图像右视图中修正后的预 估轮廓,d(qi,qj)为两个像素之间的欧拉距离。 In the formula, C r is the corrected estimated contour in the right view of the stereo image, and d(q i , q j ) is the Euler distance between two pixels.

步骤S4、根据所述计算出的立体匹配项及基于图像补丁的轮廓 项,合成轮廓概率图;Step S4, according to the stereo matching item calculated and the contour item based on the image patch, synthesize a contour probability map;

具体过程为,根据立体匹配项和基于图像补丁的轮廓项,合成轮 廓概率图CCThe specific process is to synthesize the contour probability map C C according to the stereo matching item and the contour item based on the image patch:

CC(pi,d(pi))=CS(pi,d(pi))×CO(pi,d(pi))C C (p i ,d(p i ))=C S (p i ,d(p i ))×C O (p i ,d(p i ))

式中,CS为立体匹配项,该项能精确的计算左右视图中未遮挡区 域一致性的轮廓,但是逐像素的匹配没有考虑到遮挡问题;CO为基 于图像补丁的轮廓项,该项估计了右视图中像素成为轮廓点的概率, 特别是能很好的检测遮挡区域成为轮廓点的概率。In the formula, C S is the stereo matching item, which can accurately calculate the consistent contour of the unoccluded area in the left and right views, but the pixel-by-pixel matching does not take into account the occlusion problem; C O is the contour item based on the image patch, the item The probability that the pixel in the right view becomes a contour point is estimated, especially the probability that the occluded area becomes a contour point can be detected very well.

步骤S5、采用测地线距离计算出所述合成的轮廓概率图的最优 轮廓线;Step S5, using geodesic distance to calculate the optimal contour of the synthetic contour probability map;

具体过程为,采用测地线距离计算最优轮廓图为求解如下最优化 问题:The specific process is to use the geodesic distance to calculate the optimal contour map to solve the following optimization problem:

即通过求解公式最优解 来实现;其中,Cm为左视图分割图像轮廓线,N(pi,pi-1)为轮廓连续性 项,因为封闭轮廓具有连续性,因此当像素pi-1为像素pi的8邻域时, 轮廓连续性项的值为0,否则轮廓连续性项的值为∞,由此可见,最 优化问题的求解转化为找到视差值d(pi),使得能量值E(pi,d(pi))取得 最小值;需要说明的是,在本发明方法中左视图分割图像轮廓线上的 点和需要提取的右视图分割图像轮廓线上的点允许存在一对多的对 应关系。采用测地线距离计算轮廓概率图的最优轮廓线,其包括三个 步骤:That is, by solving the formula The optimal solution is realized; among them, C m is the contour line of the left view segmented image, N(p i , p i-1 ) is the contour continuity item, because the closed contour has continuity, so when the pixel p i-1 is a pixel When there are 8 neighbors of p i , the value of the contour continuity item is 0, otherwise the value of the contour continuity item is ∞. It can be seen that the solution of the optimization problem is transformed into finding the disparity value d(p i ), so that the energy The value E(p i , d(p i )) takes the minimum value; it should be noted that in the method of the present invention, the points on the contour line of the left view segmented image and the points on the contour line of the right view segmented image that need to be extracted are allowed to exist One-to-many correspondence. The optimal contour line of the contour probability map is calculated using geodesic distance, which includes three steps:

(i)确定轮廓起止点;修正后的每个封闭轮廓曲线的质心记为 M=(ML1,ML2,…,MLm),首先沿着质心的横坐标搜索合成轮廓概率图 的最小值,该点标记为轮廓线的起点pm,为了确定轮廓的起始位置, 采用广度优先搜索策略进行搜索,该搜索策略分为四个步骤:(1)初 始化队列仅包含一个像素pt=pm;(2)搜索像素pt的8邻域中最小的 和次小的像素(不包括像素pt),分别标记为ps和pe;(3)如果三个像 素pt,ps和pe的纵坐标互不相等,即pt[y]≠ps[y]≠pe[y],其中[y]为提 取像素的纵坐标;则轮廓线的起始点pstart设置为ps,终点pend设置为 pe,程序结束;否则,删除队列队首像素,把像素ps和pe插入队列, 进入程序第4步;(4)队列队首像素出列赋值给像素pt,重复步骤2 直到队列为空;(i) Determine the starting and ending points of the contour; the centroid of each closed contour curve after correction is recorded as M=(M L1 , M L2 ,...,M Lm ), first search for the minimum value of the composite contour probability map along the abscissa of the centroid , this point is marked as the starting point p m of the contour line. In order to determine the starting position of the contour line, a breadth-first search strategy is used to search. The search strategy is divided into four steps: (1) The initialization queue only contains one pixel p t =p m ; (2) search for the smallest and second smallest pixel (not including pixel p t ) in the 8 neighborhoods of pixel p t , and mark them as p s and pe respectively; (3) if three pixels p t , p s The vertical coordinates of p and p e are not equal to each other, that is, p t [y]≠p s [y]≠p e [y], where [y] is the vertical coordinate of the extracted pixel; then the starting point p start of the contour line is set to p s , the end point p end is set to p e , and the program ends; otherwise, delete the first pixel of the queue, insert the pixels p s and p e into the queue, and enter the fourth step of the program; (4) dequeue the first pixel of the queue and assign it to the pixel p t , repeat step 2 until the queue is empty;

(ii)计算轮廓线的起点pstart到轮廓线的终点pend的测地线距离;(ii) calculate the geodesic distance from the starting point p start of the contour line to the end point p end of the contour line;

在计算轮廓线的起点pstart到轮廓线的终点pend的测地线距离之前, 轮廓概率图必须为非封闭区域,可以采用如下策略:如果像素和pt邻近且该像素的纵坐标和pt相同,则该点在轮廓概率图中的值设置为 无穷大;Before calculating the geodesic distance from the starting point p start of the contour line to the end point p end of the contour line, the contour probability map must be an unclosed area, and the following strategy can be adopted: if the pixel is adjacent to p t and the ordinate of the pixel and p t is the same, the value of the point in the contour probability map is set to infinity;

合成轮廓概率图中两个相邻像素之间的距离计算如下:The distance between two adjacent pixels in the synthetic contour probability map is calculated as follows:

dH(pi,pj)=|CC(pi,d(pi))-CC(pj,d(pj))|d H (p i ,p j )=|C C (p i ,d(p i ))-C C (p j ,d(p j ))|

沿着曲线C={p1,p2,L,pm}(前后两点位于8邻域中),轮廓线的起 点pstart到轮廓线的终点pend的测地线距离计算如下:Along the curve C={p 1 ,p 2 ,L,p m } (two points before and after are located in the neighborhood of 8), the geodesic distance from the starting point p start of the contour line to the end point p end of the contour line is calculated as follows:

测地线距离则定义为连接起点到终点所有曲线的距离的最小值Geodesic distance is defined as the minimum distance of all curves connecting the start point to the end point

因为测地点距离获得的曲线是所有连接起点和终点曲线的能量 最小值,且曲线前后两点之间是连续的,所以计算获得的测地线距离 所经过的曲线则为右视图轮廓图;轮廓点的值标记为1,其它点的值 标记为0;Because the curve obtained by the geodesic distance is the minimum energy value of all the curves connecting the starting point and the end point, and the two points before and after the curve are continuous, so the curve obtained by calculating the geodesic distance is the right-view profile; the profile The value of the point is marked as 1, and the value of other points is marked as 0;

步骤S6、扫描所述计算出的最优轮廓线,获得立体图像的右视 图分割图像。Step S6, scanning the calculated optimal contour line to obtain the right view segmented image of the stereo image.

具体过程为,依次逐行从左到右扫描,如果像素值为0,表示该 点为背景,直到扫描到像素值为1的点,表示该点为左轮廓点,继续 扫描直到扫描到像素值为0的点,表示左轮廓线结束,该点为前景, 接着继续扫描直到再次扫描到像素值为1的点且其对应下一个像素 点的像素值为0时,表示该点为右轮廓点,The specific process is to scan from left to right line by line. If the pixel value is 0, it means that the point is the background, until the point with a pixel value of 1 is scanned, it means that the point is the left contour point, and continue scanning until the pixel value is scanned. A point of 0 indicates the end of the left contour line, which is the foreground point, and then continues to scan until a point with a pixel value of 1 is scanned again and the pixel value corresponding to the next pixel point is 0, indicating that the point is a right contour point ,

将左轮廓点和右轮廓点之间的像素均填充为1,重复此步骤,直 到整个最优轮廓线填充完整,获得立体图像的右视图分割图像。Fill the pixels between the left contour point and the right contour point as 1, repeat this step until the entire optimal contour line is filled completely, and obtain the right view segmentation image of the stereo image.

如图2所示,对本发明实施例中的基于一致性轮廓提取的立体图 像分割方法工作原理做进一步说明,图2为立体图像分割技术流程框 架图,该方法可以概括为以下六个个步骤:首先用户输入立体图像左 右视图以及立体图像左视图分割图;然后根据输入的立体图像左右视 图,计算立体匹配项;以及根据输入的立体图像左右视图和左视图分 割图像,计算基于图像补丁的轮廓项;然后根据立体匹配项和基于图 像补丁的轮廓项,合成轮廓概率图;随后采用测地线距离计算最优轮 廓线;最后根据轮廓线,获得立体图像右视图分割图。As shown in Figure 2, the working principle of the stereoscopic image segmentation method based on consistent contour extraction in the embodiment of the present invention is further described. Figure 2 is a frame diagram of the stereoscopic image segmentation technical process. The method can be summarized as the following six steps: First, the user inputs the left and right views of the stereo image and the segmentation map of the left view of the stereo image; then calculates the stereo matching item according to the input left and right views of the stereo image; and calculates the contour item based on the image patch according to the input left and right views of the stereo image and the left view segmentation image ; Then according to the stereo matching item and the contour item based on the image patch, the contour probability map is synthesized; then the optimal contour line is calculated by using the geodesic distance; finally, the right view segmentation map of the stereo image is obtained according to the contour line.

图3为本发明提出的两种搜索最相似的图像补丁策略的应用场 景图,包括全局范围搜索策略和扫描线搜索策略。Fig. 3 is the application scenario figure of two kinds of search most similar image patch strategies that the present invention proposes, including global range search strategy and scan line search strategy.

为了验证本发明方法能量最优化框架中立体匹配项和基于图像 补丁的轮廓项对提取立体图像右视图分割图的影响,图4显示了本发 明采用不同代价匹配项的最优化结果。从图4中可以看出,仅仅使用 立体匹配项的最优化能量框架针对未遮挡区域,立体图像左右视图的 轮廓线能很好的匹配,但是针对遮挡区域,因为右视图中的像素在左 视图中不能找到匹配点,因此立体匹配项在遮挡区域出现了很大的偏 差,相比较本发明方法提出的基于图像补丁的轮廓项能很好的处理遮 挡区域的轮廓线提取,而合并立体匹配项和基于图像补丁的轮廓项则 更进一步提升了右视图最优轮廓线的精度,获得了高质量的立体图像 右视图分割结果。In order to verify the influence of the stereo matching item and the image patch-based contour item in the energy optimization framework of the present invention on extracting the right view segmentation map of the stereo image, Fig. 4 shows the optimization results of the present invention using different cost matching items. As can be seen from Figure 4, the optimized energy framework using only stereo matching items can match the contour lines of the left and right views of the stereo image well for the unoccluded area, but for the occluded area, because the pixels in the right view are in the left view Matching points cannot be found in the occluded area, so the stereo matching item has a large deviation in the occluded area. Compared with the contour item based on the image patch proposed by the method of the present invention, it can handle the contour line extraction of the occluded area well, and the combined stereo matching item And the contour item based on the image patch further improves the accuracy of the optimal contour line in the right view, and obtains high-quality stereo image right view segmentation results.

为了定量计算立体图像分割的效果,我们计算了本方法获得的最 优右视图分割图和标准右视图分割图之间的误判率,误判率 PUM(Percentage of unmatching pixels)的计算公式如下:In order to quantitatively calculate the effect of stereoscopic image segmentation, we calculated the misjudgment rate between the optimal right-view segmentation map obtained by this method and the standard right-view segmentation map. The calculation formula of the false positive rate PUM (Percentage of unmatching pixels) is as follows:

式中,p为像素位置索引,N为整个立体右视图像素的个数, MS(p)为本方法获得的最优右视图分割图,MG(p)为标准右视图分割 图,#为计数操作,统计满足条件|MS(p)-MG(p)|≠0的像素个数。误 判率越低,表示获得的最优右视图分割图和标准右视图分割图越接近。In the formula, p is the pixel position index, N is the number of pixels in the entire stereoscopic right view, M S (p) is the optimal right view segmentation map obtained by this method, M G (p) is the standard right view segmentation map, # For the counting operation, count the number of pixels satisfying the condition |MS (p) -M G (p)|≠0. The lower the misjudgment rate, the closer the obtained optimal right view segmentation map is to the standard right view segmentation map.

图5显示了本发明与目前主流的立体图像分割方法效果对比。轮 廓跟踪方法,例如Contour tracking(CT)目标为跟踪物体的运动,该方 法对跟踪物体轮廓的要求不需要太精确;视频分割方法,例如 Snapcut(SN)是以视频中物体运动为特性进行的物体分割;协同分割 方法,例如Co-segmentation(CO)通过两个或者多个不同视角的图片, 进行协同分割,这三种方法都是以多张照片或者视频中物体运动的特 性进行图像分割,而立体图像则具有更加严格的限制,比如极线约束。 从图5中可以看出这三种方法的误判率都较高,相比较立体图像分割 方法,例如Stereo cut(SC)和Stereo snake(SS)在立体图像分割的时候都加入了极线约束,Stereo cut(SC)是基于区域的立体分割方法,该方 法的图像分割效果依赖于概率密度函数,而Stereo snake(SS)和本发明 方法则是基于轮廓提取的立体分割方法,相比于Stereo snake(SS)使用 简单的立体图像色彩直方图来推导立体图像右视图分割物体轮廓,本 发明方法提出的基于图像补丁的轮廓项能更好的处理图像遮挡区域 的轮廓线提取。Fig. 5 shows the effect comparison between the present invention and the current mainstream stereoscopic image segmentation method. Contour tracking methods, such as Contour tracking (CT) aim to track the movement of objects, and this method does not need to be too precise in tracking the contours of objects; video segmentation methods, such as Snapcut (SN), are objects that are characterized by the movement of objects in the video Segmentation; Cooperative segmentation methods, such as Co-segmentation (CO), perform collaborative segmentation through two or more pictures of different perspectives. These three methods all perform image segmentation based on the characteristics of object motion in multiple photos or videos, and Stereoscopic images have stricter restrictions, such as epipolar constraints. It can be seen from Figure 5 that the misjudgment rates of these three methods are relatively high. Compared with stereo image segmentation methods, such as Stereo cut (SC) and Stereo snake (SS), epipolar constraints are added to stereo image segmentation. , Stereo cut (SC) is a stereoscopic segmentation method based on regions, and the image segmentation effect of this method depends on the probability density function, while Stereo snake (SS) and the method of the present invention are stereoscopic segmentation methods based on contour extraction, compared to Stereo snake (SS) uses a simple stereoscopic image color histogram to derive the contour of the right-view segmented object of the stereoscopic image. The contour item based on the image patch proposed by the method of the present invention can better handle the contour extraction of the image occlusion area.

其中,Contour tracking(CT)方法由Chih-Yuan Chung and Homer H Chen 2010年发表在IEEE Transactions on Circuits and Systems for Video Technology 20卷1期149-155页上的论文Video object extraction via mrf-based contour tracking提出。Among them, the Contour tracking (CT) method was published by Chih-Yuan Chung and Homer H Chen in 2010 in the paper Video object extraction via mrf-based contour tracking in IEEE Transactions on Circuits and Systems for Video Technology Volume 20, Issue 1, Page 149-155 propose.

Snapcut(SN)方法由Xue Bai,Jue Wang,David Simons,and Guillermo Sapiro2009年发表在ACM Trans.Graph.28卷3期 70:1-70:11页上的论文Video snapcut:Robustvideo object cutout using localized classifiers提出。The Snapcut(SN) method was published by Xue Bai, Jue Wang, David Simons, and Guillermo Sapiro in 2009 in the paper Video snapcut: Robustvideo object cutout using localized classifiers on ACM Trans.Graph.28 Volume 3 Issue 70:1-70:11 propose.

Co-segmentation(CO)方法由D.Batra,A.Kowdle,D.Parikh,J.Luo, and T.Chen2010年发表在In Computer Vision and Pattern Recognition (CVPR),2010IEEEConference on,3169-3176页上的论文icoseg: Interactive co-segmentation withintelligent scribble guidance提出。The Co-segmentation (CO) method was published by D.Batra, A.Kowdle, D.Parikh, J.Luo, and T.Chen in 2010 in In Computer Vision and Pattern Recognition (CVPR), 2010IEEEConference on, pages 3169-3176 The paper icoseg: Interactive co-segmentation with intelligent scribble guidance was proposed.

Stereo cut(SC)方法由B.L.Price and S.Cohen 2011年发表在In2011International Conference on Computer Vision,1148-1155页上的论 文Stereocut:Consistent interactive object selection in stereo image pairs 提出。The Stereo cut (SC) method was proposed by B.L.Price and S.Cohen in the paper Stereocut: Consistent interactive object selection in stereo image pairs published in In2011International Conference on Computer Vision, pages 1148-1155 in 2011.

Stereo snake(SS)方法由Ran Ju,Tongwei Ren,and Gangshan Wu 2015年发表在In The IEEE International Conference on Computer Vision(ICCV),1724-1732页上的论文Stereosnakes:Contour based consistent object extraction for stereo images提出。The Stereo snake (SS) method was proposed by Ran Ju, Tongwei Ren, and Gangshan Wu in the paper Stereosnakes: Contour based consistent object extraction for stereo images published in In The IEEE International Conference on Computer Vision (ICCV), 1724-1732 in 2015 .

为了定量的计算本发明方法和主流的立体图像分割方法的误判 率,我们测试了31组立体图像对(立体图像左右视图以及标准的立体 图像左右视图分割图来源于Middlebury数据集(参见 http://vision.middlebury.edu/stereo/)和Adobe开放数据集(参见 https://sourceforge.net/adobe/adobedatasets/home/Home/),这些数据集 具有非常强的纹理特性而且包换了不同深度的物体),定量的立体图 像分割精度见表1。In order to quantitatively calculate the misjudgment rate of the inventive method and the mainstream stereoscopic image segmentation method, we tested 31 groups of stereoscopic image pairs (the stereoscopic image left and right views and the standard stereoscopic image left and right view segmentation map come from the Middlebury data set (see http: //vision.middlebury.edu/stereo/) and Adobe open datasets (see https://sourceforge.net/adobe/adobedatasets/home/Home/), these datasets have very strong texture features and include different depth objects), the quantitative stereo image segmentation accuracy is shown in Table 1.

表1Table 1

从表1中可以看出本发明方法、Stereo cut(SC)和Stereo snake(SS) 方法由于采用了更为严格的极线约束,因此误判率要低于其它三种方 法,需要特别说明的是,本发明方法的误判率是最低的。因此本发明 方法获得了最高质量的立体图像分割效果。It can be seen from Table 1 that the method of the present invention, the Stereo cut (SC) and the Stereo snake (SS) methods adopt stricter epipolar constraints, so the misjudgment rate is lower than the other three methods, which need special explanation Yes, the misjudgment rate of the method of the present invention is the lowest. Therefore, the method of the present invention obtains the highest quality stereoscopic image segmentation effect.

图6为本发明对多组立体图像进行图像分割的结果,其中第一至 三列依次为立体图像左视图、右视图和左视图分割图像,第四列和第 五列为本发明方法获得右视图最优轮廓图和最优右视图分割图像,第 六列为标准的右视图分割图像。从图6中可以看出,本发明方法针对 遮挡区域、色彩重叠以及摄像头运动的立体图像分割均取得了高质量 的立体分割结果。Fig. 6 is the result that the present invention carries out image segmentation to multiple sets of stereoscopic images, wherein the first to third columns are the left view, right view and left view segmentation images of the stereoscopic images, the fourth column and the fifth column are the right images obtained by the method of the present invention View optimal contour map and optimal right view segmented image, the sixth column is the standard right view segmented image. As can be seen from Fig. 6, the method of the present invention achieves high-quality stereoscopic segmentation results for the stereoscopic image segmentation of occluded regions, color overlapping, and camera motion.

实施本发明实施例,具有如下有益效果:Implementing the embodiment of the present invention has the following beneficial effects:

本发明利用立体图像的极线约束,在能量最优化框架中除了加入 了立体匹配项,同时加入了基于图像补丁的轮廓项,根据这两项合成 了轮廓概率图,随后采用测地线距离计算最优轮廓线。提出的基于一 致性轮廓提取的立体图像分割方法能对立体图像进行一致性的图像 分割,特别是在图像遮挡区域,也能取得很好的立体分割效果,从而 达到既能获得高质量的立体图像分割效果,又有较快的运行效率的目 的。The present invention utilizes the epipolar constraints of the stereo image, and in the energy optimization framework, in addition to the stereo matching item, the contour item based on the image patch is added at the same time, and the contour probability map is synthesized according to these two items, and then the geodesic distance is used to calculate optimal contour. The proposed stereoscopic image segmentation method based on consistent contour extraction can perform consistent image segmentation on stereoscopic images, especially in image occlusion areas, and can also achieve good stereoscopic segmentation results, so as to achieve both high-quality stereoscopic images Segmentation effect, but also the purpose of faster operating efficiency.

本领域普通技术人员可以理解实现上述实施例方法中的全部或 部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以 存储于一计算机可读取存储介质中,所述的存储介质,如ROM/RAM、 磁盘、光盘等。Those of ordinary skill in the art can understand that all or part of the steps in the method of the above-mentioned embodiments can be completed by instructing related hardware through a program, and the program can be stored in a computer-readable storage medium, and the storage Media such as ROM/RAM, magnetic disk, optical disk, etc.

以上所述仅为本发明的较佳实施例而已,并不用以限制本发明, 凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等, 均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present invention shall be included in the protection of the present invention. within range.

Claims (9)

1.一种基于一致性轮廓提取的立体图像分割方法,其特征在于,包括步骤:1. A stereoscopic image segmentation method based on consistent contour extraction, characterized in that, comprising steps: S1、输入立体图像左视图、左视图分割图像及立体图像右视图;S1. Input the left view of the stereo image, the segmented image of the left view and the right view of the stereo image; S2、根据所述立体图像左右视图,计算立体匹配项;S2. Calculate a stereo matching item according to the left and right views of the stereo image; S3、根据所述立体图像左视图、左视图分割图像及立体图像右视图,预估立体图像右视图分割概率图,并根据所述预估的立体图像右视图分割概率图,提取及修正预估轮廓线条,且进一步根据所述修正后的预估轮廓线条,计算基于图像补丁的轮廓项;S3. Estimating the right view segmentation probability map of the stereo image according to the left view of the stereo image, the segmented image of the left view, and the right view of the stereo image, and extracting and correcting the estimate according to the estimated right view segmentation probability map of the stereo image Contour lines, and further calculate contour items based on image patches according to the revised estimated contour lines; S4、根据所述计算出的立体匹配项及基于图像补丁的轮廓项,合成轮廓概率图;S4. Synthesize a contour probability map according to the calculated stereo matching item and the contour item based on the image patch; S5、采用测地线距离计算出所述合成的轮廓概率图的最优轮廓线;S5. Calculate the optimal contour line of the synthesized contour probability map by using the geodesic distance; S6、扫描所述计算出的最优轮廓线,获得立体图像的右视图分割图像。S6. Scan the calculated optimal contour line to obtain a right-view segmented image of the stereo image. 2.如权利要求1所述的立体图像分割方法,其特征在于,所述步骤S2中立体匹配项是通过公式来实现;2. stereoscopic image segmentation method as claimed in claim 1, is characterized in that, in described step S2, stereo matching item is by formula to fulfill; 其中,pi,pj为像素位置索引,pi为左视图分割图像轮廓线上的像素,为以pi为中心的m×m的局部区域,m=5,d(pi)为像素pi的视差值,CAD为绝对误差和算法,且c表示彩色图像的三个通道。Among them, p i , p j are the pixel position index, p i is the pixel on the contour line of the left view segmented image, is the m×m local area centered on p i , m=5, d(p i ) is the parallax value of pixel p i , C AD is the absolute error sum algorithm, and c represents the three channels of the color image. 3.如权利要求2所述的立体图像分割方法,其特征在于,所述步骤S3中“预估立体图像右视图分割概率图”的具体步骤包括:3. stereoscopic image segmentation method as claimed in claim 2, is characterized in that, in described step S3, the specific step of " estimation stereoscopic image right view segmentation probability figure " comprises: 分别构建立体图像左视图和右视图立体图像的特征向量,特征向量的构建如下:式中,为以pi为中心的m×m的局部区域,为以pi为中心的m×m图像补丁,m=5;Construct the eigenvectors of the stereo image left view and the right view stereo image respectively, and the eigenvectors are constructed as follows: In the formula, is a m×m local area centered on p i , is an m×m image patch centered on p i , m=5; 采用预设的搜索策略在立体图像左视图中搜索与立体图像右视图中每个像素点最相似的图像补丁;Search for the image patch most similar to each pixel in the right view of the stereo image in the left view of the stereo image using a preset search strategy; 通过公式IRmask(pi)=ILmask(pj),预估立体图像右视图分割概率图;其中,ILmask为左视图分割图像,pi,pj为像素位置索引,pj为右视图中像素pi在立体图像左视图中找到的最相似图像补丁的中心位置。By the formula I Rmask (p i )=I Lmask (p j ), estimate the right view segmentation probability map of the stereo image; where, I Lmask is the left view segmentation image, p i , p j are the pixel position indexes, and p j is the right In-view pixel p i is the center position of the most similar image patch found in the left view of the stereo image. 4.如权利要求3所述的立体图像分割方法,其特征在于,所述采用预设的搜索策略在立体图像左视图中搜索与立体图像右视图中每个像素点最相似的图像补丁的步骤具体包括:4. stereoscopic image segmentation method as claimed in claim 3, is characterized in that, the step of searching for the image patch most similar to each pixel point in the stereoscopic image right view in the stereoscopic image left view using preset search strategy Specifically include: 在立体图像右视图中的被遮挡像素采用全局范围搜索策略搜索最相似补丁;其中,所述全局范围搜索策略是将搜索范围确定为立体图像左视图整个图像;The occluded pixel in the right view of the stereo image adopts a global range search strategy to search for the most similar patch; wherein, the global range search strategy is to determine the search range as the entire image of the left view of the stereo image; 在立体图像右视图中的非遮挡像素采用扫描线搜索策略搜索最相似补丁;其中,所述扫描线搜索策略是将搜索范围确定为立体图像左视图沿着相同水平线在[0,dmax]视差范围内;dmax为预估立体图像最大视差值;The non-occluded pixels in the right view of the stereo image use a scan line search strategy to search for the most similar patch; wherein, the scan line search strategy is to determine the search range as the left view of the stereo image along the same horizontal line at [0,d max ] parallax Within the range; d max is the maximum parallax value of the estimated stereo image; 其中,立体图像右视图中被遮挡像素和非遮挡像素的区分是通过以下步骤实现,具体包括:采用扫描线搜索策略在立体图像左视图中搜索和立体图像右视图中某一像素Pr最相似补丁的像素,搜索到的像素设为Pl,然后同样采用扫描线搜索策略在立体图像右视图中搜索和立体图像左视图中像素Pl最相似补丁的像素,搜索到的像素设为P′r;如果像素Pr和像素P′r为同一点,则像素Pr为非遮挡像素;反之,如果像素Pr和像素P′r为不同点,则像素Pr为被遮挡像素。Among them, the distinction between the occluded pixel and the non-occluded pixel in the right view of the stereo image is realized through the following steps, which specifically include: using the scan line search strategy to search in the left view of the stereo image and the most similar pixel P r in the right view of the stereo image The pixel of the patch, the searched pixel is set as P l , and then the same scan line search strategy is used to search the pixel in the right view of the stereo image and the pixel P l in the left view of the stereo image is the most similar patch pixel, and the searched pixel is set as P′ r ; if the pixel P r and the pixel P′ r are the same point, then the pixel P r is a non-occluded pixel; otherwise, if the pixel P r and the pixel P′ r are different points, then the pixel P r is an occluded pixel. 5.如权利要求3所述的立体图像分割方法,其特征在于,所述步骤S3中“根据所述预估的立体图像右视图分割概率图,提取及修正预估轮廓线条”的具体步骤包括:5. The stereoscopic image segmentation method according to claim 3, characterized in that, in the step S3, the specific steps of "according to the estimated right-view segmentation probability map of the stereoscopic image, extracting and correcting the estimated contour lines" include : 采用逐行扫描立体图像右视图分割概率图,并待扫描到立体图像右视图分割概率图中某一像素为非零且以该像素为中心的3×3区域内至少包含一个非零像素后,则保留该像素点为轮廓点,且进一步将所保留的轮廓点绘制成预估轮廓线条;Use progressive scanning of the right-view segmentation probability map of the stereoscopic image, and after scanning to a certain pixel in the right-view segmentation probability map of the stereoscopic image is non-zero and contains at least one non-zero pixel in the 3×3 area centered on this pixel, The pixel point is reserved as a contour point, and the retained contour point is further drawn into an estimated contour line; 统计预估轮廓线条中封闭曲线的个数,并按照封闭曲线围成的面积从小到大排列;Count and estimate the number of closed curves in the contour lines, and arrange them according to the area enclosed by the closed curves from small to large; 设定面积阈值Ah=βAk,移除预估轮廓线条中面积小于阈值Ah的封闭曲线,得到修正后的预估轮廓线条;其中,β为权重参数,设定为0.8,Ak为已排完序的第k个封闭曲线的面积,k为左视图分割图像中封闭曲线的个数。Set the area threshold A h = βA k , remove the closed curves whose area is smaller than the threshold A h in the estimated contour lines, and obtain the revised estimated contour lines; where β is the weight parameter, set to 0.8, and A k is The area of the kth closed curve that has been sorted, k is the number of closed curves in the left view segmented image. 6.如权利要求5所述的立体图像分割方法,其特征在于,所述步骤S3中“根据所述修正后的预估轮廓线条,计算基于图像补丁的轮廓项”是通过公式来实现;其中,dc(qi)为像素点qi到预估轮廓线条的最短距离,D为权重参数,设定为D=10,Td为阈值,设定为20;Cr为修正后的预估轮廓线条,d(qi,qj)为两个像素之间的欧拉距离。6. The stereoscopic image segmentation method according to claim 5, characterized in that, "according to the corrected estimated contour line, calculate the contour item based on the image patch" in the step S3 by the formula To achieve; wherein, d c (q i ) is the shortest distance from the pixel point q i to the estimated contour line, D is the weight parameter, set to D=10, T d is the threshold, set to 20; C r is the corrected estimated contour line, d(q i , q j ) is the Euler distance between two pixels. 7.如权利要求6所述的立体图像分割方法,其特征在于,所述步骤S4中“轮廓概率图”是通过公式CC(pi,d(pi))=CS(pi,d(pi))×CO(pi,d(pi))来实现;其中,CS为立体匹配项,CO为基于图像补丁的轮廓项。7. The stereoscopic image segmentation method according to claim 6, wherein the "contour probability map" in the step S4 is obtained by the formula C C (p i , d(p i ))=C S (p i , d(p i ))×C O (p i ,d(p i )); among them, C S is the stereo matching item, and C O is the contour item based on the image patch. 8.如权利要求6所述的立体图像分割方法,其特征在于,所述步骤S5中的“最优轮廓线”是通过求解公式最优解来实现;其中,Cm为左视图分割图像轮廓线,N(pi,pi-1)为轮廓连续性项,当像素pi-1为像素pi的8邻域时,轮廓连续性项的值为0,否则轮廓连续性项的值为∞。8. stereoscopic image segmentation method as claimed in claim 6, is characterized in that, " optimum contour line " in described step S5 is by solving formula The optimal solution is realized; where, C m is the contour line of the left view segmentation image, N(p i , p i-1 ) is the contour continuity item, when the pixel p i-1 is the 8 neighbors of the pixel p i , The value of the contour continuity item is 0, otherwise the value of the contour continuity item is ∞. 9.如权利要求6所述的立体图像分割方法,其特征在于,所述步骤S6具体包括:9. stereoscopic image segmentation method as claimed in claim 6, is characterized in that, described step S6 specifically comprises: 依次逐行从左到右扫描,如果像素值为0,表示该点为背景,直到扫描到像素值为1的点,表示该点为左轮廓点,继续扫描直到扫描到像素值为0的点,表示左轮廓线结束,该点为前景,接着继续扫描直到再次扫描到像素值为1的点且其对应下一个像素点的像素值为0时,表示该点为右轮廓点,Scan from left to right line by line. If the pixel value is 0, it means that the point is the background. Until the point with a pixel value of 1 is scanned, it means that the point is a left contour point. Continue to scan until the point with a pixel value of 0 is scanned. , indicating the end of the left contour line, the point is the foreground, and then continue to scan until a point with a pixel value of 1 is scanned again and the pixel value corresponding to the next pixel point is 0, indicating that the point is a right contour point, 将左轮廓点和右轮廓点之间的像素均填充为1,重复此步骤,直到整个最优轮廓线填充完整,获得立体图像的右视图分割图像。Fill the pixels between the left contour point and the right contour point with 1, repeat this step until the entire optimal contour line is filled completely, and obtain the right view segmentation image of the stereo image.
CN201710820101.7A 2017-09-13 2017-09-13 Stereo image segmentation method based on consistency contour extraction Active CN107578419B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710820101.7A CN107578419B (en) 2017-09-13 2017-09-13 Stereo image segmentation method based on consistency contour extraction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710820101.7A CN107578419B (en) 2017-09-13 2017-09-13 Stereo image segmentation method based on consistency contour extraction

Publications (2)

Publication Number Publication Date
CN107578419A true CN107578419A (en) 2018-01-12
CN107578419B CN107578419B (en) 2020-07-21

Family

ID=61033269

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710820101.7A Active CN107578419B (en) 2017-09-13 2017-09-13 Stereo image segmentation method based on consistency contour extraction

Country Status (1)

Country Link
CN (1) CN107578419B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110298884A (en) * 2019-05-27 2019-10-01 重庆高开清芯科技产业发展有限公司 A kind of position and orientation estimation method suitable for monocular vision camera in dynamic environment
CN110490877A (en) * 2019-07-04 2019-11-22 西安理工大学 Binocular stereo image based on Graph Cuts is to Target Segmentation method
CN112598784A (en) * 2020-12-23 2021-04-02 深兰人工智能芯片研究院(江苏)有限公司 Method and device for adjusting contour scanning result, electronic equipment and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011018328A (en) * 2009-07-08 2011-01-27 Fujitsu Ltd Method and apparatus for dividing three-dimensional model and image processing system including three-dimensional model division apparatus
CN103955945A (en) * 2014-05-23 2014-07-30 哈尔滨师范大学 Self-adaption color image segmentation method based on binocular parallax and movable outline
CN105513064A (en) * 2015-12-03 2016-04-20 浙江万里学院 Image segmentation and adaptive weighting-based stereo matching method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011018328A (en) * 2009-07-08 2011-01-27 Fujitsu Ltd Method and apparatus for dividing three-dimensional model and image processing system including three-dimensional model division apparatus
CN103955945A (en) * 2014-05-23 2014-07-30 哈尔滨师范大学 Self-adaption color image segmentation method based on binocular parallax and movable outline
CN105513064A (en) * 2015-12-03 2016-04-20 浙江万里学院 Image segmentation and adaptive weighting-based stereo matching method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
戴庆焰 等: "基于图割和视域相关性的立体图像分割", 《电信科学》 *
殷均平: "一种改进的基于图像分割的立体匹配算法研究", 《西南师范大学学报(自然科学版)》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110298884A (en) * 2019-05-27 2019-10-01 重庆高开清芯科技产业发展有限公司 A kind of position and orientation estimation method suitable for monocular vision camera in dynamic environment
CN110490877A (en) * 2019-07-04 2019-11-22 西安理工大学 Binocular stereo image based on Graph Cuts is to Target Segmentation method
CN112598784A (en) * 2020-12-23 2021-04-02 深兰人工智能芯片研究院(江苏)有限公司 Method and device for adjusting contour scanning result, electronic equipment and storage medium
CN112598784B (en) * 2020-12-23 2022-04-12 深兰智能科技(上海)有限公司 Method and device for adjusting contour scanning result, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN107578419B (en) 2020-07-21

Similar Documents

Publication Publication Date Title
CN105069808B (en) The video image depth estimation method split based on image
Concha et al. Using superpixels in monocular SLAM
CN102750711B (en) A kind of binocular video depth map calculating method based on Iamge Segmentation and estimation
CN105374019B (en) A kind of more depth map fusion methods and device
CN102903096B (en) Monocular video based object depth extraction method
CN102930530B (en) Stereo matching method of double-viewpoint image
CN104616286B (en) Quick semi-automatic multi views depth restorative procedure
CN104966289B (en) A kind of depth estimation method based on 4D light fields
US9578312B2 (en) Method of integrating binocular stereo video scenes with maintaining time consistency
CN102074020B (en) Method for performing multi-body depth recovery and segmentation on video
CN107862698A (en) Light field foreground segmentation method and device based on K mean cluster
CN102523464A (en) Depth image estimating method of binocular stereo video
CN103561258A (en) Kinect depth video spatio-temporal union restoration method
CN104966286A (en) 3D video saliency detection method
CN107507206B (en) A depth map extraction method based on saliency detection
CN103996174A (en) Method for performing hole repair on Kinect depth images
CN102368826A (en) Real time adaptive generation method from double-viewpoint video to multi-viewpoint video
CN104091336B (en) Stereoscopic image synchronous segmentation method based on dense disparity map
CN106060509B (en) Introduce the free view-point image combining method of color correction
Zhang et al. Simultaneous multi-body stereo and segmentation
CN104331890B (en) A kind of global disparity method of estimation and system
WO2013144418A1 (en) Image segmentation
CN107578419B (en) Stereo image segmentation method based on consistency contour extraction
CN102750694A (en) Local optimum belief propagation algorithm-based binocular video depth map solution method
CN104661014B (en) The gap filling method that space-time combines

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
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20180112

Assignee: Big data and Information Technology Research Institute of Wenzhou University

Assignor: Wenzhou University

Contract record no.: X2020330000098

Denomination of invention: A stereo image segmentation method based on consistent contour extraction

Granted publication date: 20200721

License type: Common License

Record date: 20201115

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20180112

Assignee: Wenzhou Yuanbian Digital Creativity Co.,Ltd.

Assignor: Wenzhou University

Contract record no.: X2022330000729

Denomination of invention: A stereo image segmentation method based on consistent contour extraction

Granted publication date: 20200721

License type: Common License

Record date: 20221206