CN108389173B - Opportunity cost-based parameter optimization method - Google Patents
Opportunity cost-based parameter optimization method Download PDFInfo
- Publication number
- CN108389173B CN108389173B CN201810248394.0A CN201810248394A CN108389173B CN 108389173 B CN108389173 B CN 108389173B CN 201810248394 A CN201810248394 A CN 201810248394A CN 108389173 B CN108389173 B CN 108389173B
- Authority
- CN
- China
- Prior art keywords
- psnr
- image
- value
- opportunity cost
- optimal
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000005457 optimization Methods 0.000 title claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 description 14
- 241000287828 Gallus gallus Species 0.000 description 5
- 235000013330 chicken meat Nutrition 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 241000272525 Anas platyrhynchos Species 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 241000282887 Suidae Species 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 241000758706 Piperaceae Species 0.000 description 1
- 241000287531 Psittacidae Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007728 cost analysis Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 244000144977 poultry Species 0.000 description 1
- 235000013594 poultry meat Nutrition 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
Abstract
The invention relates to a parameter optimization method based on opportunity cost, which comprises the steps of firstly inputting n images, then calculating PSNR values of different images under all parameter combinations, calculating the opportunity cost value of each image, then calculating the opportunity cost total value of all images, and finally selecting the parameter combination with the minimum opportunity cost total value as an optimal parameter combination. The method introduces the opportunity cost thought in economics to the image enhancement field and selects the optimal parameter combination.
Description
Technical Field
The invention relates to the field of nonlinear enhancement of images, in particular to a parameter optimization method based on opportunity cost.
Background
The opportunity cost means that when a certain choice is made, the analysis process of other possible choices has to be abandoned. When a decision is forced to be made to take either of the two, there is a cost of opportunity that is the value of the best available decision for the other. Opportunistic cost analysis methods focus on the cost of the analysis selection. In life, some opportunity costs can be measured in currency. For example, farmers may not choose to raise their chickens if they choose to raise their pigs when they have more land available, and the opportunity cost of raising their pigs is to give up the benefits of raising their chickens. But some opportunity costs are often not measured in currency, for example, to choose between reading books in a library and enjoying the enjoyment of a television show.
Opportunity cost generally refers to the greatest loss of the choice, and the opportunity cost changes with the cost of the choice, for example, when the preference or value of the rejected choice changes, the value obtained does not change. The following example is taken as a calculation example: when farmers obtain more land, if the selection of pig raising is carried out, other poultry can not be selected to be raised, and the opportunity cost of pig raising is the income of abandoning chicken raising or duck raising and the like. If pig raising can obtain 9 ten thousand yuan, chicken raising can obtain 7 ten thousand yuan, and duck raising can obtain 8 ten thousand yuan, the opportunity cost of pig raising is 8 ten thousand yuan, the opportunity cost of chicken raising is 9 ten thousand yuan, and the opportunity cost of duck raising is 9 ten thousand yuan.
In the field of information technology, for an algorithm with non-single parameters (the number of parameters is two or more), the selection of various parameter values is often a trade-off. The cost of selecting the optimal parameters is to forego the potential benefits of suboptimal selection. The theory is often used in economics and related areas. In the field of information technology, algorithms or architectures related to each link are often multi-parameter or multi-technical indexes. If an image is to be obtained after the original image is enhanced, the image G to be enhanced can be obtained by adding the corresponding high frequency component to the low frequency component of the image in the field of nonlinear enhancement of the image0If it is regarded as a low-frequency image, a corresponding higher-frequency image L can be constructed-1By superposition of G0And L-1Will obtain the ratio G0Clearer image G-1(i.e., the enhanced image) the basic technical flow is as shown in the attached figure 1 of the specification. For G0Middle high frequency component L0Performing nonlinear filtering, and extrapolating0Phase-aligned higher frequency image L-1Thereby realizing the enhancement of high-frequency components of the blurred image.
For G0Middle high frequency element L0Performing secondary derivation and quartic derivation to obtain L0And L0", so by adjusting L0And L0"" then, a more approximate ideal high-frequency component L can be obtained by the formula (3-1)-1Is/are as followsWherein k is1And k is2The method can be used to generate the high-frequency component required for amplitude regulation
Then, the obtained high frequency component can be obtainedAnd an input image G0Adding to obtain enhanced result imageNamely:
when the two-dimensional image is subjected to enhancement processing, a classical linear Laplacian high-pass filter can be used for carrying out high-frequency component L0And (5) calculating a second derivative and a fourth derivative. In the experimental simulation, the filter can be replaced by a 3 × 3Laplacian template shown in the attached fig. 2 of the specification.
In this image enhancement algorithm, the parameter k1And k is2The value of (a) will have an effect on the final enhancement effect. In the case of objectively evaluating image quality, the quality of the image is actually reflected by the size of an objective parameter related to the image, and here, a Peak Signal Noise Ratio (PSNR) value is used for description. Generally, the greater the PSNR value, the better the image will be when comparing different enhancement algorithms. The PSNR is specifically defined as follows:
wherein 2n-1 is the maximum possible pixel value in the image, n is the number of bits required for each pixel, and for an 8bit representation of the image this value is 255.
As can be seen from the formula (3-3), if gm,n≈fm,nThe PSNR value approaches infinity, and the distortion caused by the processing is small and the image quality is good.
Disclosure of Invention
In view of this, the present invention provides a method for optimizing parameters based on opportunity cost, which introduces the opportunity cost thought in economics into the image enhancement field and selects an optimal parameter combination.
The invention is realized by adopting the following scheme: a parameter optimization method based on opportunity cost specifically comprises the following steps:
step S1: given n input images, setting a parameter affecting high frequency components of the images to be k1And k is2(ii) a For each image, each group k1And k is2Can generate corresponding PSNR values by setting m groups of different k1And k is2A combination of parameters of (1);
step S2: calculate all k1And k is2Calculating the PSNR (i, j) of different images under the combination of (a) and (b); wherein i represents k1And k is2J represents the serial number of the current image, i is more than or equal to 1 and less than or equal to m, and j is more than or equal to 1 and less than or equal to n;
step S3: calculating the maximum PSNR value of each image, wherein the maximum PSNR value Optimal _ PSNR (j) of the jth image is calculated by adopting the following formula:
Optimal_PSNR(j)=Max(PSNR(1:m,j));
step S4: calculating an opportunity cost value of each image, wherein the opportunity cost value of the jth image cost (j) is calculated by the following formula:
Cost(j)=|Optimal_PSNR(j)-PSNR(i,j)|;
step S5: the total opportunity cost value Tcost for all images is calculated using the following equation:
step S6: selecting the i value corresponding to the minimum Tcost and the k corresponding to the i value1And k is2Combine and select k1And k is2The combination is the most optimal parameter combination.
Further, in step S1, m ≦ n.
Compared with the prior art, the invention has the following beneficial effects:
1. the method applies the opportunity cost in economics to the field of image enhancement, and can obtain optimal parameters in a certain sense for a multi-parameter algorithm.
2. The algorithm of the invention can realize the automatic initialization selection process of the parameters for the multi-parameter algorithm, and finally can finely adjust and optimize the parameters according to the specific algorithm.
Drawings
Fig. 1 is a basic block diagram of a non-linear image enhancement scheme in the background of the invention.
Fig. 2 is a schematic diagram of a linear high-pass Laplacian filter template M in the background art of the present invention.
Fig. 3 is a schematic diagram of the algorithm in the embodiment of the present invention.
FIG. 4 is a schematic diagram of a test image according to an embodiment of the present invention.
FIG. 5 is a schematic diagram of an application of the embodiment of the present invention.
Detailed Description
The invention is further explained below with reference to the drawings and the embodiments.
The present embodiment provides a method for parameter optimization based on opportunity cost, which is the loss caused by selecting one solution among a plurality of different solutions to meet a specific objective. For example, let X1And X2Respectively representing two images, selected for X2Will generate the opportunity cost (X) corresponding to it22) In the case of this scheme, the result is also X1Opportunity cost (X)21). Likewise, select for X1Will generate the opportunity cost (X) corresponding to it11) In the case of this scheme, the result is also X2Opportunity cost (X)12). Adding the opportunity costs of different images for the same solution can result in X for image1And X2Total opportunity cost in terms of, i.e. total opportunity cost Tcost of solution 11And the total opportunity cost Tcost of solution 22Comparing their sizes, selecting a group of smaller values to obtain the image X for enhancement1And X2In terms of an optimal solution based on opportunity cost. In the above example, the opportunity cost (X)ij) The subscripts in (a) are defined as: the first index i indicates the optimal solution for the image and the second index j indicates the corresponding image under the solution.
Given n (n > 1) input images, for each image, each group k1And k is2The parameter combination value of (2) can generate corresponding PSNR values, and the optimal parameter combination value corresponding to the maximum PSNR value can be easily found out from the PSNR values. Assuming that there are m sets of parameter combination values, it is obvious that m ≦ n, since the optimal parameter combination is the same for some images. In the m sets of parameter combinations, each image will have m opportunity cost values, the opportunity costs of different images in each set of parameter combinations are added to obtain m total opportunity cost values Tcost, the corresponding parameter combination value when the minimum value of Tcost is selected is the optimal parameter combination value, and the specific example of the algorithm is shown in fig. 3.
The method specifically comprises the following steps:
step S1: given n input images, setting a parameter affecting high frequency components of the images to be k1And k is2(ii) a For each image, each group k1And k is2Can generate corresponding PSNR values by setting m groups of different k1And k is2A combination of parameters of (1);
step S2: calculate all k1And k is2Calculating the PSNR (i, j) of different images under the combination of (a) and (b); wherein i represents k1And k is2J represents the serial number of the current image, i is more than or equal to 1 and less than or equal to m, and j is more than or equal to 1 and less than or equal to n;
step S3: calculating the maximum PSNR value of each image, wherein the maximum PSNR value Optimal _ PSNR (j) of the jth image is calculated by adopting the following formula:
Optimal_PSNR(j)=Max(PSNR(1:m,j));
step S4: calculating an opportunity cost value of each image, wherein the opportunity cost value of the jth image cost (j) is calculated by the following formula:
Cost(j)=|Optimal_PSNR(j)-PSNR(i,j)|;
step S5: the total opportunity cost value Tcost for all images is calculated using the following equation:
step S6: selecting the i value corresponding to the minimum Tcost and the k corresponding to the i value1And k is2Combine and select k1And k is2The combination is the most optimal parameter combination.
In the present embodiment, m ≦ n in step S1.
Since the optimal parameter combination is difficult to obtain intuitively, we define k empirically1And k is2The parameter value range of (2). In simulation experiments, k1And k is2Respectively has a value of k10.05,0.10,0.15,. 1 and k20.1,0.2, 0.3. 12 standard gray scale images with the size of 512 x 512 are selected, as shown in fig. 4.
TABLE 1 PSNR (dB) value analysis under different parameters of different images
TABLE 2 different k1、k2Opportunity cost and total cost corresponding to parameter combination
X1 to X2 in Table 1 and Table 2 correspond to images Aeral, Airplane, Bike, Couple, Crowd, FlintStones, Elaine, Lena, Parrot, Peppers, Stream and bridge, and Tank, respectively. Table 1 shows the PSNR values of the 12 images under the respective optimal parameter combination valuesSo as to see that the total number of the 12 selected images has 3 groups of optimal parameter combination value options, which are respectively k10.05 and k2=0.3、k10.05 and k20.4 and k10.05 and k20.5. Table 2 shows the differences of k1And k is2The opportunity cost and the total cost corresponding to the optimal parameter combination can be obviously seen10.05 and k2The total opportunity cost value Tcost corresponding to 0.4 is the minimum, so the set of parameter combinations is optimal.
TABLE 3 PSNR (dB) Performance comparison of the present enhancement algorithm
For the algorithm of the present embodiment, when using the opportunistic cost parameter optimization algorithm, it can be seen from table 3 that k is1=0.05,k20.4 is the optimum parameter combination for all test images, and the total opportunity cost value Tcost corresponding to the combination is not only the minimum (see table 2), but also the PSNR value for the test images is excellent.
The embodiment can be applied to the field of nonlinear enhancement algorithms, as shown in fig. 5.
The above description is only a preferred embodiment of the present invention, and all equivalent changes and modifications made in accordance with the claims of the present invention should be covered by the present invention.
Claims (1)
1. A parameter optimization method based on opportunity cost is characterized in that: the method comprises the following steps:
step S1: given n input images, setting a parameter affecting high frequency components of the images to be k1And k is2(ii) a For each image, each group k1And k is2Can generate corresponding PSNR values by setting m groups of different k1And k is2A combination of parameters of (1);
step S2: calculate all k1And k is2Calculating the PSNR (i, j) of different images under the combination of (a) and (b); wherein i represents k1And k is2J represents the serial number of the current image, i is more than or equal to 1 and less than or equal to m, and j is more than or equal to 1 and less than or equal to n;
step S3: calculating the maximum PSNR value of each image, wherein the maximum PSNR value Optimal _ PSNR (j) of the jth image is calculated by adopting the following formula:
Optimal_PSNR(j)=Max(PSNR(1:m,j));
step S4: calculating an opportunity cost value of each image, wherein the opportunity cost value of the jth image cost (j) is calculated by the following formula:
Cost(j)=|Optimal_PSNR(j)-PSNR(i,j)|;
step S5: the total opportunity cost value Tcost for all images is calculated using the following equation:
step S6: selecting the i value corresponding to the minimum Tcost and the k corresponding to the i value1And k is2Combine and select k1And k is2Combining the most optimal parameter combination;
in step S1, m is less than or equal to n.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810248394.0A CN108389173B (en) | 2018-03-24 | 2018-03-24 | Opportunity cost-based parameter optimization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810248394.0A CN108389173B (en) | 2018-03-24 | 2018-03-24 | Opportunity cost-based parameter optimization method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108389173A CN108389173A (en) | 2018-08-10 |
CN108389173B true CN108389173B (en) | 2021-08-31 |
Family
ID=63066788
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810248394.0A Active CN108389173B (en) | 2018-03-24 | 2018-03-24 | Opportunity cost-based parameter optimization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108389173B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104200434A (en) * | 2014-08-28 | 2014-12-10 | 哈尔滨工程大学 | Non-local mean image denoising method based on noise variance estimation |
CN104392446A (en) * | 2014-11-20 | 2015-03-04 | 江南大学 | Improved PSNR (Peak Signal to Noise Ratio)-based DCT (Discrete Cosine Transformation) domain non-reference blurred image quality evaluation method |
CN104598987A (en) * | 2014-12-16 | 2015-05-06 | 南京华苏科技股份有限公司 | Method for predicating offline tendency and probability of mobile subscriber through learning and network effects of social networking services |
US9456086B1 (en) * | 2003-03-07 | 2016-09-27 | Wai Wu | Method and system for matching entities in an auction |
CN106972488A (en) * | 2017-05-12 | 2017-07-21 | 国网江西省电力公司经济技术研究院 | A kind of life cycle management programming screening method for considering essential safety |
CN107563995A (en) * | 2017-08-14 | 2018-01-09 | 华南理工大学 | A kind of confrontation network method of more arbiter error-duration models |
-
2018
- 2018-03-24 CN CN201810248394.0A patent/CN108389173B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9456086B1 (en) * | 2003-03-07 | 2016-09-27 | Wai Wu | Method and system for matching entities in an auction |
CN104200434A (en) * | 2014-08-28 | 2014-12-10 | 哈尔滨工程大学 | Non-local mean image denoising method based on noise variance estimation |
CN104392446A (en) * | 2014-11-20 | 2015-03-04 | 江南大学 | Improved PSNR (Peak Signal to Noise Ratio)-based DCT (Discrete Cosine Transformation) domain non-reference blurred image quality evaluation method |
CN104598987A (en) * | 2014-12-16 | 2015-05-06 | 南京华苏科技股份有限公司 | Method for predicating offline tendency and probability of mobile subscriber through learning and network effects of social networking services |
CN106972488A (en) * | 2017-05-12 | 2017-07-21 | 国网江西省电力公司经济技术研究院 | A kind of life cycle management programming screening method for considering essential safety |
CN107563995A (en) * | 2017-08-14 | 2018-01-09 | 华南理工大学 | A kind of confrontation network method of more arbiter error-duration models |
Non-Patent Citations (3)
Title |
---|
Opportunity Costs and Non-Scale Free Capabilities: Profit Maximization, Corporate Scope, and Profit margins;Daniel A等;《University of Pennsylvania Scholarly Commons》;20090830;1-45 * |
农村民间借贷的形成机制和外部性影响分析及监管;简佩茹;《财政监督》;20080315(第6期);73-74 * |
基于运动矢量相位角和卷积码的大容量视频隐写算法;杨鹏;《计算机应用》;20110401;第31卷(第4期);960-962、1009 * |
Also Published As
Publication number | Publication date |
---|---|
CN108389173A (en) | 2018-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108932697B (en) | Distortion removing method and device for distorted image and electronic equipment | |
US8457398B2 (en) | Image enhancement method and apparatuses utilizing the same | |
CN108305240B (en) | Image quality detection method and device | |
CN109784358B (en) | No-reference image quality evaluation method integrating artificial features and depth features | |
CN107274379B (en) | Image quality evaluation method and system | |
US7898693B2 (en) | Fast generation of dither matrix | |
CN109410149B (en) | CNN denoising method based on parallel feature extraction | |
CN105574901B (en) | A kind of general non-reference picture quality appraisement method based on local contrast pattern | |
De Fauw et al. | Hierarchical autoregressive image models with auxiliary decoders | |
CN106169181A (en) | A kind of image processing method and system | |
CN108615231B (en) | All-reference image quality objective evaluation method based on neural network learning fusion | |
CN108596890B (en) | Full-reference image quality objective evaluation method based on vision measurement rate adaptive fusion | |
CN105894522B (en) | A kind of more distortion objective evaluation method for quality of stereo images | |
CN102790885A (en) | Image processing apparatus and image processing method, learning apparatus and learning method, program, and recording medium | |
Gautam et al. | Efficient color image contrast enhancement using range limited bi-histogram equalization with adaptive gamma correction | |
CN114842216A (en) | Indoor RGB-D image semantic segmentation method based on wavelet transformation | |
CN104616252B (en) | Digital image enhancement method based on NSCT and PCNN | |
CN111105357A (en) | Distortion removing method and device for distorted image and electronic equipment | |
CN111986275A (en) | Inverse halftoning method for multi-modal halftone image | |
CN104240208A (en) | Uncooled infrared focal plane detector image detail enhancement method | |
CN108648180B (en) | Full-reference image quality objective evaluation method based on visual multi-feature depth fusion processing | |
CN108550152B (en) | Full-reference image quality objective evaluation method based on depth feature perception inference | |
CN105635523B (en) | Image processing method, image processing apparatus and image forming apparatus | |
CN108389173B (en) | Opportunity cost-based parameter optimization method | |
CN112767385B (en) | No-reference image quality evaluation method based on significance strategy and feature fusion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |