Zero magnetic quantum sensing meridian and viscera coupling imaging method and device

Abstract

The application provides a zero-magnetic quantum sensing meridian and viscera coupling imaging method and device, the method comprising: obtaining observation data of a target object; generating an individualized meridian diagram of the target object based on a three-dimensional grid of the body surface and a body surface landmark point of the target object; constructing a viscera functional area diagram and a body source forward matrix, and combining the individualized meridian diagram to construct a meridian line source forward matrix; determining edge admittance parameters of each edge in the individualized meridian diagram based on the observation data, and constructing a time-varying meridian diagram Laplace matrix based on the same; constructing an observation equation based on the observation data, the body source forward matrix, the meridian line source forward matrix and the individualized meridian diagram, and constructing an inversion target function based on the observation data, the viscera functional area diagram and the time-varying meridian diagram Laplace matrix; jointly estimating the functional source distribution, the body surface meridian flow injection distribution and the meridian and viscera coupling operator based on the observation equation and the inversion target function, and imaging according to the estimated data.

Description

Technical Field

[0001] This invention relates to the field of computer imaging technology, and in particular to a zero-magnetic quantum sensing method and device for coupled imaging of meridians and internal organs. Background Technology

[0002] Existing meridian objectification technologies mainly include meridian resistance or conductivity detection, meridian thermal or infrared detection, and meridian optical or ultra-weak luminescence detection.

[0003] Meridian resistance or conductivity detection typically involves contacting electrodes with the skin to measure the resistance or conductivity of acupoints or meridian-related areas on the body surface. Meridian thermal or infrared detection usually infers the state of the meridians based on thermal sensitivity or infrared temperature distribution. Meridian optical or ultra-weak luminescence detection typically utilizes diffuse reflection, bio-ultra-weak luminescence, and other methods for surface characterization.

[0004] In recent years, biomagnetic imaging technologies such as zero-magnetic or near-zero-magnetic atomic magnetometers and OPM-MEG have developed rapidly, possessing fT / Hz-level sensitivity, no need for liquid helium, proximity to the body surface for placement, and dynamic imaging potential. Existing OPM / SERF biomagnetic imaging technologies are typically used for brain magnetoencephalography, cardiac magnetoencephalography, or general magnetic field reconstruction, mainly addressing noise suppression, magnetic field imaging, or general physiological discrimination issues.

[0005] Currently, the imaging methods for meridians and internal organs can be broadly categorized into two types: The first category consists of TCM meridian detection instruments and infrared visualization systems, which rely on contact conductivity measurement or thermal image fitting. Meridian information comes more from surface physical quantities and regular mapping. The second category consists of wearable magnetoencephalography (MEG) or general magnetic field imaging systems, which can measure extremely weak magnetic fields and complete organ-level reconstruction, but do not unify the topology of body surface meridians, visceral functional areas, bilateral balance, and auxiliary modal admittance changes into a solvable inversion framework.

[0006] Observations have revealed that existing solutions generally suffer from the following drawbacks: Existing meridian resistance or conductivity detection methods are easily affected by factors such as pressure, measurement time, skin moisture, and contact conditions, resulting in insufficient repeatability and standardization.

[0007] Existing meridian thermometry or infrared detection relies on temperature distribution or subjective thermal sensation. Temperature information is not magnetic field information, making it difficult to directly map deep electrical activity.

[0008] The visualization, medical interpretation, and clinical stability of existing meridian optics or ultra-weak luminescence detection methods are still insufficient.

[0009] Existing OPM / SERF biomagnetic imaging technology mainly targets classic biomagnetic tasks such as the brain and heart, and has not yet established a complete technical link between individualized meridian maps, organ coupling operators, and TCM interpretable outputs.

[0010] Existing multimodal processing methods tend to remain at the feature splicing level, lacking technical solutions that utilize auxiliary modes to dynamically adjust the side admittance of the meridian diagram and jointly invert it with the internal functional source and the meridian flow on the body surface. Summary of the Invention

[0011] To address the aforementioned technical problems, embodiments of the present invention provide a zero-magnetic quantum sensing meridian-organ coupling imaging method, comprising: The observation data of the target object is obtained, including the time series data of the three-axis magnetic field of the target object and the auxiliary modal data. The time series data of the three-axis magnetic field is acquired by a zero-magnetic quantum sensing array, and the auxiliary modal data includes the physiological modal data of the target object. Based on the three-dimensional mesh and surface markers of the target object, the preset standard meridian template is individually registered to generate an individualized meridian diagram corresponding to the target object. Construct a functional area diagram of the internal organs and a forward matrix of body sources, and construct a forward matrix of meridian sources in combination with the individualized meridian diagram; Based on the observation data, the side admittance parameters of each side in the individualized meridian diagram are determined, and the time-varying meridian diagram Laplace matrix is ​​constructed by combining the side admittance parameters. Based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram, an observation equation is constructed, and an inversion objective function is constructed by combining the observation data, the organ function area diagram, and the time-varying meridian diagram Laplace matrix. Based on the observation equation and the inversion objective function, the functional source distribution, the meridian flow distribution and the meridian-organ coupling operator in the target object are jointly estimated. Imaging is performed based on the distribution of functional sources within the body, the distribution of meridian flow on the body surface, and the coupling operator between meridians and internal organs.

[0012] In one embodiment, the step of performing individualized registration of a preset standard meridian template based on the three-dimensional mesh and surface markers of the target object to generate an individualized meridian map corresponding to the target object includes: Based on the three-dimensional mesh and surface markers of the target object, a pre-defined standard meridian template is mapped using differentiable deformation:

[0013] in, To be the optimal deformation parameters, As a standard meridian template, It is a three-dimensional mesh of the body surface. It is a set of surface landmarks. For differentiable deformation mapping, This is a projection of the standard meridian template onto the corresponding landmarks on the body surface. Controlling deformation smoothness, To constrain the length of the meridian sides from unreasonable expansion and contraction. and These are the regularization weight coefficients; Based on the differentiable deformation mapping results, an individualized meridian diagram corresponding to the target object is generated.

[0014] In one embodiment, the construction of the visceral functional area map and the body source forward matrix, and the construction of the meridian line source forward matrix in conjunction with the individualized meridian map, includes: The functional area diagram of the internal organs is constructed based on prior knowledge base data; The body source forward matrix is ​​constructed based on a pre-constructed quasi-static biomagnetic field model; The meridian source forward matrix is ​​constructed based on the individualized meridian parameters in the individualized meridian diagram.

[0015] In one embodiment, determining the marginal admittance parameters of each side in the individualized meridian diagram based on the observation data, and constructing the time-varying meridian diagram Laplace matrix in combination with the marginal admittance parameters, includes: Based on the auxiliary modal data in the observed data and the following formula, the temporal admittance parameters of each side in the individualized meridian diagram are jointly determined:

[0016] in, For the first Edge at time Admittance, The correlation matrix of the individualized meridian diagram, () is the activation function. It is the first e The transpose of the trainable parameters corresponding to the stripe edges. It is the first e The calibrable parameters corresponding to the stripe edges; Combine the aforementioned marginal admittance parameters with the following formula to construct the time-varying meridian diagram Laplace matrix:

[0017] in, Let be the Laplace matrix of the time-varying meridian diagram, and diag(w(t)) be a diagonal matrix.

[0018] In one embodiment, constructing the observation equation based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram includes: Based on the aforementioned observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram, the observation equation is constructed as follows:

[0019] in, This represents the current coefficient of the functional region within the body. The coefficient of flow of qi along the meridians on the body surface. For the sensor to observe the magnetic field vector, For the source forward matrix, For the forward matrix of the meridian source, For individualized meridian parameters in individualized meridian diagrams, For the disturbance term, This is the noise term.

[0020] In one embodiment, an inversion objective function is constructed by combining the observation data, the organ function zone diagram, and the Laplace matrix of the time-varying meridian diagram, including: Combining the aforementioned observation data, the organ function area diagram, the time-varying meridian diagram, and the Laplace matrix, an inversion objective function is constructed:

[0021] in, This represents the sensor's observed magnetic field vector after denoising and whitening. The Laplace matrix of the visceral functional area diagram. This refers to the net flow of blood at the nodes along the meridians. Surface sources and sinks induced by functional activities within the body. 'This is a set of left and right meridian pairs with the same name.' Used for directional consistency correction. For allowable delay compensation, X Indicates the distribution of functional sources in the body. Q Indicates the distribution of meridian flow on the body surface. T This represents the meridian-organ coupling operator. , , , , , These are the regularization weight coefficients. These are the weighting coefficients. It is the first e The meridian is at time t The flow intensity value, Let e ​​be the other edge in the edge pair. for t The coefficient of meridian flow on the body surface at any given time.

[0022] In one embodiment, the imaging based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the meridian-organ coupling operator includes: Imaging is performed based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the meridian-viscera coupling operator. The imaging results include at least one of the following: meridian flow map, viscera coupling map, bilateral balance index map, continuity index map, and confidence map.

[0023] In one embodiment, the method further includes: Empty field data of the measurement environment were collected under conditions without test subjects, and a baseline noise template and / or whitening matrix were established based on the empty field data; The time-series data of the body surface triaxial magnetic field are subjected to reference channel regression processing by combining the baseline noise template and / or whitening matrix; Covariance whitening, frequency band selection, and artifact suppression are performed on the time series data of the triaxial magnetic field of the body surface after reference regression processing to obtain denoised and whitened time series data of the triaxial magnetic field of the body surface.

[0024] In one embodiment, obtaining auxiliary modal data includes: Obtain at least one of the following data from the target object: thermal imaging, pulse, electrocardiogram, and respiration.

[0025] Another embodiment of the present invention also provides a zero-magnetic quantum sensing meridian-organ coupling imaging device, comprising: The acquisition module is used to acquire observation data of the target object. The observation data includes time-series data of the three-axis magnetic field on the body surface of the target object and auxiliary modal data. The time-series data of the three-axis magnetic field is acquired by a zero-magnetic quantum sensing array. The auxiliary modal data includes physiological modal data of the target object. The registration module is used to perform individualized registration of a preset standard meridian template based on the three-dimensional mesh and surface markers of the target object, and generate an individualized meridian diagram corresponding to the target object. The first construction module is used to construct the organ function area diagram and the body source forward matrix, and to construct the meridian line source forward matrix in combination with the individualized meridian diagram; The determination module is used to determine the side admittance parameters of each side in the individualized meridian diagram based on the observation data, and to construct the time-varying meridian diagram Laplace matrix in combination with the side admittance parameters; The second construction module is used to construct an observation equation based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram, and to construct an inversion objective function by combining the observation data, the organ function area diagram, and the time-varying meridian diagram Laplace matrix. The estimation module is used to jointly estimate the functional source distribution, meridian flow distribution, and meridian-organ coupling operator within the target object based on the observation equation and the inversion objective function. The imaging module is used to perform imaging based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the coupling operator of meridians and internal organs.

[0026] Based on the above-disclosed embodiments, this invention acquires time-series data of the three-axis magnetic field on the body surface and auxiliary physiological modal data through a zero-magnetic quantum sensing array. It then combines this with individualized meridian diagrams, visceral functional area diagrams, body source forward matrices, and meridian line source forward matrices generated through individualized registration to construct a time-varying meridian diagram Laplace matrix and observation equations. Finally, it estimates the distribution of functional sources within the body, the distribution of meridian flow on the body surface, and the meridian-visceral coupling operator through the observation equations and the inversion objective function. This ultimately achieves imaging based on the distribution of functional sources within the body, the distribution of meridian flow on the body surface, and the meridian-visceral coupling operator. This allows for the acquisition of high-temporal-resolution dynamic imaging results without contact with the target object, effectively suppressing environmental and physiological noise interference, and significantly improving the system robustness, data repeatability, and physical interpretability of the output results during the imaging process.

