CN101197043B - Shape-maintaining fitting algorithm in image processing - Google Patents
Shape-maintaining fitting algorithm in image processing Download PDFInfo
- Publication number
- CN101197043B CN101197043B CN 200710115786 CN200710115786A CN101197043B CN 101197043 B CN101197043 B CN 101197043B CN 200710115786 CN200710115786 CN 200710115786 CN 200710115786 A CN200710115786 A CN 200710115786A CN 101197043 B CN101197043 B CN 101197043B
- Authority
- CN
- China
- Prior art keywords
- curved surface
- image
- interpolation
- derivative
- value
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
- Editing Of Facsimile Originals (AREA)
Abstract
The invention provides a preserving type fitting algorithm in the image processing. The method firstly reverses the approximations of special points on an original curved surface according to a discrete image, and constructs a double thrice Coons interpolation curved surface on an area enclosed by four neighboring data points for the new data, wherein, a first step partial derivatives of the four corner points are obtained by estimating with a preserving type interpolation method, a second step partial derivatives of the four corner points are obtained by averaging the partial derivatives of neighboring pixel points of the four corner points. The re-sampling is done to the Coons curved surface on the curved surface according to the scaling, and we adopt a zone sampling to obtain the grey values of the re-sampling points. More pixels are obtained by increasing the sampling density, thereby achieving the magnification of a image. The new method can effectively overcome the mosaic phenomenon and preserve clearer image details.
Description
(1) technical field
The present invention relates to the image interpolation technology, belong to image processing field.
(2) background technology
Along with the development of infotech, digital picture occurs in social life in the increasing field, and having become computing machine needs extremely important class data to be processed.And the basic operation of processing as image, image amplifies the background that has a wide range of applications, most of with digital picture as the demand that all can have image to amplify in the application of data.In magic magiscan, the CT image is amplified, can bring great convenience for diagnosis.Good image magnification method can make the image that obtains have higher quality, thereby provides very large facility for further working.
It is exactly the pixel that makes new advances according to Pixel Information interpolation discrete in the original image that image amplifies.The structure interpolation curved surface has several different methods, and such as the arest neighbors interpolation, bilinear interpolation is based on the interpolation of B é zier curved surface etc.The simplest method is the zeroth order interpolation, namely the isolated ground of each pixel is amplified in proportion, and the image that produces like this has serious mosaic.The interpolation curved surface of arest neighbors interpolation is surface of discontinuity, and enlarged image has obvious mosaic; It is continuous that the bilinear interpolation curved surface can reach CO, is continuous in the pixel region inside of original image, but discontinuous on the border, transition is level and smooth not between pixel.For solving mosaic phenomenon and border continuous problem, usually for the continuous mathematical model of original image structure, according to the scaling ratio model is resampled again, to obtain the image of corresponding size.The curved surface that obtains such as the interpolation based on B é zier curved surface is that C1 is continuous in whole image space, and the enlarged image that therefore obtains is level and smooth on the whole, reasonable effect is arranged, but the outline line in the image also can be owing to smooth effect thickens.Also have a lot of similarly methods of high-order interpolation, these methods can overcome mosaic phenomenon effectively, but because the smooth effect of high-order interpolation, so that the image detail loss is serious, the objects in images soft edge is unclear.For this reason, the interpolation algorithm based on fitted dividing curve is suggested.This method is at first carried out match to the outline line in the image, for the pixel region that has outline line to pass through, for outline line the zones of different that this pixel segmentation becomes is carried out respectively interpolation, to keep the clear of outline line.It is continuous that but the boundary except outline line of the pixel region that the interpolation curved surface of the method structure passes through at outline line can not reach C1.Based on the inhomogeneous interpolation image amplification method of cubic spline, the burst continuation algorithm equal segments interpolation algorithm of image scaling is suggested in succession, also is used for the processing of outline line in the optimized image.But these piecewise interpolations have also destroyed whole interpolation curved surface in the continuity of pixel adjoiner guaranteeing that outline line place curved surface is discontinuous with when giving prominence to profile.
(3) summary of the invention
Mosaic or the serious phenomenon of loss of detail occur easily when utilizing existing method to carry out zooming, how discussion of the present invention effectively overcomes mosaic phenomenon, makes the image behind the scaling keep clearly image detail.Discrete picture is done scaling operation, and one of effective method is to obtain the former curved surface that generates discrete picture, and former curved surface is directly done the scaling operation.Structure generates the former curved surface of discrete picture, need to know the value of every bit on the former curved surface, and each known discrete picture is an area sampling value, therefore can according to the counter approximate value of asking every bit on the former curved surface of area sampling value, then construct the Proximal surface of former curved surface according to approximate value.Basic thought of the present invention can be summarized as follows.At first according to the counter approximate value of asking particular point on the former curved surface of discrete picture, the approximate value of particular point forms a quadrilateral mesh, then construct bicubic Coons interpolation curved surface on the quadrilateral mesh, all bicubic Coons interpolation curved surface amalgamations form the Proximal surface of former curved surface together.At the zone structure bicubic Coons interpolation curved surface that adjacent four pixels surround, wherein, the single order partial derivative of four angle points and second order local derviation are estimated to obtain by the shape preserving interpolation method.According to the scaling ratio Coons curved surface is resampled on this curved surface, we adopt area sampling to obtain the gray-scale value of resample points.By increasing sampling density, obtain more pixel, thereby reach the purpose that image amplifies.
The present invention can be divided into three parts, A, according to the counter approximate value of asking former Surface Data point of discrete picture; The structure of B, Coons curved surface; C, resample according to the scaling ratio.The step of then utilizing the shape preserving interpolation Coons surface to carry out the image amplification is described below:
Step1: according to the counter approximate value of asking particular point on the former curved surface of discrete picture;
Step2: according to the value of the data point of obtaining obtain data point at u to the derivative that makes progress with v;
Step3: the second order local derviation of obtaining the data point place;
Step4: the bicubic interpolation patch calculates fast.
Step5: according to the scaling ratio Coons patch is resampled, obtain more accurate pixel;
(4) description of drawings
Fig. 1, dull non-protruding situation is for P
I-2, j, P
I-1, j, P
I, jAnd P
I+1, jThe broken line of four some formation is dull non-protruding;
Fig. 2, the situation of convex polygon is for P
I-2, j, P
I-1, j, P
I, jAnd P
I+1, jFour points couple together and form a plane convex polygon;
Fig. 3, each picture element differentiate and bicubic patch organigram, the present invention be its each by adjacent four pixels, such as P
I, j, P
I+1, j, P
I, j+1, P
I+1, j+1Structure bicubic Coons interpolation curved surface is put P with adjacent image point on the zone that surrounds
I-1, j, P
I+1, jAnd P
I+1, jThree point estimation P
I, jThe u at place is to derivative (P
I, j')
u, use pixel P
I, j-1, P
I, jAnd P
I, j+1Three point estimation P
I, jThe v at place is to derivative (P
I, j')
v
Fig. 4, the area sampling synoptic diagram, wherein, the pixel of solid circle expression original image, the point that hollow circle expression resamples.For resample points 1, corresponding sample area is A, 2,3 points that resample, and its sample area is respectively B, C, and the value that resamples is determined jointly by its two adjacent Coons curved surfaces.
To be two kinds of methods amplify 2 times comparison to the CT image for Fig. 5, Fig. 6, and Fig. 5 is that bilinear interpolation method obtains, and Fig. 6 is obtained by method of the present invention.The quality of Fig. 5 and Fig. 6 two width of cloth images does not visually have significant difference.For comparing mosaic phenomenon, we amplify high power to image, see Fig. 7, Fig. 8.
To be two kinds of methods amplify 5 times comparison to the CT image for Fig. 7, Fig. 8, and Fig. 7 is that bilinear interpolation method obtains, and Fig. 8 is obtained by method of the present invention.On display, the mosaic phenomenon of Fig. 7 is very obvious, and Fig. 8 does not just have mosaic phenomenon.
(5) embodiment
The below at first illustrates the principles of science of the important step institute foundation in the above-mentioned algorithm:
The two-dimentional match that A, based on data excavate
1) mapping relations of area sampling and point sampling
Pixel value on the CT image is the area sampling value, and supposing has n * n pixel on one deck CT image, and then these pixels are the area sampling values that obtain at n * n unit area from space curved surface.Show that for example, during the scaling conversion, ideal situation is to obtain the curved surface that approaches of luv space curved surface carrying out image, then carry out resampling approaching curved surface, obtain to approach curved surface in the area sampling value of each new region.Therefore, if can obtain according to the area sampling value of CT image the more accurate point sampling value of this regional center point, just can according to these point sampling values construct precision higher approach curved surface.So, need to carry out pre-service to the area sampling value of CT image, obtain more accurate point sampling value.We adopt data mining technology to be asked by given pixel and approach curved surface.
Set up now the mapping relations of area sampling and point sampling, the situation of one dimension be discussed first:
If f (x) was n some P
i(i=0,1,2 ..., curve n-1),
(i=0,1,2 ..., n-1) be from this curve obtain corresponding to the zone [x
i, x
I+1] n area sampling value, then have following relation to set up:
By (1) formula as can be known, by setting up interpolation point P
iWith the area sampling value
(i=0,1,2 ..., relation n-1), can by
Obtain P
j(i=0,1,2 ..., n-1), namely obtain more accurate point sampling value by the area sampling value.
The situation of two dimension is discussed again:
If P (x, y) represents a curved surface, P
Ij(i, j=0,1,2 ..., n-1) be n * n point on this curved surface, obviously, P (x, y) can put approximate representation by this n * n.Again
(i, j=0,1,2 .., n-1) be from curved surface P (x, y) obtain corresponding to regional A
IjN * n area sampling value, then have following relation to set up:
Wherein, S
IjRepresent regional A
IjArea.
By (2) formula as can be known, by setting up some P
IjWith the area sampling value
(i, j=0,1,2 ..., relation n-1), can by
Obtain P
Ij(i, j=0,1,2 ..., n-1), namely obtained the approximate value of point sampling value by the area sampling value.
By the mapping relations of above-mentioned zone sampling and point sampling as can be known, by curve construction/curved surface, can be obtained by the area sampling value of known CT image comparatively accurately point sampling value.
2) calculating of point sampling value
This algorithm utilizes the mapping relations of above-mentioned area sampling and point sampling, on the basis of CT area sampling data, obtains the point sampling value of more accurate CT pixel, adopts secondary Lagrange's interpolation curved surface to approach original curved surface in the mapping process.In addition, the interior pixels point of CT image has been taked different disposal routes with the boundary pixel point, made the sampled value of the pixel that obtains more reasonable and accurate.
If
(i, j=0,1,2 ..., n-1), be n * n pixel value on the tomographic image, A
Ij(i, j=0,1,2 ..., n-1) expression is with (x
j, y
i) centered by unit square, we wish to try to achieve point (x
j, y
i) the more accurate sampled value P that locates
Ij(i, j=0,1,2 ..., n-1).
Calculation level (x
1, y
1) the sampled value P that locates
11Algorithm idea as follows:
1, by summit P
Ij(i, j=0,1,2) obtains a secondary lagrangian fit curved surface P (x, y);
2, according to the mapping relations of (2) formula, with P (x, y) respectively at A
IjCarry out integration on (i, j=0,1,2) zone, integrated value equals the known pixel value on this zone
(i, j=0,1,2) multiply by regional A
IjThe area of (i, j=0,1,2);
3, obtain nine yuan of linear function groups by the 2nd step, find the solution and to get P
IjThe value of (i, j=0,1,2).We only keep P
11Value because it is that relative accuracy is the highest in nine new values, the value of other points can adopt same process to obtain.
Specifically being calculated as follows of its each step:
Step 1: utilize summit P
IjSecondary lagrangian fit curved surface of (i, j=0,1,2) structure, surface equation is:
Step 2: the lagrangian fit curved surface P (x, y) that (3) formula is represented is at regional A
IjCarry out integration on (i, j=0,1,2), and integrated value equals
(i, j=0,1,2) and A
IjThe product of (i, j=0,1,2) area, again A
Ij(i, j=0,1,2) is unit square, so integration satisfies following relational expression:
Step 3: can obtain one with P by (3) and (4)
Ij(i, j=0,1,2) is nine yuan of linear function groups of unknown quantity, separates this system of equations and can obtain P
IjThe value of (i, j=0,1,2).Wherein, we only keep P
11Value.
In the above-mentioned steps, the integral process of (3) formula is actually Lagrangian basis function is carried out integration, is easy to obtain a general expression formula, calculates simple, convenient.In addition, the matrix of coefficients of the system of equations in the step 3 all is identical in all cases, therefore, can obtain first the inverse matrix of matrix of coefficients, then this inverse matrix is used for the solving equation group, thereby has improved the travelling speed of system, has reduced the complexity of system.Inverse matrix is as follows:
The structure of B, Coons curved surface
1) utilize the shape preserving interpolation algorithm carry out each picture element at u to the derivative that makes progress with v.Put P with adjacent image point
I-1, j, P
I, jAnd P
I+1, jThree point estimation P
I, jThe place u to derivative (P '
I, j)
u, use pixel P
I, j-1, P
I, jAnd P
I, j+1Three point estimation P
I, jThe place v to derivative (P '
I, j)
vThe below with (P '
I, j)
uBe example, carry out finding the solution of each picture element derivative.Along u to getting 3 P
I-1, j, P
I, jAnd P
I+1, jThe structure interpolation curve, owing to the step-length such as be between 3, the curve of establishing three point interpolation is P (t), then its first order derivative is calculated as follows:
(P′
i,j)
u=(P
i+1,j-P
i-1,j)/2 (5)
To border data point P
1, jAnd P
N, j, utilize following formula calculate (P '
1, j)
u=2 (P
2, j-P
1, j)-(P '
2, j)
u, (P '
N, j)
u=2 (P
N, j-P
N-1, j)-(P '
N-1, j)
u
Directly estimate that by (5) formula three Hermite curves of derivative structure do not have the shape of data point suggestion sometimes, this is unacceptable to clinical practice.Therefore need to apply to the derivative that (5) formula is estimated certain constraint, make the cubic curve of structure have the shape of data point suggestion.Theoretical analysis is known, constructs three Hermite curves by the derivative of (5) formula definition, if curve has identical monotonicity and convexity with data point, then (P '
I, j)
iNeed satisfy following condition
Wherein, Δ P
I-1, j=P
I, j-P
I-1, j, Δ P
I, j=P
I+1, j-P
I, j
(6) characteristics of the derivative of formula definition be the maximal value of three Hermite curves of constructing less than the maximal value of given data point, minimum value is greater than the minimum value of given data point.In order to improve interpolation precision and to obtain more rational curve shape, we do reduction to the constraint of (6) formula and process, and estimate derivative by (5) formula first, then it are according to circumstances adjusted.Selection mode is as follows:
If 1. for P
I-2, j, P
I-1, j, P
I, jAnd P
I+1, jThe broken line of four some formation is dull non-protruding situations, and we can directly select the derivative of (6) definition to estimate, just can satisfy guarantor's type requirement of Coons interpolation.
If 2. for P
I-2, j, P
I-1, j, P
I, jAnd P
I+1, jFour points couple together and form a plane convex polygon, and we need to do reduction to the constraint of (6) formula and process, and process as follows:
To P
I-1, j, P
I, jValue and (P '
I, j)
uCarrying out interpolation, to get curve as follows:
P=a
1u
2+b
1u+c
1 (7)
Wherein
To P
I-1, j, P
I, jValue and (P '
I-1, j)
uCarrying out interpolation, to get curve as follows:
Q=a
2u
2+b
2u+c
2 (8)
By point (P
I-2, j, u
I-2, j) and (P
I-1, j, u
I-1, j) straight line that consists of is:
P1=(P
i-1,j-P
i-2,j)(u-u
i-1,j)+P
i-1,j=A
1u+B
1 (9)
By (P
I, j, u
I, j) and (P
I+1, j, u
I+1, j) straight line that consists of is:
P2=(P
i+1,j-P
i,j)(u-u
i,j)+P
i,j=A
2u+B
2 (10)
Remember that it is (u that two straight lines intersect intersection point
Imid, j, P
Imid, j), following formula is found the solution
P-P1=a
1u
2+b
1u+c
1-(P
i,j-P
i-1,j)(u-u
i-1,j)-P
i-1,j=0 (11)
Then formula (15) has a solution at least, and namely interpolation curve and straight line P1 intersect at the end points place at least, and the end points solution is (P
I-1, j, u
I-1, j).And if have non-end points solution u
p, and satisfy u
I-1, j<u
p<u
Imid, j, prove that then interpolation curve P is not at a P
I-1, j, P
Imid, j, P
I, jIn the convex closure that surrounds, this just need to adjust (P '
I, j)
uAdjust (P '
I, j)
uAfter, as long as make interpolation curve and straight line P
1Tangent, just can make the interpolation curve P of acquisition at three some P
I-1, j, P
Imid, jAnd P
I, jConvex closure within, namely satisfy following condition:
The solving equation group can get after (12) (P '
I, j)
u, then (P '
I, j)
uCan satisfy the requirement of guarantor's type.
In like manner, following formula is found the solution:
Q-P2=a
2u
2+b
2u+c
2-(P
i+1,j-P
i,j)(u-u
i,j)-P
i,j=0 (13)
Then formula (13) is except end points solution (P
I, j, u
I, j), may there be non-end points solution u
qIf satisfy u
Imid, j<u
q<u
I, j, prove that then interpolation curve Q is not at a P
I-1, j, P
Imid, j, P
I, jIn the convex closure that surrounds, this just need to adjust (P '
I-1, j)
uAdjust (P '
I-1, j)
uAfter, as long as so that curve Q and straight line P2 tangent, Q is at a P
I-1, j, P
Imid, j, P
I, jIn the convex closure that surrounds, namely satisfy following condition:
The solving equation group can get after (14) (P '
I-1, j)
u, then (P '
I-1, j)
uCan satisfy the requirement of guarantor's type.
For (P '
I, j)
vFind the solution and adjustment and above-mentioned (P '
I, j)
uFind the solution and adjust similar, do not repeat them here.
2) the second order local derviation is estimated
The calculating of second order local derviation is usually very complicated, and in the present invention, we have selected a kind of fairly simple method of estimation.By finding the solution u to the single order local derviation, find the solution v to the single order local derviation, get it and on average then get P
I, jThe second order local derviation at place is calculated as follows:
(P″
i,j)
uv=(P
i+1,j-P
i-1,j+P
i,j+1+P
i,j-1)/4
3) structure P
I, j, P
I+1, j, P
I, j+1, P
I+1, j+1The bicubic Coons patch P that surrounds
Ij(u, v).
Bicubic Coons patch is the curved surface of being vowed, being led resultant second order local derviation information definition by the point of adjacent 4 corner points, is used for carrying out interpolation between one group of data point.For image, the Coons curved surface is applied to the interpolation of pixel gray-scale value.Get the image of a width of cloth m * n, for its each by adjacent four pixels, such as P
I, j, P
I+1, j, P
I, j+1, P
I+1, j+1Structure bicubic Coons interpolation curved surface is designated as P on the zone that surrounds
Ij(u, v).
Above-mentioned four parameter values corresponding to angle point are respectively (0,0), (1,0), (0,1), (1,1), then to above-mentioned four summit pixel value P (0,0), P (1,0), P (0,1), P (1,1), and four summits cut resultant second order local derviation totally 16 interpolation information carry out the bicubic Coons curved surface P (u, v) of interpolation, u, v ∈ [0,1] * [0,1] is defined as follows:
Wherein
P is the gray-scale value of given picture element, P
u, P
vRepresent that respectively u is to vowing P with guide v
UvThen represent the second order local derviation.In the practical application, only can obtain the gray-scale value of each pixel of piece image, by 1) and 2) part as can be known, utilize the method for shape preserving interpolation algorithm and simple method of estimation second order local derviation can obtain each picture element at u to the derivative and the second order local derviation that make progress with v, substitution formula (1) bicubic Coons patch P
Ij(u, v).
C, resample according to the scaling ratio
1) according to the scaling ratio Coons patch is resampled
If image needs amplifieroperation, we need to resample to original image according to the scaling ratio.The bicubic Coons interpolation curved surface of above-mentioned structure is applied to the gradation of image interpolation, and new sampled point is taken from this bicubic Coons interpolation curved surface sheet.
Suppose that the image enlargement factor is s, then original image I (u, v) is at u, and the v direction resamples according to the interval of l/s, namely obtains amplifying the image after s times.Digital picture I (u, the v) interpolation that is about to m * n is the image I of m ' * n ' ' (u ', v ').Wherein
m′=m*s;
n′=n*s;
Then resample points I ' (u ', v ') corresponding to the position of the pixel in the original image is:
u=u′/s;
v=v′/s;
The label of this pixel bicubic Coons patch at place in original image is i=int (u), j=int (v), and then according to this interpolation curved surface sheet (1) of hole of structure, the gray-scale value of resample points I ' (u ', v ') is P
Ij(u-i, v-j).
By the process that obtains the CT data as can be known, the CT data value of each pixel in fact is the mean value of the density of its place voxel, is an area sampling value, rather than this an accurate point sampling value.So during resampling, we also adopt area sampling to obtain the gray-scale value of resample points.Here we also adopt area sampling to obtain the gray-scale value of resample points.With the gray-scale value of the mean value on certain zonule, resample points place as resample points:
Wherein A is resample area.If on the border of former Coons patch, then the value of its area sampling should be determined by its two adjacent Coons curved surfaces resample points jointly just.That is:
By choosing different sampling densities, calculate the color value of sampled point by above-mentioned formula, just can obtain the image of any scaling ratio.
2) the quick calculating of image zooming
If integral domain is (x1, x2) ﹠amp; (y1, y2) then according to the expression-form (15) of Coons curved surface, finds the solution (16), gets (17) formula
I=(IntegralF
0(x2)-IntegralF
0(x1))*(P
i,j*(IntegralF
0(y2)-IntegralF
0(y1))+P
i,j+1*(IntegralF1(y2)-IntegralF1(y1))+(P′
i,j)
v*(IntegralG0(y2)-IntegralG0(y1))+(P′
i,j+1)
v*(IntegralG1(y2)-IntegralG1(y1)))+(IntegralF1(x2)-IntegralF1(x1))*(P
i+1,j*(IntegralF0(y2)-IntegralF0(y1))+P
i+1,j+1*(IntegralF1(y2)-IntegralF1(y1))+(P′
i+1,j)
v*(IntegralG0(y2)-IntegralG0(y1))+(P′
i+1,j+1)
v*(IntegralG1(y2)-IntegralG1(y1)))+(IntegralG0(x2)-IntegralG0(x1))*((P′
i,j)
u*(IntegralF0(y2)-IntegralF0(y1))+(P′
i,j+1)
u*(IntegraIF1(y2)-IntegralF1(y1))+(P″
i,j)
uv*(IntegralG0(y2)-IntegralG0(y1))+P″
i,j+1)
uv*(IntegralG1(y2)-IntegralG1(y1)))+(IntegralG1(x2)-IntegralG1(x1))*((P′
i+1,i)
u*(IntegralF0(y2)-IntegralF0(y1))+(P′
i+1,j+1)
u*(IntegralF1(y2)-IntegralF1(y1)+(P″
i+1,j)
uv*(IntegralG0(y2)-IntegralG0(y1))+(P″
i+1,j+1)
uv*(IntegralG1(y2)-IntegralG1(y1)))(17)
Wherein, Integral represents the integration to the function of its back, for F
0(), F
1(), G
0(), G
1The integrated form of (), IntegralF
0(t)=t
4/ 2-t
3+ t, IntegralF1 (t)=-t
4/ 2+t
3, IntegralG0 (t)=t
4/ 4-2t
3/ 3+t
2/ 2, IntegralG1 (t)=t
4/ 4-t
3/ 3
Above-mentioned expression formula is 4 order polynomial forms, and the process of each integration all will be calculated its 4 order polynomial form, and its calculated amount is too large.But according to the form that resamples as can be known, as long as enlargement factor determines, being assumed to be s, is the same for situation about resampling on all Coons patchs that are made of adjacent four pixels.That is to say, for the resample points of parameter identical (being assumed to be (uu, vv)) on the different Coons patchs, its integral domain is identical, and namely u is limited to (uu-1/2s, uu+1/2s) up and down, v is limited to (vv-1/2s, vv+1/2s) up and down.In conjunction with the form of (17) formula, if to the F in (15) formula
0(t), F
1(t), G
0(t), G
1(t) at (uu-1/2s, uu+1/2s) ﹠amp; Integration on (vv-1/2s, vv+1/2s) calculates in advance and stores, then in the process of the area sampling that calculates each resample points, and just need not be at every turn all to F
0(), F
1(), G
0(), G
1() recomputates its integrated form, thereby reduces its time complexity, accelerates computation process.For example image amplifies s doubly, for F
0(), F
1(), G
0(), G
1The integration of () is in the interval
Integration on (s is the scaling multiple)
Calculate in advance and store,
Then can search for corresponding integrated value in record when (17) formula of calculating, (17) formula of bringing into gets final product to get integral result, and experimental result shows that its computing velocity can improve 40%.
We are applied in the present invention in the CT image processing system.The CT image is amplified, can bring great convenience for diagnosis.Its realization can realize by software programming.The false code of given first the present invention programming.
In the CT image processing system, if certain width of cloth CT image is amplified, the multiple that we will need enlarged image and needs to amplify passes to the ZoomBasedCoons method.From the process of amplifying, can see that the present invention selects shape preserving interpolation to calculate the single order local derviation and the second order local derviation is constructed Coons interpolation curved surface sheet, naturally guaranteed continuity and the shape-retaining ability of the image after the amplification, also just guaranteed the flatness of image.According to the ratio of amplifying, resample at the Coons interpolation curved surface sheet of structure, with the mean value of the picture element on the zonule, the resample points place gray-scale value as resample points.Bicubic Coons interpolation curved surface sheet reaches the secondary precision, so can obtain more image detail.So this invention can make the image after the amplification that obtains have higher quality, can effectively overcome mosaic phenomenon, make the image of generation that more clearly image detail be arranged.
As an example, need in the width of cloth diagnosis process to select the image used, respectively it is amplified 2 times and 5 times.Design sketch is seen accompanying drawing Fig. 7, Fig. 8.Can obtain reasonable effect from illustrating as can be known the present invention.
Claims (6)
1. the shape-maintaining fitting method during an image is processed is characterized in that:
All pixels are taken from an original curved surface on the original image, according to the counter approximate gray-scale value of asking its corresponding point on former curved surface of each pixel gray-scale value on the original image, form a series of new datas,
To new data, at the zone structure bicubic Coons interpolation curved surface that its adjacent four data points surround, wherein,
The single order partial derivative of four angle points is estimated to obtain by the shape preserving interpolation method,
The second order local derviation of four angle points is averaged by the local derviation of its adjacent image point point and obtains,
On this bicubic Coons curved surface, resample according to the scaling ratio,
Adopt area sampling to obtain the gray-scale value of resample points, by increasing sampling density, obtain more pixel,
Wherein, the single order local derviation of four angle points has the secondary approximation accuracy, and method for solving is as follows:
Put P with adjacent image point
I-1, j, P
I, jAnd P
I+1, jThree point estimation P
I, jThe place u to derivative (P '
I, j)
u, use pixel P
I, j-1, P
I, jAnd P
I, j+1Three point estimation P
I, jThe place v to derivative (P '
I, j)
v, to (P '
I, j)
u, along u to getting 3 P
I-1, j, P
I, jAnd P
I+1, jThe structure interpolation curve, owing to the step-length such as be between 3, the curve of establishing three point interpolation is P (t), then its first order derivative is calculated as follows:
(P′
i,j)
u=(P
i+1,j-P
i-1,j)/2。
2. the shape-maintaining fitting method during image according to claim 1 is processed, it is characterized in that: the step that forms new data comprises, utilize the mapping relations of area sampling and point sampling, adopt secondary Lagrange slot curved surface to approach the point sampling value that original curved surface obtains more accurate pixel, wherein secondary Lagrange's interpolation surface equation is:
The Lagrange's interpolation curved surface P (x, y) that (1) formula is represented is at regional A
IjCarry out integration on (i, j=0,1,2), and integrated value equals
With A
IjThe product of (i, j=0,1,2) area:
Can obtain one with P by (1) formula and (2) formula
Ij(i, j=0,1,2) is nine yuan of linear function groups of unknown quantity, separates this system of equations and can obtain P
IjThe value of (i, j=0,1,2) wherein only keeps P
IjValue,
Integral process is actually Lagrangian basis function is carried out integration, and the inverse matrix of the matrix of coefficients of system of equations is as follows:
This inverse matrix is used for the solving equation group.
3. the shape-maintaining fitting method during image according to claim 1 is processed is characterized in that: in order to guarantee the shape-retaining ability of interpolation curved surface, adopt following method to P
I, jThe place u to derivative (P '
I, j)
uLimit,
If for P
I-2, j, P
I-1, j, P
I, jAnd P
I+1, jThe broken line of four some formation is dull non-protruding situations, to (P '
I, j)
uLimit as follows:
If for P
I-2, j, P
I-1, j, P
I, jAnd P
I+1, jFour points couple together and form a plane convex polygon, need to do reduction to the constraint of (3) formula and process, namely to (P '
I, j)
uAdjust, so that by P
I-1, j, P
I, jValue and (P '
I, j)
uCarry out curve and the straight line P of interpolation
I-2, jP
I-1, jTangent, just can make the interpolation curve of acquisition at three some P
I-1, j, P
Imid, jAnd P
I, jConvex closure within, then (P '
I, j)
uCan satisfy the requirement of guarantor's type, wherein, Δ u
i=u
I+1, j-u
I, j,
Δ u
I-1=u
I, j-u
I-1, j,
P
Imid, jBe straight line P
I-2, jP
I-1, jWith straight line P
I, jP
I+1, jIntersection point.
4. the shape-maintaining fitting method during image according to claim 1 is processed is characterized in that: get u to single order local derviation and v to the single order local derviation on average then P
I, jThe second order local derviation at place:
(P″
i,j)
uv=(P
i+1,j-P
i-1,j+P
i,j+1+P
i,j-1)/4。
5. the shape-maintaining fitting method during image according to claim 1 is processed is characterized in that: according to the scaling ratio Coons patch is resampled, adopt area sampling to obtain the gray-scale value of resample points:
The image enlargement factor is s, original image I (u then, v) at u, the v direction resamples according to the interval of 1/s, namely obtains amplifying the image after s times, is about to the digital picture I (u of m * n, v) interpolation is the image I of m ' * n ' ' (u ', v '), wherein resample points I ' (u ', v ') corresponding to the position of the pixel in the original image is:
u=u′/s;
v=v′/s;
The label of this pixel bicubic Coons patch at place in original image is i=int (u), j=int (v), and the employing area sampling obtains the gray-scale value of resample points, and is as follows:
If on the border of former Coons patch, then the value of its area sampling should be determined by its two adjacent bicubic Coons curved surfaces resample points jointly, that is: just
6. the shape-maintaining fitting method during image according to claim 5 is processed is characterized in that: wherein, adopt area sampling to obtain the resample points gray-scale value and comprise: utilize following method that finding the solution of integration accelerated:
(4) are found the solution, get (5) formula
I=(IntegralF
0(x2)-IntegralF
0(x1))*(P
i,j*(IntegralF
0(y2)-IntegralF
0(y1))+P
i,j+1*(IntegralF1(y2)-IntegralF1(y1))+(P′
i,j)v*(IntegralG0(y2)-IntegralG0(y1))+(P′
i,j+1)
v*(IntegralG1(y2)-IntegralG1(y1)))+(IntegralF1(x2)-IntegralF1(x1))*(P
i+1,j*(IntegralF0(y2)-IntegralF0(y1))+P
i+1,j+1*(IntegralF1(y2)-IntegralF1(y1))+(P′
i+1,j)
v*(IntegralG0(y2)-IntegralG0(y1))+(P′
i+1,j+1)
v*(IntegralG1(y2)-IntegralG1(y1)))+(IntegralG0(x2)-IntegralG0(x1))*((P′
i,j)
u*(IntegralF0(y2)-IntegralF0(y1))+(P′
i,j+1)
u*(IntegralF1(y2)-IntegralF1(y1))+(P″
i,j)
uv*(IntegralG0(y2)-IntegralG0(y1))+(P″
i,j+1)
uv*(IntegralG1(y2)-IntegralG1(y1)))+(IntegralG1(x2)-IntegralG1(x1))*((P′
i+1,j)
u*(IntegralF0(y2)-IntegralF0(y1))+(P′
i+1,j+1)
u*(IntegralF1(y2)-IntegralF1(y1))+(P″
i+1,j)
uv*(IntegralG0(y2)-IntegralG0(y1))+(P″
i+1,j+1)
uv*(IntegralG1(y2)-IntegralG1(y1)))(5)
Wherein, F
0(), F
1(), G
0(), G
1() is
Integral represents the integration to the function of its back, then F
0(), F
1(), G
0(), G
1The integrated form of () is as follows: IntegralF0 (t)=t
4/ 2-t
3+ t, IntegralF1 (t)=-t
4/ 2+t
3, IntegralG0 (t)=t
4/ 4-2t
3/ 3+t
2/ 2, IntegralG1 (t)=t
4/ 4-t
3/ 3, for F
0(), F
1(), G
0(), G
1The integration of () is in the interval
Integration on 1<=i<=s-1
1<=i<=s-1,
1<=i<=s-1,
1<=i<=s-1,
1<=i<=s-1 calculates in advance and stores, and wherein s is the scaling multiple, then when the sampled value of zoning, searches for corresponding integrated value in record, and (5) formula of bringing into gets final product to get integral result, wherein, and P
I, j, P
I+1, j, P
I, j+1, P
I+1, j+1Be four consecutive point, P
I, jThe place u to derivative be (P '
I, j)
u, v to derivative be (P '
I, j)
v, P
I, j+1The place u to derivative be (P '
I, j+1)
u, v to derivative be (P '
I, j+1)
v, P
I+1, jThe place u to derivative be (P '
I+1, j)
u, v to derivative be (P '
I+1, j)
v, P
I+1, j+1The place u to derivative be (P '
I+1, j+1)
u, v to derivative be (P '
I+1, j+1) v, x1, y1, x2, y2 are range of integration.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200710115786 CN101197043B (en) | 2007-12-19 | 2007-12-19 | Shape-maintaining fitting algorithm in image processing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200710115786 CN101197043B (en) | 2007-12-19 | 2007-12-19 | Shape-maintaining fitting algorithm in image processing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101197043A CN101197043A (en) | 2008-06-11 |
CN101197043B true CN101197043B (en) | 2013-03-13 |
Family
ID=39547425
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200710115786 Expired - Fee Related CN101197043B (en) | 2007-12-19 | 2007-12-19 | Shape-maintaining fitting algorithm in image processing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101197043B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140064590A1 (en) * | 2011-09-20 | 2014-03-06 | Toshiba Medical Systems Corporation | Image processing apparatus and medical image diagnosis apparatus |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663666B (en) * | 2012-03-27 | 2014-05-14 | 中国人民解放军国防科学技术大学 | Two-dimensional image resampling algorithm accelerator based on field-programmable gate array (FPGA) |
-
2007
- 2007-12-19 CN CN 200710115786 patent/CN101197043B/en not_active Expired - Fee Related
Non-Patent Citations (5)
Title |
---|
Cai—Ming Zhan et al..Determining Knots by Minimizing Energy.《计算机科学技术学报(英文版,2006)》.2006,第21卷(第2期),261-264. * |
张彩明,汪嘉业.可调整C2四次Bézier 插值曲线的构造.《计 算 机 学 报》.2004,第27卷(第12期),1665-1671. |
张彩明,汪嘉业.可调整C2四次Bézier 插值曲线的构造.《计 算 机 学 报》.2004,第27卷(第12期),1665-1671. * |
张爱武,张彩明.几何造型中不同目标函数特点研究.《系 统 仿 真 学 报》.2005,第17卷(第3期),674-681. * |
毕重科,张彩明.应用于CT数据的插值算法.《第一届中国图学大会》.2007,371-374. * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140064590A1 (en) * | 2011-09-20 | 2014-03-06 | Toshiba Medical Systems Corporation | Image processing apparatus and medical image diagnosis apparatus |
US9031301B2 (en) * | 2011-09-20 | 2015-05-12 | Kabushiki Kaisha Toshiba | Image processing apparatus and medical image diagnosis apparatus |
Also Published As
Publication number | Publication date |
---|---|
CN101197043A (en) | 2008-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nöllenburg et al. | Morphing polylines: A step towards continuous generalization | |
CN107256557B (en) | Error-controllable subdivision surface image vectorization method | |
CN111445579B (en) | Three-dimensional terrain model adjusting method considering vector element natural feature limitation | |
CN102521794A (en) | Image interpolation method and device based on spline surface | |
CN107330860A (en) | A kind of rational interpolation Zoom method based on CT image borders | |
CN103366342A (en) | Piecewise linear interpolation method applied to video image amplification | |
CN105957002B (en) | Image interpolation amplification method and device based on triangular mesh | |
Zeng et al. | Region-based bas-relief generation from a single image | |
CN104657978A (en) | Road extracting method based on shape characteristics of roads of remote sensing images | |
Hoschek et al. | Turbine blade design by lofted B-spline surfaces | |
CN101197043B (en) | Shape-maintaining fitting algorithm in image processing | |
CN101714258A (en) | Graphics processing systems | |
Alboul et al. | Best data-dependent triangulations | |
Berntsson et al. | Coefficient identification in PDEs applied to image inpainting | |
Akel et al. | Dense DTM generalization aided by roads extracted from LiDAR data | |
Ruprecht et al. | Deformed cross‐dissolves for image interpolation in scientific visualization | |
Fortes et al. | Filling holes with shape preserving conditions | |
CN103020936A (en) | Super-resolution reconstruction method of facial image | |
CN104915921A (en) | Triangular-mesh-deformation-based geometrical boundary image mapping method with kept content | |
Yang | Matrix weighted rational curves and surfaces | |
JP2003132353A (en) | Center line generator | |
Sarfraz | Fitting curve to planar digital data | |
CN101599169A (en) | Remote sensing image amplification method based on PDE and small echo | |
Islam et al. | Overview and challenges of different image morphing algorithms | |
Rivera et al. | Entropy controlled gauss-markov random measure field models for early vision |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130313 Termination date: 20131219 |