Diamond compact cobalt removal boundary inclination angle measuring method, device, equipment and medium

Abstract

The invention discloses a diamond compact cobalt removal boundary inclination angle measurement method, device and equipment and a medium, and the method comprises the following steps: controlling a ray emission module to emit rays to a fixed to-be-measured diamond compact at different heights, so as to obtain a plurality of corresponding attenuation images; wherein the ray emission module and the to-be-tested diamond compact are coaxially arranged; performing threshold segmentation on each attenuation image, and extracting feature point coordinates of a target boundary surface in each attenuation image based on the ray attenuation intensity difference between a cobalt-free layer and a cobalt-free layer of the diamond compact to be detected; wherein the target boundary surface is a boundary surface between the cobalt-free layer and the cobalt-free layer; carrying out plane fitting on the coordinates of all the feature points to obtain a three-dimensional model of the target boundary surface; according to the three-dimensional model, the inclination angle of the target boundary surface is obtained, and the method has the advantage that nondestructive testing of the inclination angle of the cobalt removal boundary surface of the diamond compact can be achieved.

Description

Technical Field

[0001] This application relates to the field of diamond composite sheet testing technology, and in particular to a method, apparatus, equipment and medium for measuring the decobalt boundary tilt angle of diamond composite sheets. Background Technology

[0002] Diamond composite sheets (PDCs) are core components of high-performance cutting tools, drilling bits, and other equipment. The quality of their cobalt removal treatment directly affects their performance and service life. After cobalt removal treatment, diamond composite sheets form a layered structure consisting of a matrix, an unremoved cobalt layer, and a cobalt-removed layer connected sequentially from bottom to top (e.g., ...). Figure 1 As shown), during the cobalt removal process, the inclination angle of the boundary surface formed by the uncobalt-removed layer and the cobalt-removed layer (i.e., Figure 1 The α) is a key quality parameter. Excessive deviation of the boundary surface inclination angle can lead to stress concentration, reduced wear resistance, and even cracking during operation. Existing measurement methods mostly employ destructive testing (such as section microscopy), requiring the diamond composite sheet to be cut and ground, resulting in product scrap after testing, and are characterized by low efficiency and high cost. Some non-destructive testing solutions (such as ultrasonic testing) can only determine the presence or absence of the decobalt layer, but cannot accurately calculate the boundary surface inclination angle, lacking quantitative detection capabilities. Therefore, there is an urgent need to develop a non-destructive, accurate, and efficient method for measuring the decobalt boundary surface inclination angle of diamond composite sheets. Summary of the Invention

[0003] The main purpose of this application is to provide a method, apparatus, equipment and medium for measuring the decobalt boundary inclination angle of diamond composite sheets, aiming to solve the technical problems of existing methods for measuring the decobalt boundary inclination angle of diamond composite sheets, which require damaging the product structure, have low detection efficiency and high cost.

[0004] To achieve the above objectives, this application provides a method for measuring the dip angle of the cobalt removal boundary of a diamond composite sheet, comprising the following steps: The X-ray emission module is controlled to emit X-rays at different heights onto a fixed diamond composite sheet under test to obtain multiple corresponding attenuation images; wherein the X-ray emission module is coaxially arranged with the diamond composite sheet under test. Each attenuation image is segmented by thresholding. Based on the difference in ray attenuation intensity between the uncobalt-removed layer and the cobalt-removed layer of the diamond composite sheet under test, the coordinates of feature points of the target boundary surface in each attenuation image are extracted. The target boundary surface is the boundary surface between the uncobalt-removed layer and the cobalt-removed layer. Plane fitting is performed on the coordinates of all feature points to obtain a three-dimensional model of the target boundary surface; Based on the 3D model, obtain the inclination angle of the target boundary surface.

[0005] Optionally, threshold segmentation is performed on each attenuation image, and based on the difference in ray attenuation intensity between the uncobalt-removed and cobalt-removed layers of the diamond composite sheet under test, the coordinates of feature points on the target boundary surface in each attenuation image are extracted, including: Based on the preset radiation attenuation intensity model, the average radiation intensity I1 in the region without cobalt removal and the average radiation intensity I2 in the region with cobalt removal after radiation penetrates the diamond composite sheet under test are obtained respectively. Based on the average ray intensity I1 and the average ray intensity I2, the segmentation threshold T is obtained, T=(I1+I2) / 2; Based on the segmentation threshold T, the attenuated image is segmented into regions without cobalt layer removal and regions with cobalt layer removal; Obtain the coordinates of feature points on the target boundary surface located between the region without cobalt layer removal and the region with cobalt layer removal.

[0006] Optionally, the expression for the ray attenuation intensity model is: I=I'·e -μd (1-β)+Is·β In the formula, I is the intensity of the radiation after it penetrates the diamond composite sheet under test, I' is the intensity of the radiation at the time of incident, e is the natural constant, μ is the radiation attenuation coefficient of the material, and the radiation attenuation coefficient of the uncobalt layer is greater than that of the cobalt-free layer, d is the radiation penetration path length, Is is the average intensity of the scattered radiation when there is no sample, and β is the calibrated scattering correction coefficient.

[0007] Optionally, the calibration method for the scattering correction coefficient β includes the following steps: A standard sample identical to the diamond composite sheet to be tested was selected, and the radiation intensity I after the radiation penetrated the standard sample under the same testing conditions was obtained. 实测 Among them, the thicknesses of the uncobalt-removed layer and the cobalt-removed layer of the known standard sample are d1 and d2, respectively; The theoretical attenuation intensity of X-rays penetrating the cobalt-depleted layer and the non-cobalt-depleted layer is obtained as I. 理论 Among them, I 理论 =I0·e -(μ1d1+μ2d2) I0 is the incident intensity when the radiation has no attenuation and no scattering under ideal conditions, μ1 is the radiation attenuation coefficient without the cobalt layer, and μ2 is the radiation attenuation coefficient with the cobalt layer. Obtain the scattering correction coefficient β; where β = (I 理论 -I 实测 ) / (I 理论 -Is).

[0008] Optionally, a planar fit is performed on the coordinates of all feature points to obtain a 3D model of the target boundary surface, including: Based on the least squares method, the coordinates of the feature points ( x k ,y k , z k We perform a fitting to minimize the weighted sum of squared errors S, where S is expressed as:

[0009] In the formula, k = 1, 2, 3...m, where m is the total number of feature points. w k Let be the weight adjustment parameter for the k-th feature point, and satisfy 0 < 0. w k <1, w k The magnitude of the feature point is positively correlated with the reliability of the feature point, and A, B, and C are the coefficients of the plane equation to be solved. Based on the weighted sum of squared errors S, partial derivatives with respect to A, B, and C are taken and set equal to 0 to obtain the following system of linear equations:

[0010] Solve the linear equations to obtain the plane equation coefficients A, B, and C, and then perform plane fitting to obtain a three-dimensional model of the target boundary surface.

[0011] Optionally, let the inclination angle of the target boundary surface be θ, and the expression for θ is:

[0012] In the formula, Δ x Let Δ be the angle of inclination of the bottom surface of the diamond composite sheet under test relative to the horizontal plane in the x-direction. y The angle of inclination of the bottom surface of the diamond composite sheet under test relative to the horizontal plane in the y-direction.

[0013] Optionally, weight correction parameters w k The expression is: w k =1-γ· s k / s max

[0014] In the formula, γ is the weighting adjustment coefficient. s k The extraction error of the k-th feature point is... s max The maximum value of the extraction error for all feature points. x k邻 , yk邻 This represents the coordinates of the feature points adjacent to the current feature point.

[0015] To achieve the above objectives, this application also provides a device for measuring the dip angle of the cobalt removal boundary of a diamond composite sheet, comprising: The attenuation image acquisition module is used to control the X-ray emission module to emit X-rays at different heights onto a fixed diamond composite sheet under test in order to obtain multiple corresponding attenuation images; wherein, the X-ray emission module is arranged coaxially with the diamond composite sheet under test. The feature point acquisition module is used to perform threshold segmentation on each attenuation image. Based on the difference in ray attenuation intensity between the uncobalt-removed layer and the cobalt-removed layer of the diamond composite sheet under test, the feature point coordinates of the target boundary surface in each attenuation image are extracted; where the target boundary surface is the boundary surface between the uncobalt-removed layer and the cobalt-removed layer. The plane fitting module is used to perform plane fitting on the coordinates of all feature points to obtain a 3D model of the target boundary surface. The tilt angle calculation module is used to obtain the tilt angle of the target boundary surface based on the 3D model.

[0016] To achieve the above objectives, this application also provides a computer device, which includes a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the above-described method.

[0017] To achieve the above objectives, this application also provides a computer-readable storage medium storing a computer program, on which a processor executes the computer program to implement the above-described method.

[0018] The beneficial effects that this application can achieve are as follows: This application is based on X-ray inspection. X-rays are emitted from the diamond composite sheet under test at different heights using an X-ray emission module. This is because a single-height X-ray image can only reflect the two-dimensional contour of the boundary surface at that "slice" and cannot restore the spatial tilt attitude of the boundary surface. Multi-height imaging combined with three-dimensional fitting is a key prerequisite for accurate calculation of the tilt angle. Therefore, the attenuation images obtained by emitting X-rays from different heights can reflect the two-dimensional cross-section of the boundary surface in the corresponding Z-axis coordinate. Then, each attenuation image is thresholded. Based on the difference in X-ray attenuation intensity between the uncobalt-removed layer and the cobalt-removed layer of the diamond composite sheet under test, the coordinates of the feature points of the target boundary surface in each attenuation image can be extracted. Finally, plane fitting is performed on all feature point coordinates to obtain a three-dimensional model of the target boundary surface. Based on this three-dimensional model, the tilt angle of the target boundary surface can be accurately calculated. In summary, this application employs radiographic nondestructive testing technology, which eliminates the need to damage the diamond composite sheet, enabling full inspection of finished and semi-finished products and reducing testing costs. Furthermore, through multi-view imaging and three-dimensional fitting algorithms, the tilt angle measurement accuracy is high, and the testing process is highly automated, with a single sample testing time of ≤5 minutes. The testing efficiency is significantly better than traditional destructive testing methods. Attached Figure Description

[0019] To more clearly illustrate the technical solutions in the specific embodiments of this application or the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. In all the drawings, similar elements or parts are generally identified by similar reference numerals. In the drawings, the elements or parts are not necessarily drawn to scale.

[0020] Figure 1 This is a schematic diagram of the decobalt-free layered structure of a diamond composite sheet. Figure 2 This is a flowchart illustrating the method for measuring the decobalt removal boundary inclination of a diamond composite sheet in an embodiment of this application. Figure 3 This is a schematic diagram of the frame of the diamond composite sheet decobalt boundary inclination measuring device in an embodiment of this application.

[0021] The realization of the purpose, functional features and advantages of this application will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0022] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this application, and not all of the embodiments. Based on the embodiments of this application, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of this application.

