CN100511280C - New method for restoring disrepaired image - Google Patents

Abstract

A digitalized repairing-method of damaged image includes recording image to be repaired as I, using image-dividing means to confirm out damaged region of image, generating a binary image according damaged and undamaged regions in image I, labeling undamaged region in image I as omega, using structure element E to make mathematic form corrosion operation on binary image to obtain corrosion binary image, calculating local variation rate of corrosion image, using structure element E to make expansion operation on region omega and carrying out repairing process according to preset steps.

Description

The new method that disrepaired image through digitization is repaired

(1), technical field:

Relate to method, the especially method that the image of local damage can be repaired by digitizing technique that a kind of view data is handled.

(2), background technology:

Image repair (Image Inpainting) is the digitizing recovery technique to archives, ancient painting, film: by take pictures, mode such as scanning, mould/number conversion becomes them after the digital picture, utilizes computing machine reparation again.The standard of repairing is not see the vestige of handling is arranged as far as possible.Restorative procedure master commonly used is pure manual method, and softwares such as digital picture use Photoshop are dealt with, and progressively smears damaged zone, and this processing procedure efficient is very low.Some methods of handling have automatically been studied at present in the world, known two classes that mainly contain, one class is based on the method for partial differential equation, use some smoothness assumption to separate the partial differential equation problem repairing the question resolves itself into, the shortcoming of these class methods is to repair little and elongated zone, and the damaged zone for roomy excessively level and smooth effect can occur after the repairing, leave obvious marks, poor effect.The another kind of synthetic restorative procedure of texture that is based on is searched for a suitable texture automatically in image, it is copied to damaged zone, but the process of this search is quite time-consuming, causes the reparation speed of these class methods very slow.

(3), summary of the invention:

Conventional images restorative procedure speed is slow, efficient is low in order to overcome, and the shortcoming that can't repair enlarged regions, the purpose of this invention is to provide a kind of new image repair method, this method not only can be repaired roomy damaged zone, and fast operation, than approximately fast 100 times of the methods of synthesizing based on texture.

The technical solution adopted for the present invention to solve the technical problems is as follows:

The new method that a kind of disrepaired image through digitization of the present invention is repaired, the step of this method is as follows:

[1], an image to be repaired is designated as I, and by some known image partition methods or man-machine interactively method determine damaged zone in the image.Not damaged and damaged zone according to image I generates a width of cloth bianry image, and wherein the not damaged area relative of image I partly is labeled as Ω.

[2], utilization structure element E does the computing of mathematics morphological erosion to bianry image, the bianry image ε that obtains corroding _E(Ω).

[3], calculate corrosion diagram as ε _EBorder (Ω) On the localized variation rate of being had a few, method is as follows: establish image I and have n component (I is a gray level image during n=1, and I is the RGB coloured image during n=3), the structure tensor G by following formula computed image I at first,

$G=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\▿I}_i{{\▿I}_i}^T,$

I in the formula _iI component of presentation video, Be the gradient of this component, can calculate with known central difference method.Obtain characteristic direction and the eigenwert of tensor G then by following characteristic value decomposition,

$G=\left[{\mathrm{\θ}}_{+}{\mathrm{\θ}}_{-}\right]\left[\begin{array}{cc}{\mathrm{\λ}}_{+}& 0\\ 0& {\mathrm{\λ}}_{-}\end{array}\right]\left[\begin{array}{c}{{\mathrm{\θ}}_{+}}^T\\ {{\mathrm{\θ}}_{-}}^T\end{array}\right]$

θ wherein ₊It is bigger eigenvalue ₊Pairing proper vector, it can be defined as the localized variation rate of image I at set point.

[4], regional Ω is made dilation operation, obtain regional δ with structural element E _E(Ω).Make d Ω=δ _E(Ω)-Ω represents current zone to be repaired, then for arbitrfary point x ^*∈ d Ω, repair (referring to Fig. 1) according to following steps:

(4.1), seek

Go up and x ^*The point that distance is enough approaching is promptly sought the set of point

$N\left(x^{*}\right)=\{x\∈\∂\left[{\mathrm{\ϵ}}_E\left(\mathrm{\Ω}\right)\right]:\left|\right|x-x^{*}\left|\right|\≤r\},---\left(1\right)$

Wherein r is a predefined search radius.The value of r will be considered the compromise of search speed and search precision.

(4.2), use following formula set of computations N (x ^*) in an optimum x ₀:

$x_0=\arg \min \{\frac{|{\mathrm{\θ}}_{+}\left(x\right)\·(x^{*}-x)|}{\left|\right|x^{*}-x\left|\right|}+\frac{\mathrm{\α}}{\sqrt{{\mathrm{\λ}}_{+}\left(x\right)}},x\∈N\left(x^{*}\right)\},---\left(2\right)$

Wherein α be one greater than zero parameter.The connotation of this formula is to guarantee at an x ₀The straight line x of place ^*-x ₀With characteristic direction θ ₊(x ₀) between near vertical, guarantee eigenvalue simultaneously ₊(x ₀) be maximum.

(4.3), defining point x _m=(x ^*+ x ₀)/2 use following formula to repair breaking point x ^*The pixel value at place:

$I\left(x^{*}\right)=\left\{\begin{array}{cc}2I\left(x_m\right)-I\left(x_0\right),& x_m\&Element;\mathrm{\&Omega;and}\left|\right|\mathrm{\&Delta;I}\left|\right|<\mathrm{\&epsiv;},\\ 2I\left(x_m\right),& x_m\&Element;\mathrm{\&Omega;and}\left|\right|\mathrm{\&Delta;I}\left|\right|\&GreaterEqual;\mathrm{\&epsiv;},\\ I\left(x_0\right),& x_m\&NotElement;\mathrm{\&Omega;},\end{array}\right.---\left(3\right)$

‖ Δ I ‖=‖ I (x wherein _m)-I (x ₀) ‖, Ω still represents unbroken image-region, ε is a positive parameter.

[5], make ε _E(Ω)=and Ω, Ω=δ _E(Ω), and returned for the 3rd step, all breaking points are repaired and finish in image I.

Related symbol connotation is concluded and is described as follows in the above-mentioned steps:

The image that I is to be repaired

Unbroken zone in the Ω image I

The structural element of E mathematical morphology

ε _E(Ω) regional Ω is made the result of erosion operation by structural element E

Zone ε _EBorder (Ω)

The structure tensor of G image I

I _iI the component of presentation video I

I _iThe gradient of component

λ ₊, two eigenwerts of λ _ structure tensor G, λ ₊〉=λ _

θ ₊, θ _ eigenvalue ₊And λ _ pairing proper vector

δ _E(Ω) regional Ω is made the result of dilation operation by structural element E

Zone to be repaired, i.e. δ in the d Ω current iteration _E(Ω)-Ω

x ^*Point to be repaired among the d Ω of zone

N (x ^*) border

Go up and some x ^*Enough set of approaching point of distance are defined by formula (1)

The predefined search radius of r

x ₀Set N (x ^*) in an optimum, define by formula (2)

α is greater than zero parameter

x _mPoint x ^*With an x ₀Mid point, i.e. (x ^*+ x ₀)/2

‖ Δ I ‖ pixel value I (x _m) and I (x ₀) the norm of difference, i.e. ‖ I (x _m)-I (x ₀) ‖

One of ε is on the occasion of parameter

Technique effect of the present invention is as follows:

View data for breakage provides effective digitizing restorative procedure, can repair roomy damaged zone (as Fig. 3, Fig. 4, Fig. 5 and shown in Figure 6), overcome the shortcoming that the existing image repair method method of partial differential equation (particularly based on) can only be repaired the damaged zone of " fine rule " shape.Method remediation efficiency height of the present invention, speed is fast, than general fast two orders of magnitude of the method for synthesizing based on texture.

(4), description of drawings:

Fig. 1 the present invention carries out the schematic diagram of image repair.

The software flow pattern of Fig. 2 a specific embodiment of the present invention.

Fig. 3 carries out the example of image repair with method of the present invention.

Fig. 4 carries out the example of image repair with method of the present invention.

Fig. 5 carries out the example of image repair with method of the present invention.

Fig. 6 carries out the example of image repair with method of the present invention.

(5) embodiment:

A specific embodiment of the present invention is as follows:

Consulting shown in Figure 1ly, is the schematic diagram that the present invention carries out image repair.Pixel in each grid representative image among the figure, the zone that the region representation of dotted line grid is to be repaired, wherein the grid of grey is represented the breaking point x that current needs are repaired ^*The not damaged regional Ω of the region representation of solid line grid.Part adjacent with the area to be repaired among the Ω is called the border

Set N (the x of the light gray areas representative point among the Ω ^*).

Consult the process flow diagram of Fig. 2, at first import an image I to be repaired, according to the not damaged regional Ω and the damaged zone generation bianry image of image I.

In second step, the square structure element E that uses 3 * 3 makes erosion operation to regional Ω, the bianry image ε that obtains corroding _E(Ω).

In the 3rd step, utilize aforementioned techniques computation schemes border

On the localized variation rate θ that had a few ₊

The 4th step, with structural element E regional Ω is made dilation operation, obtain δ _E(Ω).Make d Ω=δ _E(Ω)-and Ω, for the some x that does not repair among the d Ω ^*, repair with the following methods:

(a) according to the set N (x of formula (1) calculation level ^*), wherein get parameter $r=2\sqrt{5};$

(b) according to formula (2) set of computations N (x ^*) in an optimum x ₀, wherein get parameter alpha=1;

(c) repair some x according to formula (3) ^*Pixel value I (the x at place ^*), wherein get parameter ε=5.

In the 5th step, make ε _E(Ω)=and Ω, Ω=δ _E(Ω), check whether Ω comprises the entire image zone.If then repair and finish; Otherwise returning for the 3rd step continues to repair.

Seeing also shown in Fig. 3,4,5,6, is the example that method of the present invention is carried out image repair.The left figure of Fig. 3 is damaged image, and the centre is the result who repairs with Partial Differential Equation method, and right figure is the result who repairs with method of the present invention.Fig. 4 erases a bunch of flowers among the left figure automatically with method of the present invention.Fig. 5 is with the damaged zone in the method repairing retina image of the present invention.Fig. 6 erases arm and the microphone of speaker among the left figure automatically with method of the present invention.

Claims (1)

1, a kind of new method of disrepaired image through digitization reparation, it is characterized in that: this method step is as follows:

[1], an image to be repaired is designated as I, and by image partition method or man-machine interactively method determine damaged zone in the image; Not damaged and damaged zone according to image I generates a width of cloth bianry image, and wherein the not damaged area relative of image I partly is labeled as Ω;

[2], utilization structure element E does the computing of mathematics morphological erosion to bianry image, the bianry image ε that obtains corroding _E(Ω);

[3], calculate corrosion diagram as ε _EBorder (Ω)

On the localized variation rate of being had a few, method is as follows: establish image I and have n component, the structure tensor G by following formula computed image I at first,

$G=\underset{i=1}{\overset{n}{\mathrm{\&Sigma;}}}\&dtri;I_i\&dtri;{I_i}^T,$

I in the formula _iI component of presentation video,

Be the gradient of this component, calculate, obtain characteristic direction and the eigenwert of tensor G then by following characteristic value decomposition with known central difference method,

$G=\left[{\mathrm{\&theta;}}_{+}{\mathrm{\&theta;}}_{-}\right]\left[\begin{array}{cc}{\mathrm{\&lambda;}}_{+}& 0\\ 0& {\mathrm{\&lambda;}}_{-}\end{array}\right]\left[\begin{array}{c}{{\mathrm{\&theta;}}_{+}}^T\\ {{\mathrm{\&theta;}}_{-}}^T\end{array}\right]$

θ wherein ₊It is bigger eigenvalue ₊Pairing proper vector, it is defined as the localized variation rate of image I at set point;

[4], regional Ω is made dilation operation, obtain regional δ with structural element E _E(Ω); Make d Ω=δ _E(Ω)-Ω represents current zone to be repaired, then for arbitrfary point x ^*∈ d Ω, repair according to following steps:

(4.1), seek

Go up and x ^*The point that distance is enough approaching is promptly sought the set of point

$N\left(x^{*}\right)=\{x\&Element;\&PartialD;\left[{\mathrm{\&epsiv;}}_E\left(\mathrm{\&Omega;}\right)\right]:\left|\right|x-x^{*}\left|\right|\&le;r\},$

Wherein r is a predefined search radius, and the value of r will be considered the compromise of search speed and search precision;

(4.2), use following formula set of computations N (x ^*) in an optimum x ₀:

$x_0=\arg \min \{\frac{|{\mathrm{\&theta;}}_{+}\left(x\right)\&CenterDot;(x^{*}-x)|}{\left|\right|x^{*}-x\left|\right|}+\frac{\mathrm{\&alpha;}}{\sqrt{{\mathrm{\&lambda;}}_{+}\left(x\right)}},x\&Element;N\left(x^{*}\right)\},$

Wherein α be one greater than zero parameter; The connotation of this formula is to guarantee at an x ₀The straight line x of z place ^*-x ₀With characteristic direction θ ₊(x ₀) between near vertical, guarantee eigenvalue simultaneously ₊(x ₀) be maximum;

(4.3), defining point x _m=(x ^*+ x ₀)/2 use following formula to repair breaking point x ^*The pixel value at place:

$I\left(x^{*}\right)=\left\{\begin{array}{cc}2I\left(x_m\right)-I\left(x_0\right),& x_m\&Element;\mathrm{\&Omega;and}\left|\right|\mathrm{\&Delta;I}\left|\right|<\mathrm{\&epsiv;},\\ 2I\left(x_m\right),& x_m\&Element;\mathrm{\&Omega;and}\left|\right|\mathrm{\&Delta;I}\left|\right|\&GreaterEqual;\mathrm{\&epsiv;},\\ I\left(x_0\right),& x_m\&NotElement;\mathrm{\&Omega;},\end{array}\right.$

‖ Δ I ‖=‖ I (x wherein _m)-I (x ₀) ‖, Ω still represents unbroken image-region, ε is a positive parameter;

[5], make ε _E(Ω)=and Ω, Ω=δ _E(Ω), and returned for [3] step, all breaking points are repaired and finish in image I.

Images