Extremely weak magnetic determination method for traditional Chinese medicine acupuncture and moxibustion

Abstract

The invention discloses an extremely weak magnetic measurement method for traditional Chinese medicine acupuncture and moxibustion, and relates to the technical field of biological magnetic measurement. Calculating and analyzing the image data to obtain a needle body deviation factor, judging whether the needle body deviates or not, and returning an alarm if the needle body is abnormal; and finally calculating and analyzing magnetic field data, and repairing the magnetic field to be normal according to a result. Preprocessed image data are collected in real time, the pin body deviation is judged through an algorithm and a threshold value, an alarm is given, interference is reduced, and data accuracy is guaranteed; and the magnetic field data is deeply analyzed, interference factors are generated, and an intelligent decoupling repair algorithm is used for eliminating anomalies, continuing data, eliminating interference and avoiding misjudgment.

Description

Technical Field

[0001] The present invention relates to the technical field of biomagnetic measurement, in particular to an extremely weak magnetic measurement method for traditional Chinese medicine acupuncture. Background Art

[0002] The detection of biomagnetic field during acupuncture in traditional Chinese medicine is mainly through detecting the changes in weak magnetic fields of human meridians or acupoints, assisting in judging physical condition, acupoint sensitivity, and evaluating physiological reactions before and after acupuncture treatment, providing modern scientific and technological reference for accurate acupoint selection and efficacy analysis.

[0003] For example, Chinese Patent Publication No. CN118710991A provides an acupuncture parameter measurement device and an acupuncture manipulation classification method based thereon, comprising an upper platform, a lower platform, three flexible hinges, and three optical fibers; the three flexible hinges are fixed between the upper and lower platforms, spaced 120° apart and arranged on a circumference around the central axis, forming an elastic body together with the upper and lower platforms; three optical fibers are respectively arranged between every two flexible hinges, two of which are in a tightly suspended state and arranged in parallel and symmetrical manner, while the third optical fiber is only fixed to the lower platform at its lower end and is in a free state; the acupuncture parameter measurement device calculates the axial force and axial torque of the acupuncture needle installed on the elastic body by converting the grating center wavelength drift of the three optical fibers. Acupuncture manipulation classification is based on the measured axial force and axial torque, and a neural network model is trained to classify unknown acupuncture manipulations. The present invention can improve the accuracy of acupuncture parameter measurement and acupuncture manipulation classification.

[0004] However, the above patent does not take into account the patient's movement during acupuncture, which causes the human tissue to mechanically pull or squeeze the needle body, resulting in the needle body being offset at the target point of the acupuncture point, thereby affecting the collected biomagnetic field data and reducing the accuracy. Summary of the Invention

[0005] Technical problems solved

[0006] In response to the shortcomings of the existing technology, the present invention provides an extremely weak magnetic measurement method for traditional Chinese medicine acupuncture, which solves the problem that the patient's movement causes the human tissue to mechanically pull or squeeze the needle body, resulting in the displacement of the needle body target point, thereby affecting the collected biomagnetic field data and reducing the detection accuracy.

[0007] Technical Solution

[0008] To achieve the above purpose, the present invention is implemented through the following technical solutions: A very weak magnetic measurement method for traditional Chinese medicine acupuncture, comprising the following specific steps: Step 1: Start real-time acquisition of image data and magnetic field data, and perform preprocessing; Step 2: Perform comprehensive calculation on the preprocessed image data to obtain patient movement parameters and needle body dislocation parameters, the patient movement parameters are used to analyze whether the patient has moved, and the needle body dislocation parameters are used to analyze whether the needle body is dislocated, and a needle body offset factor is obtained based on the comprehensive calculation of the patient movement parameters and the needle body dislocation parameters, and the needle body offset is classified according to the needle body offset factor. Analysis, if the analysis shows that the needle body is normal, execute step three, if the analysis shows that the needle body is offset, issue an alarm and return to step one; Step three: perform comprehensive calculation on the preprocessed magnetic field data to obtain magnetic field anomaly parameters, which are used to analyze whether the magnetic field is abnormal, and perform comprehensive calculation based on the needle body misalignment parameters and the magnetic field anomaly parameters to obtain the magnetic field interference factor, and analyze whether the magnetic field is interfered with based on the magnetic field interference factor. If the analysis shows that the magnetic field is interfered with, repair the magnetic field through the biomagnetic field intelligent decoupling repair algorithm until the magnetic field is analyzed to be normal. If the analysis shows that the magnetic field is normal, end.

[0009] Furthermore, the specific method of analyzing whether the needle body is deviated based on the needle body deviation factor is as follows: a needle body deviation threshold is set through historical experiments and stored in a database as a quantitative standard for whether the needle body is deviated, and the needle body deviation factor is compared with the needle body deviation threshold. If the needle body deviation factor is less than the needle body deviation threshold, it is analyzed that the needle body is normal; if the needle body deviation factor is greater than or equal to the needle body deviation threshold, it is analyzed that the needle body is deviated.

[0010] Furthermore, the specific method of obtaining the needle body offset factor is as follows: the preprocessed image data is processed by an image contour segmentation algorithm and an image contour tracking algorithm to obtain contour image data, which includes various contours, and each contour includes the number of pixels of the contour and the coordinates of the pixel points of the contour. The patient movement parameters and the needle body misalignment parameters are normalized, and then comprehensively calculated to obtain the needle body offset factor.

[0011] Furthermore, the specific method for obtaining the patient movement parameters is as follows: various contours are screened to obtain the patient contour, a standard contour of the needle body is set and stored in a database, various contours of the contour image data are matched with the standard contour of the needle body to obtain the actual contour of the needle body, the start time of acupuncture of the needle body is determined according to the time series and the actual contour of the needle body, the time after the start time of acupuncture of the needle body is recorded as the acupuncture time, the center point of the patient contour is calculated to obtain the coordinates of the patient's center point, and only within the acupuncture time, the distance of the patient's center point coordinates at different times is calculated according to the time series and the Euclidean distance formula. The Euclidean distance formula is used to calculate the straight-line distance between the patient's center point coordinates at two different times, thereby obtaining the patient movement parameters.

[0012] Furthermore, the specific method for obtaining the patient contour is as follows: setting the acupuncture position contour, matching the pixel point coordinates within various contours with the pixel point coordinates within the acupuncture position contour to obtain the number of identical coordinates, sorting the number of identical coordinates of various contours and the acupuncture position contour in ascending order through a sorting algorithm to obtain the contour with the largest number of identical coordinates, setting a standard threshold for the human body contour area, calculating the area of ​​various contours, and obtaining the actual contour area. If the actual contour area is within the standard threshold for the human body contour area and this contour is the contour with the largest number of identical coordinates, then this contour is marked as the patient contour; if the actual contour area is outside the standard threshold for the human body contour area or this contour is not the contour with the largest number of identical coordinates, then it is not marked.

[0013] Furthermore, the specific method for obtaining the actual outline of the needle body is as follows: the standard outline of the needle body includes color and area, and the area is calculated by summing the number of pixels within the outline. The color and area of ​​the standard outline of the needle body are matched with various outlines. If the color and area of ​​a certain outline are the same as the color and area of ​​the standard outline of the needle body, it is marked as the actual outline of the needle body. If the color and area of ​​a certain outline are different from the color or area of ​​the standard outline of the needle body, it is not marked.

[0014] Furthermore, the specific method for obtaining the acupuncture time is as follows: the center point of the actual contour of the needle body is calculated to obtain the coordinates of the center point of the needle body. At a certain moment in the time series, recorded as the first moment, vector calculation is performed on the coordinates of the center point of the needle body and the coordinates of each pixel point of the patient's contour to obtain the acupuncture preview vector. The movement trend is calculated based on the continuous movement of the coordinates of the center point of the needle body in the time series to obtain the needle body action vector. If the needle body action vector is equal to a certain acupuncture preview vector, the area of ​​the patient's contour and the area of ​​the actual contour of the needle body are comprehensively calculated to obtain the acupuncture area value. If the needle body action vector is equal to a certain acupuncture preview vector and the acupuncture area value is greater than zero, the first moment is marked as the acupuncture moment, and the acupuncture moment and the moment after the acupuncture moment are summed to obtain the acupuncture time. If the needle body action vector is not equal to a certain acupuncture preview vector or the acupuncture area value is equal to zero, the first moment is not marked.