[0023] It should be noted that if the embodiments of this application involve descriptions such as "first" or "second," these descriptions are for descriptive purposes only and should not be construed as indicating or implying their relative importance or implicitly specifying the number of technical features indicated. Therefore, features defined with "first" or "second" may explicitly or implicitly include at least one of those features. Furthermore, the technical solutions of the various embodiments can be combined with each other, but this must be based on the ability of those skilled in the art to implement them. When the combination of technical solutions is contradictory or impossible to implement, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection claimed in this application.

[0024] Example 1 Reference Figure 2 This embodiment provides a method for measuring the dip angle of the decobalt removal boundary of a diamond composite sheet, including the following steps: The X-ray emission module is controlled to emit X-rays at different heights onto a fixed diamond composite sheet under test to obtain multiple corresponding attenuation images; wherein the X-ray emission module is coaxially arranged with the diamond composite sheet under test. Each attenuation image is segmented by thresholding. Based on the difference in ray attenuation intensity between the uncobalt-removed layer and the cobalt-removed layer of the diamond composite sheet under test, the coordinates of feature points of the target boundary surface in each attenuation image are extracted. The target boundary surface is the boundary surface between the uncobalt-removed layer and the cobalt-removed layer. Plane fitting is performed on the coordinates of all feature points to obtain a three-dimensional model of the target boundary surface; Based on the 3D model, obtain the inclination angle of the target boundary surface.

[0025] In this embodiment, based on the X-ray detection method, X-rays are emitted from the diamond composite sheet under test at different heights through the X-ray emission module. This is because a single-height X-ray image can only reflect the two-dimensional contour of the boundary surface at that "slice" and cannot restore the spatial tilt attitude of the boundary surface. Multi-height imaging combined with three-dimensional fitting is a key prerequisite for accurate calculation of the tilt angle. Therefore, the attenuation images obtained by emitting X-rays from different heights can reflect the two-dimensional cross section of the boundary surface in the corresponding Z-axis coordinate. Then, each attenuation image is thresholded. Based on the difference in X-ray attenuation intensity between the uncobalt-removed layer and the cobalt-removed layer of the diamond composite sheet under test, the feature point coordinates of the target boundary surface in each attenuation image can be extracted. Finally, plane fitting is performed on all feature point coordinates to obtain a three-dimensional model of the target boundary surface. Based on this three-dimensional model, the tilt angle of the target boundary surface can be accurately calculated. In summary, this embodiment employs radiographic nondestructive testing technology, which eliminates the need to damage the diamond composite sheet, enabling full inspection of finished and semi-finished products and reducing testing costs. Furthermore, through multi-view imaging and three-dimensional fitting algorithms, the tilt angle measurement accuracy is high, and the testing process is highly automated, with a single sample testing time of ≤5 minutes. The testing efficiency is significantly better than traditional destructive testing methods.

[0026] It should be noted that, in terms of hardware, the testing equipment includes a radiation emission module, a radiation receiving and imaging module, and a sample fixing platform. The radiation emission module can emit X-rays or gamma rays, the radiation receiving and imaging module is used to acquire and receive radiation attenuation images after penetrating the diamond composite sheet, and the sample fixing platform is used to fix the diamond composite sheet to be tested.

[0027] As an optional implementation, threshold segmentation is performed on each attenuation image. Based on the difference in X-ray attenuation intensity between the uncobalt-removed and cobalt-removed layers of the diamond composite sheet under test, the coordinates of feature points on the target boundary surface in each attenuation image are extracted, including: Based on the preset radiation attenuation intensity model, the average radiation intensity I1 in the region without cobalt removal and the average radiation intensity I2 in the region with cobalt removal after radiation penetrates the diamond composite sheet under test are obtained respectively. Based on the average ray intensity I1 and the average ray intensity I2, the segmentation threshold T is obtained, T=(I1+I2) / 2; Based on the segmentation threshold T, the attenuated image is segmented into regions without cobalt layer removal and regions with cobalt layer removal; Obtain the coordinates of feature points on the target boundary surface located between the region without cobalt layer removal and the region with cobalt layer removal.

[0028] In this embodiment, when rays penetrate a material, their energy attenuates due to absorption and scattering, conforming to the Lambert-Beer law. This law is derived from the physical principle that "ray intensity decreases exponentially with penetration thickness." Therefore, the difference in attenuation capabilities between the uncobalt-free and cobalt-free layers can be utilized (the uncobalt-free layer contains cobalt, has a high atomic number, and its ray attenuation coefficient is significantly greater than that of the cobalt-free layer). The intensity difference after ray penetration divides the two regions, allowing for the extraction of boundary feature points, which is the basis for boundary localization. Based on this principle and a pre-defined ray attenuation intensity model, the average ray intensity I1 of the uncobalt-free layer region and the average ray intensity I2 of the cobalt-free layer region after ray penetration of the diamond composite sheet under test can be obtained. The average of these two values ​​is then used as the segmentation threshold T. The rationale is that the uncobalt-free layer has a strong attenuation ability, resulting in a lower intensity of the transmitted rays; while the cobalt-free layer has a weak attenuation ability, resulting in a higher intensity of the transmitted rays. The intensity distribution of the two layers exhibits a "bimodal distribution." Therefore, a segmentation threshold T based on the median can accurately divide the two regions, avoiding misjudgment of regions due to intensity fluctuations and ensuring the accuracy of feature point extraction on the boundary surface. Thus, the judgment rule is: when the intensity of image pixels I(x,y) in the attenuated image is ≤T, it is judged as the region of the uncobalt-free layer; when I(x,y)>T, it is judged as the region of the cobalt-free layer, thereby accurately extracting the feature point coordinates of the target boundary surface.

[0029] As an optional implementation, the expression for the ray attenuation intensity model is: I=I'·e -μd (1-β)+Is·β In the formula, I is the intensity of the radiation after it penetrates the diamond composite sheet under test, I' is the intensity of the radiation at the time of incident, e is the natural constant, μ is the radiation attenuation coefficient of the material, and the radiation attenuation coefficient of the uncobalt layer is greater than that of the cobalt-free layer, d is the radiation penetration path length, Is is the average intensity of the scattered radiation when there is no sample, and β is the calibrated scattering correction coefficient.

[0030] In this embodiment, the ray attenuation coefficient of the uncobalt-removed layer is set to μ1, and the ray attenuation coefficient of the cobalt-removed layer is set to μ2. Substituting the corresponding parameters into the above formula, the average ray intensity I1 of the uncobalt-removed layer region and the average ray intensity I2 of the cobalt-removed layer region can be obtained respectively. In the above formula, the ray intensity I' is determined by the parameters of the ray emission module and is a constant value. The ray intensity I can be quantized based on the pixel gray value of the attenuation image. The ray penetration path length d is the ray propagation distance in the composite sheet, which is related to the thickness of the composite sheet and the tilt angle of the boundary surface. At the same time, since scattering (Compton scattering) occurs when the ray penetrates the diamond composite sheet in actual detection, the ray intensity attenuation deviates from the ideal Lambert-Beer law. Therefore, a scattering correction parameter β (β can be characterized as the scattering intensity ratio, with a value range of 0~0.15, calibrated experimentally) is also introduced here to compensate for the intensity deviation caused by scattering, thereby improving the calculation accuracy of the ray intensity I and providing an accurate and reliable data basis for the subsequent accurate fitting of the target boundary surface.

[0031] It should be noted that the above-mentioned "no sample" refers to a measurement state in which no diamond composite sheet is placed on the sample fixing platform, and only the detection parameters such as the X-ray emission module, the receiving module, and the platform position remain unchanged.

[0032] As an optional implementation method, the calibration method for the scattering correction coefficient β includes the following steps: A standard sample identical to the diamond composite sheet to be tested was selected, and the radiation intensity I after the radiation penetrated the standard sample under the same testing conditions was obtained. 实测 Among them, the thicknesses of the uncobalt-removed layer and the cobalt-removed layer of the known standard sample are d1 and d2, respectively; The theoretical attenuation intensity of X-rays penetrating the cobalt-depleted layer and the non-cobalt-depleted layer is obtained as I. 理论 Among them, I 理论 =I0·e -(μ1d1+μ2d2) I0 is the incident intensity when the radiation has no attenuation and no scattering under ideal conditions, μ1 is the radiation attenuation coefficient without the cobalt layer, and μ2 is the radiation attenuation coefficient with the cobalt layer. Obtain the scattering correction coefficient β; where β = (I 理论 -I 实测 ) / (I 理论 -Is).

[0033] In this embodiment, the core physical meaning of the scattering correction parameter β is "the proportion of scattered ray intensity to the total detected intensity." Its calibration essentially involves "comparing the ideal theoretical value of the standard sample with the actual measured value" to inversely deduce the proportion of intensity deviation caused by scattering. The core logic is as follows: the thickness of the un-cobalt-removed / cobalt-removed layer of the standard sample is known, and the theoretical ray intensity without scattering can be calculated using the ideal Lambert-Beer law; simultaneously, the measured ray intensity with the sample and the scattering background intensity Is without the sample are used, and the difference between these three values ​​is used to inversely deduce β, ensuring that the corrected ray intensity formula accurately reflects the true attenuation difference between the two layers. For example, given that the un-cobalt-removed layer thickness d1 = 2.0 mm, the cobalt-removed layer thickness d2 = 1.0 mm, and the ray attenuation coefficient of the un-cobalt-removed layer μ1 = 0.8 mm... 1 The X-ray attenuation coefficient μ2 of the cobalt-free layer is 0.3 mm. 1 (This can be pre-calibrated using material composition or by referring to publicly available data on similar materials); When measuring Is, do not place any sample and keep the positions of the X-ray emitting and receiving modules fixed; continuously collect X-ray intensity data 10 times and take the average value as Is. Example measurement result: Is = 200; Measure the X-ray intensity I with a sample. 实测 At that time, the standard sample was horizontally fixed on the sample fixing platform, ensuring that the X-ray penetration path was consistent with the sample thickness direction (i.e., the X-rays passed perpendicularly through the uncobalt-removed layer and the cobalt-removed layer, with a path length d = d1 + d2 = 3.0 mm). Under the same detection parameters, X-ray intensity data were continuously collected 10 times. After removing outliers, the average value was taken. Example measurement results: I 实测 =1380; Calculate the theoretical attenuation intensity I 理论 At that time, based on formula I 理论 =I0·e -(μ1d1+μ2d2) It can be calculated that I0 needs to be calibrated using the "ideal intensity without a sample". In the ideal state without a sample, the radiation has no attenuation and no scattering, and its intensity is I0. However, the actual measured intensity without a sample is Is (containing only scattering). Therefore, I0 needs to be supplemented by a blank experiment with a standard sample: remove the standard sample and substitute a "virtual sample" with a known attenuation coefficient μ0 = 0 (such as air). At this time, the ideal intensity I0 = the incident intensity without scattering, which can be directly read from the equipment parameters. For example, the incident intensity calibrated by the equipment is I0 = 10000. Substituting the above parameters into the formula yields: I 理论 =10000×e -(0.8×2+0.3×1) =1496, and the final calculated β is β=(1496-1380) / (1496-200)=0.0895.

[0034] As an optional implementation, a planar fit is performed on the coordinates of all feature points to obtain a three-dimensional model of the target boundary surface, including: Based on the least squares method, the coordinates of the feature points ( x k , y k , z k We perform a fitting to minimize the weighted sum of squared errors S, where S is expressed as:

[0035] In the formula, k = 1, 2, 3...m, where m is the total number of feature points. w k Let be the weight adjustment parameter for the k-th feature point, and satisfy 0 < 0. w k <1, w k The magnitude of the feature point is positively correlated with the reliability of the feature point, and A, B, and C are the coefficients of the plane equation to be solved. Based on the weighted sum of squared errors S, partial derivatives with respect to A, B, and C are taken and set equal to 0 to obtain the following system of linear equations:

[0036] Solve the linear equations to obtain the plane equation coefficients A, B, and C, and then perform plane fitting to obtain a three-dimensional model of the target boundary surface.

[0037] In this embodiment, it is assumed that the cobalt removal boundary is a plane, with its general equation being Ax + By + Cz-1 = 0. The core objective of the fitting is to minimize the sum of squared distances from all feature points to the fitting plane (i.e., minimize S), because minimizing the sum of squared distances effectively suppresses the influence of outliers on the fitting results and ensures the stability of the planar model. Furthermore, considering that feature points in the edge region are more affected by noise during the feature point extraction process, a weight correction parameter is also introduced. w k By assigning different weights to feature points with different levels of reliability, the fitting accuracy is improved. Then, based on the weighted sum of squared errors S, the partial derivatives of A, B, and C are calculated and set to 0 to obtain the above linear equation system. This linear equation system can be solved by matrix inversion or Gaussian elimination, thereby calculating the plane equation coefficients A, B, and C that determine the spatial position and orientation of the fitting plane. As an optional implementation, let the inclination angle of the target boundary surface be θ, and the expression for θ is:

[0039] In the formula, Δ x Let Δ be the angle of inclination of the bottom surface of the diamond composite sheet under test relative to the horizontal plane in the x-direction. y The angle of inclination of the bottom surface of the diamond composite sheet under test relative to the horizontal plane in the y-direction.

[0040] In this embodiment, the calculation of the tilt angle is essentially done indirectly by using the angle between the normal vectors of the two planes to calculate the angle between the cobalt removal boundary surface and the reference plane (e.g., taking the bottom surface of the composite sheet as the XY plane). The core principle is to utilize the geometric property that "the angle between two planes is equal to the angle between their normal vectors or its supplementary angle, and the acute angle is taken." Let the normal vector of the reference plane (i.e., the XY plane) be... =(0,0,1), and considering that the reference plane (bottom surface) of the diamond composite sheet may have a slight tilt (not an ideal XY plane) when the sample is fixed, a reference plane tilt correction parameter Δ is also introduced. x Δ y (These are the tilt angles of the reference plane around the X and Y axes, respectively), thereby compensating for the tilt angle calculation error caused by the fixed deviation. The tilt angle Δ of the reference plane after the sample is fixed can be measured using a level or laser interferometer. x Δ y Therefore, the corrected reference plane normal vector is =(sinΔ y , -sinΔ x cosΔ x ·cosΔ y ), and Δ x Δ y Generally small, therefore it can be simplified to sinΔ ≈Δ (Radians), cosΔ ≈1, that is =(Δ y , -Δ x ,1), and the fitted normal vector of the cobalt removal boundary surface is If the boundary surface is (A, B, C), then the inclination angle between the boundary surface and the datum surface satisfies:

[0041] The expression for θ can be derived from the above formula, and the calculation is accurate and reliable.

[0042] As an optional implementation method, the weight correction parameter w k The expression is: w k =1-γ· s k / s max

[0043] In the formula, γ is the weighting adjustment coefficient. s k The extraction error of the k-th feature point is... s max The extraction error for all feature points (denoted as ) s 1, s 2, s 3... s m The maximum value of ) x k邻 , y k邻 This represents the coordinates of the feature points adjacent to the current feature point.

[0044] In this embodiment, w k In the calculation formula, the weight adjustment coefficient γ is generally taken as 0.4-0.6, which can be adjusted according to the detection accuracy. s k It can reflect the degree of coordinate fluctuation in the "adjacent frame attenuation image", reflecting the stability of the feature point extraction result. s k The smaller the value, the more stable the position of the feature point in imaging at different heights, and the higher the extraction accuracy (stronger reliability); conversely, s k The larger the value, the greater the interference from noise, scattering, etc., and the weaker the reliability of the feature point. s k This can be calculated using the "coordinate deviation of the same feature point in adjacent frame attenuated images" (see the formula above), and s max The extraction error for all feature points (denoted as ) s 1, s 2, s 3... s m The maximum value in the set is used to normalize the extraction error of all feature points to ensure... s k / s maxThe value of is between 0 and 1 to avoid imbalance in weight allocation due to differences in the absolute error of feature points, thus accurately calculating the weight correction parameter. w k This further improved data accuracy.

[0045] Example 2 Reference Figure 3 Based on the same inventive concept as the aforementioned embodiments, this embodiment also provides a device for measuring the decobalt removal boundary tilt angle of a diamond composite sheet, comprising: The attenuation image acquisition module is used to control the X-ray emission module to emit X-rays at different heights onto a fixed diamond composite sheet under test in order to obtain multiple corresponding attenuation images; wherein, the X-ray emission module is arranged coaxially with the diamond composite sheet under test. The feature point acquisition module is used to perform threshold segmentation on each attenuation image. Based on the difference in ray attenuation intensity between the uncobalt-removed layer and the cobalt-removed layer of the diamond composite sheet under test, the feature point coordinates of the target boundary surface in each attenuation image are extracted; where the target boundary surface is the boundary surface between the uncobalt-removed layer and the cobalt-removed layer. The plane fitting module is used to perform plane fitting on the coordinates of all feature points to obtain a 3D model of the target boundary surface. The tilt angle calculation module is used to obtain the tilt angle of the target boundary surface based on the 3D model. The explanations and examples of each module in the device of this embodiment can be referred to the methods of the foregoing embodiments, and will not be repeated here.

[0046] Example 3 Based on the same inventive concept as the foregoing embodiments, this embodiment provides a computer device, which includes a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the above-described method.

[0047] Example 4 Based on the same inventive concept as the foregoing embodiments, this embodiment provides a computer-readable storage medium storing a computer program, and a processor executes the computer program to implement the above-described method.

[0048] The above are merely preferred embodiments of this application and do not limit the patent scope of this application. Any equivalent structural or procedural transformations made using the content of this application's specification and drawings, or direct or indirect applications in other related technical fields, are similarly included within the patent protection scope of this application.

Claims

1. A method for measuring the inclination angle of the decobalt removal boundary of a diamond composite sheet, characterized in that, The method comprises the following steps: controlling the ray emitting module to emit rays at different heights to the fixed diamond composite to be tested, so as to obtain a plurality of corresponding attenuation images; wherein the ray emitting module is coaxially arranged with the diamond composite to be tested; threshold segmentation is performed on each attenuation image, and feature point coordinates of a target boundary surface in each attenuation image are extracted based on the difference in ray attenuation intensity between the undecobbed layer and the decobbed layer of the diamond composite to be tested; wherein the target boundary surface is the boundary surface between the undecobbed layer and the decobbed layer; plane fitting is performed on all the feature point coordinates to obtain a three-dimensional model of the target boundary surface; an inclination angle of the target boundary surface is obtained according to the three-dimensional model.

2. The method of claim 1, wherein the method further comprises: Threshold segmentation is performed on each attenuation image, and feature point coordinates of a target boundary surface in each attenuation image are extracted based on the difference in ray attenuation intensity between the undecobbed layer and the decobbed layer of the diamond composite to be tested, comprising: based on a preset ray attenuation intensity model, the average ray intensity I1 of the undecobbed layer region and the average ray intensity I2 of the decobbed layer region after the rays penetrate the diamond composite to be tested are obtained respectively; the segmentation threshold T is obtained according to the average ray intensity I1 and the average ray intensity I2, T=(I1+I2) / 2; based on the segmentation threshold T, the attenuation image is segmented into the undecobbed layer region and the decobbed layer region; the feature point coordinates of the target boundary surface located between the undecobbed layer region and the decobbed layer region are obtained.

