Myocardial blood flow quantification method of boric acid 99mTc labeled single photon SPECT imaging agent and application of myocardial blood flow quantification method

Abstract

The invention relates to the field of heart nuclear medicine, in particular to a technical method for carrying out myocardial blood flow quantitative analysis through dynamic SPECT imaging of a boric acid 99mTc labeled single-photon SPECT imaging agent, image processing, image correction and myocardial blood flow quantitative calculation, and the technical method is used for generating a dynamic bull's-eye map and a left heart blood pool time activity curve. The method provided by the invention can be used for correcting right heart blood pool overflow interference and myocardial uptake residue interference, quantifying myocardial blood flow and improving fitting operation efficiency, can also be used for establishing a myocardial blood flow volume diagram, and is applied to more intuitively evaluating the resting blood flow, the load blood flow and the high and low areas of the blood flow reserve of the whole myocardial and blood vessel partition. The method provided by the invention can solve the technical problem of myocardial blood flow quantitative analysis by dynamic SPECT imaging of the boric acid 99mTc labeled single-photon SPECT imaging agent and the application problem of the dynamic SPECT imaging of the boric acid 99mTc labeled single-photon SPECT imaging agent.

Description

Technical Field

[0001] The present invention relates to the field of cardiac nuclear medicine, and in particular to the use of single photon emission computed tomography (SPECT) to perform boric acid-based 99m Dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents, a technical method for quantitative analysis of myocardial blood flow through image processing, image correction and kinetic model calculation, and the use of the technical method. Background Art

[0002] In the field of cardiac nuclear medicine, 99m The quantitative application of myocardial blood flow by Tc-labeled myocardial perfusion drugs is easier to popularize and apply. It is now widely used at home and abroad. 99m Tc radionuclide labeled myocardial perfusion drugs (including methoxyisobutylisocyanate and tetrofosmin). In recent years, new drugs with boric acid as the main structure 99m Tc radionuclide-labeled drugs have been developed. This new type of myocardial perfusion drug has the additional characteristics of faster and higher myocardial uptake, which helps to achieve more accurate SPECT myocardial blood flow quantitative analysis, and therefore has a broader practical application prospect.

[0003] However, this new type of myocardial perfusion drug also has certain defects. For example, after the drug enters the body and reaches the right ventricle through the systemic circulation, due to the limitation of SPECT spatial resolution, the radiation spillover from the right heart blood pool to the ventricular septum is more prominent. When the drug reaches the left ventricle through the pulmonary circulation, the radiation spillover from the left heart blood pool to the entire myocardium is also more prominent. When the drug circulates to the deep root of the aorta and enters the coronary bed to reach the myocardium, although it can be quickly taken up by the myocardium through passive diffusion, the drug is more obviously eluted from the myocardium over time within a few minutes after myocardial uptake, resulting in more complex changes in the radiation signals generated by the drug around the myocardium and in the myocardial tissue over time, especially when the actual application requires the completion of resting / load dynamic SPECT imaging on the same day. Before SPECT imaging, a low dose of drug needs to be pre-injected in advance to accurately position the heart before imaging. Before the actual resting drug injection, there is already a certain amount of residual myocardial uptake. The residual myocardial uptake formed by the pre-injection affects the accuracy of myocardial blood flow quantification. Subsequently, resting dynamic SPECT imaging is completed after the actual resting drug injection. The resting myocardial uptake before the stress dynamic SPECT also leads to residual myocardial uptake, which will also affect the accuracy of stress SPECT myocardial blood flow quantification. When the actual application requires the completion of stress / resting dynamic SPECT imaging on the same day, the same problem is encountered. Therefore, even when reconstructing the dynamic SPECT image, the dynamic SPECT raw data can be physically corrected to produce a dynamic SPECT image in a quantitative format. For boric acid 99m The distribution of Tc-labeled single-photon SPECT imaging agents around the myocardium and the characteristics of myocardial uptake are important for the effective implementation of boric acid-based 99m There are still great challenges in the quantification of myocardial blood flow using Tc-labeled single-photon SPECT imaging agents, which hinders the practical application of this type of drug in SPECT myocardial blood flow quantitative analysis.

[0004] Therefore, a solution to boric acid is provided. 99m The method and application of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents for quantitative analysis of myocardial blood flow have important practical significance. Summary of the invention

[0005] In view of this, the present invention provides a method for detecting boric acid by single photon computed tomography (SPECT). 99m Dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents, a technical method for quantitative analysis of myocardial blood flow through image processing, image correction and kinetic model calculation, and the use of the technical method.

[0006] In order to achieve the above-mentioned object of the invention, the present invention provides the following technical solutions:

[0007] In a first aspect, the present invention provides an image processing method for dynamic SPECT imaging of a boric acid nuclide labeled single-photon SPECT imaging agent, comprising the following steps:

[0008] Step 1. Generation of quantitative image format: Dynamic SPECT image reconstruction of boronate radionuclide labeled single photon SPECT imaging agent, after complete physical correction, and the pixel value of the reconstructed image is converted into the quantitative format of radiation concentration, the pixel value unit is (Bq / ml);

[0009] Step 2. Image superposition: superimposing the dynamic SPECT image of the boric acid radionuclide labeled single-photon SPECT imaging agent over a period of time to generate a static SPECT image of the boric acid radionuclide labeled single-photon SPECT imaging agent;

[0010] Step 3. Image segmentation: cutting the static SPECT image with linear planes to distinguish the myocardial area from the non-myocardial area;

[0011] Step 4. Image coordinate system conversion: convert the static SPECT image from a rectangular coordinate system to a spherical coordinate system to generate a myocardial image of the static SPECT image in a spherical coordinate system, and perform an internal difference calculation on the pixel value of each central pixel through the pixel values ​​of adjacent pixels and the corresponding spatial coordinates to ensure the correctness and continuity of the conversion of each pixel value in the spherical coordinate system;

[0012] Step 5. Generation of myocardial midline surface: using ray tracing, obtaining the radial myocardial maximum radiation concentration point from the myocardial region image in the spherical coordinate system to generate a 3D myocardial midline surface;

[0013] Step 6. Generation of the bull's eye geometric transformation matrix: Taking the apex of the heart in the spherical coordinate system, the rectangular coordinate system and the 2D bull's eye coordinate system as the central reference point, a continuous space vector transformation is performed on the 3D myocardial midline curved surface in the radial direction of the spherical coordinate system to obtain the geometric correlation of the radial myocardial maximum radiation concentration point in the 2D bull's eye coordinate system, and used to generate the bull's eye transformation geometric matrix;

