Apparatus and method for noncausal, predictive, digital image compression

Claim 1

A method of converting a noncausal, two dimensional signal representative of a digitized image into a statistically equivalent recursive model for use in processing the digital image, comprising the steps of: (i) digitizing an image to produce a two dimensional signal; (ii) filtering the two dimensional signal to arrive at data representative of an optimal set of parameters of the two dimensional signal by: a. generating data representative of a potential matrix which represents the two dimensional signal for a given set of parameters; b. filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said given set of parameters using a Riccati iteration filtering technique; c. inputting white noise to said equivalent recursive model to generate k samples of the two dimensional signal; d. obtaining an average value for χ.sub.τ.sup.(k) where χ.sup.(k).sub.τ is data representative of a correlation of lag τ computed from sample k; e. generating data representative of a gradient of a likelihood function by solving ##EQU28## f. minimizing data representative of said likelihood function utilizing said data representative of said gradient of said likelihood function to arrive at a set of parameters; g. comparing said arrived at set of parameters to a previously arrived at set of parameters to determine whether the difference therebetween is within a desired tolerance; and h. repeating steps (a) through (g) until said data representative of said optimal set of parameters is determined; (iii) generating data representative of a potential matrix which represents the two dimensional signal for said data representative of said optimal set of parameters; (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through said Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using said Riccati iteration filtering technique; and (v) processing the two dimensional signal using said recursive model.

Claim 2

A method of converting a noncausal, two dimensional signal representative of a digitized image into a statistically equivalent recursive model for use in processing the digital image, comprising the steps of: (i) digitizing an image to produce a two dimensional signal; (ii) filtering the two dimensional signal to arrive at data representative of an optimal set of parameters of the two dimensional signal by: a. generating data representative of a potential matrix which represents the two dimensional signal for a given set of parameters; b. filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said given set of parameters using a Riccati iteration filtering technique; c. inputting white noise to said equivalent recursive model to generate k samples of the two dimensional signal; d. obtaining an average value for χ.sub.τ.sup.(k) where χ.sup.(k).sub.τ is data representative of a correlation of lag τ computed from sample k; e. generating data representative of a gradient of a likelihood function by solving ##EQU28## f. minimizing data representative of said likelihood function utilizing said data representative of said gradient of said likelihood function to arrive at a set of parameters; g. comparing said arrived at set of parameters to a previously arrived at set of parameters to determine whether the difference therebetween is within a desired tolerance; and h. repeating steps (a) through (g) until said data representative of said optimal set of parameters is determined; (iii) generating data representative of a potential matrix which represents the two dimensional signal for said data representative of said optimal set of parameters; (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through said Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using said Riccati iteration filtering technique; and (v) processing the two dimensional signal using said recursive model.

Claim 6

A method of converting a noncausal, two dimensional signal representative of a digitized image expressed as a first order Gauss Markov random field (GMRF) into a statistically equivalent recursive model for use in processing the digitized image, comprising the steps of: (i) digitizing an image to produce a two dimensional signal; (ii) filtering the two dimensional signal to arrive at data representative of an optimal set of parameters of the two dimensional signal by: a. generating data representative of a gradient of a likelihood function by solving ##EQU30## b. minimizing data representative of said likelihood function utilizing said data representative of said gradient of said likelihood function to arrive at a set of parameters; c. comparing said arrived at set of parameters to a previously arrived at set of parameters to determine whether the difference therebetween is within a desired tolerance; and d. repeating steps (a) through (c) until said data representative of said optimal set of parameters is determined; (iii) generating data representative of a potential matrix which represents the two dimensional signal for said data representative of said optimal set of parameters; and (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using a Riccati iteration filtering technique and (v) processing the two dimensional signal using said recursive model.

Claim 7

