System and method for splicing images

Patent metadata

Field	Value
Publication number	US-2020202530-A1
Application number	US-202016806207-A
Country	US
Kind code	A1
Filing date	Mar 2, 2020
Priority date	Jul 14, 2016
Publication date	Jun 25, 2020
Grant date	—

Abstract

The present disclosure relates to systems and methods for image splicing. The systems and methods may acquire a first image and a second image, determine a plurality of first feature points in a first region of the first image, determine a plurality of second feature points in a second region of the second image, then match the plurality of first feature points with the plurality of second feature points to generate a plurality of point pairs. Based on the plurality of point pairs, a third region on the first image and a fourth region on the second image may be determined. Finally, a third image may be generated based on the first image and the second image, wherein the third region of the first image may overlap with the fourth region of the second image in the third image.

First claim

1 - 20 . (canceled)

21 . A method for determining a plurality of feature points of an image implemented on a computing device including at least one processor and a storage, the method comprising:
decomposing the image; generating a difference image corresponding to the decomposed image; generating a plurality of extreme points based on the difference image and grayscale values thereof; filtering out at least one extreme point on a sharp boundary of the image from the plurality of extreme points; and determining the plurality of feature points based on the filtered extreme points.

22 . The method of claim 21 , wherein the decomposing the image includes:
decomposing the image based on Gaussian pyramid decomposition.

23 . The method of claim 21 , wherein the generating a difference image corresponding to the decomposed image includes:
executing a convolution between the image and a Gaussian function to generate a plurality of smoothing images; down-sampling the plurality of smoothing images to generate the decomposed image in different scales; determining at least one difference among the plurality of down-sampled smoothing images included in the decomposed image; and generating the difference image based on the at least one difference.

24 . The method of claim 23 , the plurality of smoothing images being arranged in layers, wherein the determining at least one difference among the plurality of smoothing images includes:
determining the at least one difference between two adjacent layers of the plurality of down-sampled smoothing images.

25 . The method of claim 21 , each of the plurality of extreme points includes a pixel point that has a largest absolute value of a grayscale in a region with a preset size in the image.

26 . The method of claim 25 , the difference image comprising a plurality of layers, wherein the generating a plurality of extreme points based on the difference image and grayscale values thereof includes:
in each of the plurality of layers of the difference image, comparing absolute values of grayscales of pixel points with each other in a region with the preset size in the layer; and determining a preliminary pixel point whose absolute value of a grayscale is larger than absolute values of the grayscales of other pixel points in the region as one of the plurality of extreme points.

27 . The method of claim 26 , further comprising:
modifying positions of the pixel points in the region with the preset size in the layer; and determining the plurality of extreme points based on the pixel points with the modified positions.

28 . The method of claim 27 , wherein the modifying positions of the pixel points in the region with the preset size in the layer includes:
determining a displacement of a pixel point among the pixel points according to formula:

δ
=

-

∂

D
T

∂
X



(

∂
2


D

∂

X
2

)

-
1

,

where D denotes the difference image, X=(Δx, Δy, Δσ) T , Δx denotes a difference between horizontal coordinate x of the pixel point and horizontal coordinate x p of the modified pixel point, Δy denotes a difference between vertical coordinate y of the pixel point and vertical coordinate y p of the modified pixel point, and Δσ denotes a difference between σ and σ p ; and
modifying the position of the pixel point based on the displacement.

29 . The method of claim 26 , wherein the generating a plurality of extreme points based on the difference image and grayscale values thereof includes:
sorting the preliminary pixel points in a descending order according to the absolute values of grayscales of the preliminary pixel points; determining N pixel points with the highest absolute values of grayscales among the preliminary pixel points; modifying positions of the N pixel points; and determining the N pixel points with the modified positions as at least a part of the plurality of extreme points.

30 . The method of claim 21 , wherein the filtering out at least one extreme point on a sharp boundary of the image from the plurality of extreme points includes:
calculating a largest principal curvature and a smallest principal curvature of each of the extreme points based on a change of surrounding pixel points; and determining whether the extreme point is on the sharp boundary of the image according to a ratio of the largest principal curvature to the smallest principal curvature; if the ratio of the largest principal curvature to the smallest principal curvature is larger than a preset threshold, determining that the extreme point is on the sharp boundary of the image; and filtering out the extreme point that is determined on the sharp boundary of the image.

31 . The method of claim 30 , wherein the calculating a largest principal curvature and a smallest principal curvature of each of the extreme points based on a change of surrounding pixel points includes:
determining the largest principal curvature and a smallest principal curvature f each of the extreme points according to formula:

H
=

[

D
xx

D
xy

D
xy

D
yy

]

,

where D denotes a Gaussian difference function, D xx =D (x+1,y,σ)+D(x−1,y,σ)−2×D(x,y,σ), D yy =D(x,y+1,σ)+D(x,y−1,σ)−2×D(x,y,σ), D xy =0.25×[D+1,y+1,σ)+D(x−1,y−1,σ)−D(x−1,y+1,σ)−D(x+1,y−1,σ], and (x,y) denotes a coordinate of the extreme point.

32 . The method of claim 30 , wherein the preset threshold varies between 1 and 30.

Assignees

• Shanghai United Imaging Healthcare Co Ltd

Inventors

• Hu Yang
• Niu Jie
• Xu Liang
• Han Yecheng
• Wang Hanyu
• Bian Yueyan
• Teng Wanli
• Cui Kai
• Yang Le

Classifications

• G06T7/33 Primary

  using feature-based methods · CPC title

• G06T2207/10081

  Computed x-ray tomography [CT] · CPC title

• G06T7/174

  involving the use of two or more images · CPC title

• G06T7/0014 Primary

  using an image reference approach · CPC title

• G06T7/11

  Region-based segmentation · CPC title