[0014] Step 7. Generation of dynamic bull's eye diagram: Applying the bull's eye diagram geometric transformation matrix to boric acid 99m Each frame of dynamic SPECT image of the dynamic SPECT imaging of the Tc-labeled single-photon SPECT imaging agent is used to generate an uncorrected dynamic bull's-eye image of the myocardial area;

[0015] Step 8. Generation of the left heart blood pool time-activity curve: The 3D region of interest is automatically placed between the left atrium and the left ventricle in the static SPECT image, and the dynamic SPECT image is applied to this position to obtain the left heart blood pool time-activity curve.

[0016] In some specific embodiments of the present invention, the dynamic SPECT imaging of the boronic acid radionuclide labeled single-photon SPECT imaging agent includes dynamic SPECT imaging of the resting boronic acid radionuclide labeled single-photon SPECT imaging agent and / or dynamic SPECT imaging of the loaded boronic acid radionuclide labeled single-photon SPECT imaging agent.

[0017] In a second aspect, the present invention also provides an uncorrected resting dynamic bull's eye diagram, a resting left heart blood pool time-activity curve, a load uncorrected load dynamic bull's eye diagram and / or a load left heart blood pool time-activity curve obtained by the image processing technology method as described.

[0018] In a third aspect, the present invention also provides an image correction method for dynamic SPECT imaging of a boric acid nuclide labeled single-photon SPECT imaging agent, comprising the following steps:

[0019] Correction of right heart blood pool overflow: using the maximum radiation concentration in the left heart blood pool time-activity curve obtained by the image processing technology method as described in claim 1 or 2 as the time demarcation point for the drug to enter the left heart, measuring the radiation concentration of the left ventricular lateral wall in the uncorrected dynamic bull's eye map through the region of interest, when the pixel value in the left ventricular septal wall is higher than 1.2 to 3 times the average value of the radiation concentration of the left ventricular septal wall, it is considered to be interference from right heart blood pool overflow, and adjusted to the radiation concentration of the left ventricular lateral wall to obtain a dynamic bull's eye map of right ventricular blood pool overflow correction;

[0020] Correction of pre-injection myocardial uptake residue: for pre-injection myocardial uptake residue, the pre-injection myocardial residual amount of each frame of the dynamic SPECT image is estimated through the myocardial uptake of the first frame of the dynamic SPECT image and the myocardial washout model (the myocardial washout model can be any time-related decreasing function, and the myocardial washout model is preferably: W washout rate (t) = U first frame myocardial uptake × exp (-a × T), where a is a constant), and deducted to generate a right ventricular blood pool overflow correction and pre-injection correction dynamic bull's eye map;

[0021] Correction of resting myocardial uptake residue: For dynamic SPECT imaging of boric acid radionuclide labeled single-photon SPECT imaging agents at rest first and then under load, the myocardial residual degree is estimated from the static bull's eye diagram superimposed with the dynamic bull's eye diagram corrected for right ventricular overflow by using the time difference between the starting points of the two imaging times, the physical half-life of the radionuclide, the physical index model, and the resting myocardial washout model obtained from the time-varying characteristics of the resting dynamic bull's eye diagram, so as to generate a resting myocardial residual bull's eye diagram, and the load right ventricular overflow corrected dynamic bull's eye diagram is deducted from the load myocardial washout model obtained by considering the time-varying characteristics of the load dynamic bull's eye diagram, so as to generate a load dynamic bull's eye diagram corrected for right ventricular overflow and resting myocardial residual;

[0022] Correction of residual myocardial uptake under load: For dynamic SPECT imaging of boric acid radionuclide labeled single-photon SPECT imaging agents with load followed by rest, the residual degree of myocardial uptake under load is estimated from the superposition of the dynamic bull's eye diagram corrected for right ventricular overflow through the time difference between the starting points of the two imaging times, the physical half-life of the radionuclide, the physical index model, and the load myocardial washout model obtained from the time-varying characteristics of the dynamic bull's eye diagram under load, so as to generate a residual bull's eye diagram of myocardial uptake under load, and the dynamic bull's eye diagram corrected for right ventricular overflow is deducted from the resting myocardial washout model obtained by considering the time-varying characteristics of the resting dynamic bull's eye diagram, so as to generate a resting dynamic bull's eye diagram corrected for right ventricular overflow and myocardial residual.

[0023] Generation of fully corrected myocardial time-activity curves: A pair of resting dynamic bull's-eye diagrams corrected for right heart overflow and pre-injection, and a stress dynamic bull's-eye diagram corrected for right heart overflow and resting myocardial residual; or a pair of resting dynamic bull's-eye diagrams corrected for right heart overflow and myocardial residual, and a stress dynamic bull's-eye diagram corrected for right heart overflow and pre-injection; used to generate a set of fully corrected resting myocardial time-activity curves and fully corrected stress myocardial time-activity curves.

[0024] In a fourth aspect, the present invention further provides application of the image correction method in correcting right heart blood pool overflow interference and / or correcting myocardial uptake residual interference.

[0025] In a fifth aspect, the present invention further provides a method for quantitatively calculating myocardial blood flow in dynamic SPECT imaging of a boric acid nuclide labeled single-photon SPECT imaging agent, comprising the following steps:

[0026] Transformation of kinetic model application: A single-tissue two-chamber model and a combination coefficient FBV are used, and the myocardial radiation concentration measured at each time point is regarded as a combination of the true myocardial uptake and the interference of the left heart blood pool on the myocardial overflow, where (1-FBV) is also the partial volume coefficient caused by the limited resolution of SPECT, and FBV is also the overflow of the left heart blood pool on the myocardium, to generate a kinetic component equation with kinetic parameters (K1, k2, FBV);

[0027]

[0028] Calculation of the kinetic parameters of the resting myocardial part: firstly, a set of fixed initial values ​​are used to fit a set of corrected resting myocardial time-activity curves and resting left heart blood pool time-activity curves with the kinetic composition equation to generate the kinetic parameters of the resting whole myocardium as the variable initial values, and then the corrected resting myocardial time-activity curves and resting left heart blood pool time-activity curves of each point in the myocardial part are fitted with the kinetic composition equation to calculate the kinetic parameters of each point in the resting myocardial part, and through the resting Renkin-Crone model, the resting K1 is regarded as a function of myocardial blood flow and converted, and then the resting myocardial blood flow value is corrected by the product of heart rate and blood pressure, and displayed as a bull's eye diagram, thereby generating a resting blood flow bull's eye diagram;

[0029] Calculation of dynamic parameters of loaded myocardial sites: firstly, a set of fixed initial values ​​are used to fit a set of corrected load myocardial time-activity curves and load left heart blood pool time-activity curves with the dynamic composition equation to generate dynamic parameters of the loaded whole myocardium as the variable initial values, and then the corrected load myocardial time-activity curves and load left heart blood pool time-activity curves of each point in the myocardial site are fitted with the dynamic composition equation to calculate the dynamic parameters of each point in the loaded myocardial site; through the load Renkin-Crone model, the load K1 is regarded as a function of myocardial blood flow and converted, and then the load myocardial blood flow value is corrected by the product of heart rate and blood pressure, and displayed as a bull's eye diagram to generate a load blood flow bull's eye diagram; through the rest and load Renkin-Crone models, the rest K1 and load K1 are regarded as functions of myocardial blood flow, and converted to obtain the resting myocardial blood flow value and the load myocardial blood flow value; the load blood flow bull's eye diagram is divided by the resting blood flow bull's eye diagram to obtain the blood flow reserve bull's eye diagram.

[0030] In some specific embodiments of the present invention, in the transformation of the kinetic model application, for the rest and load fitting operations, the k2 of the whole myocardium is used as the k2 value of each point of the myocardium for fitting operations before the fitting operation to improve the efficiency of the fitting operation.

[0031] In a sixth aspect, the present invention also provides the application of the myocardial blood flow quantitative calculation method in myocardial blood flow quantification and / or improving the fitting operation efficiency.

[0032] In the seventh aspect, the present invention also provides a method for constructing a myocardial blood flow volume map, such as constructing a myocardial blood flow volume map using the resting myocardial blood flow value, the load myocardial blood flow value, the blood flow reserve value and the 2D coordinate system obtained by the myocardial blood flow quantitative calculation method.

[0033] In an eighth aspect, the present invention also provides a myocardial blood flow volume map obtained by the construction method.

[0034] In the ninth aspect, the present invention also provides a method for evaluating myocardial blood flow volume, wherein the resting blood flow bull's eye diagram, the load blood flow bull's eye diagram and the blood flow reserve bull's eye diagram obtained by the myocardial blood flow quantitative calculation method are divided into the whole myocardium and three vascular branch areas, the three vascular branch areas include the left anterior descending branch, the left circumflex branch and the right coronary branch, the average value of each indicator in the whole myocardium and each vascular branch area is calculated respectively, and then the blood flow volume area of ​​the whole myocardium or each vascular branch area is evaluated using the myocardial blood flow volume diagram as described.

[0035] In the tenth aspect, the present invention also provides the application of the evaluation method in evaluating the high and low areas of resting blood flow, load blood flow and / or blood flow reserve of the whole myocardium and vascular divisions.

[0036] In some specific embodiments of the present invention, in the image processing method, the image correction method, the myocardial blood flow quantitative calculation method, the construction method and / or the evaluation method, the nuclides include but are not limited to 99m Tc.

[0037] In the eleventh aspect, the present invention also provides a processor, which is used to run a program, wherein the program executes the image processing technology method, the image correction method, the myocardial blood flow quantitative calculation method, the construction method and / or the evaluation method when running.

[0038] In a twelfth aspect, the present invention further provides an electronic device or apparatus, the electronic device comprising: a memory and the processor;

[0039] Wherein, one or more computer programs are stored in the memory, and the one or more computer programs include instructions; when the instructions are executed by the processor, the electronic device executes the image processing technical method, the image correction method, the myocardial blood flow quantitative calculation method, the construction method and / or the evaluation method;

[0040] The device includes one or more components, and the one or more components implement the image processing technical method, the image correction method, the myocardial blood flow quantitative calculation method, the construction method and / or the evaluation method.

[0041] The present invention utilizes boric acid 99m Dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents, and the technical method of quantitative analysis of myocardial blood flow through image processing, image correction and quantitative calculation of myocardial blood flow, using: 1) boric acid 99m Image processing method for dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent, used to generate dynamic bull's eye image and left heart blood pool time activity curve; 2) Boric acid 99m Image correction method for dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents, used to correct right heart blood pool overflow interference and correct myocardial uptake residual interference; 3) Boric acid 99m A quantitative calculation method for myocardial blood flow in dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents, which is used for quantitative myocardial blood flow and improving fitting calculation efficiency; 4) Boric acid 99m The quantitative analysis method of myocardial blood flow by dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent is used to establish a myocardial blood volume map, which is used to more intuitively evaluate the resting blood flow, load blood flow and high and low areas of blood flow reserve of the whole myocardium and vascular divisions. The technical method of the present invention can solve the problem of boric acid 99m Technical difficulties and application issues in quantitative analysis of myocardial blood flow using dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art are briefly introduced below.

[0043] Figure 1 The steps in the image processing method implemented according to the present invention and the results of applying the steps to generate a dynamic bull's eye image;

[0044] Figure 2 The uncorrected left heart blood pool time activity curve according to the present invention is generated by automatically placing a sampling rectangular box between the left ventricle and the left atrium, and is used to determine the time cutoff point of SPECT drug entering the left heart;

[0045] Figure 3 1. The dynamic bull's eye diagram (A) before correction at the time cutoff point according to the present invention, the dynamic bull's eye diagram after right heart blood pool overflow correction (B), and the radiation concentration measured at the left ventricular septum;

[0046] Figure 4 The figure shows the myocardial time activity curve after the correction of the pre-injection myocardial residual according to the present invention;

[0047] Figure 5 The figure shows a resting dynamic bull's eye diagram after right heart overflow correction and myocardial residual correction according to the present invention, a group of 4 resting myocardial time-activity curves consisting of the anterior wall, lateral wall, inferior wall and septal wall; a load dynamic bull's eye diagram after right heart overflow correction and myocardial residual correction, a group of 4 load myocardial time-activity curves consisting of the anterior wall, lateral wall, inferior wall and septal wall;