[0027] Other features and advantages of this application will be set forth in the following description. The objectives and other advantages of this application can be realized and obtained through the structures particularly pointed out in the written description and drawings.

[0028] The technical solution of this application will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0029] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0030] Figure 1 This is a schematic flowchart of the zero-magnetic quantum sensing meridian-organ coupling imaging method in an embodiment of the present invention.

[0031] Figure 2 This is a flowchart illustrating the zero-magnetic quantum sensing meridian-organ coupling imaging method in another embodiment of the present invention.

[0032] Figure 3 This is a structural block diagram of the zero-magnetic quantum sensing meridian and viscera coupling imaging device in an embodiment of the present invention. Detailed Implementation

[0033] The specific embodiments of the present invention will now be described in detail with reference to the accompanying drawings, but these are not intended to limit the scope of the invention.

[0034] It should be understood that various modifications can be made to the embodiments disclosed herein. Therefore, the following description should not be considered as limiting, but merely as an example of embodiments. Other modifications within the scope of this disclosure will be apparent to those skilled in the art.

[0035] The accompanying drawings, which are included in and form part of this specification, illustrate embodiments of the present disclosure and, together with the general description of the disclosure given above and the detailed description of the embodiments given below, serve to explain the principles of the disclosure.

[0036] These and other features of the invention will become apparent from the following description of preferred forms of embodiments given as non-limiting examples, with reference to the accompanying drawings.

[0037] It should also be understood that although the invention has been described with reference to some specific examples, those skilled in the art can certainly implement many other equivalent forms of the invention.

[0038] The above and other aspects, features and advantages of this disclosure will become more apparent when taken in conjunction with the accompanying drawings and in view of the following detailed description.

[0039] Specific embodiments of the present disclosure are described thereafter with reference to the accompanying drawings; however, it should be understood that the disclosed embodiments are merely examples of the present disclosure and can be implemented in various ways. Well-known and / or repeated functions and structures are not described in detail to avoid unnecessary or redundant details that could obscure the present disclosure. Therefore, the specific structural and functional details disclosed herein are not intended to be limiting, but are merely representative of the present disclosure and are intended to teach those skilled in the art to use the present disclosure in a variety of substantially any suitable detailed structures.

[0040] This specification may use the phrases “in one embodiment,” “in another embodiment,” “in yet another embodiment,” or “in still another embodiment,” all of which may refer to one or more of the same or different embodiments according to this disclosure.

[0041] The embodiments of the present invention will now be described in detail with reference to the accompanying drawings.

[0042] In view of the shortcomings of existing solutions, the technical problems to be solved by this application include: First, how to collect the ultra-weak magnetic field on the surface of a target object in a non-contact, low-disturbance manner under zero or near-zero magnetic environment, and obtain high temporal resolution data that can be used for dynamic imaging.

[0043] Second, how to transform the meridians from traditional regular diagrams into a calculable, individualizable, and invertible directed graph structure, so that it is both constrained by surface topology and marker points, and can establish a mathematical coupling relationship with the functional areas of the internal organs.

[0044] Third, how to jointly estimate the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the coupling parameters of meridians and internal organs under the conditions of sparse, noisy, and multi-source interference magnetic field observation, instead of just generating a single magnetic field thermogram or single organ source localization results.

[0045] Fourth, how to output auditable, quantifiable, and repeatable imaging and index results, rather than directly outputting disease diagnosis conclusions or treatment recommendations.

[0046] To solve the above technical problems, such as Figure 1 As shown, this embodiment of the invention provides a zero-magnetic quantum sensing meridian-organ coupling imaging method, comprising: S1: Obtain observation data of the target object, the observation data including the time series data of the three-axis magnetic field of the target object's body surface and auxiliary modal data, the three-axis magnetic field time series data being acquired by a zero-magnetic quantum sensing array, and the auxiliary modal data including the physiological modal data of the target object; S2: Based on the three-dimensional mesh and surface markers of the target object, perform individualized registration on the preset standard meridian template to generate an individualized meridian diagram corresponding to the target object; S3: Construct a functional area diagram of the internal organs and a forward matrix of body sources, and construct a forward matrix of meridian sources in combination with the individualized meridian diagram; S4: Based on the observation data, determine the side admittance parameters of each side in the individualized meridian diagram, and construct the time-varying meridian diagram Laplace matrix in combination with the side admittance parameters; S5: Based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram, construct the observation equation, and combine it with the observation data, the organ function area diagram, and the time-varying meridian diagram Laplace matrix to construct the inversion objective function; S6: Based on the observation equation and the inversion objective function, jointly estimate the functional source distribution, the meridian flow distribution and the meridian-organ coupling operator within the target object; S7: Imaging is performed based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the coupling operator of meridians and internal organs.