A method of converting a noncausal, two dimensional signal representative of a digitized image expressed as a first order Gauss Markov random field (GMRF) into a statistically equivalent recursive model for use in processing the digitized image, comprising the steps of: (i) digitizing an image to produce a two dimensional signal; (ii) filtering the two dimensional signal to arrive at data representative of an optimal set of parameters of the two dimensional signal by: a. generating data representative of a gradient of a likelihood function by solving ##EQU30## b. minimizing data representative of said likelihood function utilizing said data representative of said gradient of said likelihood function to arrive at a set of parameters; c. comparing said arrived at set of parameters to a previously arrived at set of parameters to determine whether the difference therebetween is within a desired tolerance; and d. repeating steps (a) through (c) until said data representative of said optimal set of parameters is determined; (iii) generating data representative of a potential matrix which represents the two dimensional signal for said data representative of said optimal set of parameters; and (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using a Riccati iteration filtering technique and (v) processing the two dimensional signal using said recursive model.

Claim 8

A method of converting a noncausal, two dimensional signal representative of a digitized image expressed as a first order Gauss Markov random field (GMRF) with zero boundary conditions into a statistically equivalent recursive model for use in processing the digitized image, comprising the steps of: (i) digitizing an image to produce a two dimensional signal; (ii) filtering the two dimensional signal to arrive at data representative of an optimal set of parameters of the two dimensional signal by solving the following equations: ##EQU31## (iii) generating data representative of a potential matrix which represents the two dimensional signal for said data representative of said optimal set of parameters; and (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using a Riccati iteration filtering technique; and (v) processing the two dimensional signal using said recursive model.

Claim 9

A method of converting a noncausal, two dimensional signal representative of a digitized image expressed as a first order Gauss Markov random field (GMRF) with zero boundary conditions into a statistically equivalent recursive model for use in processing the digitized image, comprising the steps of: (i) digitizing an image to produce a two dimensional signal; (ii) filtering the two dimensional signal to arrive at data representative of an optimal set of parameters of the two dimensional signal by solving the following equations: ##EQU31## (iii) generating data representative of a potential matrix which represents the two dimensional signal for said data representative of said optimal set of parameters; and (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using a Riccati iteration filtering technique; and (v) processing the two dimensional signal using said recursive model.

Claim 10

A method of compressing a signal representative of a digital image, comprising the steps of: (i) digitizing an image to produce a digital image signal (ii) filtering the digital image signal to arrive at data representative of an optimal set of parameters of the digital image signal by: a. generating data representative of a potential matrix which represents the digital image signal for a given set of parameters; b. filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said given set of parameters using a Riccati iteration filtering technique; c. inputting white noise to said equivalent recursive model to generate k samples of the digital image signal; d. obtaining an average value for χ.sub.τ.sup.(k) where χ.sup.(k) .sub.τ is data representative of a correlation of lag τ computed from sample k; e. generating data representative of a gradient of a likelihood function by solving ##EQU32## f. minimizing data representative of said likelihood function utilizing said data representative of said gradient of said likelihood function to arrive at a set of parameters; g. comparing said arrived at set of parameters to a previously arrived at set of parameters to determine whether the difference therebetween is within a desired tolerance; and h. repeating steps (a) through (g) until said data representative of said optimal set of parameters is determined; (iii) generating data representative of a potential matrix which represents the digital image signal for said data representative of said optimal set of parameters; (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through said Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using said Riccati iteration filtering technique; and (v) compressing said digital image signal by comparing said statistically equivalent recursive model for said data representative of said optimal set of parameters to the digital image signal to establish an error image signal, said error image signal and said optimal set of parameters being an equivalent of the digital image signal but requiring fewer bits than the digital image signal for the same level of perceived quality.

Claim 11

