2 * Copyright 2007 ZXing authors
\r
4 * Licensed under the Apache License, Version 2.0 (the "License");
\r
5 * you may not use this file except in compliance with the License.
\r
6 * You may obtain a copy of the License at
\r
8 * http://www.apache.org/licenses/LICENSE-2.0
\r
10 * Unless required by applicable law or agreed to in writing, software
\r
11 * distributed under the License is distributed on an "AS IS" BASIS,
\r
12 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
\r
13 * See the License for the specific language governing permissions and
\r
14 * limitations under the License.
\r
17 package com.google.zxing.common;
\r
19 import com.google.zxing.ReaderException;
\r
22 * <p>Encapsulates logic that estimates the optimal "black point", the luminance value
\r
23 * which is the best line between "white" and "black" in a grayscale image.</p>
\r
25 * <p>For an interesting discussion of this issue, see
\r
26 * <a href="http://webdiis.unizar.es/~neira/12082/thresholding.pdf">http://webdiis.unizar.es/~neira/12082/thresholding.pdf</a>.
\r
30 * @author dswitkin@google.com (Daniel Switkin)
\r
32 public final class BlackPointEstimator {
\r
34 private BlackPointEstimator() {
\r
38 * <p>Given an array of <em>counts</em> of luminance values (i.e. a histogram), this method
\r
39 * decides which bucket of values corresponds to the black point -- which bucket contains the
\r
40 * count of the brightest luminance values that should be considered "black".</p>
\r
42 * @param histogram an array of <em>counts</em> of luminance values
\r
43 * @return index within argument of bucket corresponding to brightest values which should be
\r
44 * considered "black"
\r
45 * @throws ReaderException if "black" and "white" appear to be very close in luminance in the image
\r
47 public static int estimate(int[] histogram) throws ReaderException{
\r
49 int numBuckets = histogram.length;
\r
50 int maxBucketCount = 0;
\r
51 // Find tallest peak in histogram
\r
53 int firstPeakSize = 0;
\r
54 for (int i = 0; i < numBuckets; i++) {
\r
55 if (histogram[i] > firstPeakSize) {
\r
57 firstPeakSize = histogram[i];
\r
59 if (histogram[i] > maxBucketCount) {
\r
60 maxBucketCount = histogram[i];
\r
64 // Find second-tallest peak -- well, another peak that is tall and not
\r
65 // so close to the first one
\r
67 int secondPeakScore = 0;
\r
68 for (int i = 0; i < numBuckets; i++) {
\r
69 int distanceToBiggest = i - firstPeak;
\r
70 // Encourage more distant second peaks by multiplying by square of distance
\r
71 int score = histogram[i] * distanceToBiggest * distanceToBiggest;
\r
72 if (score > secondPeakScore) {
\r
74 secondPeakScore = score;
\r
78 // Put firstPeak first
\r
79 if (firstPeak > secondPeak) {
\r
80 int temp = firstPeak;
\r
81 firstPeak = secondPeak;
\r
85 // Kind of aribtrary; if the two peaks are very close, then we figure there is so little
\r
86 // dynamic range in the image, that discriminating black and white is too error-prone.
\r
87 // Decoding the image/line is either pointless, or may in some cases lead to a false positive
\r
88 // for 1D formats, which are relatively lenient.
\r
89 // We arbitrarily say "close" is "<= 1/16 of the total histogram buckets apart"
\r
90 if (secondPeak - firstPeak <= numBuckets >> 4) {
\r
91 throw new ReaderException("Too little dynamic range in luminance");
\r
94 // Find a valley between them that is low and closer to the white peak
\r
95 int bestValley = secondPeak - 1;
\r
96 int bestValleyScore = -1;
\r
97 for (int i = secondPeak - 1; i > firstPeak; i--) {
\r
98 int fromFirst = i - firstPeak;
\r
99 // Favor a "valley" that is not too close to either peak -- especially not the black peak --
\r
100 // and that has a low value of course
\r
101 int score = fromFirst * fromFirst * (secondPeak - i) * (maxBucketCount - histogram[i]);
\r
102 if (score > bestValleyScore) {
\r
104 bestValleyScore = score;
\r