3. The method of claim 2, wherein the method further comprises: The expression of the ray attenuation intensity model is: I = I' - e -μd (1 - β) + Is - β In the formula, I is the ray intensity after the rays penetrate the diamond composite to be tested, I' is the ray intensity when the rays are incident, e is a natural constant, μ is the ray attenuation coefficient of the material, the ray attenuation coefficient of the undecobbed layer is greater than that of the decobbed layer, d is the ray penetration path length, Is is the average intensity of scattered rays without samples, and β is the calibrated scattering correction coefficient.

4. The method of claim 3, wherein the method further comprises: The calibration method of the scattering correction coefficient β comprises the following steps: A standard sample same as the diamond composite to be tested is selected, and the ray intensity I after the ray penetrates the standard sample under the same detection condition is obtained 实测 ; wherein the thicknesses of the non-decobbed layer and the debobbed layer of the known standard sample are d1 and d2, respectively; The theoretical attenuation intensity of the rays penetrating the decobalt layer and the non-decobalt layer is I 理论 ; wherein, I 理论 = I0·e -(μ1d1+μ2d2) , I0 is the incident intensity in an ideal state without attenuation and scattering of the rays, μ1 is the ray attenuation coefficient of the non-decobalt layer, and μ2 is the ray attenuation coefficient of the decobalt layer; Obtaining a scattering correction coefficient β; wherein, β=(I 理论 -I 实测 ) / (I 理论 -Is).

5. The method of claim 1-4, wherein, Plane fitting is performed on all the feature point coordinates to obtain a three-dimensional model of the target boundary surface, comprising: Based on the least squares method, the coordinates of the feature points ( x k , y k , z k We perform a fitting to minimize the weighted sum of squared errors S, where S is expressed as: In the formula, k = 1, 2, 3...m, m is the total number of feature points, w k is the weight correction parameter of the kth feature point, and satisfies 0 < k < m, w k <1, w k The size of is positively correlated with the reliability of the feature point, and A, B, and C are respectively the plane equation coefficients to be solved. based on the weighted error sum of squares S, the partial derivatives of A, B and C are solved and set to 0 to obtain the following linear equation group: the linear equation group is solved to obtain the plane equation coefficients A, B and C, so as to complete the plane fitting to obtain the three-dimensional model of the target boundary surface.

6. The method of claim 5, wherein the method further comprises: Let the inclination angle of the target boundary surface be θ, and the expression of θ is: where Δ x is the inclination of the bottom surface of the diamond compact to be measured in the x direction relative to the horizontal plane, Δ y is the inclination of the bottom surface of the diamond compact to be measured in the y direction relative to the horizontal plane.

7. The method of claim 5, wherein the method further comprises: weight correction parameter w k The expression is: w k = 1 - γ · σ k / σ max where γ is a weight adjustment coefficient, σ k is the extraction error of the kth feature point, σ max is the maximum value of the extraction errors of all the feature points, x k邻 , y k邻 denotes the coordinates of the feature points adjacent to the current feature point.

8. A device for measuring the angle of a decobalt boundary of a diamond compact, comprising: comprising: an attenuation image acquisition module configured to control the ray emitting module to emit rays at different heights to the fixed diamond composite to be tested, so as to obtain a plurality of corresponding attenuation images; wherein the ray emitting module is coaxially arranged with the diamond composite to be tested; a feature point acquisition module configured to perform threshold segmentation on each attenuation image, and extract feature point coordinates of a target boundary surface in each attenuation image based on the difference in ray attenuation intensity between the undecobbed layer and the decobbed layer of the diamond composite to be tested; wherein the target boundary surface is the boundary surface between the undecobbed layer and the decobbed layer; a plane fitting module configured to perform plane fitting on all the feature point coordinates to obtain a three-dimensional model of the target boundary surface; an inclination angle calculation module configured to obtain an inclination angle of the target boundary surface according to the three-dimensional model.

9. A computer device, comprising: The computer device comprises a memory and a processor, the memory stores a computer program, and the processor executes the computer program to realize the method in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, The computer readable storage medium stores a computer program, and the processor executes the computer program to realize the method in any one of claims 1-7.

Images