Iterative image restoration device

Claim 2

An image restoration device comprising: an optical viewing member for providing a view of a subject, the optical viewing member responding to each point of the subject according to a response function to provide a viewed image of the subject; a data processor for receiving from the optical viewing member and storing a series of viewed images of the subject as viewed through the optical viewing member; computer processing means coupled to the data processor for responding to the viewed images and restoring the viewed images to substantially noiseless, undistorted images, the computer processing means restoring each viewed image by iteratively determining for each point in a viewed image a factor which minimizes noise and distortion at that point, the factor being iteratively determined, according to a formulation of a mathematical definition of the factor, by performance of transform calculations including a mathematical division operation in a transform domain of functions which include the response function of the optical viewing member to determine the factor, transforms of the response function being denominators in the division operation, the division operation enabling avoidance of performance of a matrix inversion operation, and with the determined factors for the points of a viewed image the computer processing means minimizing noise and distortion at each point of the image such that a substantially noiseless, undistorted image is formed; and display means coupled to the computer processing means for providing a view of the substantially noiseless, undistorted images in a non-transform domain.

Claim 3

An image restoration device as claimed in claim 1 wherein the factor is iteratively determined through a division operation in a Fourier transform domain of the response function.

Claim 4

An image restoration device as claimed in claim 1 wherein the factor is iteratively determined by the relationship c.sub.n+1 =c.sub.n - Δ where c is the factor; n is the number of iterations of c; and Δ=inverse Fourier transform of .sup.(Ψ' /.sub.(H H* +α+λ)); where Ψ' is the Fourier transform of ψ', ψ'=(M+αI) c - g; M is a matrix of viewed image points i,j, and defined as ##EQU6## D(c) is the set of viewed image points in which Σk i (x) c i >0, and k i (x) is the response function of the optical viewing member at image point i over surrounding points x; α and λ are constants (≧0); I is the Identity matrix; g is an array of measured values of the points in the image; H is the Fourier transform of h(x) which is another form of the response function k i (x); and H* is the Fourier transform of h(x) in which h(x)=h(-x).

Claim 5

An image restoration device as claimed in claim 3 wherein h(x) is the point spread function of the optical viewing member such that k.sub.i (x)=h(y.sub.i - x) where i is a point of the viewed image; and Y i is the corresponding point on the subject.

Claim 6

An image restoration device as claimed in claim 1 wherein the series of viewed images received by the data processor from the optical viewing member includes a series of 2-D images of different planes through a 3-D image of the subject.

Claim 7

An image restoration device as claimed in claim 1 where the formulation of the mathematical definition of the factor includes a related matrix which replaces a Hessian matrix of the factor such that the division operation replaces a matrix inversion operation and the factor is iteratively determinable in a transform domain of the response function of the optical viewing member.

Claim 10

A method as claimed in claim 7 wherein the factor is iteratively determined through a division operation with a Fourier transform of the response function.

Claim 11

A method as claimed in claim 7 wherein the factor is iteratively determined by the relationship c.sub.n+ 1=c.sub.n - Δ where c is the factor; n is the number of iterations of c; and Δ=inverse Fourier transform of .sup.(Ψ' /.sub.(H H* +α+λ)); where Ψ' is the Fourier transform of ψ', ψ'=(M+αI) c - g; M is a matrix of viewed image points i,j, and defined as ##EQU7## D(c) is the set of viewed image points in which Σk i (x) c i >0, and k i (x) is the response function of the optical viewing member at image point i over surrounding points x; α and λ are constants (≧0); I is the Identity matrix; g is an array of measured values of the points in the image; H is the Fourier transform of h(x) which is another form of the response function k i (x); and H* is the Fourier transform of h(x) in which h(x)=h(-x).

Claim 12

A method as claimed in claim 9 wherein h(x) is the point spread function of the optical viewing member such that k.sub.i (x)=h(y.sub.i - x) where i is a point of the viewed image; and y 1 is the corresponding point on the subject.

Claim 13

A method as claimed in claim 7 wherein the step of receiving in a data processor a series of viewed images of the subject includes receiving series of images of parallel planes through a 3-D image of the subject as viewed through the optical viewing member.

Claim 14

A method as claimed in claim 7 wherein the step of restoring each viewed image includes a formulation of a mathematical definition of the factor which replaces a Hessian matrix of the factor with a related matrix such that a division operation replaces a matrix inversion operation and the factor is iteratively determinable in a transform domain of the response function of the optical viewing member.

Claim 15

An image restoration device comprising: (i) an optical viewing member for providing a view of a subject, the optical viewing member responding to each point of the subject according to a response function to provide a viewed image of the subject; (ii) a data processor for receiving from the optical viewing member and storing a series of viewed images of the subject as viewed through the optical viewing member; (III) computer processing means associated with the data processor for restoring the viewed images to substantially noiseless, undistorted images, the computer processing means restoring each viewed image by: (a) iteratively determining for each point in a viewed image a factor which minimizes noise and distortion at that point, the factor being iteratively determined according to a formulation of a mathematical definition of the factor, by performance of transform calculations including a division operation in a transform domain of functions which include the response function of the optical viewing member to determine the factor; the factor being iteratively determined by the relationship c.sub.n+1 =c.sub.n - Δ where c is the factor; n is the number of iterations of c; and Δ=inverse Fourier transform of (.sup.Ψ' /.sub.( H H*+α+λ)); where Ψ' is the Fourier transform of ψ', ψ'=(M+αI) c - g; M is a matrix of viewed image points i, j, and defined as ##EQU8## D(c) is the set of viewed image points in which Σk i (x) c i >0, and k i (x) is the response function of the optical viewing member at image point i over surrounding points x; α and λ are constants (≧0); I is the Identity matrix; g is an array of measured values of the points in the image; H is the Fourier transform of h(x) which is another form of the response function k i (x); and H* is the Fourier transform of h(x) in which h(x)=h(-x); and (b) with the determined factors of a viewed image, minimizing noise and distortion at each point of the viewed image such that a substantially noiseless, undistorted image is formed; and (iv) display means coupled to the computer processing means for providing a view of the substantially noiseless, undistorted images in a non-transform domain.

Claim 16

An image restoration device comprising: (i) an optical viewing member for providing a view of a subject, the optical viewing member responding to each point of the subject according to a response function to provide a viewed image of the subject; (ii) a data processor for receiving from the optical viewing member and storing a series of viewed images of the subject as viewed through the optical viewing member; (III) computer processing means associated with the data processor for restoring the viewed images to substantially noiseless, undistorted images, the computer processing means restoring each viewed image by: (a) iteratively determining for each point in a viewed image a factor which minimizes noise and distortion at that point, the factor being iteratively determined according to a formulation of a mathematical definition of the factor, by performance of transform calculations including a division operation in a transform domain of functions which include the response function of the optical viewing member to determine the factor; the factor being iteratively determined by the relationship c.sub.n+1 =c.sub.n - Δ where c is the factor; n is the number of iterations of c; and Δ=inverse Fourier transform of (.sup.Ψ' /.sub.( H H*+α+λ)); where Ψ' is the Fourier transform of ψ', ψ'=(M+αI) c - g; M is a matrix of viewed image points i, j, and defined as ##EQU8## D(c) is the set of viewed image points in which Σk i (x) c i >0, and k i (x) is the response function of the optical viewing member at image point i over surrounding points x; α and λ are constants (≧0); I is the Identity matrix; g is an array of measured values of the points in the image; H is the Fourier transform of h(x) which is another form of the response function k i (x); and H* is the Fourier transform of h(x) in which h(x)=h(-x); and (b) with the determined factors of a viewed image, minimizing noise and distortion at each point of the viewed image such that a substantially noiseless, undistorted image is formed; and (iv) display means coupled to the computer processing means for providing a view of the substantially noiseless, undistorted images in a non-transform domain.

Claim 17

An image restoration device as claimed in claim 13 wherein h(x) is the point spread function of the optical viewing member such that k.sub.i (x)=h(y.sub.i - x) where is a point of the viewed image; and y i is the corresponding point on the subject.

Claim 20

A method as claimed in claim 15 wherein h(x) is the point spread function of the optical viewing member such that k.sub.i (x)=h(y.sub.i - x) where i is a point of the viewed image; and y i is the corresponding point on the subject.