[0048] Figure 6 The diagrams show the resting K1 bull's eye diagram, the resting k2 bull's eye diagram and the resting FBV bull's eye diagram implemented according to the present invention; the load K1 bull's eye diagram, the load k2 bull's eye diagram and the load FBV bull's eye diagram; and the blood flow reserve bull's eye diagram; wherein A shows the load dynamic bull's eye diagram corrected for right ventricular blood pool overflow and resting myocardial residual volume, B shows 4 complete corrected load myocardial activity curves generated by the load dynamic bull's eye diagram corrected for right ventricular blood pool overflow and resting myocardial residual volume, C shows the resting dynamic bull's eye diagram corrected for right ventricular blood pool overflow and pre-injection correction, and D shows 4 complete corrected resting myocardial activity curves of the resting dynamic bull's eye diagram corrected for right ventricular blood pool overflow and pre-injection correction;

[0049] Figure 7 1 shows a myocardial blood flow volume diagram with 14 levels of regions according to the present invention; wherein A shows a resting blood flow bull's eye diagram, a resting K1 bull's eye diagram, a resting k2 bull's eye diagram and a resting FBV bull's eye diagram; B shows a load blood flow bull's eye diagram, a load K1 bull's eye diagram, a load k2 bull's eye diagram and a resting FBV bull's eye diagram, and C shows a blood flow reserve bull's eye diagram;

[0050] Figure 8 The present invention is implemented to display a group of resting blood flow bull's eye diagrams, load blood flow bull's eye diagrams and blood flow reserve bull's eye diagrams to evaluate the myocardial blood flow volume in the vascular branch area, and obtain the myocardial blood flow volume levels in the three vascular branch areas, namely the anterior descending branch, the circumflex branch and the right coronary branch;

[0051] Fig. 9 The myocardial blood volume in the vascular branch area was evaluated by using the resting blood flow bull's-eye diagram, stress blood flow bull's-eye diagram and blood flow reserve bull's-eye diagram, and the areas of myocardial blood volume were obtained in the three vascular branch areas, namely the left anterior descending artery, left circumflex artery and right coronary artery. DETAILED DESCRIPTION

[0052] The invention discloses a method for solving the problem of boric acid 99mThe technical method and application of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent for quantitative analysis of myocardial blood flow can be achieved by those skilled in the art by referring to the content of this article and appropriately improving the process parameters. It should be particularly pointed out that all similar substitutions and modifications are obvious to those skilled in the art, and they are all considered to be included in the present invention. The method and application of the present invention have been described through preferred embodiments, and relevant personnel can obviously modify or appropriately change and combine the method and application described herein without departing from the content, spirit and scope of the present invention to realize and apply the technology of the present invention.

[0053] The present invention uses single photon emission computed tomography (SPECT) to perform boric acid 99m The dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents, the technical method for quantitative analysis of myocardial blood flow through image processing, image correction and kinetic model calculation, and the raw materials and reagents used in the use of the technical method can all be purchased on the market.

[0054] The present invention will be further described below in conjunction with embodiments:

[0055] Example 1 Boric Acids 99m Image processing technology for dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents

[0056] The present invention is applied to generate an uncorrected dynamic bull's eye diagram and a left heart blood pool time-activity curve, and the technical method comprises the following steps:

[0057] (1) Quantitative image format generation steps: The reconstruction of dynamic SPECT images requires complete physical correction (including tissue attenuation correction, scatter correction, spatial resolution correction, and image noise correction). At the same time, the pixel values ​​in the reconstructed image need to be converted into a quantitative format representing the radiation concentration. The unit of pixel value is (Bq / ml).

[0058] (2) Image overlay step: Boric acid 99m Dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents Dynamic SPECT images are superimposed over a period of time to produce boric acid-based 99m The static image of Tc-labeled single-photon SPECT imaging agent is first marked with the spatial coordinates and time coordinates of each pixel in the image during superposition. At the same time, the proportion of the time length of each frame to the total time range is used as the weighting coefficient, and the spatial coordinates of each pixel are superimposed at each time point according to the overall superposition time range to produce the boric acid-based 99m Static SPECT images of Tc-labeled single-photon SPECT imaging agents;

[0059] (3) Image segmentation step: Boric acid 99m The static image of Tc-labeled single-photon SPECT imaging agent is cut and distinguished between myocardial and non-myocardial areas by a 3D plane equation (z = ax + by + c) in a rectangular coordinate system representing the image space;

[0060] (4) Image coordinate system conversion steps: The myocardial image takes the mass center of the myocardium as the origin. 99m The static image of Tc-labeled single-photon SPECT imaging agent is transformed from the rectangular coordinate system (x, y, z) to the spherical coordinate system (r, θ, ), taking each pixel as the center during the coordinate system conversion, the pixel value of each central pixel is interpolated by the pixel values ​​of adjacent pixels and the corresponding spatial coordinates to ensure the correctness and continuity of the conversion of each pixel value in the spherical coordinate system;

[0061] (5) Steps for generating the myocardial midline surface: using boric acid transformed by spherical coordinate system 99m Static image of Tc-labeled single-photon SPECT imaging agent, with a radial ray parallel to the radius r, and each angular coordinate of the radial ray passing through The corresponding pixel values ​​are compared to obtain the maximum radial radiation concentration point, and all the maximum radial radiation concentration points are connected to form a 3D myocardial midline surface;

[0062] (6) Steps for generating the geometric transformation matrix of the bull's eye diagram: Arrange the distances of all radial maximum radiation concentration points from the origin in spherical coordinates into a 2D polar square diagram According to the position of the apex in the spherical coordinate system, the rectangular coordinate system and the 2D bull's eye diagram, the apex is taken as the central reference point to perform continuous spatial vector transformation on all radial maximum radiation concentration points to calculate the distance and azimuth (r) between all radial maximum radiation concentration points and the apex in the bull's eye diagram. p ,θ p ), thereby obtaining the geometric correlation of the maximum radiation concentration points of all radial directions of the myocardium in the 2D bull's eye coordinate system, and using it to generate the bull's eye geometric transformation matrix;

[0063] (7) Dynamic bull's eye image generation step: through the image coordinate system conversion step of step (4), the boric acid 99m The dynamic SPECT image of the Tc-labeled single-photon SPECT imaging agent is transformed into a coordinate system for each frame of the image, and then the geometric transformation matrix of the bull's eye image generated in step (6) is used to transform the boric acid 99m Each frame of the dynamic SPECT image of the Tc-labeled single-photon SPECT imaging agent is converted into a bull's-eye image, thereby generating an uncorrected dynamic bull's-eye image; Figure 1 Displaying the image processing results from step (1) to step (7);

