CN101441770B - Method for excavating optimum image based on information entropy and logarithm contrast weight sum - Google Patents
Method for excavating optimum image based on information entropy and logarithm contrast weight sum Download PDFInfo
- Publication number
- CN101441770B CN101441770B CN2008102331569A CN200810233156A CN101441770B CN 101441770 B CN101441770 B CN 101441770B CN 2008102331569 A CN2008102331569 A CN 2008102331569A CN 200810233156 A CN200810233156 A CN 200810233156A CN 101441770 B CN101441770 B CN 101441770B
- Authority
- CN
- China
- Prior art keywords
- image
- iqaf
- delta
- value
- information entropy
- 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 Analysis (AREA)
Abstract
The invention relates to a method for digging an optimal image based on information entropy and logarithm contrast weighted sum, which comprises the following steps: (1) detecting the gray scale spectrum distribution of original images, determining the existence of a descending gray scale spectrum distribution, and then calculating the gray scale spectrum continuous bandwidth BW of the original images; (2) establishing an image quality assessment function IQAF; (3) digging a bottom layer image, and calculating the value of the assessment function IQAF by using the information entropy InEn corresponding to a dug image and an average logarithm contrast lgC; and (4) using a digging parameter corresponding to the maximum value of the image quality assessment function IQAF to dig an image which is the image with the optimal quality according to the bottom layer image digging formula in the step (3). The method can quickly dig out the image with the optimal quality objectively from the images with the descending gray scale spectrum distribution of the original images, and the images shot under a condition with quite weak illumination can obtain clear optimal images with objective sceneries through the method.
Description
Technical field
The invention belongs to the computer digital image processing technology field.
Background technology
The bottom layer image mining technology becomes the vision available information with the unavailable information of human vision, can be used for excavating captured image information under the mal-condition, as escaping the information of hiding in financial institution's monitoring image etc. behind the automobile accident.Yet in mining process, just can excavate picture quality the best owing to judging with what excavation parameter; Excavate repeatedly so have in practice so that the image information of excavating to reach human vision the best in quality; Waste time and energy very much; And excavate the subjective feeling that the so-called best quality image that obtains depends on human Different Individual, be not the best quality image on the objective sense.At present, become human still unresolved good problem for how in excavation, directly obtaining image the best in quality both at home and abroad.
Summary of the invention
The purpose of this invention is to provide a kind ofly, use this method can obtain the best quality image on the objective sense apace for having source images that decline type gray scale spectrum distributes method in bottom layer image mining based on the excavating optimum image of information entropy and logarithm contrast weighted sum.
The objective of the invention is to realize that through such technical scheme this method may further comprise the steps:
(1) gray scale of detection resources image spectrum distributes, and confirms to have decline type gray scale spectrum and distributes, then the gray scale of calculation sources image spectrum continuous bandwidth BW;
(2) set up image quality evaluation function IQAF, its computing formula is:
InEn is an image information entropy,
Be the average contrast, the computing formula of image information entropy InEn is:
P in the formula (i) is the pixel count of image on i gray level, when p (i)=0, makes Log
2P (i)=0; The average contrast
Computing formula be:
M, N be image at x, the pixel count of y direction, (x is that ((x+1 y) is pixel (x+l, gray-scale value y) to G to pixel for x, gray-scale value y) y) to G;
(3) carry out bottom layer image mining; And with excavating the value of the corresponding information entropy InEn of image with average logarithm contrast
calculating evaluation function IQAF, the bottom layer image mining formula is:
Constraint condition:
TG (x, y), OG (x, y) expression is excavated image and source images at pixel (x respectively; Y) gray-scale value, Theta [0,255] and Delta [1; 255] for excavating parameter, starting point and the gray-level that excavates image are excavated in expression respectively, and K is the space telescopic factor; For bottom layer image mining, get k=255, Theta=0;
Make Theta=0; Delta increases 1 from 1 to gray scale spectrum continuous bandwidth BW at every turn changes; Source images is carried out bottom layer image mining; Use the value of information entropy InEn with average logarithm contrast
the computed image quality assessment function IQAF of the pairing excavation image of each Delta value, up to the maximal value that obtains IQAF;
(4) image of excavating according to the said bottom layer image mining formula of step (3) with the pairing excavation parameter of image quality evaluation function IQAF maximal value Delta, Theta is image the best in quality.
The detection method (this method is documented in the patent of invention that the patent No. is ZL200610054324.9) that is used for the gradation of image spectrum that the present invention is used has strong bottom layer image information excavating function; Visual precision with a pixel; General grey level histogram and logarithm grey level histogram can not be compared; The image of under low-light (level) very, taking; Compose all information that shown this image through the gray scale behind the gradually flattening, the scope that the clear image information of the width table of its spectrum exists, we are referred to as the gray scale spectrum continuous bandwidth of image and represent with BW.Used formula is:
In the formula, O (g), T (g) represent original image and the pixel count of target image on g gray level respectively; G=0,1,2 ... N-1, N=255, expression gray scale or gamut of chromaticities; M is a positive integer, is called the planarization level;
The present invention is applicable to that source images has the bottom layer image mining that decline type gray scale spectrum distributes; The image that decline type gray scale spectrum distributes is meant that image information comes into existence and maximum in 0 grayscale image information from 0 gray level, promptly is the situation that decline type gray scale spectrum distributes from 0 gray level.According to practical experience, the original image of under the very low condition of illuminance, taking has this gray scale spectrum to distribute more, and therefore when excavation has the best quality image of the original image that decline type gray scale spectrum distributes, excavating that starting point Theta fixedly installs is 0.
Find according to experimental result repeatedly, the image that Delta is each time excavated after changing, from the value of the excavation image information entropy InEn that calculates, along with the increase of Delta, image information entropy InEn is dull ascendant trend with the increase of Del ta.
According to experimental result discovery repeatedly; To the image of excavating after the Delta variation each time; From the image averaging logarithm contrast of excavating
value that calculates; Along with the increase of Delta, the average logarithm contrast of image
increases variation with Delta and is dull downtrending.
Find according to experimental result repeatedly: the image that Delta is each time excavated after changing; From the image quality evaluation function IQAF value of calculating; Increase along with Delta; The weighted sum of image information entropy InEn that excavates and average logarithm contrast
has maximal value with the increase variation of Delta, does not promptly have ginseng image quality evaluation function
maximal value is arranged; We are excavating image with the variation of Delta on one side in the experiment repeatedly; Calculate on one side when excavating image quality in images evaluation function IQAF value and find: change with Delta excavate image quality from poor to good; And then to the degradation of image; Consistent with image quality evaluation function IQAF value with the trend of the variation of Delta; And when the image quality evaluation function was got maximal value, it was the best in quality to utilize the picture quality that the described bottom layer image mining formula of above-mentioned steps (3) excavates and the image comparison that all are excavated with its pairing excavation parameter Delta and Theta.
The invention has the beneficial effects as follows: this method can promptly obtain the best quality image on the objective sense in the bottom layer image mining process of the source images with the distribution of decline type spectrum, the image of under the very low shooting condition of illumination, taking can obtain the optimized image of objective natural scene clearly through this method; This method is embedded existing camera arrangement can realize the shooting of general camera night scene, and obtain best night scene quality image, embed existing supervisory system and can realize round-the-clock monitoring, embed existing lll night vision system raising existing system function.
Description of drawings
Description of drawings of the present invention is following:
Fig. 1 is the process flow diagram of the inventive method;
Fig. 2 is for information entropy, logarithm contrast in the embodiment of the invention and do not have the variation diagram of ginseng image quality evaluation function with Delta;
Wherein I is that information entropy InEn is with the Delta change curve; II be average logarithm contrast
with the Delta change curve, III is that image quality evaluation function IQAF value is with the Delta change curve.
Fig. 3 is for excavating in the present embodiment best quality image process, and picture quality is with the variation of not having ginseng image quality evaluation functional value.
Wherein Fig. 3 (a) is source images and gray scale spectrum thereof; BW=19; Fig. 3 (b) is Theta=0; Image of excavating when corresponding image quality assessment function IQAF=2.2411 gets maximal value during Delta=3 and gray scale spectrum thereof; Fig. 3 (c) is Theta=0; During Delta=1, image of excavating during corresponding image quality assessment function IQAF=2.0052 and gray scale spectrum thereof.
Embodiment
Below in conjunction with accompanying drawing and embodiment the present invention is described further:
The process flow diagram of the inventive method is as shown in Figure 1, and concrete steps are following:
(1) gray scale of detection resources image spectrum distributes, and confirms to have decline type gray scale spectrum and distributes, then the gray scale of calculation sources image spectrum continuous bandwidth BW; The gray scale spectrum that shown in Fig. 3 (a), can obviously find out present embodiment is the distribution of decline type spectrum; Gray scale spectrum continuous bandwidth BW=19.
The detection method (this method is documented in the patent of invention that the patent No. is ZL200610054324.9) that is used for the gradation of image spectrum used in the step of the present invention (1) has strong bottom layer image information excavating function; Visual precision with a pixel; General grey level histogram and logarithm grey level histogram can not be compared; The image of under low-light (level) very, taking; Compose all information that shown this image through the gray scale behind the gradually flattening, the scope that the clear image information of the width table of its spectrum exists, we are referred to as the gray scale spectrum continuous bandwidth of image and represent with BW.Used formula is:
In the formula, O (g), T (g) represent original image and the pixel count of target image on g gray level respectively; G=0,1,2 ... N-1, N=255, expression gray scale or gamut of chromaticities; M is a positive integer, is called the planarization level;
(2) set up image quality evaluation function IQAF, its computing formula is:
InEn is an image information entropy;
is the average contrast, and the computing formula of image information entropy InEn is:
P in the formula (i) is the pixel count of image on i gray level, when p (i)=0, and order
Log
2p(i)=0;
M, N be image at x, the pixel count of y direction, (x is that ((x+1 y) is pixel (x+1, gray-scale value y) to G to pixel for x, gray-scale value y) y) to G;
(3) carry out bottom layer image mining; And with excavating the value of the corresponding information entropy InEn of image with average logarithm contrast
calculating evaluation function IQAF, the bottom layer image mining formula is:
Constraint condition:
TG (x, y), OG (x, y) expression is excavated image and source images at pixel (x respectively; Y) gray-scale value, Theta [0,255] and Delta [1; 255] for excavating parameter, starting point and the gray-level that excavates image are excavated in expression respectively, and K is the space telescopic factor; For bottom layer image mining, get k=255, Theta=0;
The present invention is applicable to that source images has the bottom layer image mining that decline type gray scale spectrum distributes; The image that decline type gray scale spectrum distributes is meant that image information comes into existence and maximum in 0 grayscale image information from 0 gray level, promptly is the situation that decline type gray scale spectrum distributes from 0 gray level.According to practical experience, the original image of under the very low condition of illuminance, taking has this gray scale spectrum to distribute more, and therefore when excavation had the best quality image of the original image that decline type gray scale spectrum distributes, it was 0 that the initial value of Theta fixedly installs.
Make Theta=0; Delta increases 1 from 1 to gray scale spectrum continuous bandwidth BW at every turn changes; Source images is carried out bottom layer image mining; Use the value of information entropy InEn with average logarithm contrast
the computed image quality assessment function IQAF of the pairing excavation image of each Delta value, up to the maximal value that obtains IQAF;
According to experimental result discovery repeatedly; To the image of excavating after the Delta variation each time; Value from the excavation image information entropy InEn that calculates; Along with the increase of Delta, image information entropy InEn is dull ascendant trend with the increase of Delta, and I as shown in Figure 2 is that information entropy InEn is with the Delta change curve.
According to experimental result discovery repeatedly; To the image of excavating after the Delta variation each time; From the image averaging logarithm contrast of excavating
value that calculates; Increase along with Delta; The average logarithm contrast of image
increases variation with Delta and is dull downtrending, and II as shown in Figure 2 is that average logarithm contrast
is with the Delta change curve.
Find according to experimental result repeatedly: the image that Delta is each time excavated after changing; From the image quality evaluation function IQAF value of calculating; Increase along with the gradation of image level Delta that excavates; The weighted sum of image information entropy InEn that excavates and average logarithm contrast
has maximal value with the increase variation of Delta; Promptly not having ginseng image quality evaluation function
, maximal value, III as shown in Figure 2 are arranged is that image quality evaluation function IQAF value is with the Delta change curve;
(4) image of excavating according to the said bottom layer image mining formula of step (3) with the pairing excavation parameter of image quality evaluation function IQAF maximal value Delta, Theta is image the best in quality.
We excavate image with the variation of Delta on one side in experiment; Calculate on one side when excavating image quality in images evaluation function IQAF value and find: increase according to Delta change each width of cloth image of excavating quality from poor to good; And then to the degradation of image; Consistent with image quality evaluation function IQAF value with the trend that Delta changes; And when the image quality evaluation function is got maximal value, utilize image that the described excavation formula of above-mentioned steps (3) excavates and its visual effect the best of image comparison that all are excavated with its pairing excavation parameter Delta and Theta.
Utilize the inventive method that the source images of present embodiment shown in Fig. 3 (a) carried out best quality image that bottom layer image mining obtained shown in Fig. 3 (b).
Following table is when excavating the present embodiment best quality image, and image information entropy InEn, average logarithm contrast
image quality evaluation function IQAF are with the situation of change of excavating parameter Delta.
Can find out from the experimental data of last table; Work as Theta=0; During Delta=3; The image quality evaluation function of present embodiment is to there being maximal value, and the (a) and (b) in the comparison diagram 3, (c) can find out: the picture quality of utilizing the corresponding excavation parameter Delta of this maximal value to excavate is best shown in Fig. 3 (b).
For the inventive method clearly is described; Last tabular has gone out present embodiment Delta when l to 9 changes; The situation of change of image quality evaluation function, information entropy, average logarithm contrast; As long as but calculating Delta in practice when getting a certain value, the value of image quality evaluation function begins to diminish, and promptly stops to the excavation of image with to the calculating of image quality evaluation functional value; And utilizing corresponding excavation parameter Delta of image quality evaluation function maximal value and Theta to excavate image, the image that is obtained is image the best in quality.For example work as in the present embodiment and calculate Delta=4; During IQAF=2.1656, the value of finding image quality evaluation function IQAF is less than working as Delta=3, during the value of IQAF=2.2411; Promptly stop the excavation of image and the calculating of image quality evaluation functional value; And the image when utilizing the described excavation formula of this method to excavate Delta=3, Theta=0, the image of being excavated is picture quality the best, shown in Fig. 3 (b).
The image information entropy InEn that can find out present embodiment from the experimental data of last table is varied to dull ascendant trend with the increase of Delta; The average logarithm contrast
of excavation image is dull downtrending with the increase variation of Delta, and image quality evaluation function IQAF has maximal value with the increase variation of Delta.In the present embodiment, when Delta=3, image quality evaluation function IQAF=2.2411 is a maximal value, utilizes the corresponding excavation parameter Theta=0 of this maximal value, and the image that Delta=3 excavates is shown in Fig. 3 (b).Corresponding excavation parameter Delta=1 for example when relatively the source images of present embodiment is got non-maximal value with image quality evaluation function IQAF shown in Fig. 3 (a); The image of excavating during Theta=0 can be found out shown in Fig. 3 (c): it is the best in quality that image quality evaluation function IQAF gets image such as Fig. 3 (b) that the maximal value correspondence excavates.
Camera of embedding the inventive method and video camera can be implemented in shoots picture rich in detail night.
The digital telescope that embeds the inventive method can be realized night-time observation.
Claims (1)
1. method based on the excavating optimum image of information entropy and logarithm contrast weighted sum is characterized in that may further comprise the steps:
(1) gray scale of detection resources image spectrum distributes, and confirms to have decline type gray scale spectrum and distributes, then the gray scale of calculation sources image spectrum continuous bandwidth BW;
(2) set up image quality evaluation function IQAF, its computing formula is:
InEn is an image information entropy,
Be the average contrast, the computing formula of image information entropy InEn is:
P in the formula (i) is the pixel count of image on i gray level, when p (i)=0, makes Log
2P (i)=0;
M, N be image at x, the pixel count of y direction, (x is that ((x+1 y) is pixel (x+1, gray-scale value y) to G to pixel for x, gray-scale value y) y) to G;
(3) carry out bottom layer image mining; And with excavating the value of the corresponding information entropy InEn of image with average logarithm contrast
calculating evaluation function IQAF, the bottom layer image mining formula is:
Constraint condition:
TG (x, y), OG (x, y) expression is excavated image and source images at pixel (x respectively; Y) gray-scale value, Theta [0,255] and Delta [1; 255] for excavating parameter, starting point and the gray-level that excavates image are excavated in expression respectively, and K is the space telescopic factor; For bottom layer image mining, get k=255, Theta=0;
Make Theta=0; Delta increases 1 from 1 to gray scale spectrum continuous bandwidth BW at every turn changes; Source images is carried out bottom layer image mining; With the information entropy InEn of the pairing excavation image of each Delta value and the value of average contrast
computed image quality assessment function IQAF, up to the maximal value that obtains IQAF;
(4) image of excavating according to the said bottom layer image mining formula of step (3) with the pairing excavation parameter of image quality evaluation function IQAF maximal value Delta, Theta is image the best in quality.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102331569A CN101441770B (en) | 2008-11-28 | 2008-11-28 | Method for excavating optimum image based on information entropy and logarithm contrast weight sum |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102331569A CN101441770B (en) | 2008-11-28 | 2008-11-28 | Method for excavating optimum image based on information entropy and logarithm contrast weight sum |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101441770A CN101441770A (en) | 2009-05-27 |
CN101441770B true CN101441770B (en) | 2012-03-21 |
Family
ID=40726189
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102331569A Expired - Fee Related CN101441770B (en) | 2008-11-28 | 2008-11-28 | Method for excavating optimum image based on information entropy and logarithm contrast weight sum |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101441770B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101620729B (en) * | 2009-07-31 | 2011-11-30 | 重庆医科大学 | Method for producing gray image with best quality |
CN102186054A (en) * | 2011-03-25 | 2011-09-14 | 谢丹玫 | Adaptive bottom video online mining system and method based on contrast resolution compensation |
CN102968800B (en) * | 2012-12-14 | 2016-01-20 | 宁波江丰生物信息技术有限公司 | A kind of evaluation method of image definition |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999053444A1 (en) * | 1998-04-15 | 1999-10-21 | Mayo Foundation For Medical Education And Research | Autocorrection of mr images for motion artifacts |
CN1866295A (en) * | 2006-05-29 | 2006-11-22 | 重庆医科大学 | Mining method for low-level image and image mining apparatus employing the same |
CN1885344A (en) * | 2006-05-25 | 2006-12-27 | 重庆医科大学 | High-resolution detection method for image gray scale/chromaticity information for base image mining |
CN101051351A (en) * | 2007-05-23 | 2007-10-10 | 重庆医科大学 | Image band parameter two-valued method and device using said method |
-
2008
- 2008-11-28 CN CN2008102331569A patent/CN101441770B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999053444A1 (en) * | 1998-04-15 | 1999-10-21 | Mayo Foundation For Medical Education And Research | Autocorrection of mr images for motion artifacts |
CN1885344A (en) * | 2006-05-25 | 2006-12-27 | 重庆医科大学 | High-resolution detection method for image gray scale/chromaticity information for base image mining |
CN1866295A (en) * | 2006-05-29 | 2006-11-22 | 重庆医科大学 | Mining method for low-level image and image mining apparatus employing the same |
CN101051351A (en) * | 2007-05-23 | 2007-10-10 | 重庆医科大学 | Image band parameter two-valued method and device using said method |
Non-Patent Citations (2)
Title |
---|
刘玉红,熊兴良,王志芳,王颖,李虹,李剑波,谢正祥.图像的带参二值化理论和技术及其应用.《计算机应用研究》.2008,第25卷(第1期), * |
谢正祥,刘玉红,王志芳,王颖,李虹,吴少波.图像的带参数二值化方法.《中国医学物理学杂志》.2007,第24卷(第2期), * |
Also Published As
Publication number | Publication date |
---|---|
CN101441770A (en) | 2009-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101419707B (en) | Best quality image acquiring method in bottom layer image mining | |
KR101742120B1 (en) | Apparatus and method for image processing | |
US20070104387A1 (en) | Method and apparatus for enhancing contrast | |
KR20100053287A (en) | The apparatus for enhancing image by generating and blending multiple image and method therefor | |
Hong et al. | An improved human visual system based reversible data hiding method using adaptive histogram modification | |
WO2012108094A1 (en) | Image processing device, image processing method, image processing program, and image pick-up device | |
EP2611173A1 (en) | Image processing device, image processing method, and program | |
CN101441770B (en) | Method for excavating optimum image based on information entropy and logarithm contrast weight sum | |
CN103327219A (en) | Video image contrast adjustment method and device | |
JP5648849B2 (en) | Image processing apparatus and image processing method | |
CN103324033A (en) | Exposure method and exposure system | |
EP1861844B1 (en) | Method and/or apparatus to improve the visual perception of an image displayed on a screen | |
CN102156960B (en) | Picture noise adding method | |
Wang et al. | Underwater image quality optimization: Researches, challenges, and future trends | |
CN105654424A (en) | Image adjusting proportion display method, display system and display device, and terminal | |
CN104021518A (en) | Reversible data hiding method for high dynamic range image | |
CN109801257A (en) | No reference DIBR generates image quality evaluating method | |
US9153038B2 (en) | Method for carrying out a dynamic range compression in traffic photography | |
JP6543787B2 (en) | Image processing apparatus and image processing method | |
US20110122284A1 (en) | Multiple illuminations automatic white balance digital cameras | |
EP1796375A3 (en) | Tone conversion apparatus, image display apparatus using the same, and tone conversion method | |
CN101957985B (en) | Automatic self-adaptive optimum compensation method of human vision contrast resolution | |
US20080187049A1 (en) | Moving image encoding device and moving image encoding method | |
JP5038190B2 (en) | Image capturing apparatus and setting method thereof | |
CN102186054A (en) | Adaptive bottom video online mining system and method based on contrast resolution compensation |
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: 20120321 Termination date: 20121128 |