[0015] Furthermore, the acupuncture area value is specifically obtained as follows: the area of ​​the actual contour of the needle body and the area of ​​the patient's contour are summed to obtain the contour experimental area; according to the time series, the contour experimental area at the previous moment and the contour experimental area at the next moment are subtracted to obtain the contour experimental area deviation value; according to the number of moments in the time series, the contour experimental area deviation values ​​are summed to obtain the acupuncture area value; Among them, ZM represents the acupuncture area value, M represents the number of moments in the time series, and SM tIndicates the contour experimental area of ​​the previous moment, SM t+1 Indicates the contour experimental area at the next moment.

[0016] Furthermore, the specific method for obtaining the magnetic field anomaly parameters is as follows: the magnetic field data includes the magnetic field amplitude and the magnetic field frequency, the magnetic field amplitude includes the maximum magnetic field amplitude and the minimum magnetic field amplitude, the maximum magnetic field amplitude and the minimum magnetic field amplitude are comprehensively calculated to obtain the magnetic field amplitude deviation value, the magnetic field frequency includes the maximum magnetic field frequency and the minimum magnetic field frequency, the maximum magnetic field frequency and the minimum magnetic field frequency are comprehensively calculated to obtain the magnetic field frequency deviation value; the magnetic field anomaly parameters are obtained by comprehensively calculating the magnetic field amplitude deviation value and the magnetic field frequency deviation value.

[0017] Furthermore, the specific steps of repairing the magnetic field through the biomagnetic field intelligent decoupling repair algorithm are: eliminating the maximum magnetic field amplitude, minimum magnetic field amplitude, maximum magnetic field frequency and minimum magnetic field frequency, and connecting the magnetic field data that are not eliminated to obtain continuous magnetic field data.

[0018] Beneficial effects

[0019] Compared with the prior art, the embodiments of the present invention have at least the following advantages or beneficial effects:

[0020] 1. By real-time acquisition and preprocessing of image data, the image contour segmentation and tracking algorithm is used to extract the contours of the patient and needle body. The patient movement parameters are calculated with the help of the Euclidean distance formula. The needle body misalignment parameters are analyzed in combination with the change value of the acupuncture area. After normalization, the needle body offset factor is generated. By comparing it with the preset threshold, it can accurately determine whether the needle body is offset. Once the offset is detected, an alarm is immediately issued, effectively avoiding the problem of needle body target displacement caused by patient movement, reducing human interference at the source of data collection, ensuring the initial accuracy of image and magnetic field data, and providing a reliable basis for subsequent analysis.

[0021] 2. Perform in-depth analysis of the pre-processed magnetic field data, calculate the magnetic field amplitude and frequency deviation values ​​to obtain magnetic field anomaly parameters, and combine the needle body misalignment parameters to generate a magnetic field interference factor to determine whether the magnetic field is interfered with. When interference is detected, the biomagnetic field intelligent decoupling and repair algorithm is used to eliminate the abnormal maximum or minimum magnetic field amplitude and maximum or minimum magnetic field frequency, and continue the normal data to form a continuous sequence. This mechanism effectively eliminates sudden interference of the needle body offset or the external environment on the magnetic field signal, avoiding misjudgment caused by data jumps or interruptions.

[0022] Of course, any product implementing the present invention does not necessarily need to achieve all of the advantages described above at the same time. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] Figure 1 This invention is a flow chart of an extremely weak magnetic measurement method for traditional Chinese medicine acupuncture. DETAILED DESCRIPTION

[0024] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0025] It should be noted that, in this document, relational terms such as first and second are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "include," "comprise," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that includes a list of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such process, method, article, or apparatus.

[0026] Example 1:

[0027] like Figure 1 As shown, the embodiment of the present invention provides an extremely weak magnetic measurement method for traditional Chinese medicine acupuncture, which includes the following specific steps:

[0028] Step 1: Start collecting image data in real time through the camera, where the image data is a two-dimensional image, and the coordinate system of the two-dimensional image is a two-dimensional coordinate system. Collect magnetic field data in real time through the fluxgate sensor. The fluxgate sensor is arranged near acupuncture points or on the meridian path to improve the accuracy of detecting the human body's biomagnetic field. The image data is cleaned to remove blurred images, and the magnetic field data is filtered and denoised to remove interference noise, which helps to improve the data quality of the image data and magnetic field data.