A method of compressing a signal representative of a digital image, comprising the steps of: (i) digitizing an image to produce a digital image signal (ii) filtering the digital image signal to arrive at data representative of an optimal set of parameters of the digital image signal by: a. generating data representative of a potential matrix which represents the digital image signal for a given set of parameters; b. filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said given set of parameters using a Riccati iteration filtering technique; c. inputting white noise to said equivalent recursive model to generate k samples of the digital image signal; d. obtaining an average value for χ.sub.τ.sup.(k) where χ.sup.(k) .sub.τ is data representative of a correlation of lag τ computed from sample k; e. generating data representative of a gradient of a likelihood function by solving ##EQU32## f. minimizing data representative of said likelihood function utilizing said data representative of said gradient of said likelihood function to arrive at a set of parameters; g. comparing said arrived at set of parameters to a previously arrived at set of parameters to determine whether the difference therebetween is within a desired tolerance; and h. repeating steps (a) through (g) until said data representative of said optimal set of parameters is determined; (iii) generating data representative of a potential matrix which represents the digital image signal for said data representative of said optimal set of parameters; (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through said Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using said Riccati iteration filtering technique; and (v) compressing said digital image signal by comparing said statistically equivalent recursive model for said data representative of said optimal set of parameters to the digital image signal to establish an error image signal, said error image signal and said optimal set of parameters being an equivalent of the digital image signal but requiring fewer bits than the digital image signal for the same level of perceived quality.

Claim 17

A method of smoothing a signal representative of a distorted digital image, comprising the steps of: (i) digitizing an image to produce a digital image signal; (ii) filtering the digital image signal to arrive at data representative of an optimal set of parameters of the digital image signal by: a. generating data representative of a potential matrix which represents the digital image signal for a given set of parameters; b. filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said given set of parameters using a Riccati iteration filtering technique; c. inputting white noise to said equivalent recursive model to generate k samples of the digital image signal; d. obtaining an average value for χ.sub.τ.sup.(k) where χ.sup.(k).sub.τ is data representative of a correlation of lag τ computed from sample k; e. generating data representative of a gradient of a likelihood function by solving ##EQU33## f. minimizing data representative of said likelihood function utilizing said data representative of said gradient of said likelihood function to arrive at a set of parameters; g. comparing said arrived at set of parameters to a previously arrived at set of parameters to determine whether the difference therebetween is within a desired tolerance; and h. repeating steps (a) through (g) until said data representative of said optimal set of parameters is determined; (iii) generating data representative of a potential matrix which represents the digital image signal for said data representative of said optimal set of parameters; (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through said Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using said Riccati iteration filtering technique; (v) subjecting said statistically equivalent recursive model to a recursive smoothing technique; and (vi) producing a signal representative of an enhanced digital image from said model.

Claim 18

A method of smoothing a signal representative of a distorted digital image, comprising the steps of: (i) digitizing an image to produce a digital image signal; (ii) filtering the digital image signal to arrive at data representative of an optimal set of parameters of the digital image signal by: a. generating data representative of a potential matrix which represents the digital image signal for a given set of parameters; b. filtering said potential matrix data to arrive at a statistically equivalent recursive model through a Cholesky factorization of said potential matrix data for said given set of parameters using a Riccati iteration filtering technique; c. inputting white noise to said equivalent recursive model to generate k samples of the digital image signal; d. obtaining an average value for χ.sub.τ.sup.(k) where χ.sup.(k).sub.τ is data representative of a correlation of lag τ computed from sample k; e. generating data representative of a gradient of a likelihood function by solving ##EQU33## f. minimizing data representative of said likelihood function utilizing said data representative of said gradient of said likelihood function to arrive at a set of parameters; g. comparing said arrived at set of parameters to a previously arrived at set of parameters to determine whether the difference therebetween is within a desired tolerance; and h. repeating steps (a) through (g) until said data representative of said optimal set of parameters is determined; (iii) generating data representative of a potential matrix which represents the digital image signal for said data representative of said optimal set of parameters; (iv) filtering said potential matrix data to arrive at a statistically equivalent recursive model through said Cholesky factorization of said potential matrix data for said data representative of said optimal set of parameters using said Riccati iteration filtering technique; (v) subjecting said statistically equivalent recursive model to a recursive smoothing technique; and (vi) producing a signal representative of an enhanced digital image from said model.