CN109408888A - A kind of roughness calculation method, computer readable storage medium and the terminal on two-dimensional cutting surface - Google Patents
A kind of roughness calculation method, computer readable storage medium and the terminal on two-dimensional cutting surface Download PDFInfo
- Publication number
- CN109408888A CN109408888A CN201811106196.7A CN201811106196A CN109408888A CN 109408888 A CN109408888 A CN 109408888A CN 201811106196 A CN201811106196 A CN 201811106196A CN 109408888 A CN109408888 A CN 109408888A
- Authority
- CN
- China
- Prior art keywords
- contour
- dimensional cutting
- cutting surface
- profile
- calculation method
- 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.)
- Granted
Links
- 238000005520 cutting process Methods 0.000 title claims abstract description 42
- 238000004364 calculation method Methods 0.000 title claims abstract description 24
- 238000003860 storage Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 26
- 230000003746 surface roughness Effects 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 13
- 230000008569 process Effects 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000000605 extraction Methods 0.000 claims abstract description 4
- 238000004590 computer program Methods 0.000 claims description 11
- 238000010586 diagram Methods 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 5
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 abstract description 10
- 229910002804 graphite Inorganic materials 0.000 abstract description 10
- 239000010439 graphite Substances 0.000 abstract description 10
- 238000004088 simulation Methods 0.000 abstract description 6
- 239000000284 extract Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000007770 graphite material Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Image Analysis (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种二维切削表面的粗糙度计算方法、计算机可读存储介质及终端，涉及石墨切削加工中表面粗糙度计算的技术领域，其方法包括以下步骤：对二维切削表面轮廓的图像进行边缘提取，通过滤波将轮廓边缘降噪平滑贴合，将轮廓边缘曲线数字化，轮廓中线与轮廓边缘曲线拟合，计算表面粗糙度。在将轮廓边缘曲线数字化的过程，先将轮廓边缘绘制在极坐标系上，进行数值化，再选取轮廓范围，最后将极坐标系转化为直角坐标系。本发明将轮廓曲线绘制在极坐标上，便于在360°范围内选取轮廓范围；将仿真图像的轮廓提取出来，再进行数值化，更便于在数学角度上的处理；本发明提供了一种高效、准确而低成本的二维切削表面的粗糙度计算方法。
A method for calculating the roughness of a two-dimensional cutting surface, a computer-readable storage medium and a terminal, which relate to the technical field of surface roughness calculation in graphite cutting, and the method includes the following steps: performing edge extraction on an image of a two-dimensional cutting surface profile , the contour edge is denoised and fitted smoothly through filtering, the contour edge curve is digitized, the contour center line is fitted with the contour edge curve, and the surface roughness is calculated. In the process of digitizing the contour edge curve, the contour edge is first drawn on the polar coordinate system, digitized, and then the contour range is selected, and finally the polar coordinate system is converted into a rectangular coordinate system. The invention draws the contour curve on polar coordinates, which is convenient to select the contour range within the range of 360°; the contour of the simulation image is extracted, and then digitized, which is more convenient for processing from a mathematical point of view; the invention provides an efficient , Accurate and low-cost two-dimensional cutting surface roughness calculation method.
Description
Technical field
The present invention relates to the fields that surface roughness in graphite machining calculates technology, and in particular to a kind of two-dimensional cutting
Roughness calculation method, computer readable storage medium and the terminal on surface.
Background technique
Graphite material is easy in the cutting process because local plastic deformation occurs, and cause the increase of surface roughness,
And then it is used for performance and generates certain negative effect.Therefore there is certain reality to the research of graphite machining roughness
Border meaning, existing research trend are processed piece surface using different cutting parameters, then in adopting experimentally
Roughness concentration is carried out to the piece surface after processing using measuring instrument, to obtain on the basis of many experiments optimal
Technological parameter.But cutting process is influenced by many variable factors simultaneously, this is excellent to graphite cutting technology parameter
Change experiment out tape and carry out very big difficulty, while this experimental method also needs to expend a large amount of money and time cost.
Summary of the invention
It is an object of the invention to provide a kind of efficient, accurate and inexpensive two dimension to overcome the shortcomings of existing technologies
Roughness calculation method, computer readable storage medium and the terminal of cutting surface, the purpose of the present invention pass through following technical side
Case is realized:
A kind of roughness calculation method on two-dimensional cutting surface, method the following steps are included:
(1) image on two-dimensional cutting surface is subjected to contour edge extraction；
(2) contour edge noise reduction is smoothly bonded by filtering；
(3) contour edge curve is digitized；
(4) profile middle line and contour edge curve matching；
(5) gauging surface roughness；
Wherein, by the digitized process of contour edge curve, first contour edge is plotted in polar coordinate system, carries out numerical value
Change, then choose profile, finally converts rectangular coordinate system for polar coordinate system.
Preferably, the process at two-dimensional cutting surface profile edge is extracted, first reading emulating image colour profile diagram, then will
Colored profile diagram is converted into profile grayscale image, then converts binary map for profile grayscale image, and segmentation contour finally extracts profile
Edge.
Preferably, contour edge process is extracted using the edge function in software Matlab, canny operator.
Preferably, the filtering of contour edge noise reduction is used into median filter method.
Preferably, profile middle line and contour edge curve fitting process, by the profile side of profile middle line and selected range
Edge curve is drawn in the same rectangular co-ordinate and fastens.
Preferably, profile middle line in Matlab by using least square method to determine its position.
Preferably, gauging surface roughness process asks the mean absolute error i.e. surface of contour curve and profile middle line thick
Rugosity.
Preferably, filtering carries out intermediate value filter according to the size of the degree selection filter window of noise jamming in data
Wave, noise jamming is strong, then the filter window selected is also big；Noise jamming is weak, then filter window is small.
Preferably, profile grayscale image is before carrying out binaryzation, and degree of comparing stretches enhancing processing, to improve image clearly
Degree,
And then it promotes roughness and calculates reliability.
A kind of computer readable storage medium is stored with computer program, when the computer program is executed by processor
The step of realizing the roughness calculation method on above-described two-dimensional cutting surface.
A kind of terminal, comprising: memory, processor and be stored on the memory and can run on the processor
Computer program, when the computer program is executed by the processor realize the process described above the step of.
Beneficial effects of the present invention: the contour curve after median filtering is plotted on polar coordinates by the present invention, is convenient for
Profile is chosen within the scope of 360 °；Contours extract in graphite cutting process simulation picture is come out, then is quantized, more
Convenient for the processing on mathematical angle；Contour curve is filtered, exceptional data point is removed, keeps calculated result more reliable
Property；
The roughness calculation method on this two-dimensional cutting surface provided by the invention, the method for passing through field experiment with the prior art
It compares, the influence of factor can not be surveyed by avoiding many other of cutting process, by directly to stabilization, reliable cutting simulation image
Roughness calculate, without devoting considerable time, energy and money under the premise of provide it is a kind of efficient, accurate and inexpensive
Graphite Cutting Process parameter optimization approach.
Detailed description of the invention
Invention is described further using attached drawing, but the embodiments in the accompanying drawings do not constitute any limitation to the present invention,
For those of ordinary skill in the art, without creative efforts, it can also be obtained according to the following drawings
Its attached drawing.
Fig. 1 is a kind of flow chart of the roughness calculation method on two-dimensional cutting surface.
Fig. 2 is graphite two-dimensional cutting grayscale image.
Fig. 3 is segmentation contour binary map.
Fig. 4 is contour edge extraction figure.
Fig. 5 is contour edge numeralization polar diagram.
Fig. 6 is that surface roughness calculates rectangular coordinate system figure.
Specific embodiment
The invention will be further described with the following Examples, a kind of roughness on two-dimensional cutting surface of the present embodiment
Calculation method calculates software by a emulating image surface roughness based on Matlab exploitation to realize.The present embodiment is answered
It is easy to user with GUI in Maltab R2015a (Graphical Us-er Interfaces, GUI) exploitation Environment Design to hand over
Mutual graphical interfaces, and the call back function by writing control, achievable graphite cutting simulation getImage is located in advance in interface
Reason, image outline extracts and digitlization, system initialization, surface roughness calculate and all operations of the detection process such as output.
And calculated result is compared with roughmeter (TR100S) actual measured results, verify the feasibility of calculation method.
Software main interface mainly includes " Loading Image ", and " calculating roughness " and " initialization " (system calibrating) function are pressed
Button." Loading Image " button, for being loaded into the emulation picture for needing gauging surface roughness；" calculating roughness " button,
It is shown on interface for calculating the surface roughness of the picture for the system that is already loaded into, and by the picture being related in the process, most
Calculated result is shown on interface afterwards.The wherein image polar coordinates range of the adjustable load of Xmin and Xmax button, it is convenient
The data of imaging surface profile are extracted, Xmin value is inputted and Xmax value is used to select the range of contour curve in polar coordinates, generally
Select top profile.Complete to click " calculate roughness " button after image load and Coordinate Adjusting, system background program automatically into
Row operation original image is converted into profile diagram, median filtering, contour curve digitlization, generates evaluation reference line, and is shown by image
Out, last to obtain calculated result.Wherein, the original image display box on interface, for showing the emulation picture of loaded system,
It is only able to display one, if the emulation picture of change loading system, this display box can also follow automatic change；Profile diagram display box,
For showing the profile successfully extracted；Profile polar diagram display box, is plotted on polar coordinates with by the profile of specified range
And it here it is shown that is shown in frame；Profile rectangular plots display box is plotted in rectangular co-ordinate with by the profile of specified range
Above and here it is shown that being shown in frame.
As shown in Figure 1, its calculation method the following steps are included:
(1) emulating image colour profile diagram is read first；As shown in Fig. 2, converting profile grayscale image for the colored profile diagram；
As shown in figure 3, converting binary map for the profile grayscale image again, and profile in the binary map is split；Such as Fig. 4 institute
Show, extracts contour edge；
(2) the contour edge noise reduction is smoothly bonded by filtering, removes exceptional data point, keeps calculated result more reliable
Property；
(3) it as shown in figure 5, the contour edge is plotted in polar coordinate system, quantizes, chooses profile, it will
The polar coordinate system is converted into rectangular coordinate system；
(4) as shown in fig. 6, the contour edge curve of profile middle line and the selected range is drawn in the same rectangular coordinate system
On；
(5) the contour edge curve of the selected range and the mean absolute error of profile middle line are calculated.
Preferably, profile grayscale image is before carrying out binaryzation, and degree of comparing stretches enhancing processing, to improve image clearly
Degree,
And then it promotes roughness and calculates reliability.
Preferably, the contour edge noise reduction is smoothly bonded using the method for median filtering.
Preferably, filtering carries out intermediate value filter according to the size of the degree selection filter window of noise jamming in data
Wave, noise jamming is strong, then the filter window selected is also big；Noise jamming is weak, then filter window is small.
Preferably, the contour edge process of extracting is using the edge function in software Matlab, canny operator.
Preferably, the profile middle line in Matlab by using least square method to determine its position.
Mean absolute error, that is, surface roughness of gained contour curve and least squares.
The present embodiment comes out the contours extract in graphite cutting process simulation picture, then quantizes, and is more convenient for
Processing on mathematical angle；Contour curve after median filtering is plotted on polar coordinates, convenient for choosing wheel within the scope of 360 °
Wide range；
The roughness calculation method on this two-dimensional cutting surface provided in this embodiment, avoid cutting process it is many other can not
The influence of survey factor is calculated by the roughness directly to stabilization, reliable cutting simulation image, without when putting into a large amount of
Between, a kind of efficient, accurate and inexpensive graphite Cutting Process parameter optimization approach is provided under the premise of energy and money.
The present embodiment also provides a kind of computer readable storage medium, is stored with computer program, the computer program
The step of roughness calculation method on above-described two-dimensional cutting surface is realized when being executed by processor.
The present embodiment also provides a kind of terminal, comprising: memory, processor and is stored on the memory and can be in institute
The computer program run on processor is stated, the computer program realizes the process described above when being executed by the processor
The step of.
Finally it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than the present invention is protected
The limitation of range is protected, although explaining in detail referring to preferred embodiment to the present invention, those skilled in the art are answered
Work as understanding, it can be with modification or equivalent replacement of the technical solution of the present invention are made, without departing from the reality of technical solution of the present invention
Matter and range.
Claims (10)
1. a kind of roughness calculation method on two-dimensional cutting surface, it is characterised in that: the described method comprises the following steps:
Contour edge extraction is carried out to the image on two-dimensional cutting surface；
The contour edge noise reduction is smoothly bonded by filtering；
The contour edge curve is digitized；
Profile middle line and the contour edge curve matching；
Gauging surface roughness；
Wherein, the step (3) includes following two sub-step:
Step A: first the contour edge is plotted in polar coordinate system, is quantized；
Step B: then choosing the profile of the contour edge, finally converts rectangular coordinate system for the polar coordinate system.
2. a kind of roughness calculation method on two-dimensional cutting surface as described in claim 1, it is characterised in that: the step
(1) include following sub-step:
Step C: reading emulating image colour profile diagram first；
Step D: profile grayscale image is converted by the colored profile diagram again；
Step E: and then binary map, segmentation contour are converted by the profile grayscale image；
Step F: contour edge is finally extracted.
3. a kind of roughness calculation method on two-dimensional cutting surface as claimed in claim 2, it is characterised in that: the step E
In, for profile grayscale image before carrying out binaryzation, degree of comparing stretches enhancing processing.
4. a kind of roughness calculation method on two-dimensional cutting surface as described in claim 1, it is characterised in that: the step
(1) in, contour edge process is extracted using the edge function in software Matlab, canny operator；In the step (2), it will take turns
The filtering of wide edge noise reduction uses median filter method.
5. a kind of roughness calculation method on two-dimensional cutting surface as claimed in claim 4, it is characterised in that: the step
(2) in filtering, the degree of data noise interference and the size positive correlation of filter window.
6. a kind of roughness calculation method on two-dimensional cutting two-dimensional cutting surface as described in claim 1, it is characterised in that: institute
Stating step (4) is specifically: the contour edge curve of profile middle line and the selected range is drawn in the same rectangular coordinate system
On.
7. a kind of roughness calculation method on two-dimensional cutting surface as described in claim 1, it is characterised in that: the step
(4) the profile middle line in is by using least square method to determine its position in Matlab.
8. a kind of roughness calculation method on two-dimensional cutting surface as described in claim 1, it is characterised in that: the step
(5) method that surface roughness is calculated in is specifically: seeking the mean absolute error i.e. rough surface of contour curve Yu profile middle line
Degree.
9. a kind of computer readable storage medium, which is characterized in that be stored with computer on the computer readable storage medium
Program realizes such as two-dimensional cutting surface described in any item of the claim 1 to 8 when the computer program is executed by processor
Roughness calculation method the step of.
10. a kind of terminal characterized by comprising memory, processor and be stored on the memory and can be at the place
The computer program run on reason device is realized when the computer program is executed by the processor as appointed in claim 1 to 8
The step of method described in one.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811106196.7A CN109408888B (en) | 2018-11-27 | 2018-11-27 | Roughness calculation method of two-dimensional cutting surface, computer-readable storage medium and terminal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811106196.7A CN109408888B (en) | 2018-11-27 | 2018-11-27 | Roughness calculation method of two-dimensional cutting surface, computer-readable storage medium and terminal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109408888A true CN109408888A (en) | 2019-03-01 |
| CN109408888B CN109408888B (en) | 2024-01-05 |
Family
ID=65466100
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811106196.7A Active CN109408888B (en) | 2018-11-27 | 2018-11-27 | Roughness calculation method of two-dimensional cutting surface, computer-readable storage medium and terminal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109408888B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023173570A1 (en) * | 2022-03-17 | 2023-09-21 | 中钢集团马鞍山矿山研究总院股份有限公司 | Method for determining interface roughness of block stone based on block stone contour |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102954756A (en) * | 2012-11-27 | 2013-03-06 | 青岛理工大学 | Method and device for predicting roughness of nano particle jet flow minimal quantity lubrication grinding surface |
| CN103136520A (en) * | 2013-03-25 | 2013-06-05 | 苏州大学 | Shape matching and target recognition method based on PCA-SC algorithm |
| CN104697476A (en) * | 2015-03-19 | 2015-06-10 | 北京时代之峰科技有限公司 | Automatic detection method and device for roughness light-incision profile curve |
| CN106767564A (en) * | 2016-11-03 | 2017-05-31 | 广东工业大学 | A kind of detection method for being applied to phone housing surface roughness |
| CN108303045A (en) * | 2018-02-01 | 2018-07-20 | 北京科技大学 | A kind of surface roughness measuring method and device |
-
2018
- 2018-11-27 CN CN201811106196.7A patent/CN109408888B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102954756A (en) * | 2012-11-27 | 2013-03-06 | 青岛理工大学 | Method and device for predicting roughness of nano particle jet flow minimal quantity lubrication grinding surface |
| CN103136520A (en) * | 2013-03-25 | 2013-06-05 | 苏州大学 | Shape matching and target recognition method based on PCA-SC algorithm |
| CN104697476A (en) * | 2015-03-19 | 2015-06-10 | 北京时代之峰科技有限公司 | Automatic detection method and device for roughness light-incision profile curve |
| CN106767564A (en) * | 2016-11-03 | 2017-05-31 | 广东工业大学 | A kind of detection method for being applied to phone housing surface roughness |
| CN108303045A (en) * | 2018-02-01 | 2018-07-20 | 北京科技大学 | A kind of surface roughness measuring method and device |
Non-Patent Citations (2)
| Title |
|---|
| 黄强等: ""基于虚拟加工的主轴跳动对表面粗糙度的作用规律分析"", 《重庆理工大学学报（自然科学）》 * |
| 黄强等: ""车床主轴轴颈精度的可视化分析方法"", 《机床与液压》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023173570A1 (en) * | 2022-03-17 | 2023-09-21 | 中钢集团马鞍山矿山研究总院股份有限公司 | Method for determining interface roughness of block stone based on block stone contour |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109408888B (en) | 2024-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Spontón et al. | A review of classic edge detectors | |
| Lindeberg | Principles for automatic scale selection | |
| Juneja et al. | Performance evaluation of edge detection techniques for images in spatial domain | |
| CN104359403B (en) | Planar part dimension measurement method based on sub-pixel edge algorithm | |
| Marie et al. | The delta medial axis: a fast and robust algorithm for filtered skeleton extraction | |
| WO2016062159A1 (en) | Image matching method and platform for testing of mobile phone applications | |
| CN107153848A (en) | Instrument image automatic identifying method based on OpenCV | |
| KR101032446B1 (en) | Apparatus and method for detecting vertices of images | |
| CN107607053B (en) | A measurement method for standing tree diameter at breast height based on machine vision and 3D reconstruction technology | |
| CN105930858A (en) | Fast high-precision geometric template matching method enabling rotation and scaling functions | |
| CN103136525B (en) | High-precision positioning method for special-shaped extended target by utilizing generalized Hough transformation | |
| Jeon et al. | Ring difference filter for fast and noise robust depth from focus | |
| CN107192716A (en) | A kind of workpiece, defect quick determination method based on contour feature | |
| CN115096206B (en) | High-precision part size measurement method based on machine vision | |
| CN111160477B (en) | An image template matching method based on feature point detection | |
| CN113705564B (en) | Pointer type instrument identification reading method | |
| Gupta et al. | Image edge detection: a review | |
| CN113159103B (en) | Image matching method, device, electronic equipment and storage medium | |
| CN106529548A (en) | Sub-pixel level multi-scale Harris corner detection algorithm | |
| Cai | Robust filtering-based thinning algorithm for pattern recognition | |
| CN111260564A (en) | Image processing method and device and computer storage medium | |
| CN109408888A (en) | A kind of roughness calculation method, computer readable storage medium and the terminal on two-dimensional cutting surface | |
| Liu et al. | Feature-preserving mesh denoising based on guided normal filtering | |
| KR102398260B1 (en) | Method and apparatus for camera calibration | |
| CN109254320A (en) | Seismic properties optimization and sand body Overlay District prediction technique based on image processing method |
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 |