CN105310712A - Method and system for conducting image measurement to obtain tumor standardized uptake value through nuclear medicine single photons - Google Patents
Method and system for conducting image measurement to obtain tumor standardized uptake value through nuclear medicine single photons Download PDFInfo
- Publication number
- CN105310712A CN105310712A CN201410257648.7A CN201410257648A CN105310712A CN 105310712 A CN105310712 A CN 105310712A CN 201410257648 A CN201410257648 A CN 201410257648A CN 105310712 A CN105310712 A CN 105310712A
- Authority
- CN
- China
- Prior art keywords
- image
- single photon
- computer tomography
- patient
- uptake value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 206010028980 Neoplasm Diseases 0.000 title claims abstract description 57
- 238000005259 measurement Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000009206 nuclear medicine Methods 0.000 title abstract 2
- 238000002591 computed tomography Methods 0.000 claims abstract description 86
- 238000012937 correction Methods 0.000 claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000005855 radiation Effects 0.000 claims description 14
- 230000000155 isotopic effect Effects 0.000 claims description 13
- 238000013016 damping Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 9
- 239000000523 sample Substances 0.000 claims description 9
- 208000031229 Cardiomyopathies Diseases 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 5
- 210000000577 adipose tissue Anatomy 0.000 claims description 4
- 230000037396 body weight Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 238000011084 recovery Methods 0.000 abstract description 3
- 238000002603 single-photon emission computed tomography Methods 0.000 abstract 2
- 238000010586 diagram Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 9
- 238000003384 imaging method Methods 0.000 description 4
- 238000011002 quantification Methods 0.000 description 4
- 238000002601 radiography Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 239000003814 drug Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- PNDPGZBMCMUPRI-HVTJNCQCSA-N 10043-66-0 Chemical compound [131I][131I] PNDPGZBMCMUPRI-HVTJNCQCSA-N 0.000 description 1
- 208000018084 Bone neoplasm Diseases 0.000 description 1
- GYHNNYVSQQEPJS-OIOBTWANSA-N Gallium-67 Chemical compound [67Ga] GYHNNYVSQQEPJS-OIOBTWANSA-N 0.000 description 1
- 208000001132 Osteoporosis Diseases 0.000 description 1
- GKLVYJBZJHMRIY-OUBTZVSYSA-N Technetium-99 Chemical compound [99Tc] GKLVYJBZJHMRIY-OUBTZVSYSA-N 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940006110 gallium-67 Drugs 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 229940056501 technetium 99m Drugs 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Landscapes
- Nuclear Medicine (AREA)
Abstract
The invention relates to a method and system for conducting image measurement to obtain the tumor standardized uptake value through nuclear medicine single photons. Image measurement is conducted to obtain the tumor standardized uptake value through single-photon emission computed tomography (SPECT) or obtain the tumor standardized uptake value through single-photon emission computer tomography/computer tomography (SPECT/CT). The method includes the photography step, the scattering correction step, the attenuation correction step, the image resolution ratio recovery step, the noise removal step and the uptake value calculation step. The method solves the problems of physical interference on dynamic single-photon emission computer tomography and patient movement, and therefore the capacity for quantifying the tumor standardized uptake value of patients can be achieved.
Description
Technical field
The present invention relates to the design of a kind of image measuring method and system, particularly about a kind of method and system of core doctor single photon radiographic measurement tumor standard uptake value.
Background technology
Clinically at present, in order to measure the standard uptake value (standardizeduptakevalue of tumor for medicine, SUV), first need first to set up as by positron radiography (positronemissiontomography, or the image of equal attribute that obtains of positron radiography/computer tomography (PET/ComputedTomography, PET/CT) PET).Although but PET or PET/CT can accurate measurement tumor standard uptake value, the problems such as high born by the popularization degree and the inspection funds that are limited to radiography medicine, PET and PET/CT cannot widely use in each level hospital.
The image mode of SPECT and SPECT/CT has the advantage of economic par compared to PET and PET/CT, can be widely used in each level hospital, but traditional SPECT and SPECT/CT radiography is due to Physical Interference and the problem such as spatial resolution is low, diagnosis often suffers from image artifacts (artifacts) and the inadequate puzzlement of accuracy, and cannot be used to quantize tumor standard uptake value, therefore limit to the one Eco-power clinical tool of SPECT and SPECT/CT as diagnosing tumour disease.In view of the above, real needs one Eco-power quantification imaging system, for the quantification realizing tumor standard uptake value.
Summary of the invention
Edge this, namely object of the present invention is to provide the method and system of a kind of core doctor single photon radiographic measurement tumor standard uptake value, to improve the problem of prior art, reaches the effect of economy and accurate measurement.
The technological means that the problem that the present invention is solution prior art adopts is the method for a kind of core doctor single photon radiographic measurement tumor standard uptake value, utilize Single Photon Emission computer tomography or utilize Single Photon Emission computer tomography/computer tomography to carry out radiographic measurement tumor standard uptake value, method comprises: a photography step, by Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography to take an image of a patient; One scatter correction step, deducts the scattering component of scattering component and correcting image by the original projection deriving from a photopeak energy window; One correction for attenuation step, calculates the decay number of each irradiation image pixel of the image through scatter correction by coordinate conversion computer tomography image and lschemic cardiomyopathy; One image resolution recovering step, by setting up a point spread function matrix from a point source traverse measurement and point spread function matrix number being given iteration integration reconstruction and again recover resolution; One noise remove step, analyzes noise filter by the image through image reconstruction via one and fore-telling Herba Orostachyos (Poisson) emulator gives iteration integration reconstruction, to remove the noise in image; And a uptake values calculation procedure, according to a concentration image rays relation, calculate the isotopes concentration corresponding to a region of interest of the image through removing noise, and according to the body weight of patient and the isotopic dosage putting on patient, calculate the tumor standard uptake value of patient.
In one embodiment of this invention, in scatter correction step, more decay according to the rotation time of a lonizing radiation probe and the isotope of isotopic half-life compensation image.
In one embodiment of this invention, in correction for attenuation step, use derives from computer tomography or radionuclide penetrates the decay number of each irradiation image pixel of the decay number calculating image of image and creates a damping matrix, and integrates damping matrix with correcting attenuation in an iterative approximation.
In one embodiment of this invention, in image resolution recovering step, by setting up the point spread function matrix number of depth dependence to measure the characteristic of a collimator (collimator).
In one embodiment of this invention, more comprise patient's shift calibrating step, the patient that the patient comprised in correct scan moves between (Intra-scanpatientmotion) and correct scan moves (Inter-scanpatientmotion), patient wherein in correct scan moves the angle comprising and pop one's head in for lonizing radiation measures patient's movement displacement in iterative projection mode, and the patient between correct scan moves to comprise and manually again aims at computer tomography image and lschemic cardiomyopathy.
In one embodiment of this invention, in uptake values calculation procedure, more according to the whole body quality of patient, body fat rate and/or total body surface area and regular tumor standard uptake value.
In one embodiment of this invention, in uptake values calculation procedure, concentration image rays relation according to flow into a prosthese isotopic concentration, calculated by the image of prosthese captured by Single Photon Emission computer tomography and the intensity of radiation that is injected with isotopic prosthese and obtain.
In one embodiment of this invention, in uptake values calculation procedure, wherein prosthese is a symmetry prosthese.
The present invention is a kind of system utilizing the core doctor single photon radiographic measurement tumor standard uptake value of Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography for another technological means that the problem solving prior art adopts, and comprises a Single Photon Emission computer tomography/computed tomography apparatus, a computer and a program product.Single Photon Emission computer tomography/computed tomography apparatus is for taking the image of a Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography.Computer has a processing unit.Program product comprises a computer readable program code, and the processing unit via computer loads and performs the following step: the image receiving Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography; The scattering of correcting image; The decay of correcting image; Recover the spatial resolution of image; Remove the image noise of image; And the tumor standard uptake value corresponding to calculating image.
Via the technology used in the present invention means, by the problem such as the Physical Interference and spatial resolution that solve SPECT and SPECT/CT be low, the image artifacts avoiding tumor disease to diagnose and the low erroneous judgement problem that may cause of accuracy, and as quantizing the Eco-power clinical tool of one of tumor standard uptake value, reach the effect of economy and accurate measurement.
Accompanying drawing explanation
Specific embodiment of the present invention, by by following embodiment and to be attachedly further described in graphic.
The display of Fig. 1 system is according to the flow chart of the method for the core doctor single photon radiographic measurement tumor standard uptake value of one embodiment of the invention.
The display of Fig. 2 system is according to the equipment drawing of the system of the core doctor single photon radiographic measurement tumor standard uptake value of one embodiment of the invention.
The display of Fig. 3 system is according to the block chart of the system of the core doctor single photon radiographic measurement tumor standard uptake value of one embodiment of the invention.
The system of the core doctor single photon radiographic measurement tumor standard uptake value of Fig. 4 A system display foundation one embodiment of the invention takes the schematic diagram of patient image to left back oblique side from right front oblique side.
The system of the core doctor single photon radiographic measurement tumor standard uptake value of Fig. 4 B system display foundation one embodiment of the invention takes the schematic diagram of patient image to rear from front.
The display of Fig. 4 C system is according to the schematic diagram of system extremely right back oblique side's shooting patient image from LAO side of the core doctor single photon radiographic measurement tumor standard uptake value of one embodiment of the invention.
The system of the core doctor single photon radiographic measurement tumor standard uptake value of Fig. 4 D system display foundation one embodiment of the invention takes the schematic diagram of patient image to front from rear.
The display of Fig. 5 A system is according to the flow chart of patient's shift calibrating step of one embodiment of the invention.
The display of Fig. 5 B system is according to the schematic diagram of patient image before resolution is recovered of one embodiment of the invention.
The display of Fig. 5 C system is according to the schematic diagram of patient image after resolution is recovered of one embodiment of the invention.
Fig. 6 A to 6C is that display is according to the schematic diagram of patient image of the present invention before SPECT and CT para-position.
Fig. 7 A to 7C is that display is according to the schematic diagram of patient image of the present invention after SPECT and CT para-position.
The display of Fig. 8 system is according to the schematic diagram of the damping matrix of the patient image of one embodiment of the invention.
The patient image of Fig. 9 A system display foundation one embodiment of the invention is in the schematic diagram without timing.
Fig. 9 B to 9F is that the patient image of display according to one embodiment of the invention is in the schematic diagram of every timing.
Figure 10 A to 10C is CT, SPECT, SPECT/CT image schematic diagram in oncotherapy before of display according to patient of the present invention.
Figure 11 A to 11C is CT, SPECT, SPECT/CT image schematic diagram in oncotherapy after of display according to patient of the present invention.
Wherein:
The system of 100 core doctor single photon radiographic measurement tumor standard uptake values
1 Single Photon Emission computer tomography/computed tomography apparatus
10 lonizing radiation probes
12 lonizing radiation probes
14 computer tomography cameras
16 imaging tables
18 frames
2 computers
3 program products
30 computing modules
31 image receiver modules
32 scatter correction modules
33 fall-off correction module
34 resolution recover module
35 noise remove modules
36 patient's shift calibrating modules
A, b, c, d, e, f, g reference point
M1 damping matrix
P1 processing unit
U patient
Detailed description of the invention
The invention provides quantification Single Photon Emission computer tomography (SPECT) reconstructing system, carry out the quantification of tumor standard uptake value with the dynamic imaging of SPECT or Single Photon Emission computer tomography/computer tomography (SPECT/CT).Refer to shown in Fig. 2 and Fig. 3, according to the system 100 of the core doctor single photon radiographic measurement tumor standard uptake value of one embodiment of the invention, system utilizes the radiographic measurement tumor standard uptake value of Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography, and system 100 comprises Single Photon Emission computer tomography/computed tomography apparatus 1, computer 2 and a program product 3 (loading on computer 2).
Program product 3 can comprise computer readable program code, and machine 2 loads and performs some treatment steps as calculated.Computer 2 can including but not limited to the computer of any type, and it comprises desktop computer, notebook computer and intelligent device, such as, and tablet PC and intelligent mobile phone.Program product 2 can comprise software, and it can be loaded on computer 2 or accessed by computer 2.Software through loading can comprise the program loading on computer 2 hard disk, or an application program of intelligent device.In certain embodiments, a web browser is used to access this software by computer 2.Computer 2 can utilize the Internet, external network, internal network, host server (hostserver) and network high in the clouds (internetcloud) etc. and access this software.The present invention can be exchanged into the computer software of medical imaging in medical environment, and patient is carried out to the research of SPECT or SPECT/CT.
Single Photon Emission computer tomography/computed tomography apparatus 1 is for taking the image of a Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography for an example patient U (coordinating Fig. 4 A to Fig. 4 B).In order to take this image, example patient U can be placed on an imaging table (imagingtable) 16 of Single Photon Emission computer tomography/computed tomography apparatus 1.Single Photon Emission computer tomography/computed tomography apparatus 1 can have a lonizing radiation probe 10,12.One computer tomography camera 14 can be placed in frame 18.Frame 18 can comprise a column scanner, to take this image.
Computer 2 has a processing unit P1.Program product 3 comprises a computer readable program code, and the processing unit P1 via computer 2 loads and performs the following step: the image receiving Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography; The scattering of correcting image; The decay of correcting image; Recover the spatial resolution of image; Remove the image noise of image; And the tumor standard uptake value corresponding to calculating image.As shown in Figure 3, in the present embodiment, program product 3 can comprise image receiver module 31, scatter correction module 32, fall-off correction module 33, resolution and recover module 34, noise remove module 35, patient's shift calibrating module 36 and a computing module 30, sequentially performs above-mentioned steps respectively.
Refer to shown in Fig. 1, according to the method for the core doctor single photon radiographic measurement tumor standard uptake value of one embodiment of the invention, be utilize Single Photon Emission computer tomography or utilize Single Photon Emission computer tomography/computer tomography to carry out radiographic measurement tumor standard uptake value, the method of core doctor single photon radiographic measurement tumor standard uptake value comprises a photography step (step S1), one scatter correction step (step S2), one correction for attenuation step (step S3), one image resolution recovering step (step S4), one noise remove step (step S5), and a uptake values calculation procedure (step S6).Below and coordinate the system 100 of the core doctor single photon radiographic measurement tumor standard uptake value shown in Fig. 2, the method for core doctor single photon radiographic measurement tumor standard uptake value of the present invention is described in detail as follows:
First, photography step (step S1), by Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography to take the image of patient U.The acquisition of image is to be detoured the image data form obtained in the patient U injected after isotope with 180 degree or 360 degree by lonizing radiation probe 10,12.Take turns in rotation at each, image data system is represented by the two indices of dynamic framework (frame) and angle.Isotope Ke Wei Technetium-99m, iodine-131, iodo-123, the medicine such as gallium-67, Value linear, but be not limited thereto.As shown in Fig. 4 A to Fig. 4 D, lonizing radiation probe 10,12 can from a right front oblique side to left front rear (Fig. 4 A), front to rear (Fig. 4 B), LAO side and right back oblique side (Fig. 4 C), rear to front (Fig. 4 D) direction of rotation take the image of this patient U, or in contrast to the direction shooting patient U image of this direction of rotation.The present invention can comprise the medical digital image transmission agreement (digitalimagingandcommunicationsinmedicine with a standard, DICOM) form or with machine proprietary formats shooting and receive image data, for follow-up image procossing.Those images can be from the image captured by SPECT or SPECT/CT, and derive from the DICOM form CT image of SPECT/CT or CT, or derive from the data of the machine proprietary formats of SPECT or the lschemic cardiomyopathy of DICOM form.
Then, scatter correction step (step S2), deducts the scattering component (step S21) of scattering component and correcting image by the original projection deriving from a photopeak energy window.This scatter correction by readjust derive from such as Compton scattering (Comptonscatter) can the counting of original projection of window, and deduct derive from a photopeak energy window original projection (measuring projection) in scattering component and calculated.In addition, in preferred embodiment, in scatter correction step (step S2), more decay (step S22) according to the lonizing radiation probe rotation time of 10,12 and the isotope of isotopic half-life compensation image.In certain embodiments, isotope decay can be performed to image to correct.By readjusting the counting in the original projection of the framework in capturing corresponding to SPECT and angle, and the isotope decay that can perform the original projection of SPECT corrects.When SPECT frame 18 is surrounded on patient U, this readjusts factor by possessing the exponential disintegration module (exponentialdecaymodel) that is defined in the angle of corresponding time point and the time factor of framework and being calculated.
Thereafter, correction for attenuation step (step S3), calculates the decay number of each irradiation image pixel of the image through scatter correction (step S2) by coordinate conversion computer tomography image and lschemic cardiomyopathy.Specifically, in the present embodiment, system's use derives from computer tomography or radionuclide penetrates the decay number of each irradiation image pixel of the decay number calculating image of image and creates a damping matrix M1, and in an iterative approximation, integrates damping matrix M1 with correcting attenuation.Decay derives from computer tomography or 140 kilo electron volts (keV) by using or is corrected corresponding to the radionuclide of the isotopic particular energy of other SPECT penetrates the attenuation quotient (attenuationcoefficients) of image.Computerized tomography image or radionuclide penetrate image can be converted into damping matrix M1, and damping matrix M1 can record and correspond to the amplitude of the photon attenuation of each pixel in an image of five mark (x, y, z, angle, framework) institute labellings.Image para-position is carried out by SPECT/CT and CT before correction for attenuation, improve the accuracy of correction for attenuation, after the para-position of SPECT/CT and CT image, (Fig. 7 A to 7B) has the position data of correct image compared to (Fig. 6 A to 6B) before the para-position of SPECT/CT and CT image.
Then, image resolution recovering step (step S4), integrate and rebuild by setting up a point spread function matrix from a point source traverse measurement and point spread function matrix number (pointspreadfunction, PSF) being given iteration and again recover the resolution of image.Specifically, in the present embodiment, be integrate and rebuild by setting up a point spread function matrix from a point source traverse measurement and point spread function matrix number being given iteration and again recover the resolution of image.By this, can further by setting up the point spread function matrix number of depth dependence to measure the characteristic of a collimator (not shown) of Single Photon Emission computer tomography/computed tomography apparatus 1.One group of PSF can be measured for a SPECT scanner type with a particular type collimator.PSF modeling by analytical function (analyticfunctions), and then produce one group of complete PSF for having to each image pixel (indicated by six labellings (x, y, z, distance, angle, framework)) of this collimator surfaces known distance.Then this PSF can be applied to this irradiation image before the front projection step in iterative approximation.In addition, above-mentioned damping matrix M1 also can be applied in iterative approximation front projection step before irradiation image.
In order to avoid the interference of patient U movement, preferably, more comprise patient's shift calibrating step (step S15), the patient U that the patient U comprised in correct scan moves between (Intra-scanpatientmotion) and correct scan moves (Inter-scanpatientmotion), system is performed by patient's shift calibrating module 36, coordinating consults shown in Fig. 5 A, patient U wherein in correct scan moves and comprises for lonizing radiation probe 10, the angle of 12 measures the displacement of patient U movement in iterative projection mode, patient U between correct scan moves to comprise and manually again aims at computer tomography image and lschemic cardiomyopathy.Wherein, patient U in correct scan moves (Intra-scanpatientmotion) must in scatter correction step (step S2) front execution, and the patient U between correct scan moves (Inter-scanpatientmotion) as long as perform correction for attenuation step (step S3) is front.
It is that a kind of patient occurring in dynamic SPECT scan period moves kenel that patient in scanning moves (Intra-scanpatientmotion), and the movement of the patient U in this scanning may cause artifact.The mobile system of the patient U in the scanning in the projection through measuring of each dynamic framework is individually evaluated to correct.In certain embodiments, this correction is by being performed for projection through measuring of angle and framework displacement iteratively.For the displacement of the vertical of angle and framework and horizontal direction, by maximize to compare at vertical and horizontal direction indivedual cross correlation between the two have before projection projection this projection through measuring and define.Until each projection all reaches stable maximum cross correlation, identical process can perform iteratively for multiple angle and framework.
It is the kenel penetrating a kind of patient's movement caused between scanning at dynamic SPECT scanner uni CT or radionuclide that patient between scanning moves (Inter-scanpatientmotion), and the patient between scanning moves and can cause correction for attenuation artifact.Between SPECT and CT, or SPECT and the patient penetrating between image between scanning move, by manually again aiming at x, y, z direction and σ, δ, θ angle individually corrects.After the para-position of SPECT/CT and CT image, (Fig. 7 A to 7B) has the position data of correct image compared to (Fig. 6 A to 6B) before the para-position of SPECT/CT and CT image.
After image resolution recovering step (step S4), each reference point a of Fig. 4 A to Fig. 4 B, the image resolution of b, c, d, e, f, g in recovery after (Fig. 5 C) compared to recovery before (Fig. 5 B) tool improve significantly.
Subsequently, noise remove step (step S5), image through image reconstruction (step S4) is foretold Herba Orostachyos (Poisson) emulator via an analysis noise filter and one and gives iteration integration reconstruction, to remove the noise in image.The noise of image is by applying the data for projection of a controllable analysis filter to front projection, and produce the homogenize noise profile of the data for projection be measured and removed, this homogenize noise profile through the data for projection of measurement also compares in the step of two data sets by equivalent analysis wave filter (equivalentanalyticfilter) and is filtered in iterative approximation.As a substituting embodiment, noise in the projection of front projection by the Random process of Poisson distribution (Poissondistribution), to produce the homogenize noise profile of the data for projection through measuring in the step comparing two data sets at iterative approximation.
Finally, uptake values calculation procedure (step S6), according to a concentration image rays relation, calculate the isotopes concentration corresponding to a region of interest of the image through removing noise (step S5), and according to the body weight of patient U and the isotopic dosage putting on patient U, calculate the tumor standard uptake value of patient U.Wherein, in order to the radioactivity penetrance for human body corrects, in the present embodiment, concentration image rays relation according to flow into a prosthese isotopic concentration, calculated by the image of prosthese captured by Single Photon Emission computer tomography and the intensity of radiation that is injected with isotopic prosthese and obtain.Conveniently measure, preferably, prosthese is a symmetry prosthese.But, due to the specificity of biology, in the present embodiment, more according to the whole body quality of patient U, body fat rate and total body surface area and regular tumor standard uptake value.Again in other embodiment, also can only according to the whole body quality of patient U, body fat rate, total body surface area one of them or wherein two and regular tumor standard uptake value.Certainly, the present invention is not limited thereto, also can according to other physiological parameters of patient U regular tumor standard uptake value, such as, utilize patient U osteoporosis degree normalization bone tumor tumor standard uptake value.
Fig. 9 A to 9F is the not calibrated image (Fig. 9 A) representing patient U respectively, through the image (Fig. 9 B) of noise remove, through noise remove, the image (Fig. 9 C) of patient's shift calibrating, through noise remove, patient's shift calibrating, the image (Fig. 9 D) of correction for attenuation, through noise remove, patient's shift calibrating, correction for attenuation, the image (Fig. 9 E) of scatter correction, through noise remove, patient's shift calibrating, correction for attenuation, scatter correction, the image (Fig. 9 F) of image reconstruction.Can obviously observe and define knub position and image rays intensity by the contrast of Fig. 9 A to 9F.
Figure 10 A to 10C is CT, SPECT, SPECT/CT image schematic diagram in oncotherapy before of display according to patient of the present invention.Figure 11 A to 11C is CT, SPECT, SPECT/CT image schematic diagram in oncotherapy after of display according to patient of the present invention.Provable compared to prior art, via technological means of the present invention, the difference of tumor before and after treatment can more obviously be observed out by naked eyes.
Above computer basic data processing System and method for only for illustrative object, and is embodied as the form of any computer system, programming or processing environment, or separately or a computer program of combined with hardware.The present invention is also embodied as the software be stored on a computer fetch medium, and performs as general service or a computer program of special purpose computer.For simplicity's sake, only describe the aspect of those and the closely-related system of the present invention herein, and omit in product details well known in the art.Based on same cause, also this computer hardware is not further described.Therefore, should be understood that, the present invention is not limited to any specific computer language, program or computer.What can further contemplate that is, the present invention can operate on an independently computer system or can in a computer server system running, and this computer server system can be interconnected and access or can pass through the Internet and accessed by client through an internal network by multiple client computer system.In addition, many embodiments of the present invention are applicable in the industry of broad range.The method that the present invention discloses a system, this system is implemented and to be stored on computer fetch medium and to perform the software of the method in the computer of general service or special-purpose as a computer program, all in scope of the present invention.Further, disclose a method in the present invention, implement the apparatus system of the method once setting, they all comprise within the scope of the present invention.
Above describing is only preferred embodiment of the present invention and illustrates, is allly skillful in this those skilled in the art when can do other all improvement according to above-mentioned explanation, but these change and still belong in the spiritual and the scope of the claims that defines of invention of the present invention.
Claims (9)
1. the method for a core doctor single photon radiographic measurement tumor standard uptake value, the party's genealogy of law utilizes Single Photon Emission computer tomography or utilizes Single Photon Emission computer tomography/computer tomography to carry out radiographic measurement tumor standard uptake value, and the method comprises:
One photography step, by Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography to take an image of a patient;
One scatter correction step, deducts scattering component by the original projection deriving from a photopeak energy window and corrects the scattering component of this image;
One correction for attenuation step, calculates the attenuation quotient of each irradiation image pixel of this image through scatter correction by coordinate conversion computer tomography image and lschemic cardiomyopathy;
One image resolution recovering step, by setting up a point spread function matrix from a point source traverse measurement and this point spread function matrix number being given iteration integration reconstruction and again recover resolution;
One noise remove step, analyzes noise filter by the image through image reconstruction via one and fore-telling Herba Orostachyos (Poisson) emulator gives iteration integration reconstruction, to remove the noise in this image; And
One uptake values calculation procedure, according to a concentration image rays relation, calculate this isotopes concentration corresponding to a region of interest through removing the image of noise, and according to the body weight of this patient and the isotopic dosage putting on this patient, calculate the tumor standard uptake value of this patient.
2. the method for core doctor single photon radiographic measurement tumor standard uptake value as claimed in claim 1, is characterized in that, in this scatter correction step, more decay according to the rotation time of a lonizing radiation probe and the isotope of this isotopic this image of half-life compensation.
3. the method for core doctor single photon radiographic measurement tumor standard uptake value as claimed in claim 1, it is characterized in that, in this correction for attenuation step, system uses and derives from the attenuation quotient that computer tomography or radionuclide penetrate image and calculate the attenuation quotient of each irradiation image pixel of this image and create a damping matrix, and in an iterative approximation, integrates this damping matrix with correcting attenuation.
4. the method for core doctor single photon radiographic measurement tumor standard uptake value as claimed in claim 1, is characterized in that, in this image resolution recovering step, is that this point spread function matrix number by setting up depth dependence is to measure the characteristic of a collimator.
5. the method for core doctor single photon radiographic measurement tumor standard uptake value as claimed in claim 1, it is characterized in that, more comprise patient's shift calibrating step, the patient that the system patient comprised in correct scan moves between (Intra-scanpatientmotion) and correct scan moves (Inter-scanpatientmotion), patient wherein in this correct scan moves the angle comprising and pop one's head in for lonizing radiation measures this patient's movement displacement in iterative projection mode, patient between this correct scan moves to comprise and manually again aims at this computer tomography image and this lschemic cardiomyopathy.
6. the method for core doctor single photon radiographic measurement tumor standard uptake value as claimed in claim 1, it is characterized in that, in this uptake values calculation procedure, more this tumor standard uptake value of normalization according to the whole body quality of this patient, body fat rate and/or total body surface area.
7. the method for core doctor single photon radiographic measurement tumor standard uptake value as claimed in claim 1, it is characterized in that, in this uptake values calculation procedure, this concentration image rays relation according to flow into a prosthese this isotopic concentration, calculated by the image of this prosthese captured by this Single Photon Emission computer tomography and the intensity of radiation that is injected with this this prosthese isotopic and obtain.
8. the method for core doctor single photon radiographic measurement tumor standard uptake value as claimed in claim 7, it is characterized in that, in this uptake values calculation procedure, wherein this prosthese is a symmetry prosthese.
9. utilize a system for the core doctor single photon radiographic measurement tumor standard uptake value of Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography, comprise:
One Single Photon Emission computer tomography/computed tomography apparatus, for taking the image of a Single Photon Emission computer tomography or Single Photon Emission computer tomography/computer tomography;
One computer, has a processing unit; And
One program product, comprises a computer readable program code, and the processing unit via this computer loads and performs the following step:
Receive the image of this Single Photon Emission computer tomography or this Single Photon Emission computer tomography/computer tomography;
Correct the movement of a patient;
Correct the scattering of this image;
Correct the decay of this image;
Recover the spatial resolution of this image;
Remove the image noise of this image; And
Calculate the tumor standard uptake value corresponding to this image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410257648.7A CN105310712A (en) | 2014-06-11 | 2014-06-11 | Method and system for conducting image measurement to obtain tumor standardized uptake value through nuclear medicine single photons |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410257648.7A CN105310712A (en) | 2014-06-11 | 2014-06-11 | Method and system for conducting image measurement to obtain tumor standardized uptake value through nuclear medicine single photons |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105310712A true CN105310712A (en) | 2016-02-10 |
Family
ID=55239718
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410257648.7A Pending CN105310712A (en) | 2014-06-11 | 2014-06-11 | Method and system for conducting image measurement to obtain tumor standardized uptake value through nuclear medicine single photons |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105310712A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105816197A (en) * | 2016-05-25 | 2016-08-03 | 刘丽 | SPECT dopamine development quantitative analysis technique and application thereof in brain dopamine assessment |
| CN110599472A (en) * | 2019-09-03 | 2019-12-20 | 佛山原子医疗设备有限公司 | Method and system for calculating SUV normalization coefficient in SPECT quantitative tomographic image |
| CN111882499A (en) * | 2020-07-15 | 2020-11-03 | 上海联影医疗科技有限公司 | PET image noise reduction method and device and computer equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101884545A (en) * | 2009-05-13 | 2010-11-17 | 株式会社东芝 | Nuclear medical imaging apparatus, image processing apparatus and image processing method |
| JP2013044615A (en) * | 2011-08-23 | 2013-03-04 | Toshiba Corp | Medical image diagnostic device, image information display device and control program |
| CN103099634A (en) * | 2011-10-14 | 2013-05-15 | 美国西门子医疗解决公司 | Method and apparatus for correcting medical imaging data |
| WO2013077963A1 (en) * | 2011-11-23 | 2013-05-30 | General Electric Company | Systems and methods for communicating dose calibration information |
-
2014
- 2014-06-11 CN CN201410257648.7A patent/CN105310712A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101884545A (en) * | 2009-05-13 | 2010-11-17 | 株式会社东芝 | Nuclear medical imaging apparatus, image processing apparatus and image processing method |
| JP2013044615A (en) * | 2011-08-23 | 2013-03-04 | Toshiba Corp | Medical image diagnostic device, image information display device and control program |
| CN103099634A (en) * | 2011-10-14 | 2013-05-15 | 美国西门子医疗解决公司 | Method and apparatus for correcting medical imaging data |
| WO2013077963A1 (en) * | 2011-11-23 | 2013-05-30 | General Electric Company | Systems and methods for communicating dose calibration information |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105816197A (en) * | 2016-05-25 | 2016-08-03 | 刘丽 | SPECT dopamine development quantitative analysis technique and application thereof in brain dopamine assessment |
| CN110599472A (en) * | 2019-09-03 | 2019-12-20 | 佛山原子医疗设备有限公司 | Method and system for calculating SUV normalization coefficient in SPECT quantitative tomographic image |
| CN111882499A (en) * | 2020-07-15 | 2020-11-03 | 上海联影医疗科技有限公司 | PET image noise reduction method and device and computer equipment |
| CN111882499B (en) * | 2020-07-15 | 2024-04-16 | 上海联影医疗科技股份有限公司 | PET image noise reduction method and device and computer equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11576628B2 (en) | Full dose PET image estimation from low-dose PET imaging using deep learning | |
| CN104161534B (en) | Dynamic single photon Emission Computed Tomography or Single Photon Emission computerized tomography/computerized tomography myocardial blood flow quantify system and method | |
| Jin et al. | List-mode reconstruction for the Biograph mCT with physics modeling and event-by-event motion correction | |
| CN106491151B (en) | PET image acquisition method and system | |
| Oddstig et al. | Reduced administered activity, reduced acquisition time, and preserved image quality for the new CZT camera | |
| Llacer et al. | Statistical analysis of maximum likelihood estimator images of human brain FDG PET studies | |
| CN105147312A (en) | PET image acquiring method and system | |
| Schwartz et al. | Repeatability of SUV measurements in serial PET | |
| US10489940B2 (en) | System and computer-implemented method for improving image quality | |
| CN104408758A (en) | Low-dose processing method of energy spectrum CT image | |
| US20160078646A1 (en) | Prca-based method and system for dynamically re-establishing pet image | |
| Park et al. | Metal artifact reduction in CT by identifying missing data hidden in metals | |
| CN104306011A (en) | Quantitative analysis technology for tumor uptake by virtue of SPECT (Single Photon Emission Computed Tomography) imaging and application to tumor evaluation | |
| CN105476653A (en) | Quantitative analysis technology of SPECT bone images and application of quantitative analysis technology to bone assessment | |
| Li et al. | Multienergy cone-beam computed tomography reconstruction with a spatial spectral nonlocal means algorithm | |
| CN105310712A (en) | Method and system for conducting image measurement to obtain tumor standardized uptake value through nuclear medicine single photons | |
| TWI547266B (en) | Method and system to quantify standardized uptake value in tumors with single photon emission computed tomography (spect) and spect/computed tomography (spect/ct) | |
| Knoll et al. | Improved scatter correction with factor analysis for planar and SPECT imaging | |
| Boisson et al. | NEMA NU 4-2008 validation and applications of the PET-SORTEO Monte Carlo simulations platform for the geometry of the Inveon PET preclinical scanner | |
| Häggström et al. | Compartment modeling of dynamic brain PET—The impact of scatter corrections on parameter errors | |
| Vaska et al. | Quantitative Imaging with the Micro‐pet Small‐Animal Pet Tomograph | |
| Thomas et al. | A dual modality approach to quantitative quality control in emission tomography | |
| Bayerlein et al. | Development of a Monte Carlo-based scatter correction method for total-body PET using the uEXPLORER PET/CT scanner | |
| Hughes et al. | A multi-center phantom study comparing image resolution from three state-of-the-art SPECT-CT systems | |
| Jian et al. | Applications of the line-of-response probability density function resolution model in PET list mode reconstruction |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20160823 Address after: Sao Zhen GUI Zhu Cun tower of 414202 cities in Hunan province Yueyang city Huarong County three Group No. 012 Applicant after: Liu Li Applicant after: Xu Bailing Address before: Taiwan Taitung County Chinese chishang Fuxing Village 8 o No. 113 Applicant before: Xu Bailing |
|
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160210 |