[0029] Step 2: The cleaned image data is processed by an image contour segmentation algorithm and an image contour tracking algorithm. An image contour segmentation algorithm, such as the Sobel algorithm, first converts the image data into a grayscale image to simplify calculations. Then, it is convolved with 3×3 convolution kernels in the horizontal and vertical directions to obtain horizontal and vertical gradient matrices respectively. The gradient amplitude and direction are then calculated. The former determines the edge, and the latter determines the edge extension direction. Finally, through threshold binarization, pixels above the threshold are marked as contours, and contour image data is finally obtained. The contour image data includes various contours, but it is initially uncertain which contour is the patient's contour. Various contours include the number of contour pixels and the coordinates of the contour pixels. An image contour tracking algorithm, such as the Lucas-Kanade optical flow method, establishes an optical flow constraint equation based on the assumption that the brightness of pixels in adjacent frames is constant. The spatial and temporal gradients are calculated using the Sobel algorithm and frame difference. Assuming that the neighborhood motion is consistent, an overdetermined set of equations is constructed and the optical flow velocity is solved using the least squares method. Finally, feature points such as corner points are selected for tracking and position update, thereby tracking the contours in the contour image data.

[0030] Standardizing the contour image data helps eliminate the dimension of the contour image data and convert the values ​​of the contour image data of different orders of magnitude into a unified numerical range. Then, a comprehensive calculation is performed to obtain the patient movement parameter and the needle body misalignment parameter. The patient movement parameter and the needle body misalignment parameter are normalized to further eliminate the dimension of the patient movement parameter and the needle body misalignment parameter and convert the values ​​of the patient movement parameter and the needle body misalignment parameter of different orders of magnitude into a unified numerical range. Then, a comprehensive calculation is performed to obtain the needle body offset factor, ZY=YC+WC; wherein ZY represents the needle body offset factor, which reflects whether the needle body is offset, YC represents the patient movement parameter, which reflects whether the patient has moved, and WC represents the needle body misalignment parameter, which reflects whether the needle body is misaligned. Whether the needle body is offset is analyzed based on the needle body offset factor. If the analysis shows that the needle body is normal, step three is executed. If the analysis shows that the needle body is offset, an alarm is issued to remind the patient not to move to avoid mechanical traction or compression of the needle body by human tissue, resulting in displacement of the needle body target point, thereby affecting the collected biomagnetic field data and reducing the detection accuracy, and returning to step one.

[0031] The specific method of analyzing whether the needle body is deflected according to the needle body deflection factor is as follows:

[0032] A needle body deviation threshold is set through historical experiments and stored in a database. This serves as a quantitative standard for determining whether the needle body is deviated. The needle body deviation factor is compared with the needle body deviation threshold. If the needle body deviation factor is less than the needle body deviation threshold, the needle is considered normal. If the needle body deviation factor is greater than or equal to the needle body deviation threshold, the needle is considered deviated.

[0033] The specific method of obtaining the patient movement parameters is as follows:

[0034] Since the contour image data includes various contours, the various contours are screened to obtain the patient contour, and the standard contour of the needle body is set and stored in the database. The various contours of the contour image data are matched with the standard contour of the needle body to obtain the actual contour of the needle body. The start time of acupuncture of the needle body is determined according to the time series and the actual contour of the needle body. The time after the start time of acupuncture of the needle body is recorded as the acupuncture time. The center point of the patient contour is calculated to obtain the coordinates of the patient's center point. Only within the acupuncture time, the distance between the coordinates of the patient's center point at different times is calculated according to the time series and the Euclidean distance formula. The Euclidean distance formula is used to calculate the straight-line distance between the coordinates of the patient's center point at two different times, thereby obtaining the patient's movement parameters.

[0035]

[0036] Among them, YC represents the patient movement parameter, which reflects the degree of patient movement. The larger the patient movement parameter, the more the patient moves. j+1 The horizontal coordinate of the patient's center point at the j+1th moment, x j The horizontal coordinate of the patient's center point at the jth moment, y j+1 The ordinate of the patient's center point at the j+1th moment, y j The ordinate represents the coordinate of the patient's center point at the jth moment.

[0037] The specific method of obtaining the patient profile is as follows:

[0038] An acupuncture position contour is set to constrain the patient's position range during acupuncture. Since the contour image data includes various contours, the pixel coordinates within the various contours are matched with the pixel coordinates within the acupuncture position contour to obtain the number of identical coordinates. This is used to analyze which contour is closest to the acupuncture position contour. The closer the contour is to the acupuncture position contour, the greater the probability that it is the patient's contour. A sorting algorithm is used to sort the various contours in ascending order of the number of identical coordinates with the acupuncture position contour. The sorting algorithm, such as bubble sort, compares and exchanges adjacent contours with the same number of coordinates, gradually moving the contour with the largest number of identical coordinates to the end of the array, and finally achieving a sort from small to large order to obtain the contour with the largest number of identical coordinates, i.e., the contour with the highest probability of being the patient's contour. A standard threshold for the area of ​​human body contours is set, and the areas of various contours are calculated, i.e., the number of pixels within various contours is summed to obtain the actual contour area. If the actual contour area is within the standard threshold for the area of ​​human body contours and this contour has the largest number of identical coordinates, then this contour is marked as the patient's contour. If the actual contour area is outside the standard threshold for the area of ​​human body contours or this contour is not the contour with the largest number of identical coordinates, then this contour is not the patient's contour and is not marked.

[0039] The specific method of obtaining the actual outline of the needle body is as follows:

[0040] The standard needle body contour includes color and area. The area is calculated by summing the number of pixels within the contour. The color and area of ​​the standard needle body contour are matched with various contours. If the color and area of ​​a contour are the same as those of the standard needle body contour, it is marked as the actual needle body contour. If the color and area of ​​a contour are different from the color or area of ​​the standard needle body contour, it is not marked.

[0041] The specific method of obtaining acupuncture time is as follows:

[0042] The center point of the actual contour of the needle body is calculated to obtain the coordinates of the center point of the needle body. At a certain moment in the time series, recorded as the first moment, vector calculation is performed on the coordinates of the center point of the needle body and the coordinates of each pixel point of the patient's contour to obtain the acupuncture preview vector. The movement trend is calculated based on the continuous movement of the coordinates of the center point of the needle body in the time series to obtain the needle body action vector, that is, the needle body's movement trend is analyzed to determine whether it is piercing the patient's contour. If the needle body action vector is equal to a certain acupuncture preview vector, it means that the needle body is piercing the patient, and the area of ​​the patient's contour and the area of ​​the actual needle body contour are comprehensively calculated to obtain the acupuncture area value. If the needle body action vector is equal to a certain acupuncture preview vector and the acupuncture area value is greater than zero, the first moment is marked as the acupuncture moment, and the acupuncture moment and the moment after the acupuncture moment are summed to obtain the acupuncture time. If the needle body action vector is not equal to a certain acupuncture preview vector or the acupuncture area value is zero, the first moment is not marked.

[0043] The specific method of obtaining the acupuncture preview vector is as follows:

[0044] The coordinate of the needle center is marked as (x0, y0), and the coordinate within the patient's outline is marked as (x i ,y i ), i is not 0, so in Represents the acupuncture preview vector.

[0045] The specific method of obtaining the needle action vector is as follows:

[0046] Vector calculation is performed on the coordinates of the center point of the needle at the first moment and at each moment after the first moment to obtain an actual needle movement vector group, and each actual needle movement vector in the actual needle movement vector group is averaged to obtain the needle action vector.

[0047] The specific method of obtaining the acupuncture area value is as follows:

[0048] The area of ​​the actual needle contour and the area of ​​the patient contour are summed to obtain the contour experimental area. Based on the time series, the difference between the contour experimental area at the previous moment and the contour experimental area at the next moment is calculated to obtain the contour experimental area deviation value, which reflects that the needle has penetrated into the patient's body. Therefore, the contour experimental area at the next moment is smaller than the contour experimental area at the previous moment. The contour experimental area deviation values ​​are summed according to the number of moments in the time series to obtain the acupuncture area value, which reflects the degree of needle penetration into the patient.

[0049]

[0050] Among them, ZM represents the acupuncture area value, M represents the number of moments in the time series, and SM t Indicates the contour experimental area of ​​the previous moment, SM t+1 Indicates the contour experimental area at the next moment.

