A medical image parameter reconstruction method, system, device and medium
By combining ultrasound and optical probes with a composite probe, an optical scattering model was constructed and the absorption and scattering coefficients were solved using a regression algorithm. This solved the problem of insufficient lesion localization accuracy in DOT technology and achieved high-precision detection of breast tissue lesions.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- SYMBOW MEDICAL TECH
- Filing Date
- 2025-12-17
- Publication Date
- 2026-04-10
AI Technical Summary
The optical spatial resolution defects in DOT technology result in poor lesion localization accuracy, while ultrasound imaging technology has low soft contrast of tissues, making it difficult to distinguish subtle differences in tissues.
By employing a composite probe that combines an ultrasound probe and an optical probe, the location and structural information of suspicious lesions in the human tissue under test are obtained, an optical scattering model is constructed, and the absorption coefficient and scattering coefficient are solved by using a forward analytical model and a regression algorithm to improve the positioning accuracy.
By employing dual positioning with ultrasound and optics, the accuracy of lesion detection in human tissues under test is significantly improved, especially the accuracy of lesion localization in breast tissue.
Smart Images

Figure CN121313119B_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of medical detection, in particular to a medical image parameter reconstruction method, system, device and medium. BACKGROUND
[0002] Diffuse optical tomography (DOT) is a technology for tissue imaging using optical methods. DOT technology has high soft contrast and can provide optical parameter information of tissue, such as absorption coefficient and scattering coefficient. Specifically, it uses the spectral absorption difference of breast tissue under different waveband light irradiation, calculates the optical amplitude and phase data received by the detector after light scattering in the tissue according to the radiation transfer equation, and inversely solves the optical information inside the tissue by constructing the spatial distribution of the breast tissue.
[0003] However, the positioning accuracy of the lesion itself is poor due to the defect of optical spatial resolution in DOT technology, and the soft contrast of the tissue is low in ultrasound imaging technology, which is difficult to distinguish the subtle differences of the tissue. SUMMARY
[0004] Therefore, the present application aims to provide a medical image parameter reconstruction method, system, device and medium to solve the technical problems in the prior art.
[0005] In order to achieve the above-mentioned purpose, the technical scheme adopted by the present application is as follows:
[0006] In a first aspect, the present application provides a medical image parameter reconstruction method, comprising:
[0007] acquiring the suspicious lesion position and the structure information of the suspicious lesion of the human tissue to be measured by the ultrasonic probe in the compound probe, guiding the light source end and the optical detection end of the optical probe in the compound probe to collect optical signals according to the structure information of the suspicious lesion;
[0008] constructing an optical scattering model according to the structure information and the optical signals collected by the light source end and the optical detection end;
[0009] constructing a forward analytical model using the optical scattering model, and obtaining an analytical expression corresponding to the absorption coefficient and the scattering coefficient of the human tissue to be measured using the forward analytical model;
[0010] using a regression algorithm to inversely solve the analytical expression to obtain a solution parameter of the analytical expression;
[0011] solving a target parameter according to a solution relationship, and solving the absorption coefficient and the scattering coefficient of the human tissue to be measured using the target parameter.
[0012] In optional embodiments, the human tissue to be measured is breast tissue, and the step of constructing an optical scattering model according to the optical information collected from the structure information to be measured, the light source end and the optical detection end comprises:
[0013] constructing a spherical wave solution expression of the breast tissue according to the characteristics of the breast tissue;
[0014]
[0015] wherein, is the photon density, is the speed of light, is the point source intensity, is the diffusion coefficient, is the spatial position vector, is the angular modulation frequency of the light source, is the imaginary unit, is the time;
[0016] setting an extrapolated boundary for the breast tissue, and transforming the spherical wave solution expression according to the extrapolated boundary to obtain a Helmholtz equation:
[0017]
[0018] wherein, is the point source intensity, is the distance from the source point to the observation point, is the distance from the mirror source point to the observation point;
[0019] constructing the optical scattering model according to the Helmholtz equation.
[0020] In optional embodiments, the step of constructing the optical scattering model according to the Helmholtz equation comprises:
[0021] constructing a first optical information expression from the light source end to the optical detection end according to the Helmholtz equation:
[0022]
[0023] wherein, is the distance from the light source end to the optical detection end, is the distance from the mirror light source end to the optical detection end;
[0024] constructing a second optical information expression from the light source end to the breast tissue according to the Helmholtz equation:
[0025]
[0026] wherein, is the distance from the light source end to the breast tissue, is the distance from the mirror light source end to the breast tissue;
[0027] constructing the optical scattering model according to the first optical information expression and the second optical information expression:
[0028]
[0029] wherein, is the Green function.
[0030] In an optional embodiment, the method further comprises:
[0031] constructing a weight matrix:
[0032]
[0033] normalizing the weight matrix to obtain a normalized matrix:
[0034]
[0035] constructing a conversion parameter:
[0036]
[0037] converting the optical scattering model according to the normalized matrix and the conversion parameter to obtain a target equation group:
[0038]
[0039] wherein, is the absorption coefficient, is an optical parameter.
[0040] In an optional embodiment, the step of inversely solving the analytical expression by using a regression algorithm to obtain a solution parameter of the analytical expression comprises:
[0041] determining a total error E of the target equation group according to the regression algorithm;
[0042] defining a parameter , a parameter ;
[0043] converting the target equation group according to the parameter A and the parameter q to obtain a simplified equation:
[0044]
[0045] defining a parameter B=AQ, F=EQ, y=Q -1q, transforming the simplified equation according to the parameters B, F and y to obtain the analytical expression:
[0046]
[0047] The parameter corresponding to the minimum value of the analytical expression is taken as the solution parameter.
[0048] In an optional implementation, the step of taking the parameter corresponding to the minimum value of the analytical expression as the solution parameter comprises:
[0049] obtaining a quadratic functional of the analytical expression;
[0050] obtaining a stationary point of the quadratic functional, determining the parameter corresponding to the minimum value of the analytical expression according to the stationary point of the quadratic functional, and taking the parameter corresponding to the minimum value as the solution parameter.
[0051] In a second aspect, an embodiment of the present application provides a medical image parameter reconstruction system, comprising: an electronic device and a composite probe, the electronic device being in communication connection with the composite probe;
[0052] The electronic device is configured to acquire a suspicious lesion position and structural information of a suspicious lesion of a human body tissue to be measured by an ultrasonic probe in the composite probe, and guide a light source end and an optical detection end of an optical probe in the composite probe to collect optical signals according to the structural information of the suspicious lesion.
[0053] An optical scattering model is constructed according to the structural information and the optical signals collected by the light source end and the optical detection end.
[0054] A forward analytical model is constructed by using the optical scattering model, and an analytical expression corresponding to absorption coefficients and scattering coefficients of the human body tissue to be measured is obtained by using the forward analytical model.
[0055] A regression algorithm is used to inversely solve the analytical expression, so as to obtain a solution parameter of the analytical expression.
[0056] A target parameter is solved according to a solution relationship, and the absorption coefficients and the scattering coefficients of the human body tissue to be measured are solved by using the target parameter.
[0057] In a third aspect, an embodiment of the present application provides an electronic device, comprising a processor and a memory, the memory storing machine executable instructions capable of being executed by the processor, and the processor being capable of executing the machine executable instructions to implement the medical image parameter reconstruction method of the first aspect.
[0058] In a fourth aspect, the embodiments of the present application provide a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the medical image parameter reconstruction method in the first aspect.
[0059] The medical image parameter reconstruction method, system, device and medium provided by the embodiments of the present application improve the positioning accuracy of the human tissue to be measured through dual positioning of ultrasound and optics, then construct an optical scattering model according to the positioning result, and calculate the absorption coefficient of the human tissue to be measured through an algorithm, thereby greatly improving the detection accuracy of the lesion in the human tissue.
[0060] In order to make the above objectives, characteristics and advantages of the present application more apparent, clear and easy to understand, the following preferred embodiments are specifically described below, and the accompanying drawings are referred to for detailed description. BRIEF DESCRIPTION OF DRAWINGS
[0061] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the embodiments, and it should be understood that the following drawings only show some embodiments of the present application, and therefore should not be regarded as a limitation to the scope, and for those skilled in the art, other related drawings can also be obtained without creative labor on the basis of these drawings.
[0062] Figure 1 A block schematic diagram of an electronic device provided by the embodiments of the present application is shown;
[0063] Figure 2 A flowchart of a medical image parameter reconstruction method provided by the embodiments of the present application is shown;
[0064] Figure 3 An extrapolation boundary diagram provided by the embodiments of the present application is shown;
[0065] Figure 4 A functional module diagram of a medical image parameter reconstruction system provided by the embodiments of the present application is shown.
[0066] Icon: 100-electronic device; 110-memory; 120-processor; 130-communication module; 400-medical image parameter reconstruction system; 200-composite probe; 300-ultrasound probe; 500-optical probe; 600-light source end; 700-optical detection end. DETAILED DESCRIPTION
[0067] The technical solutions in the embodiments of the present application will be described clearly and completely below in combination with the drawings in the embodiments of the present application, and obviously, the described embodiments are only some of the embodiments of the present application, but not all the embodiments. The components of the embodiments of the present application described and shown in the drawings herein can be arranged and designed in various different configurations.
[0068] Therefore, the following detailed description of the embodiments of the application provided in the accompanying drawings is not intended to limit the scope of the application claimed, but merely represents selected embodiments of the application. Based on the embodiments of the application, all other embodiments obtained by those skilled in the art without creative work fall within the scope of the application.
[0069] It should be noted that the relational terms herein such as first and second and the like are used solely to distinguish one from another entity or action, without necessarily requiring or implying any actual relationship or order between such entities or actions. Moreover, the terms "comprises", "comprising", or any other variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but can include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without more limitations, an element preceded by "comprises... " does not, without more limitations, foreclose the existence of additional identical elements in the process, method, article, or apparatus that comprises the stated elements.
[0070] Please refer to Figure 1 is a block diagram of the electronic device 100. The electronic device 100 includes a memory 110, a processor 120, and a communication module 130. The memory 110, the processor 120, and the communication module 130 are electrically connected to each other directly or indirectly to realize data transmission or interaction. For example, these elements can be electrically connected to each other through one or more communication buses or signal lines.
[0071] The memory 110 is configured to store programs or data. The memory 110 can be, but is not limited to, a random access memory (RAM), a read only memory (ROM), a programmable read only memory (PROM), an erasable programmable read only memory (EPROM), an electrically erasable programmable read only memory (EEPROM), and the like.
[0072] The processor 120 is configured to read / write data or programs stored in the memory 110 and perform corresponding functions.
[0073] The communication module 130 is used to establish a communication connection between the electronic device 100 and other communication terminals through the network, and to send and receive data through the network.
[0074] It should be understood that, Figure 1 The structure shown is only a schematic diagram of the electronic device 100. The electronic device 100 may also include components that are larger than... Figure 1 The more or fewer components shown, or having the same Figure 1 The different configurations shown. Figure 1 The components shown can be implemented using hardware, software, or a combination thereof.
[0075] Please refer to Figure 2 , Figure 2 This is a flowchart illustrating a medical image parameter reconstruction method provided in this embodiment. The method includes:
[0076] S101. Obtain the location and structural information of suspicious lesions in the human tissue to be tested through the ultrasound probe in the composite probe, and guide the light source end and optical detection end of the optical probe in the composite probe to collect optical signals based on the structural information.
[0077] To improve the accuracy of human tissue detection, a composite probe can be used for human tissue localization and detection. The composite probe includes both an ultrasound probe and an optical probe.
[0078] The ultrasound probe is a high-frequency linear array ultrasound probe, and the optical probe consists of a multi-channel optical fiber and a multi-channel bundled optical guide, forming the optical emission path at the light source end and the reception path at the optical detection end. The ultrasound probe and optical probe are arranged in a specific configuration; for example, the ultrasound probe is positioned at the center, and the optical probes can be arranged around it, forming an ultrasound and optical acquisition matrix for acquiring ultrasound location and structural information as well as optical signals from the human tissue. In practical operation, the light source end of the optical probe emits a laser of a specific wavelength, which is reflected by the human tissue and captured by the optical detection end.
[0079] S102. Construct an optical scattering model based on the structural information and the optical signals collected by the light source and the optical detection end.
[0080] One can arbitrarily select a test point in the human tissue to be tested, and then construct a corresponding optical scattering model based on the principle of optical scattering and the spatial relationship between the test point, the light source end and the optical detection end.
[0081] S103. Construct a forward analytical model using the optical scattering model, and obtain analytical expressions for the absorption coefficient and scattering coefficient of the human tissue to be tested using the forward analytical model.
[0082] The forward analysis is that the optical scattering model is forward mode automatic differentiation, the local derivative is forward propagated while the evaluation is performed, so that the gradient in one direction is obtained, so that the analytical expression corresponding to the absorption coefficient and the scattering coefficient is obtained.
[0083] In S104, a solving parameter of the analytical expression is obtained by inversely solving the analytical expression by using a regression algorithm.
[0084] The parameters in the analytical expression of the absorption coefficient and the scattering coefficient are unknown, and cannot be obtained by forward reasoning, so that the regression algorithm is used to inversely solve the analytical expression, so that the solving parameter of the analytical expression is obtained. The regression algorithm includes but is not limited to the least square method, the ridge regression, the LASSO regression and the robust regression.
[0085] In S105, a target parameter is solved according to a solving relationship, and the absorption coefficient and the scattering coefficient of the human tissue to be measured are solved by using the target parameter.
[0086] The solving parameter is iteratively calculated by using iteration, and when the iteration number meets a preset condition or the accuracy meets a preset condition, the parameter at this time is taken as the target parameter, the analytical expression is solved by using the target parameter, and the absorption coefficient and the scattering coefficient of the human tissue to be measured are obtained. The absorption coefficient and the scattering coefficient of the human tissue represent the absorption parameter and the scattering parameter of the optical signal of the optical probe, and according to the difference between the absorption parameter and the scattering parameter, it can be determined whether the human tissue is normal tissue or abnormal tissue.
[0087] The embodiment improves the positioning accuracy of the human tissue to be measured by dual positioning of ultrasound and light, then constructs an optical scattering model according to the positioning result, and calculates the absorption coefficient and the scattering coefficient of the human tissue to be measured by using an algorithm, so that the detection accuracy of the lesion in the human tissue is greatly improved.
[0088] In an embodiment, the human tissue to be measured is breast tissue, and the step of constructing an optical scattering model according to the optical information collected by the light source end and the optical detection end according to the structure information to be measured includes:
[0089] The spherical wave analytical expression of the breast tissue is constructed according to the characteristics of the breast tissue.
[0090] The traditional optical scattering model is:
[0091]
[0092] In the formula, is a diffusion coefficient, is a photon density, is an absorption coefficient, a spatial position vector, a time, a light speed, a light source term, a diffusion term.
[0093] The human breast tissue can be regarded as an infinite homogeneous medium, in which case the traditional optical scattering model corresponds to a spherical wave solution expression:
[0094]
[0095] wherein, a photon density, a point source intensity, a light source angular modulation frequency, a unit imaginary number.
[0096] Please refer to Figure 3 , Figure 3 a schematic diagram of an extrapolated boundary provided by the embodiment.
[0097] The light flux on the extrapolated boundary is zero, and the extrapolated boundary is set for the breast tissue, at this time the spherical wave solution expression is transformed according to the extrapolated boundary to obtain a Helmholtz equation:
[0098]
[0099] wherein, a point source intensity, a distance from a source point to an observation point, a distance from a mirror source point to the observation point; , .
[0100] The optical scattering model is constructed according to the Helmholtz equation.
[0101] According to the geometric relationship in Figure 3 , the Helmholtz equation can be further arranged as:
[0102]
[0103] wherein, a source intensity amplitude, a radial coordinate, a decay constant, a phase constant.
[0104] The amplitude and phase thereof can be expressed by the following linear equation:
[0105]
[0106] is an amplitude, is a phase.
[0107] In an embodiment, the step of constructing the optical scattering model according to the Helmholtz equation comprises:
[0108] constructing a first optical information expression from the light source end to the optical detection end according to the Helmholtz equation:
[0109]
[0110] wherein, is a distance from the light source end to the optical detection end, is a distance from the mirror image light source end to the optical detection end;
[0111] constructing a second optical information expression from the light source end to the breast tissue according to the Helmholtz equation:
[0112]
[0113] wherein, is a distance from the light source end to the breast tissue, is a distance from the mirror image light source end to the breast tissue;
[0114] constructing the optical scattering model according to the first optical information expression and the second optical information expression:
[0115]
[0116] wherein, is a Green function.
[0117] In an embodiment, the method further comprises:
[0118] constructing a weight matrix:
[0119]
[0120] normalizing the weight matrix to obtain a normalized matrix:
[0121]
[0122] constructing a conversion parameter:
[0123]
[0124] converting the optical scattering model according to the normalized matrix and the conversion parameter to obtain a target equation group:
[0125]
[0126] wherein, is the absorption coefficient, is an optical parameter, which can be calculated according to the amplitude and phase in the optical parameter.
[0127] In order to solve the above objective equation group, a regression algorithm can be used for solving, and there are two errors in the regression algorithm solving, one is the measurement error, and the other is the error of the weight matrix W, therefore, the solving of the regression algorithm can be converted into the minimization of the two errors.
[0128] In an embodiment, the step of solving the analytical expression by using the regression algorithm to obtain the solving parameter of the analytical expression comprises:
[0129] determining the total error E of the objective equation group according to the regression algorithm;
[0130] defining a parameter , and a parameter ;
[0131] converting the objective equation group according to the parameter A and the parameter q to obtain a simplified equation:
[0132]
[0133] defining a parameter B=AQ, F=EQ, y=Q -1 converting the simplified equation according to the parameters B, F and y to obtain the analytical expression:
[0134]
[0135] taking the parameter corresponding to the minimum value of the analytical expression as the solving parameter.
[0136] In an embodiment, the step of taking the parameter corresponding to the minimum value of the analytical expression as the solving parameter comprises:
[0137] obtaining a quadratic functional of the analytical expression;
[0138] obtaining a stationary point of the quadratic functional, determining the parameter corresponding to the minimum value of the analytical expression according to the stationary point of the quadratic functional, and taking the parameter corresponding to the minimum value as the solving parameter.
[0139] After obtaining the analytical expression, the solving of the above regression algorithm can be converted into solving the minimum singular value of the matrix satisfying the expression of the to-be-solved parameter .
[0140] taking a quadratic functional:
[0141]
[0142] To find the stationary point of the above functional, we only need to let That is, from the above formula:
[0143]
[0144] Corresponding to:
[0145]
[0146] Using rank(F) = 1, we can give:
[0147]
[0148] Therefore, the to-be-solved
[0149]
[0150] Matrix theory shows that the above formula is equivalent to finding the minimum value of its Rayleigh quotient F(q), where:
[0151]
[0152] Solving minF(q) has the mathematical meaning of finding the eigenvector q corresponding to the minimum eigenvalue of the Hermite matrix ( .
[0153] In order to perform the corresponding steps in the above embodiment and each possible way, an implementation of a medical image parameter reconstruction system is given below, please refer to Figure 4 , Figure 4 A functional module diagram of a medical image parameter reconstruction system provided by an embodiment of the present application. It should be noted that the medical image parameter reconstruction system provided by the present embodiment has the same basic principles and technical effects as the above-mentioned embodiments. For brief description, the part not mentioned in the present embodiment can refer to the corresponding content in the above-mentioned embodiments. The medical image parameter reconstruction system 400 comprises:
[0154] An electronic device 100 and a composite probe 200, the electronic device 100 is in communication connection with the composite probe 200;
[0155] The electronic device 100 is configured to acquire suspicious lesion position and structure information of the suspicious lesion of the human body tissue to be measured by an ultrasonic probe 300 in the composite probe 200, and guide the light source end 600 and the optical detection end 700 of the optical probe 500 in the composite probe 200 to collect optical signals according to the structure information;
[0156] constructing an optical scattering model according to the structure information and optical signals collected by the light source end 600 and the optical detection end 700;
[0157] constructing a forward analytical model by using the optical scattering model, and obtaining an analytical expression corresponding to the absorption coefficient and the scattering coefficient of the human tissue to be measured by using the forward analytical model;
[0158] solving the analytical expression by using a regression algorithm to obtain a solving parameter of the analytical expression;
[0159] solving a target parameter according to a solving relationship, and solving the absorption coefficient and the scattering coefficient of the human tissue to be measured by using the target parameter.
[0160] Optionally, the modules described above can be stored in the memory shown in the form of software or firmware (Firmware) or solidified in an operating system (Operating System, OS) of the electronic device, and can be executed by the processor in the electronic device. Figure 1 The data, program codes and the like required for executing the modules described above can be stored in the memory. Figure 1
[0161] In several embodiments provided in the present application, it should be understood that the disclosed apparatus and method can also be implemented by other manners. The apparatus embodiments described above are only schematic, for example, the flowcharts and block diagrams in the drawings show the possible implementation architectures, functions and operations of the apparatus, method and computer program product according to the embodiments of the present application. In this regard, each block in the flowcharts or block diagrams can represent a module, a program segment or a part of code, which contains one or more executable instructions for implementing the specified logic function. It should also be noted that, in some alternative implementation manners, the functions noted in the blocks can also occur in different orders from those noted in the drawings. For example, two consecutive blocks can actually be executed in parallel, and sometimes they can be executed in reverse order, depending on the functions involved. It should also be noted that each block in the block diagrams and / or flowcharts, and the combination of blocks in the block diagrams and / or flowcharts, can be implemented by a dedicated hardware-based system for implementing the specified functions or actions, or can be implemented by a combination of dedicated hardware and computer instructions.
[0162] In addition, the functional modules in each embodiment of the present application can be integrated together to form an independent part, or each module can exist independently, or two or more modules can be integrated to form an independent part.
[0163] If the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer readable storage medium. Based on this understanding, the technical solutions of the present application or the parts of the technical solutions that essentially contribute to the prior art or the parts of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the method described in the various embodiments of the present application. The aforementioned storage medium includes a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.
[0164] The above only describes the preferred embodiments of the present application and is not intended to limit the present application. For those skilled in the art, the present application can have various modifications and changes. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.
Claims
1. A medical image parameter reconstruction method, characterized by, The method comprises the following steps: obtaining suspicious lesion position and structure information of the suspicious lesion of the human body tissue to be measured by an ultrasonic probe in a composite probe, guiding a light source end and an optical detection end of an optical probe in the composite probe to collect optical signals according to the structure information of the suspicious lesion; constructing an optical scattering model according to the structure information and the optical signals collected by the light source end and the optical detection end; constructing a forward analytical model by using the optical scattering model, and obtaining an analytical expression corresponding to absorption coefficients and scattering coefficients of the human body tissue to be measured by using the forward analytical model; solving the analytical expression by using an iterative algorithm to obtain a solving parameter of the analytical expression; solving a target parameter according to a solving relationship, and solving the absorption coefficients and the scattering coefficients of the human body tissue to be measured by using the target parameter; the human body tissue to be measured is breast tissue, and the step of constructing the optical scattering model according to the structure information and the optical signals collected by the light source end and the optical detection end comprises the following steps: constructing a spherical wave analytical expression of the breast tissue according to characteristics of the breast tissue: wherein is the photon density, is the speed of light, is the point source strength in an infinite homogeneous medium, is the diffusion coefficient, is the spatial position vector, is the angular modulation frequency of the light source, is the imaginary unit, is time; setting an extrapolation boundary for the breast tissue, and transforming the spherical wave analytical expression to obtain a Helmholtz equation according to the extrapolation boundary: wherein, is the equivalent point source strength after extrapolating the boundary, is the distance from the source point to the observation point, is the distance from the mirror source point to the observation point; constructing the optical scattering model according to the Helmholtz equation.
2. The medical image parameter reconstruction method of claim 1, wherein, The step of constructing the optical scattering model according to the Helmholtz equation comprises the following steps: constructing a first optical information expression from the light source end to the optical detection end according to the Helmholtz equation: wherein, is the distance from the light source end to the optical detection end, is the distance from the mirror image light source end to the optical detection end; constructing a second optical information expression from the light source end to the breast tissue according to the Helmholtz equation: wherein, is the distance from the light source end to the breast tissue, is the distance from the mirror light source end to the breast tissue; constructing the optical scattering model according to the first optical information expression and the second optical information expression: wherein G is the Green function.
3. The medical image parameter reconstruction method of claim 2, wherein, The method further comprises the following steps: constructing a weight matrix: normalizing the weight matrix to obtain a normalized matrix: constructing a conversion parameter: converting the optical scattering model according to the normalized matrix and the conversion parameter to obtain a target equation group: wherein is the absorption coefficient, is the optical parameter.
4. The medical image parameter reconstruction method of claim 3, wherein, The step of solving the analytical expression by using an iterative algorithm to obtain a solving parameter of the analytical expression comprises the following steps: determining a total error E of the target equation group according to the iterative algorithm; define parameter , parameter ; converting the target equation group to obtain a simplified equation according to the parameter A and the parameter q: defining parameters B = AQ, F = EQ, y = Q -1 q, transforming said simplified equation according to said parameters B, F and y, obtaining said analytical expression: taking a parameter corresponding to a minimum value of the analytical expression as the solving parameter.
5. The medical image parameter reconstruction method of claim 4, wherein, The step of taking a parameter corresponding to a minimum value of the analytical expression as the solving parameter comprises the following steps: obtaining a quadratic functional of the analytical expression; obtaining a stationary point of the quadratic functional, determining a parameter corresponding to a minimum value of the analytical expression according to the stationary point of the quadratic functional, and taking the parameter corresponding to the minimum value as the solving parameter.
6. An electronic device, comprising: The computer program is executed by a processor to implement the medical image parameter reconstruction method according to any one of claims 1-5.
7. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by a processor to implement the medical image parameter reconstruction method according to any one of claims 1-5.
Citation Information
Patent Citations
Method of medical imaging using combined near infrared diffusive light and ultrasound
US20040215072A1