KR101667145B1 - Method for correction of shaking of gamma camera - Google Patents
Method for correction of shaking of gamma camera Download PDFInfo
- Publication number
- KR101667145B1 KR101667145B1 KR1020150055190A KR20150055190A KR101667145B1 KR 101667145 B1 KR101667145 B1 KR 101667145B1 KR 1020150055190 A KR1020150055190 A KR 1020150055190A KR 20150055190 A KR20150055190 A KR 20150055190A KR 101667145 B1 KR101667145 B1 KR 101667145B1
- Authority
- KR
- South Korea
- Prior art keywords
- photon
- correction reference
- gamma camera
- reference plane
- gamma
- Prior art date
Links
- 238000012937 correction Methods 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000003384 imaging method Methods 0.000 claims abstract description 20
- 238000003475 lamination Methods 0.000 claims abstract description 8
- 238000010030 laminating Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000005251 gamma ray Effects 0.000 description 13
- 230000003287 optical effect Effects 0.000 description 8
- 239000013078 crystal Substances 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 230000005855 radiation Effects 0.000 description 6
- 206010028980 Neoplasm Diseases 0.000 description 5
- 201000011510 cancer Diseases 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 210000001165 lymph node Anatomy 0.000 description 4
- 206010027476 Metastases Diseases 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009401 metastasis Effects 0.000 description 2
- 238000009206 nuclear medicine Methods 0.000 description 2
- 230000005693 optoelectronics Effects 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000019522 cellular metabolic process Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 239000012857 radioactive material Substances 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 239000010907 stover Substances 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/164—Scintigraphy
- G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
- G01T1/1642—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras using a scintillation crystal and position sensing photodetector arrays, e.g. ANGER cameras
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4258—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector for detecting non x-ray radiation, e.g. gamma radiation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Optics & Photonics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Biophysics (AREA)
- General Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine (AREA)
Abstract
The present invention relates to a shake correction method of a gamma camera, comprising: (a) setting a correction reference position value for the gamma camera; (b) a photon event generated by the gamma camera during a predetermined photographing time and a sensor position value sensed by a position sensor installed in the gamma camera are collected; (c) setting a correction reference plane on an x-y plane based on the correction reference position value, setting a three-dimensional inverse projection space having the correction reference plane as one side and the z axis direction as a depth; (d) imaging each of the photon events at a corresponding location on the correction reference plane and comparing the correction reference location value with the sensor location value corresponding to each photon event, Projecting back from the correction reference plane to the three-dimensional inverse projection space; (e) generating a gamma image by projecting the photon event back-projected in the three-dimensional inverse projection space back in the z-axis direction on a lamination plane facing the correction reference plane. Accordingly, the blur phenomenon of the gamma image caused by the operation noise such as the shake of the gamma camera can be eliminated, so that the gamma image of higher resolution and higher quality can be obtained.
Description
The present invention relates to a shake correction method of a gamma camera, and more particularly, to a shake correction method of a gamma camera capable of correcting a shake of a gamma camera applied to a gamma ray imaging apparatus.
Nuclear medicine imaging technology is a technique of reconstructing images by injecting and diffusing a small amount of radioactive material into the human body to detect gamma rays emitted by radioactive isotopes (tracers) in specific organs, bones and tissues.
Positron Emission Tomography (PET) and Single Photon Emission Computed Tomography (SPECT) are two types of gamma-ray imaging techniques. They are CT, MRI, In addition, biochemical functional images such as blood flow and cell metabolism are provided, so that cancer and disease can be diagnosed early.
The gamma ray imaging apparatus to which the gamma ray imaging technique is applied is provided with a gamma camera 300 (or gamma probe) for detecting gamma ray photons. 1 is a view showing a configuration of a
The
The radiation passing through the
The
In the nuclear medicine imaging technique using the
For example, in the radiographic image processing apparatus disclosed in Japanese Patent Application Laid-Open No. 2011-027601, a look-up-table generated in advance and a photon event are selectively obtained in comparison with the detected photon events Discloses a technique capable of quickly processing a large number of photon events without reducing the resolution and quantification.
However, in the case of the
In the case of a conventional optical camera, an optical correction method is often used to remove motion noise such as a camera shake. A sensor shift method and a lens shift method are typically used.
2 (a) conceptually illustrates the sensor shift method, and a method of moving the
In the case of the
In addition, in the case of the shake correction method used in a general optical camera image processing technique, since a method of correcting the operation noise by moving the entire shot to be photographed is used, the photon is photographed in water every few seconds, It is not suitable to apply the general shake correction method of the optical camera to the
SUMMARY OF THE INVENTION Accordingly, the present invention has been made in order to solve the above-mentioned problems, and it is an object of the present invention to provide a gamma camera which removes operational noise generated by shaking of a gamma camera, And to provide a camera shake correction method.
According to another aspect of the present invention, there is provided a shake correction method of a gamma camera, including: (a) setting a correction reference position value for the gamma camera; (b) a photon event generated by the gamma camera during a predetermined photographing time and a sensor position value sensed by a position sensor installed in the gamma camera are collected; (c) setting a correction reference plane on an x-y plane based on the correction reference position value, setting a three-dimensional inverse projection space having the correction reference plane as one side and the z axis direction as a depth; (d) imaging each of the photon events at a corresponding location on the correction reference plane and comparing the correction reference location value with the sensor location value corresponding to each photon event, Projecting back from the correction reference plane to the three-dimensional inverse projection space; (e) laminating the photon events, which are projected back to the three-dimensional inverse projection space, in the z-axis direction on a lamination plane facing the correction reference plane to generate a gamma image. Is achieved by a shake correction method.
Here, the photon event includes information on the sensing time of the photon, the energy of the photon, and the sensing position of the photon in the gamma camera; The sensor position value may include a coordinate value for the x-axis, the y-axis, and the z-axis, and an angle value for each of the x-axis, the y-axis, and the z-axis.
In the step (c), the correction reference plane is set in the form of an MxN matrix, and the three-dimensional back projection space is an M × N × K three-dimensional cubic shape having a K layer in the z- Is set; In the step (d), the energy of the photons is imaged on a row and a column on the correction reference plane corresponding to the photon sensing position of the photon event, and the sensing time of the photon of the photon event Axis, the y-axis and the z-axis direction are determined according to a deviation between the sensor position value and the correction reference position value at a corresponding time, The energy of the photon can be projected back to the layer of the three-dimensional back projection space at the back projection angle.
In the step (d), the energy of the photon is projected back to each unit cube forming the three-dimensional inverse projection space in the form of M × N × K three-dimensional cubes, and is projected back to each of the unit cubes Wherein the energy of the photons is summed in each of the unit cube units to form a unit cube photon value; In the step (e), the cube photon values of the unit cubes located in the same row and column may be summed and stacked on the corresponding rows and columns on the layer plane.
In the step (a), the correction reference position value may be set to any one of the sensor position values collected during the photographing time in the step (b).
Further, the laminating plane may be set as the correction reference plane.
According to the present invention, according to the present invention, a blur phenomenon of a gamma image caused by a motion noise such as a shake of a gamma camera is eliminated, and a gamma image of higher resolution and higher quality can be obtained .
This makes it possible to improve the accuracy of the surveillance lymph node search of the cancer. Through this, it is possible to more accurately detect the cancer metastasis, and the incision of the lymph node during surgery can be minimized.
1 is a view showing a configuration of a general gamma camera,
2 is a view for explaining examples of a method of correcting shaking in a general optical camera,
FIG. 3 is a view showing a configuration of a gamma ray imaging apparatus according to the present invention,
4 is a control flowchart for explaining a shake correction method of a gamma camera according to the present invention,
5 is a view showing an example of a three-dimensional back projection space in the shake correction method of the gamma camera according to the present invention,
6 to 10 are views for explaining the principle of a shake correction method of a gamma camera according to the present invention.
Hereinafter, the present invention will be described in detail with reference to the accompanying drawings showing embodiments of the present invention.
FIG. 3 is a view showing a configuration of a gamma
The
The
The gamma imaging processor receives the photon event from the
Hereinafter, a method of correcting the shake of the
First, the gamma
Then, photon events are collected from the
When all the photon events and the
More specifically, the gamma
When the correction reference plane (CRP) is set, a three-dimensional inverse projection space (BPA) in which the correction reference plane (CRP) is set as one side and the z axis direction is set as a depth is set. 5 is a diagram showing an example of a three-dimensional back projection space (BPA) set in the shake correction method of the
5, the correction reference plane CRP is set in the form of an M × N matrix on the xy plane, and the three-dimensional back projection space BPA is set as M X N x K three-dimensional cubes.
At this time, the number of rows and columns of the correction reference plane (CRP) forming the three-dimensional back projection space BPA is set corresponding to the resolution of the
When the three-dimensional back projection space BPA is set as described above, the photon event is imaged at a corresponding position on the correction reference plane CRP in units of photon events (S45), and the correction reference position value and the sensor position value The photon events imaged on the correction reference plane CRP are reflected from the correction reference plane CRP into the three-dimensional back projection space BPA (S46).
Hereinafter, a method of reversely projecting a photon event to a three-dimensional inverse projection space (BPA) in the shake correction method of the
First, energy of a photon is imaged on a row and column on a correction reference plane (CRP) corresponding to a sensing position of a photon event sensed at a time t. Figures 6 through 8 illustrate schematically the first and second figures from which the energy of the photons is imaged in rows and columns on the correction reference plane (CRP).
Then, the backprojection angle in the x-axis, y-axis, and z-axis directions is determined according to the deviation between the sensor position value sensed at time t and the correction reference position value. FIG. 6 is a diagram illustrating a case where there is no deviation between the sensor position value and the correction reference position value and the back projection angle of the corresponding photon is equal to the correction reference plane (CRP). FIGS. 7 and 8 show the relationship between the sensor position value and the correction reference position value And the back projection angle of the photon is changed from the correction reference plane (CRP). The angle of the correction reference plane (CRP) is changed to schematically show the change of the back projection angle.
As described above, when the back projection angle is determined, the energy of the photons imaged on the row and column on the correction reference plane (CRP) is converted to the z-axis Lt; RTI ID = 0.0 > direction. ≪ / RTI >
At this time, the energy of the photon is reversely projected to each unit cube UC forming the three-dimensional back projection space BPA in the form of M × N × K three-dimensional cubes. As shown in FIGS. 6 to 8 When the back projection for all the photon events collected during the imaging time is completed, the magnitudes of the photons reflected back to each unit cube UC are summed in units of a unit cube (UC).
That is, when the back projection for the first photon event is completed, the energy of the photon is projected back to the unit cube UC in the z-axis direction as shown in the third figure from above in Fig. 7, when the back projection for the second photon event is completed, in the unit cube (UC) where the first photon event and the second photon event are both reversely projected, the photon intensity of the two photon events Summed.
Similarly, as shown in Fig. 8, when the back projection for the third photon event is completed, the unit photon event, the second photon event, and the third photon event are all projected backward in the unit cube (UC) The intensity of the photons of the photodetector becomes summed. On the other hand, in a unit cube (UC) in which one or two photon events are reversely projected, only the intensity of the photons of one or two photon events is reflected. Hereinafter, the energy of the summed photons in one unit cube (UC) will be described as unit cube photon values.
When the back projection to the three-dimensional back projection space (BPA) is completed for all the photon events collected from the gamma camera 300 (S47), the back projection (BPA) And the generated photon events are laminated in a z-axis direction on a lamination plane opposite to the correction reference plane (CRP) to generate a two-dimensional gamma image (S48).
9, when a photon event is backprojected to each unit cube UC on a three-dimensional back projection space BPA and a photon event is added to each unit cube UC, The unit cube photon values of the unit cubes UC located in the same row and column in the back projection space BPA are summed and stacked on the corresponding rows and columns on the lamination plane.
In other words, the unit cubic photon values of the unit cubes (UC) of the corresponding row and column in the z-axis direction in the three-dimensional back projection space (BPA) having the three-dimensional cubic shape are summed, And the unit cube photon values added on the plane of the stack form a corresponding row and column on the plane of the stack, that is, a photon value at the corresponding pixel, , Gamma images can be generated. In FIG. 9, the correction reference plane (CRP) is set as the lamination plane.
According to the above configuration, the blur phenomenon of the gamma image caused by the operation noise such as the shake of the
FIG. 10 is a graphical representation of a gamma image when the shake correction method according to the present invention is not applied in the state of FIGS. 6 to 8. FIG. As shown in FIG. 10, blurring occurs in the gamma image due to shaking of the
In this manner, the blur phenomenon of the gamma image is removed to obtain a high-quality gamma image with high resolution, so that the accuracy of the surveillance lymph node search of the cancer can be improved. Through this, it is possible to more accurately detect the cancer metastasis, and the incision of the lymph node during surgery can be minimized.
Although several embodiments of the present invention have been shown and described, those skilled in the art will readily appreciate that many modifications may be made without departing from the spirit or scope of the invention . The scope of the invention will be determined by the appended claims and their equivalents.
100: gamma ray imaging device 300: gamma camera
310: collimator 320: detection crystal
330: optoelectronic amplifier array 340: position logic circuit
400: position sensor 500: gamma image processor
PA: 3D inverse projection space CRP: Correction plane
UC: Unit Cube
Claims (6)
(a) setting a correction reference position value for the gamma camera;
(b) a photon event generated by the gamma camera during a predetermined photographing time and a sensor position value sensed by a position sensor installed in the gamma camera are collected;
(c) setting a correction reference plane on the xy plane based on the correction reference position value, setting a three-dimensional inverse projection space having the correction reference plane as one side and the z axis direction as a depth;
(d) imaging each of the photon events at a corresponding location on the correction reference plane and comparing the correction reference location value with the sensor location value corresponding to each photon event, Projecting back from the correction reference plane to the three-dimensional inverse projection space;
(e) laminating the photon events, which are projected back to the three-dimensional inverse projection space, in the z-axis direction on a lamination plane facing the correction reference plane to generate a gamma image. Shake correction method.
Wherein the photon event comprises information on the sensing time of the photon, the energy of the photon and the sensing position of the photon in the gamma camera;
Wherein the sensor position value includes a coordinate value for the x-axis, the y-axis, and the z-axis, and an angle value for each of the x-axis, the y-axis, and the z-axis.
In the step (c), the correction reference plane is set in the form of an M × N matrix, and the three-dimensional back projection space is set in an M × N × K three-dimensional cubic form having a K layer in the z-axis direction ;
In the step (d)
Energy of the photons is imaged in a row and a column on the correction reference plane corresponding to a sensing position of the photon of the photon event,
Axis, the y-axis and the z-axis direction are determined according to a deviation between the sensor position value and the correction reference position value in a time corresponding to the sensing time of the photon of the photon event,
Wherein the energy of the photons imaged in the rows and columns on the correction reference plane is projected back to the layer of the three-dimensional back projection space at the back projection angle.
In the step (d), the energy of the photon is projected back to each unit cube forming the three-dimensional inverse projection space in the form of M × N × K three-dimensional cubes, The energy of each unit cube being summed to form a unit cube photon value;
In the step (e)
Wherein the cube photon values of the unit cubes located in the same row and column are summed and stacked on corresponding rows and columns on the lamination plane. .
Wherein the correction reference position value in the step (a) is set to any one of the sensor position values collected during the photographing time in the step (b).
And the lamination plane is set as the correction reference plane.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020150055190A KR101667145B1 (en) | 2015-04-20 | 2015-04-20 | Method for correction of shaking of gamma camera |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020150055190A KR101667145B1 (en) | 2015-04-20 | 2015-04-20 | Method for correction of shaking of gamma camera |
Publications (1)
Publication Number | Publication Date |
---|---|
KR101667145B1 true KR101667145B1 (en) | 2016-10-17 |
Family
ID=57250083
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020150055190A KR101667145B1 (en) | 2015-04-20 | 2015-04-20 | Method for correction of shaking of gamma camera |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101667145B1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108776322A (en) * | 2018-02-26 | 2018-11-09 | 奕瑞新材料科技(太仓)有限公司 | Radiographic source localization method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05180941A (en) * | 1991-06-07 | 1993-07-23 | Digital Scintigraphics Inc | Method and apparatus for position analysis and radioactivity nuclide camera apparatus |
JP2008245695A (en) * | 2007-03-29 | 2008-10-16 | Natl Inst Of Radiological Sciences | Image reconstruction method of tomograph, failure diagnosis method, tomograph, and management program of system matrix |
JP2011027601A (en) * | 2009-07-27 | 2011-02-10 | Dainippon Printing Co Ltd | Radiation image processor, radiation image processing method, and radiation image processing program |
KR20130121317A (en) * | 2012-04-27 | 2013-11-06 | 가천대학교 산학협력단 | A real-time motion tracking of the subject and medical imaging correction method |
-
2015
- 2015-04-20 KR KR1020150055190A patent/KR101667145B1/en active IP Right Grant
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05180941A (en) * | 1991-06-07 | 1993-07-23 | Digital Scintigraphics Inc | Method and apparatus for position analysis and radioactivity nuclide camera apparatus |
JP2008245695A (en) * | 2007-03-29 | 2008-10-16 | Natl Inst Of Radiological Sciences | Image reconstruction method of tomograph, failure diagnosis method, tomograph, and management program of system matrix |
JP2011027601A (en) * | 2009-07-27 | 2011-02-10 | Dainippon Printing Co Ltd | Radiation image processor, radiation image processing method, and radiation image processing program |
KR20130121317A (en) * | 2012-04-27 | 2013-11-06 | 가천대학교 산학협력단 | A real-time motion tracking of the subject and medical imaging correction method |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108776322A (en) * | 2018-02-26 | 2018-11-09 | 奕瑞新材料科技(太仓)有限公司 | Radiographic source localization method |
CN108776322B (en) * | 2018-02-26 | 2020-10-16 | 张岚 | Ray source positioning method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5172347B2 (en) | Reconstruction of 2D planar images of nuclear medicine by iterative constraint deconvolution | |
US7778452B2 (en) | Image reconstruction method for structuring two-dimensional planar imaging into three-dimension imaging | |
US9753146B2 (en) | Pixel identification for small pitch scintillation crystal arrays | |
JP6009755B2 (en) | Diagnostic imaging apparatus and method | |
US8294109B2 (en) | Extracting location information using difference images from a non-parallel hole collimator | |
US9453804B2 (en) | Method and apparatus for generating a representation of an internal structure of an object | |
JP5459397B2 (en) | Positron CT apparatus and timing correction method | |
JP2011503535A (en) | Indirect radiation detector | |
KR102541717B1 (en) | Gap Resolution for Linear Detector Arrays | |
JP7254322B2 (en) | Collimator for medical image processing system and its image reconstruction method | |
US9658349B2 (en) | Apparatuses for and methods of generating images | |
WO2014049667A1 (en) | Digital image processing method and imaging device | |
EP0844498B1 (en) | Radiation imaging apparatus and method | |
US6140649A (en) | Imaging attenuation correction employing simultaneous transmission/emission scanning | |
US20140319360A1 (en) | Gantry-free spect system | |
KR101667145B1 (en) | Method for correction of shaking of gamma camera | |
Zhu et al. | Fast and accurate flow measurement through dual-camera light field particle image velocimetry and ordered-subset algorithm | |
JP2007132800A (en) | Nuclear medicine diagnosis device | |
KR102020531B1 (en) | Method or Apparatus for generating a high-resolution PET(positron emission tomography) image using line gamma-ray source | |
JP6595211B2 (en) | Nuclear medicine diagnostic apparatus and nuclear medicine image processing apparatus | |
US9633423B2 (en) | Method of reduction of septal shadows for thick septa collimators | |
JP2019529874A (en) | X-ray detector with low dose efficiency in the surrounding area | |
KR102193662B1 (en) | Apparatus and method for finding radioactive hotspot using gamma camera | |
TWI742891B (en) | Planar positron emission tomography device for brain | |
JP2005114739A (en) | Nuclear medicine diagnostic equipment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20190904 Year of fee payment: 4 |