[0051] The specific method of obtaining the needle body misalignment parameters is as follows:

[0052] During the acupuncture time, the difference between the acupuncture area value at the previous moment and the acupuncture area value at the next moment is calculated to obtain the acupuncture area change value. If the acupuncture area change value is less than zero, it means that the needle body has moved away from the patient. Therefore, this acupuncture area change value is recorded as the needle body misalignment parameter. If the acupuncture area change value is greater than or equal to zero, the acupuncture area change value is not marked. The reason why the coordinate change of the needle body center point coordinate is not used to analyze the movement of the needle body contour is that when the patient contour moves, although the actual needle body contour does not move relative to the patient, it moves relative to the acupuncture position contour, which causes calculation ambiguity. In addition, the area change of the needle body contour reflects whether the needle body has moved out of the patient contour.

[0053] The number of changes in the acupuncture area change value is counted to obtain the number of acupuncture movement changes. The number of acupuncture movement changes reflects the movement of acupuncture, so as to avoid the acupuncture area change value moving at time A and recovering at time B, and avoid the acupuncture area change value not changing between the two moments.

[0054] The needle body dislocation parameter is obtained by comprehensive calculation based on the change value of acupuncture area and the number of acupuncture movement changes;

[0055] WC=ZC-MB;

[0056] Among them, WC represents the needle body dislocation parameter, ZC represents the number of acupuncture movement changes, and the more the number of acupuncture movement changes, the larger the needle body dislocation parameter. The smaller the acupuncture area change value, the deeper the degree of acupuncture separation from the patient, and the larger the needle body dislocation parameter.

[0057] The difference between Example 2 and Example 1 is that:

[0058] Step 3: Standardizing the magnetic field data helps to eliminate the dimension of the magnetic field data, and converting the different orders of magnitude values ​​of the magnetic field data into a unified numerical range, and then performing comprehensive calculations to obtain the magnetic field anomaly parameters, normalizing the needle body misalignment parameters and the magnetic field anomaly parameters, and further eliminating the dimensions of the needle body misalignment parameters and the magnetic field anomaly parameters, and converting the different orders of magnitude values ​​of the needle body misalignment parameters and the magnetic field anomaly parameters into a unified numerical range, and then performing comprehensive calculations to obtain the magnetic field interference factor, CY=WC+YS; wherein CY represents the magnetic field interference factor, reflecting whether the magnetic field is disturbed, WC represents the needle body misalignment parameter, and YS represents the magnetic field anomaly parameter, reflecting whether the magnetic field is abnormal. Whether the magnetic field is disturbed is analyzed based on the magnetic field interference factor. If it is analyzed that the magnetic field is disturbed, the magnetic field is repaired by the biomagnetic field intelligent decoupling repair algorithm until the magnetic field is analyzed to be normal. If the magnetic field is analyzed to be normal, the process ends.

[0059] The specific method of analyzing whether the magnetic field is disturbed according to the magnetic field interference factor is as follows:

[0060] The magnetic field interference threshold is set through historical experiments and stored in the database as a quantitative standard for whether the magnetic field is disturbed. The magnetic field interference factor is compared with the magnetic field interference threshold. If the magnetic field interference factor is greater than the magnetic field interference threshold, it is analyzed that the magnetic field is disturbed. If the magnetic field interference factor is less than or equal to the magnetic field interference threshold, it is analyzed that the magnetic field is normal.

[0061] The specific method of obtaining magnetic field anomaly parameters is as follows:

[0062] The magnetic field data includes the magnetic field amplitude and the magnetic field frequency. The magnetic field amplitude includes the maximum magnetic field amplitude and the minimum magnetic field amplitude. The maximum magnetic field amplitude and the minimum magnetic field amplitude are comprehensively calculated to obtain the magnetic field amplitude deviation value. The magnetic field frequency includes the maximum magnetic field frequency and the minimum magnetic field frequency. The maximum magnetic field frequency and the minimum magnetic field frequency are comprehensively calculated to obtain the magnetic field frequency deviation value.

[0063] The magnetic field anomaly parameters are obtained by performing comprehensive calculation based on the magnetic field amplitude deviation value and the magnetic field frequency deviation value;

[0064] YS=FP+PP;