[0047] For example, the system first obtains observational data of the target object (user). This data includes time-series data of the three-dimensional magnetic field around the target object's body surface, collected by a quantum sensing array with zero or near-zero magnetic flux, and at least one auxiliary modal data, which is the physiological modal data of the facial object, such as heart rate, body temperature, etc. Next, based on the pre-configured three-dimensional mesh and surface markers of the target object, a pre-defined standard meridian template is individually registered to obtain a personalized meridian map for the target object. Simultaneously, the system constructs a functional area map of internal organs and a forward matrix of body sources. Specifically, this can be achieved by combining standard data and historical detection data of human internal organ functions to construct the functional area map of internal organs and the forward matrix of body sources, as well as by combining the personalized meridian map to construct the forward matrix of meridian lines. Then, based on the observational data, the marginal admittance parameters of each side in the personalized meridian map are determined, and the time-varying meridian map Laplace matrix is ​​constructed using these marginal admittance parameters. After constructing the aforementioned data, the system can build an observation equation based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram. It then combines the observation data, the organ function area diagram, and the Laplace matrix of the time-varying meridian diagram to construct an inversion objective function. Subsequently, the system jointly estimates the functional source distribution within the target object, the meridian flow distribution on the body surface, and the meridian-organ coupling operator based on the observation equation and the inversion objective function. Finally, the system can perform imaging processing based on the solved functional source distribution within the target object, the meridian flow distribution on the body surface, and the meridian-organ coupling operator to obtain the required meridian-related images.

[0048] The method of this embodiment can acquire high temporal resolution dynamic imaging results without contacting the target object, effectively suppress environmental and physiological noise interference, and significantly improve the system robustness, data repeatability and physical interpretability of the output results of the imaging process.

[0049] In one embodiment, the system acquires the time-series data of the three-axis magnetic field on the surface of the target object collected by a multi-channel quantum magnetic sensing array under zero or near-zero magnetic conditions, and simultaneously acquires at least one auxiliary mode data.

[0050] Among them, obtaining auxiliary modal data includes: S101: Obtain at least one of the following data from the target object: thermal image, pulse, electrocardiogram, and respiration.

[0051] In application, data acquisition can be performed using, but is not limited to, 16 to 64 triaxial quantum magnetic sensing channels with a magnetic field sampling rate of 500 Hz to 2000 Hz. The sampler for the auxiliary modes is not specifically limited, but its sampling rate can be, but is not limited to, 30 Hz to 1000 Hz; for PPG / ECG, a sampling rate of 256 Hz to 1000 Hz is preferred, and for thermal imaging, a sampling rate of 5 Hz to 30 Hz is preferred.

[0052] After obtaining the observation data, such as Figure 2 As shown, the method further includes: S8: Collect field data of the measurement environment without test subjects, and establish a baseline noise template and / or whitening matrix based on the field data; S9: Perform reference channel regression processing on the time series data of the body surface triaxial magnetic field by combining the baseline noise template and / or whitening matrix; S10: Perform covariance whitening, frequency band selection, and artifact suppression on the time series data of the triaxial magnetic field of the body surface after reference regression processing to obtain the denoised and whitened time series data of the triaxial magnetic field of the body surface.

[0053] For example, after obtaining the observation data, the system will also perform noise modeling and preprocessing on the triaxial magnetic field time series data of the body surface, including: (1) collecting the empty field data of the measurement environment under the condition of no subject, and establishing a baseline noise template and / or whitening matrix accordingly; (2) performing reference channel regression on the formally collected triaxial magnetic field time series data to suppress residual magnetic field, respiration, heartbeat or other reference disturbance terms; (3) performing covariance whitening, frequency band selection and artifact suppression on the magnetic field data after reference regression to obtain the denoised and whitened observation data.

[0054] In one embodiment, the step of performing individualized registration of a preset standard meridian template based on the three-dimensional mesh and surface markers of the target object to generate an individualized meridian map corresponding to the target object includes: S201: Based on the three-dimensional mesh and surface markers of the target object, perform differentiable deformation mapping on the preset standard meridian template:

[0055] in, To be the optimal deformation parameters, As a standard meridian template, It is a three-dimensional mesh of the body surface. It is a set of surface landmarks. For differentiable deformation mapping, This is a projection of the standard meridian template onto the corresponding landmarks on the body surface. Controlling deformation smoothness, To prevent unreasonable expansion and contraction of the meridian sides. and These are the regularization weight coefficients; S202: Generate an individualized meridian diagram corresponding to the target object based on the differentiable deformation mapping result.

[0056] In other words, in this embodiment, the system completes the individualized alignment of the three-dimensional mesh of the body surface, the body surface landmarks, and the standard meridian template through the aforementioned differential deformation, thereby generating an individualized meridian diagram.

[0057] After completing the construction of the individualized meridian map, the system will construct a functional area map of the internal organs and a forward matrix of body sources, and combine the individualized meridian map to construct a forward matrix of meridian line sources, including: S301: Construct the functional area diagram of the internal organs based on prior knowledge base data; S302: Construct the body source forward matrix based on the pre-constructed quasi-static biomagnetic field model; S303: Construct the meridian source forward matrix based on the individualized meridian parameters in the individualized meridian diagram.

[0058] The quasi-static biomagnetic field model is a physical-mathematical mapping relationship used to describe how the magnetic fields generated by bioelectric current sources (internal organs) and equivalent line sources (meridians) on the body surface propagate to the body surface and are measured by a quantum magnetic sensor array. For example, the quasi-static biomagnetic field model can be a forward physical model that uses the Biot-Savart law to linearly map the spatial distribution of "internal organ sources" and "meridian line sources" to "body surface sensor measurements." This model does not depend on any physiological signals from the target object and can be pre-calculated offline and embedded in the system in one go.

[0059] In one embodiment, determining the marginal admittance parameters of each side in the individualized meridian diagram based on the observation data, and constructing the time-varying meridian diagram Laplace matrix in combination with the marginal admittance parameters, includes: S401: Based on the auxiliary modal data in the observed data and the following formula, jointly determine the temporal admittance parameters of each side in the individualized meridian diagram:

[0060] in, For the first Edge at time Admittance, The correlation matrix of the individualized meridian diagram, () is the activation function. It is the first e The transpose of the trainable parameters corresponding to the stripe edges. It is the first e The calibrable parameters corresponding to the stripe edges; S402: Combine the aforementioned side admittance parameters with the following formula to construct the time-varying meridian diagram Laplace matrix:

[0061] in, Let be the Laplace matrix of the time-varying meridian diagram, and diag(w(t)) be a diagonal matrix.

[0062] In this step, auxiliary modes are used to change the propagation conditions of the meridian diagram, thereby avoiding scale inconsistencies and spurious correlations between modes.

[0063] Furthermore, the construction of the observation equation based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram includes: S501: Construct observation equations based on the aforementioned observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram:

[0064] in, This represents the current coefficient of the functional region within the body. The coefficient of flow of qi along the meridians on the body surface. For the sensor to observe the magnetic field vector, For the source forward matrix, For the forward matrix of the meridian source, For individualized meridian parameters in individualized meridian diagrams, For the disturbance term, This is the noise term.

[0065] The core of the above observation equation calculation is not to treat meridians as mysterious entities, but to abstract meridian paths into a mailbox coefficient line flow graph constrained by the topology of the body surface and marker points.

[0066] When the continuous source is approximated by discrete volume elements and discrete meridian edges, the above observation equations can be obtained from the low-frequency biomagnetic quasi-static approximation:

[0067] In the above formula This refers to the body's origin region. For the first Along the meridian, Let be the tangent vector of this edge. For the first i A quantum magnetic sensor at time t The magnetic field strength vector measured at a spatial location, i.e. The physical prototype. The permeability of free space, It is the current density vector. The coordinates of the source point within the body, It is a distance vector. The total number of meridian sides. To limit the injection intensity, Let be the arc length parameter along the meridian edge. The coordinates of the source point on the body surface. It is a distance vector. For the first iThe total noise at each sensor. The preferred boundary condition is: the volume source satisfies... Meridian flow satisfies approximate conservation at internal nodes; terminal nodes are only allowed sparse source and sink.

[0068] Based on the above formula, the system will combine the observed data, the organ function area diagram, the time-varying meridian diagram, and the Laplace matrix to construct the inversion objective function, including: S502: Construct the inversion objective function by combining the aforementioned observation data, the organ function area diagram, the time-varying meridian diagram, and the Laplace matrix of the meridian diagram:

[0069] in, This represents the sensor's observed magnetic field vector after denoising and whitening. The Laplace matrix of the visceral functional area diagram. This refers to the net flow of blood at the nodes along the meridians. Surface sources and sinks induced by functional activities within the body. 'This is a set of left and right meridian pairs with the same name.' Used for directional consistency correction. For allowable delay compensation, X Indicates the distribution of functional sources in the body. Q Indicates the distribution of meridian flow on the body surface. T This represents the meridian-organ coupling operator. , , , , , These are the regularization weight coefficients. These are the weighting coefficients. It is the first e The meridian is at time t The flow intensity value, Let e ​​be the other edge in the edge pair. for t The coefficient of meridian flow on the body surface at any given time.

[0070] Furthermore, the regularization weight coefficients, specifically the suggested generation methods, include: It can be estimated using L-curve or generalized cross-validation; Press Initialization yields the result; The larger value is determined for the body surface area where there is strong left-right consistency; The search is conducted within 0–3 analysis frames to determine the outcome. sparsity is determined by Control, preferably keeping each organ's functional area coupled with only a few meridian branches.

[0071] in, Specifically, it is the penalty weight for the change in the distribution of meridian flow on the body surface between two adjacent time points. It increases when the target object is in a resting or slow-change state and decreases when it is in an induced state.

[0072] Furthermore, the data fidelity term in the above objective function corresponds to the maximum likelihood under Gaussian noise; and Ensure local smoothing on both the organ diagram and the meridian diagram; The objective function is forced to have a small number of sources and sinks, improving coupling interpretability; a two-sided balancing term is used to suppress mirror inconsistencies caused by accidental noise; and a total time variation term is used to capture local abrupt changes under slow variations. With this structure, the objective function can be solved offline with high precision, or online with approximate solutions.

[0073] In one embodiment, after obtaining the observation equation and the objective function, the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the coupling operator of meridians and internal organs can be obtained by jointly solving the two. Specifically, the coordinate descent method can be used for iterative solution.

[0074] If a hybrid solution using sparse forward matrix and conjugate gradient / ADMM is employed, the complexity of a single iteration can be approximated as:

[0075] Used to describe the asymptotic upper bound of the computational complexity of an algorithm. The number of non-zero elements. The number of discrete elements in the functional areas of the internal organs. For coupling operators The target sparsity. The online version can be fixed. Narrowing the time window and pre-calculating and This further compresses the complexity. Because the algorithm reduces repeated transmission of the original full window and prioritizes the transmission of sparse representations and exponential results, it also has the additional technical effect of "improving the internal performance of the computer system."

[0076] In one embodiment, the imaging based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the meridian-organ coupling operator includes: S701: Imaging is performed based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the meridian-viscera coupling operator. The imaging results include at least one of the following: meridian flow map, viscera coupling map, bilateral balance index map, continuity index map, and confidence map.

[0077] In practical applications, the execution process of this application's solution includes: 1. Start the quantum sensing array in an active magnetic compensation or magnetic shielding environment to complete the empty field acquisition and establish a noise template.

[0078] 2. Collect three-dimensional grids and landmarks on the body surface of the test subjects, and automatically generate individualized meridian maps by combining them with standard meridian templates.

[0079] 3. Acquire zero or near-zero magnetic multi-channel magnetic fields, and simultaneously acquire selectable auxiliary modes.

[0080] 4. Write the multimodal data into HDF5 or Zarr according to a unified timeline and a unified packet format.

[0081] 5. Perform reference regression, whitening, frequency band selection, and artifact suppression on the magnetic field.

[0082] 6. Construct a time-varying meridian Laplace matrix by estimating the dynamic graph admittance using auxiliary modal estimation.

[0083] 7. Run the joint inversion of the observation equation and the inversion objective function to obtain the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the coupling operator of meridians and viscera.

[0084] 8. Calculate the continuity index, the two-sided balance index, the coupling index, and the confidence plot.

[0085] 9. Generate 3D visualizations and structured reports; if used for scientific research, export audit logs and model versions.

[0086] 10. If used for clinical software translation, freeze the model version and proceed with the validation or registration process.

[0087] Some of the steps can be replaced by the following alternatives: Alternative Option 1: Use only the zero magnetic field without auxiliary modes. This method results in a simpler system that better highlights the quantum biomagnetic core, but the dynamic plot admittance is unobservable, and the inversion stability is relatively weak.

[0088] Alternative Option 2: Combining zero magnetic field with thermal imaging. This method provides richer surface information, which is beneficial for admittance modeling and visualization, but the thermal imaging frame rate is low and easily affected by ambient temperature.

[0089] Alternative Option 3: Use a combination of zero magnetic field with PPG, ECG, or respiration. This method offers better dynamic synchronization and facilitates the construction of time-varying admittance, but it provides less local spatial information on the body surface compared to thermal imaging.

[0090] Alternative Option 4: Use a magnetically shielded room. This method has lower noise levels, but is more expensive and requires heavier deployment.

[0091] Alternative Option 5: Employ an active compensation approach. This method is easier to implement clinically and has lower costs, but the environmental model is more complex and requires more sophisticated algorithms.

[0092] The above alternatives reflect the differences in implementation methods / means. Even if the above alternatives are adopted, the core of the overall method remains unchanged, which is still "zero magnetic quantum sensing + individualized meridian map + dual-map admittance coupling inversion + technical image output".

[0093] like Figure 3 As shown, another embodiment of the present invention also provides a zero-magnetic quantum sensing meridian-organ coupling imaging device, comprising: The acquisition module is used to acquire observation data of the target object. The observation data includes time-series data of the three-axis magnetic field on the body surface of the target object and auxiliary modal data. The time-series data of the three-axis magnetic field is acquired by a zero-magnetic quantum sensing array. The auxiliary modal data includes physiological modal data of the target object. The registration module is used to perform individualized registration of a preset standard meridian template based on the three-dimensional mesh and surface markers of the target object, and generate an individualized meridian diagram corresponding to the target object. The first construction module is used to construct the organ function area diagram and the body source forward matrix, and to construct the meridian line source forward matrix in combination with the individualized meridian diagram; The determination module is used to determine the side admittance parameters of each side in the individualized meridian diagram based on the observation data, and to construct the time-varying meridian diagram Laplace matrix in combination with the side admittance parameters; The second construction module is used to construct an observation equation based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram, and to construct an inversion objective function by combining the observation data, the organ function area diagram, and the time-varying meridian diagram Laplace matrix. The estimation module is used to jointly estimate the functional source distribution, meridian flow distribution, and meridian-organ coupling operator within the target object based on the observation equation and the inversion objective function. The imaging module is used to perform imaging based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the coupling operator of meridians and internal organs.

[0094] In one embodiment, the step of performing individualized registration of a preset standard meridian template based on the three-dimensional mesh and surface markers of the target object to generate an individualized meridian map corresponding to the target object includes: Based on the three-dimensional mesh and surface markers of the target object, a pre-defined standard meridian template is mapped using differentiable deformation:

[0095] in, To be the optimal deformation parameters, As a standard meridian template, It is a three-dimensional mesh of the body surface. It is a set of surface landmarks. For differentiable deformation mapping, This is a projection of the standard meridian template onto the corresponding landmarks on the body surface. Controlling deformation smoothness, To constrain the length of the meridian sides from unreasonable expansion and contraction. and These are the regularization weight coefficients; Based on the differentiable deformation mapping results, an individualized meridian diagram corresponding to the target object is generated.

[0096] In one embodiment, the construction of the visceral functional area map and the body source forward matrix, and the construction of the meridian line source forward matrix in conjunction with the individualized meridian map, includes: The functional area diagram of the internal organs is constructed based on prior knowledge base data; The body source forward matrix is ​​constructed based on a pre-constructed quasi-static biomagnetic field model; The meridian source forward matrix is ​​constructed based on the individualized meridian parameters in the individualized meridian diagram.

[0097] In one embodiment, determining the marginal admittance parameters of each side in the individualized meridian diagram based on the observation data, and constructing the time-varying meridian diagram Laplace matrix in combination with the marginal admittance parameters, includes: Based on the auxiliary modal data in the observed data and the following formula, the temporal admittance parameters of each side in the individualized meridian diagram are jointly determined:

[0098] in, For the first Edge at time Admittance, The correlation matrix of the individualized meridian diagram, () is the activation function. It is the first e The transpose of the trainable parameters corresponding to the stripe edges. It is the first e The calibrable parameters corresponding to the stripe edges; Combine the aforementioned marginal admittance parameters with the following formula to construct the time-varying meridian diagram Laplace matrix:

[0099] in, Let be the Laplace matrix of the time-varying meridian diagram, and diag(w(t)) be a diagonal matrix.

[0100] In one embodiment, constructing the observation equation based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram includes: Based on the aforementioned observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram, the observation equation is constructed as follows:

[0101] in, This represents the current coefficient of the functional region within the body. The coefficient of flow of qi along the meridians on the body surface. For the sensor to observe the magnetic field vector, For the source forward matrix, For the forward matrix of the meridian source, For individualized meridian parameters in individualized meridian diagrams, For the disturbance term, This is the noise term.

[0102] In one embodiment, an inversion objective function is constructed by combining the observation data, the organ function zone diagram, and the Laplace matrix of the time-varying meridian diagram, including: Combining the aforementioned observation data, the organ function area diagram, the time-varying meridian diagram, and the Laplace matrix, an inversion objective function is constructed:

[0103] in, This represents the sensor's observed magnetic field vector after denoising and whitening. The Laplace matrix of the visceral functional area diagram. This refers to the net flow of blood at the nodes along the meridians. Surface sources and sinks induced by functional activities within the body. 'This is a set of left and right meridian pairs with the same name.' Used for directional consistency correction. For allowable delay compensation, X Indicates the distribution of functional sources in the body. Q Indicates the distribution of meridian flow on the body surface. T This represents the meridian-organ coupling operator. , , , , , These are the regularization weight coefficients. These are the weighting coefficients. It is the first e The meridian is at time t The flow intensity value, Let e ​​be the other edge in the edge pair. for tThe coefficient of meridian flow on the body surface at any given time.

[0104] In one embodiment, the imaging based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the meridian-organ coupling operator includes: Imaging is performed based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the meridian-viscera coupling operator. The imaging results include at least one of the following: meridian flow map, viscera coupling map, bilateral balance index map, continuity index map, and confidence map.

[0105] In one embodiment, the device further includes: The acquisition module is used to acquire empty field data of the measurement environment under the condition of no test subject, and to establish a baseline noise template and / or whitening matrix based on the empty field data; The first processing module is used to perform reference channel regression processing on the body surface triaxial magnetic field time series data by combining the baseline noise template and / or whitening matrix; The second processing module is used to perform covariance whitening, frequency band selection, and artifact suppression on the time series data of the body surface triaxial magnetic field after reference regression processing, so as to obtain the denoised and whitened time series data of the body surface triaxial magnetic field.

[0106] In one embodiment, obtaining auxiliary modal data includes: Obtain at least one of the following data from the target object: thermal imaging, pulse, electrocardiogram, and respiration.

[0107] Another embodiment of the present invention also provides an electronic device, comprising: One or more processors; Memory, configured to store one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors implement the zero-magnetic quantum sensing meridian-visceral coupling imaging method as described above.

[0108] Furthermore, one embodiment of the present invention also provides a storage medium storing a computer program, which, when executed by a processor, implements the zero-magnetic quantum sensing meridian-organ coupling imaging method as described above. It should be understood that the various solutions in this embodiment have the corresponding technical effects in the above-described method embodiments, and will not be repeated here.

[0109] Furthermore, embodiments of the present invention also provide a computer program product, which is tangibly stored on a computer-readable medium and includes computer-readable instructions that, when executed, cause at least one processor to perform a zero-magnetic quantum sensing meridian-visceral coupling imaging method as described in the embodiments above.

[0110] It should be noted that the computer storage medium of the present invention can be a computer-readable signal medium or a computer-readable storage medium, or any combination thereof. A computer-readable medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, system, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, a random access storage medium (RAM), a read-only storage medium (ROM), an erasable programmable read-only storage medium (EPROM or flash memory), an optical fiber, a portable compact disk read-only storage medium (CD-ROM), an optical storage medium, a magnetic storage medium, or any suitable combination thereof. In the present invention, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, system, or device. In the present invention, a computer-readable signal medium can include a data signal propagated in baseband or as part of a carrier wave, carrying computer-readable program code. Such propagated data signals can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. Computer-readable signal media can also be any computer-readable medium other than computer-readable storage media, which can send, propagate, or transmit a program configured for use by or in connection with an instruction execution system, system, or device. The program code contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to: wireless, antenna, optical fiber, RF, etc., or any suitable combination thereof.

[0111] Furthermore, those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage and optical storage) containing computer-usable program code.

[0112] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1One or more processes and / or boxes Figure 1 A system that specifies functions in one or more boxes.

[0113] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including an instruction set implemented in a process. Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0114] Those skilled in the art should understand that the discussion of any of the above embodiments is merely exemplary and is not intended to imply that the scope of protection of this application is limited to these examples; within the framework of this application, the technical features of the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other variations of different aspects of one or more embodiments of this application as described above, which are not provided in detail for the sake of brevity.

Claims

1. A zero-magnetic quantum sensing method for coupled imaging of meridians and internal organs, characterized in that, include: The observation data of the target object is obtained, including the time series data of the three-axis magnetic field of the target object and the auxiliary modal data. The time series data of the three-axis magnetic field is acquired by a zero-magnetic quantum sensing array, and the auxiliary modal data includes the physiological modal data of the target object. Based on the three-dimensional mesh and surface markers of the target object, the preset standard meridian template is individually registered to generate an individualized meridian diagram corresponding to the target object. Construct a functional area diagram of the internal organs and a forward matrix of body sources, and construct a forward matrix of meridian sources in combination with the individualized meridian diagram; Based on the observation data, the side admittance parameters of each side in the individualized meridian diagram are determined, and the time-varying meridian diagram Laplace matrix is ​​constructed by combining the side admittance parameters. Based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram, an observation equation is constructed, and an inversion objective function is constructed by combining the observation data, the organ function area diagram, and the time-varying meridian diagram Laplace matrix. Based on the observation equation and the inversion objective function, the functional source distribution, the meridian flow distribution and the meridian-organ coupling operator in the target object are jointly estimated. Imaging is performed based on the distribution of functional sources within the body, the distribution of meridian flow on the body surface, and the coupling operator between meridians and internal organs.

2. The zero-magnetic quantum sensing meridian-organ coupling imaging method according to claim 1, characterized in that, The step of individualizing the preset standard meridian template based on the three-dimensional mesh and surface markers of the target object to generate an individualized meridian map corresponding to the target object includes: Based on the three-dimensional mesh and surface markers of the target object, a pre-defined standard meridian template is mapped using differentiable deformation: in, To be the optimal deformation parameters, As a standard meridian template, It is a three-dimensional mesh of the body surface. It is a set of surface landmarks. For differentiable deformation mapping, This is a projection of the standard meridian template onto the corresponding landmarks on the body surface. Controlling deformation smoothness, To constrain the length of the meridian sides from unreasonable expansion and contraction. and These are the regularization weight coefficients; Based on the differentiable deformation mapping results, an individualized meridian diagram corresponding to the target object is generated.

