CN110473252A - The method and module of arbitrary shape contact area are calculated based on matrix operation - Google Patents
The method and module of arbitrary shape contact area are calculated based on matrix operation Download PDFInfo
- Publication number
- CN110473252A CN110473252A CN201910804534.2A CN201910804534A CN110473252A CN 110473252 A CN110473252 A CN 110473252A CN 201910804534 A CN201910804534 A CN 201910804534A CN 110473252 A CN110473252 A CN 110473252A
- Authority
- CN
- China
- Prior art keywords
- image
- contact area
- contact
- binary map
- pixel
- 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.)
- Pending
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 50
- 239000000284 extract Substances 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims description 18
- 238000001914 filtration Methods 0.000 claims description 17
- 238000000605 extraction Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 230000001788 irregular Effects 0.000 abstract description 8
- 230000008676 import Effects 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 20
- 230000008569 process Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000004205 dimethyl polysiloxane Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 229920000435 poly(dimethylsiloxane) Polymers 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- -1 Polydimethylsiloxane, dimethyl Polymers 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003733 optic disk Anatomy 0.000 description 1
- 239000005304 optical glass Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 229920005573 silicon-containing polymer Polymers 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G06T5/70—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
Abstract
The invention discloses a kind of methods and module that arbitrary shape contact area is calculated based on matrix operation, wherein method is the following steps are included: import the contact image of shooting, and contact image is changed into binary map using maximum variance between clusters；Contact area edge, and filling cavity are extracted according to binary map, to obtain image overall region；The two values matrix of image and original image after being handled according to filling determines image contact area with matrix operation, and calculates image contact area pixel and mark, and obtains contact area.This method automatically determines binarization threshold using function, random error caused by being artificially arranged is avoided, experimental precision is improved, can accurately extract contact zone outline and calculating contact area, extractable arbitrary graphic, the profile including rule and irregular shape simultaneously calculate its area.
Description
Technical field
The present invention relates to technical field of image processing, in particular to a kind of to calculate arbitrary shape contact surface based on matrix operation
Long-pending method and module.
Background technique
Contact electrification characteristic and the real contact area of texture and substrate are closely related, test contact image and test system
For principle of uniting as shown in Figure 1, a figure is contact area image, b figure is contact image formation basic theory.When measurement, texturing PDMS
(Polydimethylsiloxane, dimethyl silicone polymer) film is contacted with quartz glass, due to entering for texture contact area
It penetrates light and contact area measuring system is reflected back with low-angle, and be not in contact region with wide-angle reflection, texture contact area
Reflective light intensity be greater than and be not in contact the reflective light intensity in region, therefore practical contact portion is the area that texture becomes clear in contact image
Domain.When carrying out image procossing calculating contact area, the relevant technologies are as follows:
(1) based on the rgb value of image, artificial threshold value, computation rule image area；(2) image processing software is utilized
The profile of contact area required by manual extraction simultaneously calculates its pixel, measures contact area by pixel ratio；(3) it is mentioned using MATLAB
Graph outline is taken, and calculates outermost contour area, is i.e. the texture gross area.
However, the relevant technologies have following defects that (1) artificial threshold value, artificial random error is introduced, and can only calculate
Regular figure area；(2) manual extraction region contour causes human error, experimental result inaccuracy；(3) image can only be calculated
Outer profile area can not accurately calculate real contact area.
Summary of the invention
The present invention is directed to solve at least some of the technical problems in related technologies.
For this purpose, an object of the present invention is to provide a kind of sides for calculating arbitrary shape contact area based on matrix operation
Method, this method can accurately extract contact zone outline and calculating contact area, and can extract arbitrary graphic, simple easily to realize.
It is another object of the present invention to propose a kind of module that arbitrary shape contact area is calculated based on matrix operation.
In order to achieve the above objectives, one aspect of the present invention embodiment is proposed one kind and is connect based on matrix operation calculating arbitrary shape
The method of contacting surface product, comprising the following steps: import the contact image of shooting, and use maximum variance between clusters by the hookup
As being changed into binary map；Contact area edge, and filling cavity are extracted according to the binary map, to obtain image overall region；
According to image and the two values matrix of described image after filling processing, image contact area is determined with matrix operation, and
It calculates described image contact area pixel and marks, obtain contact area.
The method for calculating arbitrary shape contact area based on matrix operation of the embodiment of the present invention, is automatically determined using function
Binarization threshold avoids random error caused by being artificially arranged, improves experimental precision, can accurately extract contact zone outline and count
Contact area is calculated, can extract arbitrary graphic, the profile including rule and irregular shape simultaneously calculates its area.
In addition, the method according to the above embodiment of the present invention for calculating arbitrary shape contact area based on matrix operation may be used also
With following additional technical characteristic:
Further, in one embodiment of the invention, described to use maximum variance between clusters by the contact image
It is changed into binary map, comprising: entire data are divided by two classes according to preset threshold, when the variance between the two classes is maximum
When, determine optimal threshold, and the grayscale image of the contact image is become by binary map according to the optimal threshold.
Further, in one embodiment of the invention, after being changed into the binary map, further includes: using small
The method of wave threshold denoising and median filtering carries out de-noising filtering processing to the binary map.
Further, in one embodiment of the invention, described that contact area edge, packet are extracted according to the binary map
It includes: detecting the marginal point of the binary map by gradient direction and dual-threshold voltage using canny operator, extract the weak side of image
Edge obtains the contact area edge.
Further, in one embodiment of the invention, it the calculating described image contact area pixel and marks, wraps
It includes: two values matrix the operations such as being subtracted each other to and take absolute value and determines real contact area, finds the connected component of its each contact area,
The pixel summation for determining each connected region, using pixel ratio calculating contact area.
In order to achieve the above objectives, another aspect of the present invention embodiment proposes a kind of based on matrix operation calculating arbitrary shape
The module of contact area, comprising: importing and conversion unit for importing the contact image of shooting, and use maximum between-cluster variance
The contact image is changed into binary map by method；Extraction unit for extracting contact area edge according to the binary map, and is filled out
Cavity is filled, to obtain image overall region；Computing unit, for according to the two of image after filling processing and described image
Value matrix determines image contact area with matrix operation, and calculates described image contact area pixel and mark, and is contacted
Area.
The module that arbitrary shape contact area is calculated based on matrix operation of the embodiment of the present invention, is automatically determined using function
Binarization threshold avoids random error caused by being artificially arranged, improves experimental precision, can accurately extract contact zone outline and count
Contact area is calculated, can extract arbitrary graphic, the profile including rule and irregular shape simultaneously calculates its area.
In addition, the module according to the above embodiment of the present invention for calculating arbitrary shape contact area based on matrix operation may be used also
With following additional technical characteristic:
Further, in one embodiment of the invention, the importing is further used for conversion unit according to default
Entire data are divided into two classes by threshold value, when the variance maximum between the two classes, determine optimal threshold, and according to it is described most
The grayscale image of the contact image is become binary map by good threshold value.
Further, in one embodiment of the invention, further includes: processing unit, for being changed into the two-value
After figure, de-noising filtering processing is carried out to the binary map using the method for wavelet threshold denoising and median filtering.
Further, in one embodiment of the invention, the extraction unit is further used for logical using canny operator
It crosses gradient direction and dual-threshold voltage detects the marginal point of the binary map, extract the weak edge of image, obtain the contact area
Edge.
Further, in one embodiment of the invention, the computing unit is further used for subtracting each other two values matrix
And the equal operations that take absolute value determine real contact area, find the connected component of its each contact area, determine each connected region
Pixel summation, using pixel ratio reference area.
The additional aspect of the present invention and advantage will be set forth in part in the description, and will partially become from the following description
Obviously, or practice through the invention is recognized.
Detailed description of the invention
Above-mentioned and/or additional aspect and advantage of the invention will become from the following description of the accompanying drawings of embodiments
Obviously and it is readily appreciated that, in which:
Fig. 1 is contact area image and contact area image formation basic theory schematic diagram；
Fig. 2 is the process according to the method for calculating arbitrary shape contact area based on matrix operation of the embodiment of the present invention
Figure；
Fig. 3 is the contact area test macro schematic diagram according to the embodiment of the present invention；
Fig. 4 is the stream according to the method for calculating arbitrary shape contact area based on matrix operation of one embodiment of the invention
Cheng Tu；
Fig. 5 is according to each of the method based on matrix operation calculating arbitrary shape contact area of one embodiment of the invention
A step process result schematic diagram；
Fig. 6 is the two values matrix schematic diagram according to the embodiment of the present invention；
Fig. 7 is the irregular roundness ring-shaped schematic diagram according to the embodiment of the present invention；
Fig. 8 is to be shown according to the structure of the module for calculating arbitrary shape contact area based on matrix operation of the embodiment of the present invention
It is intended to.
Specific embodiment
The embodiment of the present invention is described below in detail, examples of the embodiments are shown in the accompanying drawings, wherein from beginning to end
Same or similar label indicates same or similar element or element with the same or similar functions.Below with reference to attached
The embodiment of figure description is exemplary, it is intended to is used to explain the present invention, and is not considered as limiting the invention.
Describe to propose according to embodiments of the present invention with reference to the accompanying drawings calculates arbitrary shape contact surface based on matrix operation
Long-pending method and module, describe to propose according to embodiments of the present invention first with reference to the accompanying drawings calculates arbitrary shape based on matrix operation
The method of shape contact area.
Fig. 2 is the process of the method that arbitrary shape contact area is calculated based on matrix operation of one embodiment of the invention
Figure.
As shown in Fig. 2, should based on matrix operation calculate arbitrary shape contact area method the following steps are included:
In step s 201, the contact image of shooting is imported, and is changed into contact image using maximum variance between clusters
Binary map.
It should be noted that reading micron order texture contact image first before the contact image for importing shooting.Such as figure
Shown in 3, the embodiment of the present invention can read contact image by the measuring system of micro- texture real contact area, wherein this hair
The measuring system hardware of micro- texture real contact area of bright embodiment mainly includes that linear motor, optical glass, single-pass are micro-
Mirror, object lens, light source, high-definition camera and computer.Reading contact image, detailed process is as follows: linear motor pushes micro- knit
Structure PDMS film is contacted with quartz glass, is imaged in the sensor plane of high-definition camera, is taken the photograph through object lens and monocular microscope
As light intensity signal is changed into electrical signal by head and is uploaded to computer, to read micron order texture contact image.
Further, in one embodiment of the invention, contact image is changed into two using maximum variance between clusters
Value figure, comprising: entire data are divided by two classes according to preset threshold, when the variance maximum between the two classes, are determined most
Good threshold value, and the grayscale image of contact image is become by binary map according to optimal threshold.
Specifically, as shown in figure 4, binary conversion treatment: the contact image (such as Fig. 5 a) for importing shooting is grayscale image, in order to
Grayscale image is changed into two-value using maximum variance between clusters by subsequent image processing and reduction calculation amount, the embodiment of the present invention
Figure, wherein maximum variance between clusters, which refer to, is divided into two classes for entire data using a threshold value, when the variance between the two
When maximum, then the threshold value is exactly optimal threshold value, can automatically determine optimal threshold using this method and grayscale image is become binary map.
Further, in one embodiment of the invention, after being changed into binary map, further includes: use small echo threshold
The method of value denoising and median filtering carries out de-noising filtering processing to binary map.
Specifically, as shown in figure 4, de-noising filters: being generated due to having noise in image taking and treatment process, be
Accurate extraction contact area edge, needs de-noising to be filtered.The method that the embodiment of the present invention uses wavelet threshold denoising, makes to disappear
Image after making an uproar is more smooth and can be reduced distortion；The method that filtering uses median filtering, to eliminate speckle noise and the spiced salt
Noise (impulsive noise).
In step S202, contact area edge, and filling cavity are extracted according to binary map, to obtain image entirety area
Domain.
Wherein, in one embodiment of the invention, contact area edge is extracted according to binary map, comprising: use canny
Operator detects the marginal point of binary map by gradient direction and dual-threshold voltage, extracts the weak edge of image, obtains contact area side
Edge.
Specifically, as shown in figure 4, step S202 is specifically included:
(1) it extracts image border: in order to accurately identify characteristics of image and reference area, extracting its edge contour, the present invention
Embodiment, by gradient direction and dual-threshold voltage detected edge points, extracts the weak edge of image using canny operator.
(2) filling cavity: the image (such as Fig. 5 d) after edge extracting can not calculate its area, need filling cavity.Filling
Image (such as Fig. 5 e) afterwards is texture area, i.e. image overall region.
In step S203, the two values matrix of image and original image after being handled according to filling is determined with matrix operation
Image contact area, and calculate image contact area pixel and mark, obtain contact area.
Wherein, in one embodiment of the invention, it calculates image contact area pixel and marks, comprising: by two-value square
Battle array, which subtracts each other and takes absolute value, waits operations to determine real contact area, finds the connected component of its each contact area, determines each company
The pixel summation in logical region, using pixel ratio reference area.
Specifically, as shown in figure 4, step S203 is specifically included:
(1) contact area is determined；By comparing Fig. 5 b and Fig. 5 e it is found that white overall region includes black in Fig. 5 b in Fig. 5 e
The contact area not in contact with region and white of color, different zones color is different in bianry image, and corresponding two values matrix is just not
Together, and the element of two values matrix only has " 0 " and " 1 ", then can determine image contact area by matrix operation.In binary map
As in, white is indicated with " 1 ", " 0 " indicates black；The two values matrix of Fig. 5 b is denoted as A (such as Fig. 6 a) after two-value processing, wherein base
Bottom " 1 ", not in contact with region " 0 ", contact area " 1 ", and the two values matrix of Fig. 5 e is denoted as B (such as Fig. 6 b) after filling cavity, including
Two values matrix A and two values matrix B are subtracted each other and are taken absolutely to obtain accurate acquisition contact area by substrate " 0 ", texture " 1 "
Value, obtains two values matrix C (such as Fig. 6 c), and image is shown in Fig. 5 f, there is noise at this time, and de-noising filters again, removes noise jamming,
Finally obtain contact area (such as Fig. 5 g).
(2) calculate and contact area pixel and mark: this algorithm uses pixel ratio reference area, and there are multiple contacts in image
Region needs to find the connected component of its each contact area, determines the pixel summation of each connected region；But image procossing mistake
Cheng Zhonghui introduces systematic error, as shown in Fig. 5 c and Fig. 5 g red circle picture.Due to the pixel of error and smaller, thus select it is suitable
Pixel size filters out contact area, and calculates its area.
In addition, the method for the embodiment of the present invention can be realized under MATLAB translation and compiling environment.Detailed process is as follows for it: (1)
Import image；(2) grayscale image is become using im2bw function by binary map according to the threshold value acquired, and calculates the total pixel of image；
(3) by the binaryzation hookup (5b) of acquisition, carry out wavelet noise using wdencmp function, recall medfilt2 function into
Row median filtering eliminates speckle noise, obtains Fig. 5 c；(3) canny operator extraction image border is utilized, as a result such as Fig. 5 d；(4)
It is empty (instruction is filled into outermost edge region such as Fig. 5 e automatically) that imfill function fills edge；(5) two values matrix subtracts each other simultaneously
It takes absolute value, while de-noising filtering finds out contact area, such as Fig. 5 g；(6) pixel of regionprops function statistical regions is called
Size, ismember function eliminate the error that edge introduces, and finally select specified contact area pixel such as Fig. 5 h.
Wherein, wavelet noise uses sym4 small echo, and number of plies N is 2；Median filtering calls medfilt2 function, in the function
Neighborhood block size be 5 × 5；The pixel size that error is screened in this program module is set as 80.
To sum up, the embodiment of the present invention is handled image based on image processing techniques, automatically extracts the figure of contact portion
Shape profile, and accurately calculate its contact area.In addition, the embodiment of the present invention can also can be extended in other Experiments of Optics arbitrarily
Image (rule and irregular figure) outline identification and areal calculation.Irregular roundness ring-shaped as shown in Figure 7a, the prior art
The error interference of black outer toroid can not be eliminated automatically, and the embodiment of the present invention can eliminate automatically black outer toroid and accurately count
The area of white inner circle is calculated, as a result as shown in Figure 7b.
The method for calculating arbitrary shape contact area based on matrix operation proposed according to embodiments of the present invention, utilizes function
Binarization threshold is automatically determined, random error caused by being artificially arranged is avoided, improves experimental precision, can accurately extract contact area
Profile and calculating contact area can extract arbitrary graphic, and the profile including rule and irregular shape simultaneously calculates its area.
Arbitrary shape contact surface is calculated based on matrix operation referring next to what attached drawing description proposed according to embodiments of the present invention
Long-pending module.
Fig. 8 is that the structure of the module for calculating arbitrary shape contact area based on matrix operation of one embodiment of the invention is shown
It is intended to.
As shown in figure 8, the module 10 that should calculate arbitrary shape contact area based on matrix operation includes: that importing and transformation are single
Member 100, extraction unit 200 and computing unit 300.
Wherein, the contact image for being used to import shooting with conversion unit 100 is imported, and will be connect using maximum variance between clusters
Touching image is changed into binary map.Extraction unit 200 is used to extract contact area edge, and filling cavity according to binary map, with
To image overall region.Computing unit 300 is for the two values matrix of image and original image after being handled according to filling, with matrix
Operation determines image contact area, and calculates image contact area pixel and mark, and obtains contact area.The embodiment of the present invention
Module 10 can accurately extract contact zone outline and calculating contact area, and can extract arbitrary graphic, simple easily to realize.
Further, in one embodiment of the invention, it imports and is further used for conversion unit 100 according to default threshold
Entire data are divided into two classes by value, when the variance maximum between the two classes, determine optimal threshold, and according to optimal threshold
The grayscale image of contact image is become into binary map.
Further, in one embodiment of the invention, the module 10 of the embodiment of the present invention further include: processing unit.
Wherein, processing unit is used for after being changed into binary map, using the method for wavelet threshold denoising and median filtering to two-value
Figure carries out de-noising filtering processing.
Further, in one embodiment of the invention, extraction unit 200 is further used for logical using canny operator
The marginal point for crossing gradient direction and dual-threshold voltage detection binary map, extracts the weak edge of image, obtains contact area edge.
Further, in one embodiment of the invention, computing unit 300 is further used for subtracting each other two values matrix simultaneously
The equal operations that take absolute value determine real contact area, find the connected component of its each contact area, determine each connected region
Pixel summation, using pixel ratio reference area.
It should be noted that the aforementioned explanation to the embodiment of the method for calculating arbitrary shape contact area based on matrix operation
Illustrate the module that arbitrary shape contact area is calculated based on matrix operation for being also applied for the embodiment, details are not described herein again.
What is proposed according to embodiments of the present invention calculates the module of arbitrary shape contact area based on matrix operation, utilizes function
Binarization threshold is automatically determined, random error caused by being artificially arranged is avoided, improves experimental precision, can accurately extract contact area
Profile and calculating contact area can extract arbitrary graphic, and the profile including rule and irregular shape simultaneously calculates its area.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or spy described in conjunction with this embodiment or example
Point is included at least one embodiment or example of the invention.In the present specification, schematic expression of the above terms are not
It must be directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described can be in office
It can be combined in any suitable manner in one or N number of embodiment or example.In addition, without conflicting with each other, the skill of this field
Art personnel can tie the feature of different embodiments or examples described in this specification and different embodiments or examples
It closes and combines.
In addition, term " first ", " second " are used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance
Or implicitly indicate the quantity of indicated technical characteristic.Define " first " as a result, the feature of " second " can be expressed or
Implicitly include at least one this feature.In the description of the present invention, " N number of " is meant that at least two, such as two, three
Deng unless otherwise specifically defined.
Any process described otherwise above or method description are construed as in flow chart or herein, and expression includes
One or it is more N number of for realizing custom logic function or process the step of executable instruction code module, segment or portion
Point, and the range of the preferred embodiment of the present invention includes other realization, wherein can not press shown or discussed suitable
Sequence, including according to related function by it is basic simultaneously in the way of or in the opposite order, Lai Zhihang function, this should be of the invention
Embodiment person of ordinary skill in the field understood.
Expression or logic and/or step described otherwise above herein in flow charts, for example, being considered use
In the order list for the executable instruction for realizing logic function, may be embodied in any computer-readable medium, for
Instruction execution system, device or equipment (such as computer based system, including the system of processor or other can be held from instruction
The instruction fetch of row system, device or equipment and the system executed instruction) it uses, or combine these instruction execution systems, device or set
It is standby and use.For the purpose of this specification, " computer-readable medium ", which can be, any may include, stores, communicates, propagates or pass
Defeated program is for instruction execution system, device or equipment or the dress used in conjunction with these instruction execution systems, device or equipment
It sets.The more specific example (non-exhaustive list) of computer-readable medium include the following: being electrically connected with one or N number of wiring
Socket part (electronic device), portable computer diskette box (magnetic device), random access memory (RAM), read-only memory (ROM),
Erasable edit read-only storage (EPROM or flash memory), fiber device and portable optic disk read-only storage
(CDROM)。
It should be appreciated that each section of the invention can be realized with hardware, software, firmware or their combination.Above-mentioned
In embodiment, software that N number of step or method can be executed in memory and by suitable instruction execution system with storage
Or firmware is realized.Such as, if realized with hardware in another embodiment, following skill well known in the art can be used
Any one of art or their combination are realized: have for data-signal is realized the logic gates of logic function from
Logic circuit is dissipated, the specific integrated circuit with suitable combinational logic gate circuit, programmable gate array (PGA), scene can compile
Journey gate array (FPGA) etc..
Those skilled in the art are understood that realize all or part of step that above-described embodiment method carries
It suddenly is that relevant hardware can be instructed to complete by program, the program can store in a kind of computer-readable storage medium
In matter, which when being executed, includes the steps that one or a combination set of embodiment of the method.
It, can also be in addition, each functional unit in each embodiment of the present invention can integrate in a processing module
It is that each unit physically exists alone, can also be integrated in two or more units in a module.Above-mentioned integrated mould
Block both can take the form of hardware realization, can also be realized in the form of software function module.The integrated module is such as
Fruit is realized and when sold or used as an independent product in the form of software function module, also can store in a computer
In read/write memory medium.
Storage medium mentioned above can be read-only memory, disk or CD etc..Although having been shown and retouching above
State the embodiment of the present invention, it is to be understood that above-described embodiment is exemplary, and is not considered as limiting the invention,
Those skilled in the art can make changes, modifications, alterations, and variations to the above described embodiments within the scope of the invention.
Claims (10)
1. a kind of method for calculating arbitrary shape contact area based on matrix operation, which comprises the following steps:
The contact image of shooting is imported, and the contact image is changed into binary map using maximum variance between clusters；
Contact area edge, and filling cavity are extracted according to the binary map, to obtain image overall region；
According to the two values matrix of image and original image after filling processing, image contact area is determined with matrix operation,
And calculate described image contact area pixel and mark, obtain contact area.
2. the method according to claim 1, wherein described use maximum variance between clusters by the contact image
It is changed into binary map, comprising:
Entire data are divided into two classes according to preset threshold, when the variance maximum between the two classes, determine optimal threshold,
And the grayscale image of the contact image is become by binary map according to the optimal threshold.
3. the method according to claim 1, wherein after being changed into the binary map, further includes:
De-noising filtering processing is carried out to the binary map using the method for wavelet threshold denoising and median filtering.
4. the method according to claim 1, wherein it is described according to the binary map extract contact area edge,
Include:
The marginal point for detecting the binary map by gradient direction and dual-threshold voltage using canny operator extracts the weak side of image
Edge obtains the contact area edge.
5. the method according to claim 1, wherein the calculating described image contact area pixel and mark,
Include:
The two values matrix is subtracted each other and the operation that takes absolute value determines actual contact area, finds the interconnecting part of its each contact area
Point, the pixel summation of each connected region is determined, using pixel ratio reference area.
6. a kind of module for calculating arbitrary shape contact area based on matrix operation characterized by comprising
Importing and conversion unit for importing the contact image of shooting, and use maximum variance between clusters by the contact image
It is changed into binary map；
Extraction unit, for extracting contact area edge, and filling cavity according to the binary map, to obtain image entirety area
Domain；
Computing unit, it is true with matrix operation for image and the two values matrix of described image after being handled according to the filling
Determine image contact area, and calculate described image contact area pixel and mark, obtains contact area.
7. module according to claim 6, which is characterized in that the importing is further used for conversion unit according to default
Entire data are divided into two classes by threshold value, when the variance maximum between the two classes, determine optimal threshold, and according to it is described most
The grayscale image of the contact image is become binary map by good threshold value.
8. module according to claim 6, which is characterized in that further include:
Processing unit, for after being changed into the binary map, using the method pair of wavelet threshold denoising and median filtering
The binary map carries out de-noising filtering processing.
9. module according to claim 6, which is characterized in that the extraction unit is further used for using canny operator
The marginal point that the binary map is detected by gradient direction and dual-threshold voltage extracts the weak edge of image, obtains the contact surface
Product edge.
10. module according to claim 6, which is characterized in that the computing unit is further used for the two-value square
Battle array is subtracted each other and the operation that takes absolute value determines real contact area, finds the connected component of its each contact area, determines each connection
The pixel summation in region, using pixel ratio reference area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910804534.2A CN110473252A (en) | 2019-08-28 | 2019-08-28 | The method and module of arbitrary shape contact area are calculated based on matrix operation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910804534.2A CN110473252A (en) | 2019-08-28 | 2019-08-28 | The method and module of arbitrary shape contact area are calculated based on matrix operation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110473252A true CN110473252A (en) | 2019-11-19 |
Family
ID=68513845
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910804534.2A Pending CN110473252A (en) | 2019-08-28 | 2019-08-28 | The method and module of arbitrary shape contact area are calculated based on matrix operation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110473252A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113642437A (en) * | 2021-08-03 | 2021-11-12 | 中国地质大学（北京） | Quantitative calculation method for content and radius of different components in coal |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2902969A1 (en) * | 2012-09-27 | 2015-08-05 | FUJIFILM Corporation | Image processing device, image processing method, and image processing program |
| JPWO2016139770A1 (en) * | 2015-03-04 | 2017-12-14 | オリンパス株式会社 | Image processing device |
| CN107862677A (en) * | 2017-10-16 | 2018-03-30 | 中铁第四勘察设计院集团有限公司 | The Tunnel Lining Cracks recognition methods of thresholding algorithm and system between a kind of class based on gradient |
-
2019
- 2019-08-28 CN CN201910804534.2A patent/CN110473252A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2902969A1 (en) * | 2012-09-27 | 2015-08-05 | FUJIFILM Corporation | Image processing device, image processing method, and image processing program |
| JPWO2016139770A1 (en) * | 2015-03-04 | 2017-12-14 | オリンパス株式会社 | Image processing device |
| CN107862677A (en) * | 2017-10-16 | 2018-03-30 | 中铁第四勘察设计院集团有限公司 | The Tunnel Lining Cracks recognition methods of thresholding algorithm and system between a kind of class based on gradient |
Non-Patent Citations (2)
| Title |
|---|
| 施丽莲等: "气液两相流气泡图像的形态学分割方法", 《自动化仪表》 * |
| 江丹丹: "图像跟踪分析在海水养殖试验中的应用技术研究", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113642437A (en) * | 2021-08-03 | 2021-11-12 | 中国地质大学（北京） | Quantitative calculation method for content and radius of different components in coal |
| CN113642437B (en) * | 2021-08-03 | 2022-05-31 | 中国地质大学（北京） | Quantitative calculation method for content and radius of different components in coal |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105517677B (en) | The post-processing approach and device of depth map/disparity map | |
| Trujillo-Pino et al. | Accurate subpixel edge location based on partial area effect | |
| KR100889014B1 (en) | Extraction method of tongue region using graph-based approach | |
| CN111340789A (en) | Method, device, equipment and storage medium for identifying and quantifying eye fundus retinal blood vessels | |
| JP5926937B2 (en) | Image processing apparatus, image processing method, and image processing program | |
| Golabbakhsh et al. | Vessel‐based registration of fundus and optical coherence tomography projection images of retina using a quadratic registration model | |
| Odstrcilik et al. | Thickness related textural properties of retinal nerve fiber layer in color fundus images | |
| CN105160652A (en) | Handset casing testing apparatus and method based on computer vision | |
| US20200258300A1 (en) | Method and apparatus for generating a 3d reconstruction of an object | |
| Zimmerman-Moreno et al. | Automatic detection of specular reflections in uterine cervix images | |
| CN106991688A (en) | Human body tracing method, human body tracking device and electronic installation | |
| JP2018191087A (en) | Adhesive matter detection device and adhesive matter detection method | |
| FR3067242A1 (en) | METHOD FOR EVALUATING AN ORTHODONTIC GUTTER | |
| CN109300104A (en) | Angular-point detection method, device | |
| US11439306B2 (en) | Apparatus and method for measuring blood flow direction using a fluorophore | |
| CN107016348A (en) | With reference to the method for detecting human face of depth information, detection means and electronic installation | |
| CN110378944A (en) | Depth map processing method, device and electronic equipment | |
| Muhimmah et al. | Automated cervical cell nuclei segmentation using morphological operation and watershed transformation | |
| CN110473252A (en) | The method and module of arbitrary shape contact area are calculated based on matrix operation | |
| WO2013175390A1 (en) | Detection of disease-related retinal nerve fiber layer thinning | |
| Zhao et al. | Deep Learning‐Based Denoising in High‐Speed Portable Reflectance Confocal Microscopy | |
| US7502525B2 (en) | System and method for edge detection of an image | |
| US8824768B1 (en) | Fertilization test method and apparatus | |
| CN106991379A (en) | Human body skin recognition methods and device and electronic installation with reference to depth information | |
| CN109982011A (en) | Imaging method based on infrared structure light |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191119 |
|
| RJ01 | Rejection of invention patent application after publication |