[0065] Among them, YS represents the magnetic field anomaly parameter, FP represents the magnetic field amplitude deviation value, which reflects the degree of magnetic field amplitude deviation, and PP represents the magnetic field frequency deviation value, which reflects the degree of magnetic field frequency deviation.

[0066] The specific steps for repairing the magnetic field through the biomagnetic field intelligent decoupling repair algorithm are as follows:

[0067] Since the maximum magnetic field amplitude, minimum magnetic field amplitude, maximum magnetic field frequency and minimum magnetic field frequency reflect the highest probability of needle movement, the maximum magnetic field amplitude, minimum magnetic field amplitude, maximum magnetic field frequency and minimum magnetic field frequency are eliminated, and the magnetic field data that are not eliminated are continued, that is, the magnetic field values ​​at both ends of the interruption are forced to be equal to obtain continuous magnetic field data.

[0068] The preferred embodiments of the present invention disclosed above are intended only to help illustrate the present invention. These preferred embodiments do not exhaustively describe all details, nor do they limit the present invention to the specific embodiments described. Obviously, many modifications and variations are possible based on the contents of this specification. These embodiments are selected and described in detail in this specification to better explain the principles and practical applications of the present invention, thereby enabling those skilled in the art to better understand and utilize the present invention. The present invention is limited only by the claims and their full scope and equivalents.

Claims

1. A very weak magnetic measurement method for traditional Chinese medicine acupuncture, characterized by: The specific steps include: Step 1: Start real-time acquisition of image data and magnetic field data, and perform preprocessing; Step 2: Perform comprehensive calculation on the pre-processed image data to obtain patient movement parameters and needle misalignment parameters. The patient movement parameters are used to analyze whether the patient has moved, and the needle misalignment parameters are used to analyze whether the needle has been misaligned. Comprehensive calculation is performed based on the patient movement parameters and the needle misalignment parameters to obtain a needle offset factor. Based on the needle offset factor, analyze whether the needle has been misaligned. If the analysis shows that the needle is normal, execute step 3. If the analysis shows that the needle has been misaligned, issue an alarm and return to step 1. Step 3: Perform comprehensive calculation on the pre-processed magnetic field data to obtain magnetic field anomaly parameters. The magnetic field anomaly parameters are used to analyze whether the magnetic field is abnormal. Comprehensive calculation is performed based on the needle body misalignment parameters and the magnetic field anomaly parameters to obtain the magnetic field interference factor. Based on the magnetic field interference factor, analyze whether the magnetic field is disturbed. If it is analyzed that the magnetic field is disturbed, the magnetic field is repaired by the biomagnetic field intelligent decoupling repair algorithm until the magnetic field is analyzed to be normal. If the magnetic field is analyzed to be normal, the process ends.

2. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture according to claim 1, characterized in that: The specific method of analyzing whether the needle body is deflected according to the needle body deflection factor is as follows: A needle body deviation threshold is set through historical experiments and stored in a database. This serves as a quantitative standard for determining whether the needle body is deviated. The needle body deviation factor is compared with the needle body deviation threshold. If the needle body deviation factor is less than the needle body deviation threshold, the needle is considered normal. If the needle body deviation factor is greater than or equal to the needle body deviation threshold, the needle is considered deviated.

3. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture according to claim 2, characterized in that: The specific method of obtaining the needle body offset factor is as follows: The preprocessed image data is processed by the image contour segmentation algorithm and the image contour tracking algorithm to obtain contour image data. The contour image data includes various contours, and each contour includes the number of pixels and the coordinates of the pixel points of the contour. The patient movement parameters and the needle body dislocation parameters are normalized and then comprehensively calculated to obtain the needle body offset factor.

4. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture according to claim 2, characterized in that: The specific method of obtaining the patient movement parameters is as follows: Various contours are screened to obtain the patient contour, and the standard contour of the needle body is set and stored in the database. Various contours of the contour image data are matched with the standard contour of the needle body to obtain the actual contour of the needle body. The start time of acupuncture of the needle body is determined according to the time series and the actual contour of the needle body. The time after the start time of acupuncture of the needle body is recorded as the acupuncture time. The center point of the patient contour is calculated to obtain the coordinates of the patient's center point. Only within the acupuncture time, the distance between the coordinates of the patient's center point at different times is calculated according to the time series and the Euclidean distance formula. The Euclidean distance formula is used to calculate the straight-line distance between the coordinates of the patient's center point at two different times, thereby obtaining the patient's movement parameters.

5. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture according to claim 4, characterized in that: The specific method of obtaining the patient profile is as follows: Set the acupuncture position contour, match the pixel coordinates within various contours with the pixel coordinates within the acupuncture position contour to obtain the number of identical coordinates, sort the various contours with the acupuncture position contour in ascending order of the number of identical coordinates through a sorting algorithm, and obtain the contour with the largest number of identical coordinates. Set the standard threshold for the human body contour area, calculate the area of ​​various contours, and obtain the actual contour area. If the actual contour area is within the standard threshold for the human body contour area and this contour is the contour with the largest number of identical coordinates, then this contour is marked as the patient contour. If the actual contour area is outside the standard threshold for the human body contour area or this contour is not the contour with the largest number of identical coordinates, then it is not marked.

6. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture according to claim 4, characterized in that: The specific method of obtaining the actual outline of the needle body is as follows: The standard needle body contour includes color and area. The area is calculated by summing the number of pixels within the contour. The color and area of ​​the standard needle body contour are matched with various contours. If the color and area of ​​a contour are the same as those of the standard needle body contour, it is marked as the actual needle body contour. If the color and area of ​​a contour are different from the color or area of ​​the standard needle body contour, it is not marked.

7. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture according to claim 4, characterized in that: The specific method of obtaining the acupuncture time is as follows: The center point of the actual contour of the needle body is calculated to obtain the coordinates of the center point of the needle body. At a certain moment in the time series, recorded as the first moment, vector calculation is performed on the coordinates of the center point of the needle body and the coordinates of each pixel point of the patient's contour to obtain the acupuncture preview vector. The movement trend is calculated based on the continuous movement of the coordinates of the center point of the needle body in the time series to obtain the needle body action vector. If the needle body action vector is equal to a certain acupuncture preview vector, the area of ​​the patient's contour and the area of ​​the actual needle body contour are comprehensively calculated to obtain the acupuncture area value. If the needle body action vector is equal to a certain acupuncture preview vector and the acupuncture area value is greater than zero, the first moment is marked as the acupuncture moment, and the acupuncture moment and the moment after the acupuncture moment are summed to obtain the acupuncture time. If the needle body action vector is not equal to a certain acupuncture preview vector or the acupuncture area value is equal to zero, the first moment is not marked.

8. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture according to claim 7, characterized in that: The specific method of obtaining the acupuncture area value is as follows: The area of ​​the actual needle body contour and the area of ​​the patient's contour are summed to obtain the contour experimental area. According to the time series, the difference between the contour experimental area at the previous moment and the contour experimental area at the next moment is calculated to obtain the contour experimental area deviation value. The contour experimental area deviation values ​​are summed according to the number of moments in the time series to obtain the acupuncture area value. Among them, ZM represents the acupuncture area value, M represents the number of moments in the time series, and SM t Indicates the contour experimental area of ​​the previous moment, SM t+1 Indicates the contour experimental area at the next moment.

9. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture according to claim 1, characterized in that: The specific method of obtaining the magnetic field anomaly parameters is as follows: The magnetic field data includes the magnetic field amplitude and the magnetic field frequency. The magnetic field amplitude includes the maximum magnetic field amplitude and the minimum magnetic field amplitude. The maximum magnetic field amplitude and the minimum magnetic field amplitude are comprehensively calculated to obtain the magnetic field amplitude deviation value. The magnetic field frequency includes the maximum magnetic field frequency and the minimum magnetic field frequency. The maximum magnetic field frequency and the minimum magnetic field frequency are comprehensively calculated to obtain the magnetic field frequency deviation value. The magnetic field anomaly parameters are obtained by comprehensively calculating the magnetic field amplitude deviation value and the magnetic field frequency deviation value.

10. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture according to claim 1, characterized in that: The specific steps of repairing the magnetic field by the biomagnetic field intelligent decoupling repair algorithm are as follows: The maximum magnetic field amplitude, the minimum magnetic field amplitude, the maximum magnetic field frequency and the minimum magnetic field frequency are eliminated, and the magnetic field data that are not eliminated are connected to obtain continuous magnetic field data.

Images