3. The zero-magnetic quantum sensing meridian-organ coupling imaging method according to claim 1, characterized in that, The construction of the visceral functional area map and the body source forward matrix, and the construction of the meridian line source forward matrix in combination with the individualized meridian map, includes: The functional area diagram of the internal organs is constructed based on prior knowledge base data; The body source forward matrix is ​​constructed based on a pre-constructed quasi-static biomagnetic field model; The meridian source forward matrix is ​​constructed based on the individualized meridian parameters in the individualized meridian diagram.

4. The zero-magnetic quantum sensing meridian-organ coupling imaging method according to claim 1, characterized in that, The process of determining the marginal admittance parameters of each side in the individualized meridian diagram based on the observed data, and constructing the time-varying meridian diagram Laplace matrix using the marginal admittance parameters, includes: Based on the auxiliary modal data in the observed data and the following formula, the temporal admittance parameters of each side in the individualized meridian diagram are jointly determined: in, For the first Edge at time Admittance, The correlation matrix of the individualized meridian diagram, () is the activation function. It is the first e The transpose of the trainable parameters corresponding to the stripe edges. It is the first e The calibrable parameters corresponding to the stripe edges; Combine the aforementioned marginal admittance parameters with the following formula to construct the time-varying meridian diagram Laplace matrix: in, Let be the Laplace matrix of the time-varying meridian diagram, and diag(w(t)) be a diagonal matrix.

5. The zero-magnetic quantum sensing meridian-organ coupling imaging method according to claim 1, characterized in that, The construction of the observation equation based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram includes: Based on the aforementioned observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram, the observation equation is constructed as follows: in, This represents the current coefficient of the functional region within the body. The coefficient of flow of qi along the meridians on the body surface. For the sensor to observe the magnetic field vector, For the source forward matrix, For the forward matrix of the meridian source, For individualized meridian parameters in individualized meridian diagrams, For the disturbance term, This is the noise term.

6. The zero-magnetic quantum sensing meridian-organ coupling imaging method according to claim 1, characterized in that, Combining the aforementioned observation data, the organ function area diagram, and the Laplace matrix of the time-varying meridian diagram, an inversion objective function is constructed, including: Combining the aforementioned observation data, the organ function area diagram, the time-varying meridian diagram, and the Laplace matrix, an inversion objective function is constructed: in, This represents the sensor's observed magnetic field vector after denoising and whitening. The Laplace matrix of the visceral functional area diagram. This refers to the net flow of blood at the nodes along the meridians. Surface sources and sinks induced by functional activities within the body. 'This is a set of left and right meridian pairs with the same name.' Used for directional consistency correction. For allowable delay compensation, X Indicates the distribution of functional sources in the body. Q Indicates the distribution of meridian flow on the body surface. T This represents the meridian-organ coupling operator. , , , , , These are the regularization weight coefficients. These are the weighting coefficients. It is the first e The meridian is at time t The flow intensity value, Let e ​​be the other edge in the edge pair. for t The coefficient of meridian flow on the body surface at any given time.

7. The zero-magnetic quantum sensing meridian-organ coupling imaging method according to claim 6, characterized in that, The imaging based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the coupling operator of meridians and internal organs includes: Imaging is performed based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the meridian-viscera coupling operator. The imaging results include at least one of the following: meridian flow map, viscera coupling map, bilateral balance index map, continuity index map, and confidence map.

8. The zero-magnetic quantum sensing meridian-organ coupling imaging method according to claim 6, characterized in that, The method further includes: Empty field data of the measurement environment were collected under conditions without test subjects, and a baseline noise template and / or whitening matrix were established based on the empty field data; The time-series data of the body surface triaxial magnetic field are subjected to reference channel regression processing by combining the baseline noise template and / or whitening matrix; Covariance whitening, frequency band selection, and artifact suppression are performed on the time series data of the triaxial magnetic field of the body surface after reference regression processing to obtain denoised and whitened time series data of the triaxial magnetic field of the body surface.

9. The zero-magnetic quantum sensing meridian-organ coupling imaging method according to claim 1, characterized in that, Obtain auxiliary modal data, including: Obtain at least one of the following data from the target object: thermal imaging, pulse, electrocardiogram, and respiration.

10. A zero-magnetic quantum sensing meridian-organ coupling imaging device, characterized in that, include: The acquisition module is used to acquire observation data of the target object. The observation data includes time-series data of the three-axis magnetic field on the body surface of the target object and auxiliary modal data. The time-series data of the three-axis magnetic field is acquired by a zero-magnetic quantum sensing array. The auxiliary modal data includes physiological modal data of the target object. The registration module is used to perform individualized registration of a preset standard meridian template based on the three-dimensional mesh and surface markers of the target object, and generate an individualized meridian diagram corresponding to the target object. The first construction module is used to construct the organ function area diagram and the body source forward matrix, and to construct the meridian line source forward matrix in combination with the individualized meridian diagram; The determination module is used to determine the side admittance parameters of each side in the individualized meridian diagram based on the observation data, and to construct the time-varying meridian diagram Laplace matrix in combination with the side admittance parameters; The second construction module is used to construct an observation equation based on the observation data, the body source forward matrix, the meridian line source forward matrix, and the individualized meridian diagram, and to construct an inversion objective function by combining the observation data, the organ function area diagram, and the time-varying meridian diagram Laplace matrix. The estimation module is used to jointly estimate the functional source distribution, meridian flow distribution, and meridian-organ coupling operator within the target object based on the observation equation and the inversion objective function. The imaging module is used to perform imaging based on the distribution of functional sources in the body, the distribution of meridian flow on the body surface, and the coupling operator of meridians and internal organs.

Images