[0064] (8) Steps for generating the time-activity curve of the left heart blood pool: In rectangular coordinates, a sampling rectangular box is automatically placed on the boric acid 99m Tc-labeled single-photon SPECT imaging agent between the left atrium and left ventricle in the dynamic SPECT image, and the sampling rectangle box is applied to all boric acid 99m Each frame of the dynamic SPECT image of the Tc-labeled single-photon SPECT imaging agent is used to obtain the time-activity curve of the left heart blood pool. Figure 2 As shown; Steps (1) to (8) can also be applied to resting boric acid 99m Dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents and loaded boric acid 99m Image processing of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents.

[0065] Example 2 Boric Acids 99m Image correction technology for dynamic SPECT imaging using Tc-labeled single-photon SPECT imaging agents

[0066] The method is applied to correct right heart blood pool overflow interference and correct myocardial uptake residual interference, and the method comprises the following steps:

[0067] (1) Correction steps for right heart blood pool overflow: The maximum radiation concentration in the left heart blood pool time-activity curve is used as the time demarcation point for the single-photon SPECT drug to enter the left heart, and a region of interest is used to measure the side wall average value (in terms of Myo 侧壁 The radiation concentration in the left ventricular septal wall pixels of each frame (expressed as Myo 室间隔壁 Myo 侧壁 Compare and obtain the ratio Ratio=(Myo 室间隔壁 / Myo 侧壁 ), when the Ratio value is higher than a certain multiple, it is interference from right heart blood pool overflow and is adjusted to Myo 侧壁 By adjusting the myocardial uptake of the ventricular septum, the interference of the right ventricular blood pool radiation on the left ventricular septum is corrected, thereby obtaining a dynamic bull's eye diagram of the right ventricular blood pool overflow correction. Figure 3 A shows the uncorrected dynamic bull's eye graph before the time cutoff point. Figure 3 B shows the dynamic bull's eye image corrected for right heart blood pool overflow before the time cutoff point; step (1) can be applied to both resting and stress boric acid 99mDynamic SPECT images of dynamic SPECT imaging of a Tc-labeled single-photon SPECT imaging agent, thereby generating a resting dynamic bull's-eye diagram corrected for right ventricular blood pool overflow and a stress dynamic bull's-eye diagram corrected for right ventricular blood pool overflow;

[0068] (2) Pre-injection myocardial uptake residual correction step: The initial pre-injection myocardial residual is obtained by correcting the myocardial uptake of the first frame in the dynamic bull's eye image through right ventricular blood pool overflow, and the pre-injection myocardial residual in each frame is calculated by the myocardial uptake model. The myocardial washout model is the degree of myocardial uptake decline from the peak of myocardial uptake after formal drug injection to the end of dynamic imaging. After normalizing each frame with the peak myocardial uptake, a time-varying myocardial washout function W is generated. 心肌洗脱 (t), and the first frame bull’s eye image × W 心肌洗脱 (t) estimating the residual myocardial volume in the bull's-eye image of each frame, deducting the residual myocardial volume in the bull's-eye image of each frame from the bull's-eye image of each frame, thereby obtaining a right ventricular blood pool overflow correction and pre-injection correction dynamic bull's-eye image, and applying it to generate a pre-injection correction myocardial time-activity curve, Figure 4 Display the myocardial time-activity curve without pre-injection correction, the myocardial residual time-activity curve before injection, and the myocardial time-activity curve corrected before injection;

[0069] (3) Resting myocardial uptake residual correction steps: If it is resting boric acid 99m Dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents in front, loaded with boric acid 99m The dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent is performed in the following order, and the resting dynamic bull's eye image corrected by right ventricular blood pool overflow correction and pre-injection correction in step (1) and step (2) is close to or equal to the bull's eye image at the most late time point as the resting myocardial uptake bull's eye image, and the resting boric acid is used to 99m Starting time and loading boric acid of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent 99m The time difference of the starting time of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent, 99m The physical half-life of Tc nuclide, the physical exponential model and the resting myocardial elution function obtained in step (2) and the extrapolation method are used to estimate the resting myocardial uptake in the presence of boric acid 99m The resting myocardial residual volume at the starting point of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent Residual=exp(-0.693*Δt / T 物理半衰期 )×W 静息心肌洗脱 (Δt , ), thereby generating a resting myocardial residual bull's eye diagram of the load dynamic bull's eye diagram, and obtaining the load myocardial washout function W by the method described in step (2) 负荷心肌洗脱(t) Calculating the resting myocardial residual volume of each frame of the stress dynamic bull's eye diagram, subtracting the resting myocardial residual volume of each frame from the right ventricular blood pool overflow correction and pre-injection correction dynamic bull's eye diagram, thereby generating a stress dynamic bull's eye diagram corrected for right ventricular blood pool overflow and resting myocardial residual volume, and applying it to generate a myocardial time-activity curve, Figure 5 Display the resting myocardial time-activity curve, the stress myocardial time-activity curve, the resting residual myocardial time-activity curve and the resting residual corrected stress myocardial time-activity curve;

[0070] (4) Correction steps for residual myocardial uptake during stress: If the load is boric acid 99m Dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents in front, static boric acid 99m The dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent is performed in the following order, and the stress dynamic bull's eye diagram corrected by right ventricular blood pool overflow correction and pre-injection correction through step (1) and step (2) is close to or equal to the most late time point as the stress myocardial uptake bull's eye diagram, and the stress boric acid is used to 99m The starting time of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents and the static boric acid 99m The time difference of the starting time of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent, 99m The physical half-life of Tc nuclide, the physical exponential model and the stress myocardial elution function obtained in step (2) and the extrapolation method are used to estimate the stress myocardial uptake in resting boric acid 99m The residual volume of myocardial load at the starting time of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent Residual=exp(-0.693*Δt / T 物理半衰期 )×W 负荷心肌洗脱 (Δt , ), thereby generating a resting myocardial residual bull's eye diagram of a resting dynamic bull's eye diagram, and obtaining a resting myocardial washout function W by the method described in step (2) 负荷心肌洗脱 (t) Calculate the residual amount of myocardial load in each frame of the resting dynamic bull's-eye diagram, and deduct the residual amount of myocardial load in each frame from the dynamic bull's-eye diagram corrected for right ventricular blood pool overflow and pre-injection correction, thereby generating a residual bull's-eye diagram of myocardial load in the resting dynamic bull's-eye diagram, and obtain the resting myocardial washout function W by the method described in step (2) 静息心肌洗脱 (t) estimating the residual volume of stressed myocardium in each frame of the resting dynamic bull's-eye diagram, and deducting the residual volume of stressed myocardium in each frame from the dynamic bull's-eye diagram corrected for right ventricular blood pool overflow and pre-injection, thereby generating a resting dynamic bull's-eye diagram corrected for right ventricular blood pool overflow and stressed myocardium residual volume;

[0071] (5) Steps for generating a complete corrected myocardial time activity curve: Based on the resting boric acid99m Dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agents and loaded boric acid 99m The order of dynamic SPECT imaging of Tc-labeled single-photon SPECT imaging agent is to draw a line on the time axis with each pixel in the bull's eye diagram for a group (resting dynamic bull's eye diagram with right ventricular blood pool overflow correction and pre-injection correction, and resting dynamic bull's eye diagram with right ventricular blood pool overflow correction and resting myocardial residual volume correction), or a group (resting dynamic bull's eye diagram with right ventricular blood pool overflow correction and pre-injection correction, and resting dynamic bull's eye diagram with right ventricular blood pool overflow correction and load myocardial residual volume correction). When the fully corrected resting myocardial time-activity curves of each point of the myocardium are averaged, a fully corrected overall resting myocardial time-activity curve can be obtained. When the fully corrected load myocardial time-activity curves of each point of the myocardium are averaged, a fully corrected overall load myocardial time-activity curve can be obtained. Figure 6 A shows the stress dynamic bull's eye diagram with correction for right ventricular blood pool overflow and resting myocardial residual volume. Figure 6 B shows 4 complete corrected stress myocardial activity curves generated by the stress dynamic bull's eye diagram corrected for right ventricular blood pool overflow and resting myocardial residual volume. Figure 6 C shows the resting dynamic bull's eye diagram with right ventricular blood pool overflow correction and pre-injection correction. Figure 6 D shows the four complete corrected resting myocardial activity curves of the resting dynamic bull's eye plot corrected for right ventricular blood pool overflow and pre-injection correction.

[0072] Example 3 Boric Acids 99m A technical method for quantitative calculation of myocardial blood flow in dynamic SPECT imaging using Tc-labeled single-photon SPECT imaging agents

[0073] The technical method includes the following steps:

[0074] (1) Kinetic model application conversion steps: A single-tissue two-chamber model was used, with K1 representing the rate at which the SPECT drug entered the myocardial cells (unit: ml / min / g), k2 representing the elution rate of the SPECT drug from the myocardial cells after entering the cells (unit: 1 / min), and the complete corrected resting whole myocardial time-activity curve, the complete corrected resting whole myocardial time-activity curve, the complete corrected stress whole myocardial time-activity curve, and the complete corrected stress myocardial time-activity curve representing the measured values ​​of myocardial radiation concentration (in Myo 测量 The measured values ​​at each time are regarded as the real myocardial uptake (expressed as Myo 真实 and left heart blood pool interference with myocardial overflow (expressed as Myo 左心血池溢出干扰 ) to produce Myo 测量=(1-FBV)*Myo 真实 +C*Myo 左心血池溢出干扰 The kinetic composition equation of Myo 真实 is a time function with parameters (K1, k2, FBV), Myo 左心血池溢出干扰 is the time-activity curve of the left heart blood pool, FBV is a combination coefficient, and also represents the interference degree coefficient of the left heart blood pool on myocardial overflow, (1-FBV) represents Myo 真实 to Myo 测量 The contribution coefficient of , and also represents the partial volume effect of the myocardium;

[0075] (2) Calculation steps of resting myocardial kinetic parameters: a set of (completely corrected whole resting myocardial time-activity curves and resting left heart blood pool time-activity curves) is first fitted with a set of fixed initial values ​​by using a minimum variance iterative approximation algorithm or a maximum expectation value iterative algorithm to generate resting whole myocardial kinetic parameters including (K1, k2 and FBV). Then, the resting whole myocardial kinetic parameters are used as variable initial values. A set of resting complete corrected myocardial time-activity curves and resting myocardial overflow corrected left heart blood pool time-activity curves are fitted with the kinetic component equation by using the same minimum variance iterative approximation algorithm or the maximum expectation value iterative algorithm to calculate the kinetic parameters of each point in the resting myocardial site. The resting K1 is regarded as a function of resting myocardial blood flow by using the resting Renkin-Crone model. The resting K1 is converted into resting myocardial blood flow and displayed as a bull's eye diagram to generate a resting blood flow bull's eye diagram. Figure 7 A shows the resting blood flow bull's eye diagram, resting K1 bull's eye diagram, resting k2 bull's eye diagram, and resting FBV bull's eye diagram;

[0076] (3) Steps for calculating the dynamic parameters of the loaded myocardium: a set of (completely corrected overall loaded myocardium time-activity curves and loaded left heart blood pool time-activity curves) is first fitted with a set of fixed initial values ​​using a minimum variance iterative approximation algorithm or a maximum expected value iterative algorithm to generate the dynamic parameters of the loaded overall myocardium including (K1, k2 and FBV). The dynamic parameters of the loaded overall myocardium are then used as variable initial values. A set of complete loaded myocardial time-activity curves and loaded myocardial overflow corrected left heart blood pool time-activity curves are fitted with the dynamic component equation using the same minimum variance iterative approximation algorithm or a maximum expected value iterative algorithm to calculate the dynamic parameters of each point in the loaded myocardium. The load K1 is regarded as a function of the loaded myocardial blood flow through the load Renkin-Crone model. The load K1 is converted into the loaded myocardial blood flow and displayed as a bull's eye diagram to generate a load blood flow bull's eye diagram. Figure 7B shows the stress blood flow bull's eye diagram, stress K1 bull's eye diagram, stress K2 bull's eye diagram and resting FBV bull's eye diagram. Figure 7 C shows the blood flow reserve bull's eye diagram;

[0077] (4) Steps for improving fitting operation efficiency: In order to improve fitting operation efficiency, the kinetic model in step (1) is transformed and the kinetic component equations can be simplified by reducing parameters. First, the k2 generated by steps (2) and (3) is used as a fixed value. After the k2 value of step (1) is fixed, the fitting operation of each point of the myocardium is performed to obtain (K1, k2 fixed, FBV) kinetic parameters to improve fitting operation efficiency. The following table shows the comparison of kinetic parameter values ​​and fitting operation efficiency of the dual-tissue dual-chamber model and the simplified model when the fitting operation is performed on a workstation with an Intel quad-core computing unit:

[0078] Table 1

[0079]

[0080] Example 4 Boric Acids 99m The purpose of the quantitative analysis of myocardial blood flow in dynamic SPECT imaging with Tc-labeled single-photon SPECT imaging agents is to establish a myocardial blood volume map, which is used to more intuitively evaluate the resting blood flow, load blood flow and blood flow reserve of the entire myocardium and vascular divisions.

[0081] The method comprises the following steps:

[0082] (1) Steps for establishing a myocardial blood volume map: A 2D coordinate system is used based on three indicators, namely, resting blood flow, load blood flow and blood flow reserve, with the horizontal axis being the load blood flow indicator, the vertical axis being the blood flow reserve indicator, and the diagonal line being the resting blood flow indicator. The 2D coordinate system is further divided into multiple levels according to the indicator value range, thereby establishing a myocardial blood volume map. Figure 8 Displays a 14-level area myocardial blood volume map;

[0083] (2) Myocardial blood volume assessment steps: The resting blood flow bull's-eye diagram, stress blood flow bull's-eye diagram, and blood flow reserve bull's-eye diagram are divided into the whole myocardium and three vascular branch areas (left anterior descending artery, left circumflex artery, and right coronary artery). The average values ​​of the resting blood flow index, stress blood flow index, and blood flow reserve index in the whole myocardium and / or vascular branch areas are calculated respectively. The myocardial blood volume diagram is then used to assess the myocardial blood volume of the whole myocardium or each vascular branch area. Fig. 9 It shows a group of myocardial blood flow volume assessments in vascular branch areas using resting blood flow bull's-eye diagram, stress blood flow bull's-eye diagram and blood flow reserve bull's-eye diagram, and obtains the myocardial blood flow volume areas in the three vascular branch areas, namely the left anterior descending artery, left circumflex artery and right coronary artery.

[0084] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principle of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. An image processing method for dynamic SPECT imaging of boric acid radionuclide labeled single-photon SPECT imaging agents, characterized in that: The steps include: Step 1. Generation of quantitative image format: Dynamic SPECT image reconstruction of boronate radionuclide labeled single-photon SPECT imaging agent, after complete physical correction, and the pixel value of the reconstructed image is converted into the quantitative format of radiation concentration, the pixel value unit is Bq / ml; Step 2. Image superposition: superimposing the dynamic SPECT image of the boric acid radionuclide labeled single-photon SPECT imaging agent over a period of time to generate a static SPECT image of the boric acid radionuclide labeled single-photon SPECT imaging agent; Step 3. Image segmentation: cutting the static SPECT image with linear planes to distinguish the myocardial area from the non-myocardial area; Step 4. Image coordinate system conversion: convert the static SPECT image from a rectangular coordinate system to a spherical coordinate system to generate a myocardial image of the static SPECT image in a spherical coordinate system, and perform an internal difference calculation on the pixel value of each central pixel through the pixel values ​​of adjacent pixels and the corresponding spatial coordinates to ensure the correctness and continuity of the conversion of each pixel value in the spherical coordinate system; Step 5. Generation of myocardial midline surface: using ray tracing, obtaining the radial myocardial maximum radiation concentration point from the myocardial region image in the spherical coordinate system to generate a 3D myocardial midline surface; Step 6. Generation of the bull's eye geometric transformation matrix: Taking the apex of the heart in the spherical coordinate system, the rectangular coordinate system and the 2D bull's eye coordinate system as the central reference point, a continuous space vector transformation is performed on the 3D myocardial midline curved surface in the radial direction of the spherical coordinate system to obtain the geometric correlation of the radial myocardial maximum radiation concentration point in the 2D bull's eye coordinate system, and used to generate the bull's eye transformation geometric matrix; Step 7. Generation of dynamic bull's eye diagram: Applying the bull's eye diagram geometric transformation matrix to boric acid 99m Each frame of dynamic SPECT image of the dynamic SPECT imaging of the Tc-labeled single-photon SPECT imaging agent is used to generate an uncorrected dynamic bull's-eye image of the myocardial area; Step 8. Generation of the left heart blood pool time-activity curve: The 3D region of interest is automatically placed between the left atrium and the left ventricle in the static SPECT image, and the dynamic SPECT image is applied to this position to obtain the left heart blood pool time-activity curve.

2. The image processing method according to claim 1, characterized in that: The dynamic SPECT imaging of the boron acid nuclide labeled single-photon SPECT imaging agent includes the dynamic SPECT imaging of the resting boron acid nuclide labeled single-photon SPECT imaging agent and / or the dynamic SPECT imaging of the loaded boron acid nuclide labeled single-photon SPECT imaging agent.

3. An uncorrected resting dynamic bull's eye diagram, a resting left heart blood pool time-activity curve, a load uncorrected load dynamic bull's eye diagram and / or a load left heart blood pool time-activity curve obtained by the image processing technology method as described in claim 1 or 2.

4. An image correction method for dynamic SPECT imaging of boric acid nuclide labeled single-photon SPECT imaging agents, characterized in that: The steps include: Correction of right heart blood pool overflow: using the maximum radiation concentration in the left heart blood pool time-activity curve obtained by the image processing technology method as described in claim 1 or 2 as the time demarcation point for the drug to enter the left heart, measuring the radiation concentration of the left ventricular lateral wall in the uncorrected dynamic bull's eye map through the region of interest, when the pixel value in the left ventricular septal wall is higher than 1.2 to 3 times the average value of the radiation concentration of the left ventricular septal wall, it is considered to be interference from right heart blood pool overflow, and adjusted to the radiation concentration of the left ventricular lateral wall to obtain a dynamic bull's eye map of right ventricular blood pool overflow correction; Correction of pre-injection myocardial uptake residue: for pre-injection myocardial uptake residue, the pre-injection myocardial residual amount of each frame of the dynamic SPECT image is estimated through the myocardial uptake and myocardial washout model of the first frame of the dynamic SPECT image, and deducted to generate a right ventricular blood pool overflow correction and pre-injection correction dynamic bull's eye map; Correction of resting myocardial uptake residue: For dynamic SPECT imaging of boric acid radionuclide labeled single-photon SPECT imaging agents at rest first and then under load, the myocardial residual degree is estimated from the static bull's eye diagram superimposed with the dynamic bull's eye diagram corrected for right ventricular overflow by using the time difference between the starting points of the two imaging times, the physical half-life of the radionuclide, the physical index model, and the resting myocardial washout model obtained from the time-varying characteristics of the resting dynamic bull's eye diagram, so as to generate a resting myocardial residual bull's eye diagram, and the load right ventricular overflow corrected dynamic bull's eye diagram is deducted from the load myocardial washout model obtained by considering the time-varying characteristics of the load dynamic bull's eye diagram, so as to generate a load dynamic bull's eye diagram corrected for right ventricular overflow and resting myocardial residual; Correction of residual myocardial uptake under load: For dynamic SPECT imaging of boric acid radionuclide labeled single-photon SPECT imaging agents with load followed by rest, the residual degree of myocardial uptake under load is estimated from the superposition of the dynamic bull's eye diagram corrected for right ventricular overflow through the time difference between the starting points of the two imaging times, the physical half-life of the radionuclide, the physical index model, and the load myocardial washout model obtained from the time-varying characteristics of the dynamic bull's eye diagram under load, so as to generate a residual bull's eye diagram of myocardial uptake under load, and the dynamic bull's eye diagram corrected for right ventricular overflow is deducted from the resting myocardial washout model obtained by considering the time-varying characteristics of the resting dynamic bull's eye diagram, so as to generate a resting dynamic bull's eye diagram corrected for right ventricular overflow and myocardial residual. Generation of fully corrected myocardial time-activity curves: A pair of resting dynamic bull's-eye diagrams corrected for right heart overflow and pre-injection, and a stress dynamic bull's-eye diagram corrected for right heart overflow and resting myocardial residual; or a pair of resting dynamic bull's-eye diagrams corrected for right heart overflow and myocardial residual, and a stress dynamic bull's-eye diagram corrected for right heart overflow and pre-injection; used to generate a set of fully corrected resting myocardial time-activity curves and fully corrected stress myocardial time-activity curves.

5. Application of the image correction method as claimed in claim 4 in correcting right heart blood pool overflow interference and / or correcting myocardial uptake residual interference.

6. A method for quantitative calculation of myocardial blood flow in dynamic SPECT imaging using boric acid radionuclide labeled single-photon SPECT imaging agents, characterized in that: The following steps are involved: Transformation of kinetic model application: A single-tissue two-chamber model and a combination coefficient FBV are used, and the myocardial radiation concentration measured at each time point is regarded as a combination of the true myocardial uptake and the interference of the left heart blood pool on the myocardial overflow, where (1-FBV) is also the partial volume coefficient caused by the limited resolution of SPECT, and FBV is also the overflow of the left heart blood pool on the myocardium, to generate a kinetic component equation with kinetic parameters (K1, k2, FBV); RC=(1-FBV) Calculation of the kinetic parameters of the resting myocardial part: firstly, a set of fixed initial values ​​are used to fit a set of corrected resting myocardial time-activity curves and resting left heart blood pool time-activity curves with the kinetic composition equation to generate the kinetic parameters of the resting whole myocardium as the variable initial values, and then the corrected resting myocardial time-activity curves and resting left heart blood pool time-activity curves of each point in the myocardial part are fitted with the kinetic composition equation to calculate the kinetic parameters of each point in the resting myocardial part, and through the resting Renkin-Crone model, the resting K1 is regarded as a function of myocardial blood flow and converted, and then the resting myocardial blood flow value is corrected by the product of heart rate and blood pressure, and displayed as a bull's eye diagram, thereby generating a resting blood flow bull's eye diagram; Calculation of dynamic parameters of loaded myocardial parts: firstly, a set of fixed initial values ​​are used to fit a set of corrected loaded myocardial time-activity curves and loaded left heart blood pool time-activity curves with the dynamic composition equations to generate dynamic parameters of the loaded whole myocardium as the variable initial values, and then the corrected loaded myocardial time-activity curves and loaded left heart blood pool time-activity curves of each point in the myocardial part are fitted with the dynamic composition equations to calculate the dynamic parameters of each point in the loaded myocardial part, and through the loaded Renkin-Crone model, the load K1 is regarded as a function of myocardial blood flow, and converted and displayed as a bull's eye diagram to generate a loaded blood flow bull's eye diagram; The flow reserve bull's-eye is obtained by dividing the stress flow bull's-eye by the resting flow bull's-eye.

7. The method for quantitatively calculating myocardial blood flow according to claim 6, characterized in that: In the transformation of the dynamic model application, for the rest and load fitting operations, the k2 of the entire myocardium is used as the k2 value of each point of the myocardium for fitting operations before the fitting operations to improve the fitting operation efficiency.

8. A method for constructing a myocardial blood volume map, characterized in that: The resting myocardial blood flow value, the load myocardial blood flow value, the blood flow reserve value and the 2D coordinate system obtained by the quantitative calculation method of myocardial blood flow as described in claim 6 or 7 are used to construct a myocardial blood flow volume map.

9. A method for evaluating myocardial blood volume, characterized in that: The resting blood flow bull's eye diagram, the load blood flow bull's eye diagram and the blood flow reserve bull's eye diagram obtained by the quantitative calculation method of myocardial blood flow as described in claim 6 or 7 are divided into the whole myocardium and three vascular branch areas, the three vascular branch areas include the left anterior descending artery, the left circumflex artery and the right coronary artery, and the average value of each index in the whole myocardium and each vascular branch area is calculated respectively; The evaluation method includes evaluating the high and low areas of resting blood flow, load blood flow and / or blood flow reserve of the whole myocardium and vascular divisions.

10. The image processing method according to claim 1 or 2, the image correction method according to claim 4, the quantitative calculation method of myocardial blood flow according to claim 6 or 7, the construction method according to claim 8 and / or the evaluation method according to claim 9, characterized in that: The nuclides include but are not limited to 99m Tc.

11. A processor, electronic device or apparatus, characterized in that: The processor is used to run a program, and when the program is run, the image processing method according to claim 1 or 2, the image correction method according to claim 4, the myocardial blood flow quantitative calculation method according to claim 6 or 7, the construction method according to claim 8 and / or the evaluation method according to claim 9 are executed; The electronic device comprises: a memory and the processor; Wherein, one or more computer programs are stored in the memory, and the one or more computer programs include instructions; when the instructions are executed by the processor, the electronic device executes the image processing method according to claim 1 or 2, the image correction method according to claim 4, the myocardial blood flow quantitative calculation method according to claim 6 or 7, the construction method according to claim 8 and / or the evaluation method according to claim 9; The device includes one or more components, which implement the image processing technology method as described in claim 1 or 2, the image correction method as described in claim 4, the myocardial blood flow quantitative calculation method as described in claim 6 or 7, the construction method as described in claim 8 and / or the evaluation method as described in claim 9.

Images