CN101082987B - Column diagram comparability measurement method based on average difference between windows - Google Patents
Column diagram comparability measurement method based on average difference between windows Download PDFInfo
- Publication number
- CN101082987B CN101082987B CN2007100429163A CN200710042916A CN101082987B CN 101082987 B CN101082987 B CN 101082987B CN 2007100429163 A CN2007100429163 A CN 2007100429163A CN 200710042916 A CN200710042916 A CN 200710042916A CN 101082987 B CN101082987 B CN 101082987B
- Authority
- CN
- China
- Prior art keywords
- histogram
- window
- width
- cloth
- value
- 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)
- Image Processing (AREA)
Abstract
The invention discloses a histogram similarity gauging method based on average difference among windows in the image analyzing technical domain, which comprises the following steps: calculating average difference after making the window distance as power value to powering through statisticallizing value difference between two same histogram windows and the value difference between two different distant windows; adopting the average difference as the standard to estimate the similarity of two histograms; displaying superior gauging property to traditional gauging method under deviation condition for gray histogram; providing reliable reference to index the image based on histogram.
Description
Technical field
The invention belongs to the image analysis technology field, be specifically related to a kind of histogram method for measuring similarity based on mean difference between window.
Technical background
Histogram is one of key character of digital picture, the population distribution information of characteristics such as the gray scale of reflection piece image (or wherein a part is regional), color, texture.Histogram feature is widely used in aspects [1] [2] such as image retrieval, target following.In concrete the application, the similarity between two width of cloth histograms is very important evaluation criterion.Existing histogram similarity measurement mainly is based on Euclidean distance or some criterion distance after the conversion between two width of cloth histograms, like L
1Distance, L
pDistance, χ
2Distance etc.These methods are when estimating the histogram similarity; All be that value through calculating two histogram corresponding positions is as standard; Calculate simpler; But histogram is had the situation of certain deformation (skew, compression, stretching etc.), and for example image receives factor affecting such as noise, illumination variation, and some just can not reflect its similarity well.When adopting these modules; Piece image and its receive the histogram similarity between the image behind certain noise, the illumination effect may be significantly less than the histogram similarity between this width of cloth image and another the incoherent image, and this is very disadvantageous in concrete the application.The L of some improved methods such as cumulative histogram
1The relation between the histogram diverse location has been considered in [3] [4] such as distance, but shortcoming is to have lost histogrammic local message.Other some histogram comparative approach have good effect, but on ubiquity, have certain deficiency like histogram peak matching method [5], each rank square relative method parameterized methods such as [6] of histogram under some application-specific.And to these not enough EMD that propose (EarthMover ' s Distance) method [7] [8] situation such as histogram skew can better be solved; But because coupling will solve an optimization problem between two width of cloth histograms; Complexity is higher, in the application such as image retrieval of big data quantity, comes with some shortcomings.The present invention proposes a kind of new histogram method for measuring similarity based on mean difference between window; Through calculating the value difference and the weighted sum of different windows distance between two width of cloth histograms; Obtaining between a window mean difference measures and estimates two histogrammic similaritys; Algorithm is simple, general, has solved some problems that previous method exists preferably.
Summary of the invention
The object of the present invention is to provide the histogram method for measuring similarity that a kind of computation complexity is lower, average error rate is also low.
The histogram method for measuring similarity that the present invention proposes is a kind of measure based on the window mean difference.
Introduce existing histogram similarity evaluation standard below earlier
Existing histogram similarity evaluation mainly is based on the comparison of the value at two width of cloth histogram uniform window places, like L
1Distance, L
pDistance, χ
2Distance etc.These distances and all can regard the Euclidean distance under the certain parameter as through the distance after some simple transformation (like the cosine distance etc.), more greatly then the histogram similarity is more little for distance.
Simple and be without loss of generality, consider the one dimension histogram here, and with L
1Distance, L
pDistance, χ
2Distance is an example.If two width of cloth histograms are respectively G={g (j) | j=0,1 ... R} and H={h (k) | k=0,1 ... R}, wherein j, k represent position of window in the histogram, r is that (j, k just represent gray-scale value to the maximum label of window for the grey level histogram of image; The r value is 255), g (j) expression histogram G is in the value at window j place, and h (k) expression histogram H is in the value at window k place.Two histogrammic all kinds of distance definitions are:
The greatest problem that these measures exist is exactly the relation of only having calculated between the value at two width of cloth histogram uniform window places, does not consider the correlativity between the different windows.Explain as shown in Figure 1 with a simple example.Shown the histogram distribution situation that three width of cloth are different among the figure; Wherein the Y histogram is the certain compression of X histogram process result's (for example gray level image being made compressed transform between a gray area) afterwards, and subjective judgement can draw these two histogrammic similaritys of X and Y obviously will be higher than the similarity between X and the Z.But owing to calculate L
1Distance, L
pDistance, χ
2Apart from the time only consider can draw the difference of the corresponding value in uniform window position these two histogrammic similaritys of X and Y and be less than the similarity between X and the Z.The example of another real image is following: the even number (between 0 to 254) that each point gray-scale value in the piece image is all become the minimum that is not more than former gray-scale value; And the gray-scale value of another width of cloth image each point to be first width of cloth image corresponding position gray-scale value add 1 (promptly complete be odd number); At this moment the grey level histogram of two width of cloth images is not 0 in the value of even number value part and odd number value part only respectively; Has no common factor; Therefore the above three analogous column diagrams distance of being obtained is all very big, can draw the complete dissimilar conclusion of this two width of cloth figure.And in fact this two width of cloth figure does not almost have difference, and histogram shape is also just the same, just translation a window.Certainly this problem is enlarging window size or can obtain solution to a certain degree through histogram after level and smooth, but such loss that can cause statistical information again, and still possibly have problems to the situation of some other conversion.These examples can explain because the value at each window place of histogram has certain correlativity, and the difference of only considering uniform window place value is estimated the histogram similarity and had certain limitation.
The mean difference measure is following between the histogram window that the present invention proposes:
Consider between the value at each window place of histogram correlativity, establish two width of cloth histograms and be respectively G={g (j) | j=0,1;, r} and H={h (k) | k=0,1;, r}, wherein j, k represent position of window in the histogram; R is the maximum label of window, and g (j) expression histogram G is at the statistical value at window j place, and h (k) expression histogram H is at the statistical value at window k place.If G, H satisfy the equal condition of always counting, i.e. ∑ (g (j))=∑ (h (k)), define among two width of cloth histogram G and the H between different windows apart from d=|j-k|, then mean difference tolerance is tried to achieve according to following method between histogram window:
(a) set a histogram maximized window apart from d
Max, as the loop termination sign, greater than d
MaxWindow distance all regard d as
Max
(b) get d=0, obtain n
0=∑ (min (g (i), h (i))) is promptly added up the common factor at two width of cloth histogram uniform window places.Statistical value with two each window place of width of cloth histogram all deducts corresponding common factor part then, generates the new histogram of two width of cloth, i.e. g (i) '=g (i)-min (g (i), h (i)), h (i) '=h (i)-min (g (i), h (i)).This moment this two width of cloth histogram uniform window place to correct to rare 1 be 0.
(c) get d=1; The histogram new to two width of cloth; M=0 begins from the window's position; Search m=r-d from small to large: if on the width of cloth histogram on m window and another width of cloth histogram the value of m+d window (promptly corresponding window is apart from being the window of d) all be not 0, then establishing less in two values one is p
m, the value at two window places is respectively deducted p
m, generate the histogram that two width of cloth upgrade.After accomplishing all search, the p of value will be arranged
mAdd up the statistical value n that obtains counting
1, it is defined as window apart from the histogram value difference that is at 1 o'clock.
(d) to d from 1 to d
MaxCarry out operation successively, just can obtain between two width of cloth histograms different windows apart from the value difference n of d like step (c)
1To n
Max, in conjunction with before the n that tries to achieve
0, histogrammic mean difference D between different windows is defined as:
Wherein, molecule is the weighted sum (is weights with the window) of the histogram value difference of corresponding variant window distance, and denominator is always counting in the width of cloth histogram.
D has represented mean difference between the window between two width of cloth histograms, and the similarity between more little expression two width of cloth histograms of D value is big more.
The acquiring method of bright D for instance below.Two width of cloth histogram A and B (as shown in Figure 2) always count and are 30, owing to have only 6 windows altogether, maximized window is apart from d
MaxBe made as 5.
A:(7,2,8,2,9,2)
B:(0,2,4,6,8,10)
(A, B)=(0,2,4,2,8,2), d always counts to obtain the common factor min of A, B
0=18; A, B obtain after deducting common factor respectively:
A1:(7,0,4,0,1,0)
B1:(0,0,0,4,0,8)
From left to right search for window among A1, the B1 apart from r be 1 and window place value all be not 0 point, and deduct one less in two values respectively: promptly A1, B1 become (7,0,0,0,1,0) and (0,0,0,0,0,8) at first respectively; Become (7,0,0,0,0,0) and (0,0,0,0,0,7) then respectively; D is at last always counted
1=4+1=5;
A2:(7,0,0,0,0,0)
B2:(0,0,0,0,0,7)
Increase deviation r and continue search, obtain d successively
2=0; d
3=0; d
4=0; d
5=7;
Calculate D=1.333 according to formula (4)
The comparison of several kinds of histogram method for measuring similarity performances
Chosen 10 width of cloth natural images (as shown in Figure 3) in the standard picture storehouse here, the image size is 128 pixels * 128 pixels (always count is 16384).Obtain they reach each other and separately image receive the similarity measurement between the grey level histogram after noise, grayscale shift etc. influence, with mean difference tolerance (d between the window that proposes
MaxGet 40) and classical L
1Distance, L
pDistance (p gets 2) and χ
2Distance metric compares.Ask according to (1), (2), (3) formula these three apart from the time, the distance of trying to achieve divided by total window number (256), is obtained the result's (because relative size relatively just, so the normalization operation is to not influence of result) here after the equalization.The conversion of image being carried out for emulating image influenced by noise, grayscale shift etc. comprises: 1, add the 15dB Gaussian noise; 2, the each point gray scale becomes original 90%; 3, the each point gray scale becomes original 110%.The influence of these transfer pair images itself and histogram shape thereof is not very big; No matter to the subjective feeling of image or histogram shape still from the requirement of application such as image retrieval, piece image and its histogram similarity through the image after the above-mentioned conversion generally should be less than the histogram similaritys between this image and another width of cloth image.And can find out from actual specific result, adopt these distance metrics to estimate the histogram similarity under the certain situation, can obtain opposite conclusion.
The result that (being altogether 45 groups of situation) reaches separately and four kinds of similarity measurements of employing are estimated between the image after its conversion between any two adds up and can know to 10 width of cloth images, adopts L
1Distance, L
pDistance (p gets 2) and χ
2The error rate of distance metric is higher, is respectively 7 groups, 14 groups and 18 groups.And mean difference metrology error rate is minimum between the window that employing the present invention proposes, and only goes wrong under 3 groups of situation therein, and under these 3 groups of situation, adopts other three kinds of distance metrics also all can not correctly estimate.Thus it is clear that, with the histogram method for evaluating similarity (L of classics
1Distance, L
pDistance and χ
2Distance) compare, mean difference tolerance receives factor affecting such as noise, illumination condition and takes place under the situation of histogram skew better effect is arranged between window at image.
Description of drawings
Fig. 1, histogram similarity example.
Mean difference metric calculation method example between Fig. 2, window.
Fig. 3, be used for the gray-scale map image set (yardstick: 128 pixels * 128 pixels) of comparison.
Fig. 4, the gray level image that is used for comparison and histogram thereof.First row from left to right is followed successively by this two width of cloth image and first width of cloth image and adds that 15dB Gaussian noise, gray scale become original 90% and gray scale and become the image after original 110%, and the next line correspondence position then is the grey level histogram of every width of cloth figure.
Embodiment
Contrast with the image of ' woman2 ' and ' crowd ' by name in the selected image is example (Fig. 4).Fig. 4 first row from left to right is respectively this two width of cloth image and ' woman2 ' image and adds that 15dB Gaussian noise, gray scale become original 90% and gray scale and become the image after original 110%, and the next line correspondence position then is the grey level histogram of every width of cloth figure.According to formula (1), (2), (3) and (4), adopt L
1Distance, L
pDistance (p gets 2), χ
2Mean difference tolerance (is represented d between distance and window with AD in the hereinafter table
MaxGet 40) to (a) with (b), (a) with (c), (a) with (d) with (a) compare with (e) histogram similarity, the result is as shown in table 1.
Wherein, (a) and the L (b)
1Distance (32.2) basically and (c) near (a), (a) and (d), (a) and (e) between L
1Distance (being respectively 32.3,29.3,29.1) is not easy to judge the size of similarity each other.(a) and the L (b)
pThe distance (47.5) then be significantly less than (a) and (c), (a) and (d), (a) and (e) between L
pDistance (being respectively 61.0,71.7,69.0).Similarly, (a) and the χ (b)
2The distance (10.9) be significantly less than (a) and (c), (a) and (d), (a) and (e) between χ
2Distance (being respectively 32.6,43.7,34.9).Adopt these two kinds of distance metrics can draw (a) and (b) the bigger conclusion of histogram similarity, this and actual observation to image and histogram situation all do not meet.Thus it is clear that, if utilize L
1Distance, L
pDistance (p gets 2) or χ
2The measured histogram similarity of distance is carried out image retrieval to (a), then possible errors ground (b) figure as result for retrieval, but not (c), (d) or (e).
By contrast; Adopt mean difference tolerance between the window that the present invention proposes then can draw mean difference (7.9) between (a) and histogram window (b) obviously and (c) greater than (a), (a) and (d), (a) and (e) between histogram window between mean difference (being respectively 4.9,3.6,3.4); Promptly (a) and the less conclusion of histogram similarity (b) meet truth more.
Table 1, Fig. 4 adopt distinct methods to carry out the result of histogram similarity measurement
In sum, mean difference is measured the histogram similarity between the histogram window that this paper proposes, and receives factor affecting such as noise, illumination condition at image and takes place that performance is superior to classical histogram method for measuring similarity under the situation of histogram skew.This method is calculated simple and clear, and has good versatility, is expected to be applied to the various fields based on the histogram similarity, like image retrieval, target following etc.
List of references
[1]Belongie,S.,Carson,C.,Greenspan,H.,and?Malik,J.Color?and?texture-based?imagesegmentation?using?EM?and?its?application?to?content-based?image?retrieval.IEEE?InternationalConference?on?Computer?Vision,1998,pp.675-682.
[2]Hafner,J.,Sawhney,H.S.,Equitz,W.,Flickner,M.,and?Niblack,W.Efficient?color?histogramindexing?for?quadratic?form?distance?functions.IEEE?Transactions?on?Pattern?Analysis?and?MachineIntelligence,17,(1995):729-735.
[3]Shen,H.C.and?Wong,A.K.C.Generalized?texture?representation?and?metric.Computer,Vision,Graphics,and?Image?Processing,23,(1983):187-206.
[4]Werman,M.,Peleg,S.,and?Rosenfeld,A.A?distance?metric?for?multi-dimensional?histograms.Computer,Vision,Graphics,and?Image?Processing,32,(1985):328-336.
[5]Das,M.,Riseman,E.M.,and?Draper,B.A.FOCUS:Searching?for?multi-colored?objects?in?adiverse?image?database.IEEE?Conference?on?Computer?Vision?and?Pattern?Recognition,1997,pp.756-761.
[6]Stricker,M.and?Orengo,M.Similarity?of?color?images.SPIE?Conference?on?Storage?andRetrieval?for?Image?and?Video?Databases?III,1995,pp.381-392.
[7]Rubner?Y.,Tomasi?C.,and?Guibas?L.J.The?Earth?Mover′s?Distance?as?a?metric?for?image?retrieval.International?Journal?of?Computer?Vision,40,(2000):99-121.
[8]Levina?E.,Bickel?P.?The?Earth?Mover’s?Distance?is?the?Mallows?Distance:some?insights?fromstatistics.IEEE?International?Conference?on?Computer?Vision,2001,pp.251-256.
Claims (1)
1. histogram method for measuring similarity based on mean difference between window; It is characterized in that adding up value difference and the value difference between the different distance window between two width of cloth histogram uniform windows; Obtaining with the window distance is the mean difference after weights carry out weighting; As the standard of estimating this two width of cloth histogram similarity, concrete steps are following with it:
If two width of cloth histograms are respectively G={g (j) | j=0,1 ... R} and H={h (k) | k=0,1 ... R}, wherein j, k represent position of window in the histogram, r is the maximum label of window; G (j) expression histogram G is at the statistical value at window j place, and h (k) expression histogram H is at the statistical value at window k place; If G, H satisfy the equal condition of always counting, i.e. ∑ (g (j))=∑ (h (k)), define among two width of cloth histogram G and the H between different windows apart from d=|j-k|, abbreviate the window distance as, then mean difference tolerance is tried to achieve according to following method between histogram window:
(a) set a histogram maximized window apart from d
Max, as the loop termination sign, greater than d
MaxWindow distance all regard d as
Max
(b) get d=0, add up the common factor at two width of cloth histogram uniform window places; Obtain n
0=∑ (min (g (i), h (i))); Statistical value with two each window place of width of cloth histogram all deducts corresponding common factor part then; Generate the new histogram of two width of cloth: i.e. g (i) '=g (i)-min (g (i); H (i)); H (i) '=h (i)-min (g (i), h (i)), this moment this two width of cloth histogram g (i) ' and h (i) ' uniform window place to correct to rare 1 be 0;
(c) get d=1, to the new histogram of described two width of cloth, m=0 begins from the window's position, searches m=r-d from small to large: if on the width of cloth histogram on m window and another width of cloth histogram value of m+d window all be not 0, then establishing less in two values one is p
m, the value at two window places is respectively deducted p
m, generate the histogram that two width of cloth upgrade; After accomplishing all search, the p of value will be arranged
mAdd up the statistical value n that obtains counting
1, it is defined as window apart from the histogram value difference that is at 1 o'clock;
(d) to d from 2 to d
MaxCarry out successively searching the concrete operations of m=r-d from small to large, obtain between two width of cloth histograms different windows apart from the value difference n of d as beginning from the window's position m=0 in the step (c)
2To n
Max, in conjunction with before the n that tries to achieve
0And n
1, between different windows histogrammic mean difference D be:
Wherein, molecule is the weighted sum of the histogram value difference of corresponding variant window distance, and denominator is always counting in the width of cloth histogram; Similarity between more little expression two width of cloth histograms of D value is big more.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100429163A CN101082987B (en) | 2007-06-28 | 2007-06-28 | Column diagram comparability measurement method based on average difference between windows |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100429163A CN101082987B (en) | 2007-06-28 | 2007-06-28 | Column diagram comparability measurement method based on average difference between windows |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101082987A CN101082987A (en) | 2007-12-05 |
CN101082987B true CN101082987B (en) | 2012-01-18 |
Family
ID=38912543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007100429163A Expired - Fee Related CN101082987B (en) | 2007-06-28 | 2007-06-28 | Column diagram comparability measurement method based on average difference between windows |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101082987B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108364072A (en) * | 2018-01-11 | 2018-08-03 | 华南理工大学 | A kind of transfer learning Activity recognition method based on metric |
CN108845998B (en) * | 2018-04-03 | 2022-04-29 | 南昌奇眸科技有限公司 | Trademark image retrieval and matching method |
CN108984601B (en) * | 2018-06-05 | 2020-11-03 | 北京工业职业技术学院 | Image retrieval method and system based on probability histogram area similarity |
CN109800787B (en) * | 2018-12-14 | 2020-12-29 | 西安交通大学 | Image template matching method based on relative feature distance error measurement |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1489111A (en) * | 2003-08-21 | 2004-04-14 | 上海交通大学 | Remote-sensing image mixing method based on local statistical property and colour space transformation |
CN1811820A (en) * | 2006-03-02 | 2006-08-02 | 复旦大学 | A picture-chasing algorithm based on compound correlation similarity |
-
2007
- 2007-06-28 CN CN2007100429163A patent/CN101082987B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1489111A (en) * | 2003-08-21 | 2004-04-14 | 上海交通大学 | Remote-sensing image mixing method based on local statistical property and colour space transformation |
CN1811820A (en) * | 2006-03-02 | 2006-08-02 | 复旦大学 | A picture-chasing algorithm based on compound correlation similarity |
Non-Patent Citations (2)
Title |
---|
胡敏,丁益洪,平西建.《基于属性关系直方图统计的线性纹理图像检索方法》.《计算机辅助设计与图形学学报》.2003,第15卷(第1期),48-52. * |
黄祥林,沈兰荪.《基于内容的图像检索技术研究》.《电子学报》.2002,第30卷(第7期),1065-1071. * |
Also Published As
Publication number | Publication date |
---|---|
CN101082987A (en) | 2007-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Park et al. | Multi-scale tensor voting for feature extraction from unstructured point clouds | |
Kassim et al. | A comparative study of efficient generalised Hough transform techniques | |
CN1985275B (en) | Method and system for hybrid rigid registration based on joint correspondences between scale-invariant salient region features | |
US20030146913A1 (en) | Object-correspondence identification without full volume registration | |
Dibeklioglu et al. | 3D facial landmarking under expression, pose, and occlusion variations | |
Domke | Parameter learning with truncated message-passing | |
US20030053696A1 (en) | System and method for performing edge detection in an image | |
Liu et al. | Optimal matching problem in detection and recognition performance evaluation | |
CN101082987B (en) | Column diagram comparability measurement method based on average difference between windows | |
Debayle et al. | Rigid image registration by general adaptive neighborhood matching | |
Rieger et al. | On curvature estimation of ISO surfaces in 3D gray-value images and the computation of shape descriptors | |
Saleem et al. | On visual complexity of 3D shapes | |
Pan et al. | An adaptive multifeature method for semiautomatic road extraction from high-resolution stereo mapping satellite images | |
Cheng et al. | U-Net cascaded with dilated convolution for medical image registration | |
Paulhac et al. | A solid texture database for segmentation and classification experiments | |
Hadfield et al. | Stereo reconstruction using top-down cues | |
Daudt et al. | Learning to understand earth observation images with weak and unreliable ground truth | |
CN111160300A (en) | Deep learning hyperspectral image saliency detection algorithm combined with global prior | |
Broadhurst et al. | Histogram statistics of local model-relative image regions | |
CN106709921B (en) | Color image segmentation method based on space Dirichlet mixed model | |
Hirsa et al. | Robust Rolling PCA: Managing Time Series and Multiple Dimensions | |
CN114862868A (en) | Cerebral stroke final infarct area segmentation method based on CT perfusion source data | |
CN104050486A (en) | Polarimetric SAR image classification method based on maps and Wishart distance | |
Zambanini et al. | 3D acquisition of historical coins and its application area in numismatics | |
Lloyd | Theory & Methods: Semi‐parametric estimation of ROC curves based on binomial regression modelling |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120118 Termination date: 20140628 |
|
EXPY | Termination of patent right or utility model |