CA2405592A1 - Hand held camera with tomograhic capability - Google Patents
Hand held camera with tomograhic capability Download PDFInfo
- Publication number
- CA2405592A1 CA2405592A1 CA002405592A CA2405592A CA2405592A1 CA 2405592 A1 CA2405592 A1 CA 2405592A1 CA 002405592 A CA002405592 A CA 002405592A CA 2405592 A CA2405592 A CA 2405592A CA 2405592 A1 CA2405592 A1 CA 2405592A1
- Authority
- CA
- Canada
- Prior art keywords
- detector
- radiation
- gamma
- lesion
- camera
- 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.)
- Abandoned
Links
- 238000000034 method Methods 0.000 claims abstract description 51
- 230000005855 radiation Effects 0.000 claims abstract description 48
- 230000003902 lesion Effects 0.000 claims abstract description 28
- 230000005251 gamma ray Effects 0.000 claims abstract description 25
- 238000003384 imaging method Methods 0.000 claims abstract description 20
- 230000010354 integration Effects 0.000 claims abstract description 11
- 239000000463 material Substances 0.000 claims description 11
- 239000002699 waste material Substances 0.000 claims description 9
- 230000005670 electromagnetic radiation Effects 0.000 claims description 5
- 239000002901 radioactive waste Substances 0.000 claims description 4
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims 2
- 239000000523 sample Substances 0.000 description 11
- 238000012546 transfer Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 10
- 238000009607 mammography Methods 0.000 description 9
- 239000011159 matrix material Substances 0.000 description 9
- 238000012937 correction Methods 0.000 description 8
- 206010028980 Neoplasm Diseases 0.000 description 7
- 238000001356 surgical procedure Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 210000002307 prostate Anatomy 0.000 description 4
- 230000002285 radioactive effect Effects 0.000 description 4
- 239000000941 radioactive substance Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 230000001926 lymphatic effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- 102000012406 Carcinoembryonic Antigen Human genes 0.000 description 2
- 108010022366 Carcinoembryonic Antigen Proteins 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000001574 biopsy Methods 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 238000000151 deposition Methods 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 230000005672 electromagnetic field Effects 0.000 description 2
- 239000002927 high level radioactive waste Substances 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000012857 radioactive material Substances 0.000 description 2
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 2
- 229910052721 tungsten Inorganic materials 0.000 description 2
- 239000010937 tungsten Substances 0.000 description 2
- 210000003932 urinary bladder Anatomy 0.000 description 2
- 206010009944 Colon cancer Diseases 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 206010033128 Ovarian cancer Diseases 0.000 description 1
- 206010061535 Ovarian neoplasm Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000003190 augmentative effect Effects 0.000 description 1
- 210000001099 axilla Anatomy 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 229910052793 cadmium Inorganic materials 0.000 description 1
- BDOSMKKIYDKNTQ-UHFFFAOYSA-N cadmium atom Chemical compound [Cd] BDOSMKKIYDKNTQ-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 208000029742 colonic neoplasm Diseases 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000003414 extremity Anatomy 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- APFVFJFRJDLVQX-AHCXROLUSA-N indium-111 Chemical compound [111In] APFVFJFRJDLVQX-AHCXROLUSA-N 0.000 description 1
- 229940055742 indium-111 Drugs 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000002925 low-level radioactive waste Substances 0.000 description 1
- BOYZAERJCXIRAX-UHFFFAOYSA-N lutetium(3+);trisilicate Chemical compound [Lu+3].[Lu+3].[Lu+3].[Lu+3].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-] BOYZAERJCXIRAX-UHFFFAOYSA-N 0.000 description 1
- 210000001165 lymph node Anatomy 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 201000001441 melanoma Diseases 0.000 description 1
- 210000003739 neck Anatomy 0.000 description 1
- 201000011519 neuroendocrine tumor Diseases 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 210000004197 pelvis Anatomy 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 210000004061 pubic symphysis Anatomy 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 239000002915 spent fuel radioactive waste Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 210000004881 tumor cell Anatomy 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 210000003708 urethra Anatomy 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
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
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
A tomographic imaging system includes a moveable detector (1) or detectors capable of detecting gamma radiation (3); one or more position sensor (2) fo r determining the position and angulation of the detector(s) in relation to a gamma ray emitting source; and a computational device (9) for integrating th e position and angulation of the detector(s) with information as to the energy and distribution of gamma rays detected by the detector and deriving a three - dimensional representation of the source based on the integration. A method of imaging a radiation emitting lesion located in a volume of interest also is disclosed.
Description
TITZE OF THE INVENTION
HAND HEZD CAMERA WITH TOMOGRAPHIC CAPABILITY
BACKGROUND OF THE INVENTION
[0001 The present invention relates to a tomographic imaging system and method.
[00~2~ Intraoperative visualization of target lesions and overlying tissues can reduce the time and invasiveness of surgical procedures, which results in cost savings and reductions in surgical complications.
[0003 Currently available gamma-ray surgical guidance tools include gamma-ray sensitive nonimaging "probes'°. These non-imaging gamma-ray probes resemble classic Geiger counters in appearance. Most modern nonimaging gamma-ray probes have enhanced directional responses (unlike Geiger counters) so that the surgeon can point to structures of interest, and feature a user interface that generates squawks and whistles instead of clicks. Gamma-ray probes are utilized in surgical procedures in which patients are administered radioactive substances prior to surgery.
The radioactive substances can be injected systemically, as in the case of tumor-seeking radiotracers (e. g., carcinoembryonic antigen (CEA) analogues for ovarian cancer surgery). In the case of systemic injection, the surgeon°s goal is to detect occult nests of cancer cells for removal to increase the chances for complete tumor kill during chemotherapy. The radioactive substances can also be injected locally, in order to delineate lymphatic drainage patterns (i.e., sentinel node procedure) .
[00041 Gamma-ray surgical guidance has been attempted for several tumor types. For example, neuroendocrine tumors have been detected intraoperatively with nonimaging probes, even when the tumors were initially missed on magnetic resonance images ("MRI") and computer tomography ("CT") scans.
Colon cancer deposits also have been detected with intraoperative nonimaging probes.
[0005] Once cancer has been identified, lymphatic drainage patterns can be studied to stage the patient's disease. For this application, the radioactive substances are injected locally near the site of a known primary cancer, so that the drainage patterns to local lymph nodes can be ascertained.
According to the "sentinel node" theory, a single node stands at the entryway to more distant sites.
Subscribers to this theory attempt to identify and remove this sentinel node. By examining whether the sentinel node contains tumor cells, pathologists aim to predict whether the tumor is likely to have spread to distant locations. Sampling of the sentinel node is preferable to the traditional surgical practice of removing entire blocks of nodes, because of the reduced levels of complications following node removal.
[0006] Gamma-ray probes have become the standard of care for surgical procedures involving melanoma in the extremities (i.e., legs and hands).
However, for some surgical sites involving complex anatomic regions (i.e., head & neck, axilla, abdomen) the lack of depth and size information provided by simple nonimaging probes can reduce the efficacy of surgical guidance with gamma-ray probes. For example, lymphatic drainage patterns can in some areas of the body be quite complex, and often vary widely from patient to patient.
HAND HEZD CAMERA WITH TOMOGRAPHIC CAPABILITY
BACKGROUND OF THE INVENTION
[0001 The present invention relates to a tomographic imaging system and method.
[00~2~ Intraoperative visualization of target lesions and overlying tissues can reduce the time and invasiveness of surgical procedures, which results in cost savings and reductions in surgical complications.
[0003 Currently available gamma-ray surgical guidance tools include gamma-ray sensitive nonimaging "probes'°. These non-imaging gamma-ray probes resemble classic Geiger counters in appearance. Most modern nonimaging gamma-ray probes have enhanced directional responses (unlike Geiger counters) so that the surgeon can point to structures of interest, and feature a user interface that generates squawks and whistles instead of clicks. Gamma-ray probes are utilized in surgical procedures in which patients are administered radioactive substances prior to surgery.
The radioactive substances can be injected systemically, as in the case of tumor-seeking radiotracers (e. g., carcinoembryonic antigen (CEA) analogues for ovarian cancer surgery). In the case of systemic injection, the surgeon°s goal is to detect occult nests of cancer cells for removal to increase the chances for complete tumor kill during chemotherapy. The radioactive substances can also be injected locally, in order to delineate lymphatic drainage patterns (i.e., sentinel node procedure) .
[00041 Gamma-ray surgical guidance has been attempted for several tumor types. For example, neuroendocrine tumors have been detected intraoperatively with nonimaging probes, even when the tumors were initially missed on magnetic resonance images ("MRI") and computer tomography ("CT") scans.
Colon cancer deposits also have been detected with intraoperative nonimaging probes.
[0005] Once cancer has been identified, lymphatic drainage patterns can be studied to stage the patient's disease. For this application, the radioactive substances are injected locally near the site of a known primary cancer, so that the drainage patterns to local lymph nodes can be ascertained.
According to the "sentinel node" theory, a single node stands at the entryway to more distant sites.
Subscribers to this theory attempt to identify and remove this sentinel node. By examining whether the sentinel node contains tumor cells, pathologists aim to predict whether the tumor is likely to have spread to distant locations. Sampling of the sentinel node is preferable to the traditional surgical practice of removing entire blocks of nodes, because of the reduced levels of complications following node removal.
[0006] Gamma-ray probes have become the standard of care for surgical procedures involving melanoma in the extremities (i.e., legs and hands).
However, for some surgical sites involving complex anatomic regions (i.e., head & neck, axilla, abdomen) the lack of depth and size information provided by simple nonimaging probes can reduce the efficacy of surgical guidance with gamma-ray probes. For example, lymphatic drainage patterns can in some areas of the body be quite complex, and often vary widely from patient to patient.
MISSING AT THE TIME OF PUBLICATION
[0013] The present invention also contemplates a method of imaging a radiation emitting lesion located in a volume of,interest comprising the steps of positioning a detector that is sensitive to radiation at multiple locations in relation to the lesion to obtain information as to the energy and distribution of the radiation detected by the detector;
recording the positions and angulations of the detector in relation to the lesion; integrating the positions and angulations of the detector with the energy and distribution information; and deriving a three dimensional representation of the lesion based on the integration.
[0014] The present invention also contemplates a method of imaging a gamma ray emitting lesion in a patient comprising the steps of positioning a gamma camera that is sensitive to gamma radiation at multiple locations in relation to the lesion to obtain information as to the energy and distribution of gamma rays detected by the detector; recording the positions and angulations of the camera in relation to the lesion;
integrating the positions and angulations of the camera with the energy and distribution information; and deriving a three dimensional representation of the lesion based on the integration.
[0015] The present invention also contemplates a method of radioactive waste surveillance comprising the steps of positioning a detector at multiple locations in relation to the waste, the detector being sensitive to radiation emitted by the waste and providing information as to the energy and distribution of the radiation detected by the detector;
recording the positions and angulations of the detector in relation to the waste; integrating the positions and angulations of the detector with the source and distribution information; and deriving a three dimensional representation of the waste based on the integration.
BRTEF DESCRIPTION OF THE DRAWINGS
[0016] The benefits and advantages of the present invention will become more readily apparent to those of ordinary skill in the relevant art after reviewing the following detailed description and accompanying drawings, whereine [0017] FIG. 1 shows a schematic of the inventive imaging system using a single hand-held pinhole gamma camera.
[0018] FIG. 2 shows a comparison of the geometry for positron emission mammography backprojection and hand-held pinhole gamma camera backprojection.
[0019] FIG. 3 shows an example of a graphic user interface supplying a user with feedback about the location of a hand-held pinhole gamma camera.
[0020] FIG. 4 illustrates a comparison of simulated backprojection images derived from a hand-held pinhole gamma camera. Figure 4B is inverted with respect to Figure 4C.
[0021] FIG. 5 illustrates an efficiency correction calculation for a hand-held pinhole gamma camera.
[0022] FIG. 6 compares positron emission mammography reconstruction geometry and hand-held pinhole gamma camera backprojection geometry.
[0023] FIG. 7 illustrates an alternative embodiment of the imaging system using a multiple gamma-ray sensitive cameras to define a line of backprojection between two of the cameras.
[0024] FIG. 8 illustrates an alternative embodiment utilizing a Compton scattering approach to define a line of backprojection.
[0025] FIG. 9 illustrates an alternative embodiment of the imaging system using a multiple cameras to define a line of backprojection, in which the cameras are of different sizes.
[0026]' FIG. 10 illustrates an alternative embodiment in which two gamma-ray sensitive cameras can obtain an image of the prostate.
DETAILED DESCRIPTION OF THE INVENTION
[0027] Referring now to Figure 1, it is seen that in one embodiment of the inventive imaging system, a moveable detector 1 that is sensitive to radiation 3 emitted by a source 4 in a volume of interest 5 is provided.
[0028] In a preferred embodiment, the detector can be a hand held gamma camera that provides a two dimensional image of radiation that enters the camera through an aperture 7 and strikes material on a backplane 6, which material is sensitive to the deposition of energy from incident gamma rays.
[0029) Affixed rigidly to the camera body is a position sensor 2, or some other device for recording the position and angulation of the camera in time with respect to the source 4 and volume of interest 5. Information regarding the camera's position and angulation, and the detected radiation are transmitted to a data acquisition device 8 and sent digitally to a computer 9 or other computational device with a display also sometimes referred to as a graphical user interface (an example of which is shown in Figure 3). ' [0030] The camera 1 may contain shielding material to reduce the number of events detected on the backplane that do not traverse the aperture. The aperture may be a single hole (i.e, "pinhole") or multiple pinholes ("i.e., coded aperture"), or many pinholes in a grid ("parallel hole collimator"). The pinhole grid pattern may converge ("converging hole collimator"), diverge ("diverging hole collimator"), or slant ("slant hole collimator").
[0031] In one embodiment, a one-inch square field-of-view portable gamma camera was built, using a Hamamatsu 85900 C8 position-sensitive photomultiplier.
The camera head weighed less than two pounds, and included a tungsten housing and pinhole assembly. A
pinhole size of 2.5mm diameter was selected to provide good spatial resolution. A GSO (Gd~Si05) crystal array (2mm pitch) was affixed to the photomultiplier face.
[0032] The camera head is placed on a plastic handle in which a Polhemus electromagnetic position sensor is disposed. The Polhemus system includes an electromagnetic transmitter and one or more receivers. The receiver is about one cubic centimeter in size. The Polhemus transmitter and the camera are attached via cables to an electronics controller assembly, which includes analog-to-digital conversion and power circuitry. The Polhemus circuitry interfaces with the computer via a serial port, using ASCII
commands that are interpreted by the controller circuitry. The Polhemus Euler angle specification system is similar to the standard "Goldstein" type angles.
[0033] The inventive imaging system and method integrates information as to the energy and distribution of gamma rays detected by the detector backplane 6 as well as the position and angulation of the detector backplane 6 in time. Integrating positions and angulations of the detector backplane with source information solves several problems. For example, the angulation information provided by the imaging system is needed to perform a tomographic backprojection or reconstruction with the data generated from the detector. The angulation information also is needed to determine the depth of a source or lesion. ;Further, data is not discarded when the detector is moved to different positions and angulations and the field of view of the hand-held camera is determined by the position of the camera in relation to the lesion of interest, so that the size of the camera is not a critical determinant of the field-of-view. A small camera that can fit in close quarters can be moved by the user to provide information relating to a much large field-of-view.
[0034] Referring now to Figures 2A and 2B, it is seen that a "backprojection method" can be used to determine the position of the source and to provide a three dimensional mapping or representation of the source, which for example can be a lesion in a patient that has been injected with a gamma ray emitting radiotracer. In a backprojection, a line is drawn (or calculated) from the detector in which energy was deposited, back towards the location where the gamma ray could have originated.
[0035] Referring now to Figure 2A, it is seen that in a typical positron emission mammography geometry, a positron annihilation 14 results in rays 16 and 18 from coincident events detected on two parallel detector planes 10 and 12, which are projected onto r virtual planes (plane 1-plane 6) between the two detectors. Pixel values are incremented at the locations on each virtual plane intersected by the rays.
A coincident angle a is defined by the angle of incidence of the ray 18 on the detector plane 22.
[0036] An example of a hand-held camera backprojection is shown in Figure 2B. In a typical sequence of data and image processing as shown in Figure 2B, each event recorded by the detector 20, which includes a scintillator in an array mounted on a position-sensitive photomultiplier, or 'PMT', is assigned to a particular crystal in the detector array, ' hereafter sometimes referred to as the -'event origin".
Both the location of the event origin 28 and the location of the center of the pinhole 24 are dependent on the position and angulation of the receiver. Because the pinhole, detector array, and position sensor are all on a single rigid body, the Euler angle formalism allows transformation from the rotating and translating camera frame of reference to the stationary frame of reference as for example is described in Goldstein, Classical Mechanics, Addison-Wesley Publishing Co., Copyright 1950, pages 107-109. An algorithm implemented in a computer program utilizes the angular and position information for each detected event to calculate the location of the event origin and center of the pinhole in the stationary frame of reference.
[0037] As shown in Figure 2B, ray 26 originating at event origin 28 is projected from the detector plane 20 through the center of a pinhole 24 disposed on plane 22 and onto a set of stationary virtual planes (plane 1-plane 6). Again, pixel values are incremented at the locations on each stationary virtual plane intersected by the rays. In the'hand-held case, unlike in typical positron emission mammography as for example is shown in Figure 2A, the detector plane 20 is free to move with respect to the stationary virtual planes. The output of the Polhemus receiver that is attached to the hand-held camera can be used to determine the ray angle I3.
[0038 Referring now to Figure 3 it is seen that a graphical user interface supplies a user with a virtual image of the hand-held camera 31, Polhemus transmitter 32 and a plane 40 in the backprojection volume. A line 33 is drawn between the Polhemus transmitter 32 to the Polhemus receiver 30, which is located at the base of the hand-held camera 31.
Although the origin of the stationary frame of reference can be placed in many different locations, for clarity we are placing it at the location of the Polhemus transmitter in Figure 3. A virtual ray 38 is generated arising from the event origin 35, where the energy of the gamma ray was deposited in the backplane 34 of the camera 31. A second line 37 shows the location of the rigid body hand held camera with respect to the Polhemus receiver 30. The direction of the virtual ray is determined using the two point formula for a line.
These two points (in the stationary reference frame) include the event location 35, and the center of the pinhole 36.
[0039) The virtual ray 38 is allowed to project onto one or more virtual planes, for example plane 40. At every point on each virtual plane that the virtual ray intersects (e.g., 42 in Figure 3), the virtual plane s pixel value is incremented by one. In an alternative embodiment the pixel value may be incremented by a number different from one (e.g., in excess of one) and may factor in such weighting factors as efficiency and dwell time.
[0040] In Figures 2B and 3, the virtual planes (Planes 1-6 in Figure 2B) or plane (40 in Figure 3) are set at varying distances X from the origin of the stationary frame of reference. The origin of the state of reference is shown as 32 in Figure 3, which in this figure is also the location of the Polhemus transmitter.
Typically a physician performing a biopsy will select one coordinate to be referred to as depth, and in this case "X" is the coordinate associated with the depth.
In an alternative embodiment, Y or Z can be the coordinates designated for depth. Position sensing can be implemented by utilizing a transmitter (set at the origin of the stationary frame of reference 32 in Figure 3) that creates an electromagnetic field, and one or more receivers (31 in Figure 3) that send information about the strength and direction of the electromagnetic field in order to tell the computer where the receiver is located. Alternatively, position sensing can be implemented with optical detectors or with potentiometers.
[0041] The depth information presented by this backprojection method is similar to a tomographic image in which a structure appears sharpest in the plane in which the structure is located. The positron emission mammography backprojection method reproduces this behavior, since if a source is actually present at the location along X of a particular virtual plane, rays passing through the source will converge in the same source virtual plane and form a sharp, bright point there. For positron emission mammography, the intersection of source rays form a blur in virtual planes far from the point source's plane, and are sharpest and brightest in the source's actual plane.
[0042] The behavior of a moving handheld camera in the present invention is different from the typical behavior of a stationary gamma camera. If the handheld camera were stationary during the entire data acquisition, rays from a source would form a disk on each virtual plane, with the size of the disk depending on the distance to the detector plane of the pinhole camera. This stationary handheld camera is equivalent to a simple stationary pinhole camera projection image, and would not provide any depth information.
[0043] However, when the handheld camera is moved, the motion of the handheld camera causes a blurring in each virtual plane that does not actually contain the source. In the actual source plane, rays still converge to a point. Backprojecting from a moving pinhole camera thus has the salutary effect of distinguishing overlapping sources from one another as is shown in Figures 4A-4C.
[0044] Figure 4A depicts simulation conditions in which three coplanar sources (a,b,c) are placed in plane X=100, and one source (d) is placed in plane X=0.
[0045] Figure 4B depicts an image as seen by a pinhole camera when no camera motion occurs. The sources appear inverted because of pinhole geometry.
All four sources (a,b,c,d) contribute to the image.
Because of overlap from source d, images of sources a and b cannot be individually resolved.
[0046] Figure 4C shows a backprojected image at plane X=200 obtained with a simulation of a moving pinhole camera. Figure 4C is not inverted because the software has corrected for inversion in order to present more intuitive information to the surgeon using the moving pinhole camera. In Figure 4C, backprojection removes the image of source d (at X=0) from the images of coplanar sources a, b, c.
[0047] Because the hand-held camera by definition has no predetermined orbit, the efficiency of data collection must be calculated for each actual position and angle in the orbit. The implemented efficiency correction involves backprojecting onto virtual planes the influence of all possible events that could have been detected on the detector plane for the time that the handheld detector was in each particular position. The actual event backprojection is then divided by this efficiency backprojection to get an efficiency-corrected backprojection data set. An approximate efficiency correction for the handheld detector was developed by calculating a simulated efficiency backprojection that was dependent on the positions and angulations encountered along the hand-held camera's orbit. This efficiency correction is shown schematically in Figures 5A and 5B.
[0048] Referring now to Figures 5A and 5B, it is seen that a first-order correction for inhomogeneity in efficiency has been implemented, by calculating the efficiency of the acquisition orbit in collecting counts from a simulated uniformly distributed source filling the specified volume of interest. For each simulated event with position and angulation encountered per unit time by the hand-held camera's orbit, each element of the detector backplane projects a ray through a pinhole.
[0049] Figure 5A shows an example of three detector elements 51, 53 and 55 disposed on moving detector backplane 50 contributing projections 58, 60 and 62, which pass through moving pinhole 52. A
mathematical artifice is employed to ensure that only the set of virtual rays traversing the volume of interest 59 will enter into the image formation process.
The artifice includes drawing the boundaries 54 and 56 of the volume of interest 59, backprojecting the virtual rays 58, 60 and 62 through the volume of interest 59, and determining mathematically where the virtual rays cross the boundaries of the volume of interest. The two locations on the boundaries 54 and 56 of the volume of interest 59 are then considered to have generated new lines. These new lines can be considered to be like the lines of response generated by coincident events in a positron emission tomography ("PET") scanner. These are called "pseudo-coincident" events.
[0050 Figure 5B shows an illustration of multiple simulated projections on boundary planes and reconstruction volume. It is seen that the accumulated pseudo-coincident events from all rays that intersect the boundary planes on both sides of the reconstruction volume 64 are binned into a sinogram, just as would be the case in a PET scanner. Because the orbit of the hand-held camera is not fixed, it is necessary to correct the sinogram for the effect of variations in collection efficiency along the orbit. This can be done as follows: the efficiency of the collection of radiation by the hand-held camera is estimated by simulating events in the volume of interest numerically (e.g., through Monte-Carlo statistical methods, or through analytic approaches), and using a computer to calculate the events that would have been detected by the hand-held camera as the camera ran through the actual orbit of the hand-held camera around the volume of interest. This numerical result yields an efficiency sinogram. The sinogram from actual gamma-ray events can be divided by the efficiency sinogram before proceeding to backprojection or reconstruction.' The same correction strategy follows for multiple camera heads as well, in which simulations of events arising in the volume of interest are combined with the known orbit of the camera in order to obtain an efficiency correction sinogram or other correction factor or matrix. The simulations can be updated as the camera moves to different positions.
[0051] As shown in Figures 6A and 6B, it is possible to cast the handheld pinhole camera configuration (Figure 6B) into a geometry that is similar to positron emission mammography (Figure 6A).
This is accomplished by representing the boundaries of the volume of reconstruction prescribed by the user as if there were detectors at these boundary locations.
[0052] Figure 6A shows geometry far reconstruction of positron emission mammography data sets. A positron annihilation 74 results in coincident gamma rays 71 and 73 that are detected at coordinates (Y1;Z1) on detector 70 and at coordinates (Y2,Z2) on detector 72. A reconstruction volume 76 is defined between the two detector planes 70 and 72. A sinogram of pixels on the two detector planes is submitted for reconstruction by an iterative Monte Carlo algorithm.
[0053] Figure 6B shows a formalism that converts rays from a moving pinhole into the equivalent of the dual detector geometry depicted in Figure 6A.
Boundary planes 88 and 96 are defined, which are traversed by rays passing through the pinhole 84. The detector backplane 82 is beyond the volume of reconstruction. The event origin on the detector backplane 82 is backprojected onto the boundary 88 of the reconstruction volume 90. The location of intersection 92 with this boundary 88 is then treated as if the ray actually arose from this location 92 on the boundary volume. Thus even though the physical locations of the detectors are beyond the boundaries of the volume of reconstruction, the event data is treated as if it was from lines which originated and terminated in "virtual" detectors 88 and 96 placed at the boundary locations 92 and 94.
[0054] The pixels on the boundary planes 88 and 96 that are struck by backprojected rays are considered to arise from two "pseudo-coincident " events.
The event can be considered to be equivalent to coincident gamma rays detected at coordinates (Y1,Z1) on planar detector 88 and coordinates(Y2,Z2) on planar detector 96. A sinogram is then formed from coincident events in these planar or "pseudo" detectors 88 and 96, and the sinogram is reconstructed as in the positron emission mammography case.
[0055] Creating virtual detectors on the boundary planes, which fit the handheld camera data into a positron emission mammography-like geometry (as in Figure 2), results in a loss of some angular information, because rays that do not hit both of the virtual boundary planes are discarded. In an alternative embodiment, it is possible to use boundary surfaces of other shapes and locations than parallel planes to reduce the amount of angular information that is lost.
For example, the boundary surfaces could form a six-sided cube enclosing the reconstruction volume.
[0056] Figure 7 shows an alternative embodiment where multiple detectors 100 and 116 are used. For example, both detectors can be handheld gamma cameras sensitive to radiation of 511 keV energy emitted coincidently by positron annihilation. Alternatively, each gamma camera could be sensitive to prompt radiation emitted in two directions by Indium-111. Alternatively, one of the detectors can be a handheld gamma camera and the other detector can be a stationary detector.
[0057] A positron annihilation occurring in source 104 within a volume of interest 106 emits coincident radiation 108 that is deposited simultaneously in two gamma-ray detectors at positions 112 and 118. Volume of interest 106 can be prescribed by the user of the system, or can be defined on the basis of other mathematical reasons. Position sensors 102 and 114 are affixed to the respective detectors~100 and 116.
The detectors 100 and 116 are connected to a data acquisition device 110 and a computer 120 with display.
The data acquisition device 110 gets input from both cameras as well as from position sensors 102 and 114 mounted on both cameras, and can determine whether coincidence occurred in hardware (via an AND gate) or software (by inspection of a time-stamped list mode file) or both.
[0058] Where one gamma camera and one stationary detector are used, the position-sensing receiver can be affixed to the hand-held component alone. In that case, since the other detector components) would not move, the position of these detector components can be determined through a rigid-body relationship to the Polhemus transmitter or other position-sensing apparatus.
[0059] When the lines from many energy depositions are accumulated, they generally cross at the location of the actual source (backprojection). This process can be modeled numerically as in Figure 2B, in which a set of virtual planes are imagined, and every time a plane is struck by a backprojected line, the 1s pixel on the plane which is struck is incremented by a value. The net effect is that in the planes that contain the actual source object (i.e., source plane), the pixel at the source's location in that plane has a high value.
[0060] Tt also is possible to simultaneously detect "prompt" coincidences in the two detector heads rather than annihilation coincidences, to reduce the effect of background radiation. For annihilation coincidences, gamma rays are formed from the annihilation of positions and electrons, while for prompt coincidences gamma rays are emitted essentially simultaneously but are not approximately 180 degrees apart as for annihilation radiation.
[0061] Referring now to Figure 8, it is seen that a "Compton camera" alternatively can be used to detect radiation from a source 124. A gamma ray 126 is emitted by a source 124 located in a volume of interest 125. Energy from the gamma ray 126 is deposited into a first detector plane 128 and then into a second detector plane 134 of a Compton camera 136 with position sensor 130, taking into account angular deviation 132 due to Compton interaction. Also included in the system are a data acquisition device 138 and a computer 140 with display. As in the case of the pinhole camera (Figure 1), a set of virtual rays can be generated for each event which deposits energy in both detector planes and each of these virtual rays can be backprojected to derive an image of the source 124.
[0062] Alternatively, the camera can be a mobile nonimaging detector or probe that records the amount of radiation detected at each point in time. The term " camera" is used loosely in reference to a nonimaging detector, because no image is formed at the backplane. In the case of a nonimaging detector, the number of counts received by the backplane can be integrated with the position and angulation data to backproject rays into a volume, thereby creating an image of the source. The backprojection method is similar to the method shown in Figure 2B and subsequent figures, except that the "event location" 28 of Figure 2B is the location of the detector element of the nonimaging detector, rather than a single pixel on a detector plane as represented in Figure 2B.
[0063] Alternatively, one component of the camera may emit radiation that passes through an area of interest and is then recorded by another component of the detector.
[0064] Alternatively it is possible to register data from anatomic imaging modalities such as ultrasound or magnetic resonance imaging to demonstrate needle tips or other landmarks of interest. The registration can be accomplished using position sensors that are touched by the user to fiduciary markers visible in the anatomic imaging modality or to the naked eye, and having these position sensors relay information to the computer that is displaying the images derived from data generated by the handheld scanner.
[0065] Using the principles of the invention as discussed herein, it is possible for the single camera head embodiments of Figures 1 and 8 to be placed within a body cavity or to be used intraoperatively.
[0066] It also is possible for one of the two gamma camera heads shown in Figure 7 to be miniaturized so that one or many of the camera heads can be in a body cavity and other camera heads be outside the body cavity, or for all camera heads to be inside the body cavity. For example, one of the camera heads of Figure 7 can be small, and can be inserted via a rectal probe to examine a patient's prostate, while the second detector head of Figure 7 can be larger and can be placed anterior to the patient's pelvis. An example of such an asymmetric configuration is illustrated in Figure 9. Referring now to Figure 9, it is seen that a first camera 156, which can inserted into a body cavity or orifice, includes a position sensor 158 and a backplane 160. A second camera 142 is disposed outside the patient°s body. The second camera includes a second position sensor 146 and a second backplane 144. The cameras are connected to a data acquisition device 162 and a computer 164 with a display. A gamma-ray emitting source 152 is located in a volume of interest 150 between the two cameras. The cameras detect emitted coincident or prompt gamma rays 148 and 154.
[0067] Figure 10 illustrates the use of the camera arrangement of Figure 9 to obtain an image of the prostate. Referring now to Figure 10, it is seen that a first camera 166 is disposed anterior to the urinary bladder 170 and pubic symphysis 168. A second smaller camera 172 is disposed inside the rectal vault near the prostate 180. In some cases it may be preferable to decrease the amount of background activity detected by the camera by, for example, draining the bladder with a catheter inserted through the urethra. The second camera 172 includes a handle 176 with an aperture 178 through which an optional biopsy or intervention needle 174 can be inserted. The needle, which is thus rigidly attached to camera 172, can be used to help the user localize objects of interest in space. By using rigid-body mechanics it is possible to display a "virtual needle tip" in the same graphical user interface that is showing the detector positions and the image of the object of interest.
[0066] Alternatively, the needle need not be inserted through the handle of one of the detectors and instead can be separate from the~detector. In such a case, a position sensor can be attached to the needle or to a holder for the needle.
[0069] Alternatively the needle 174 can be replaced by some other form of focal intervention such as for example, a focused ultrasound beam, radiotherapy beam, heating element, or a light guide.
[0070] Another way of deriving an image is through reconstruction rather than backprojection. In a backprojection, when the lines cross planes that do not contain the actual source (out-of-source planes), pixels on the planes are incremented. This has the effect of making it appear that background activity is present in the out-of-source planes (i.e., that do not actually contain a source). This false background activity in each plane from sources in other planes can obscure the effect of true sources in each plane. In a reconstruction, the contributions of backprojected lines from source planes to out-of-source planes are minimized or eliminated. Thus there is less background activity outside the source plane, thereby increasing detectability of sources.
[0071] ~ The classic method of performing reconstructions is called filtered backprojection, and is employed in x-ray computed tomography. In this method projection sets are decomposed into Fourier components, which are then modified with frequency-dependent filters. The filtered components are reverse transformed and backprojected.
[0072 Because a handheld gamma camera generally does not give a full set of angular ' information about the source distribution, filtered backprojections are not as useful as they would be for x-ray computed tomography where the detector goes around the source in a complete rotation. Therefore, for limited angle projection sets, which for example are provided by the inventive imaging system when a hand held gamma camera is used, it is preferable to use an iterative reconstruction, which is generally more flexible than a filtered backprojection approach.
[0073] Iterative reconstructions take an initial estimate as to the source distribution (i.e., "estimated source distribution"), and perform a projection of this estimate onto the known positions of the detectors. The projection of the estimated source distribution is then compared to the actual projection from the source, and the estimated source distribution is modified based on this comparison. In most cases, assuming that the estimate and the method of modifying the estimated source distribution were reasonable, an estimated source distribution is arrived at that is very close to the actual source distribution.
[0074) There are several methods of modifying the estimated source distribution. One such method is the maximum likelihood estimation method ("MZE"), as for example is described in Worstell et al., "Monte Carlo-based Implementation of the Mh-EM Algorithm for 3-D PET Reconstruction,°° Proceedings IEEE Nucl. Sci.
Symp. 1998. But other methods also are contemplated.
[0075, Another way of representing the process in which the estimated source distribution would cause a particular set of projections is with a matrix.
This matrix is typically referred to as a system transfer function, or as a transition matrix. The system transfer matrix can be thought of as a method of establishing which detectors would be excited by a source at which location. The system transfer function can be assembled analytically for cameras with detectors at known fixed positions. For rings of detectors, as in conventional PET scanners or SPECT scanner heads with.
fixed orbits, it is possible to solve for system transfer functions using analytical formulas, and to apply symmetries to reduce the size of the transfer functions (so-called "ordered-subset" methods).
[0076 In the case where one or more of the detectors are hand-held, the system transfer function is not readily calculated using analytical formulas, and there are few symmetries that can be used to reduce the matrix size. Accordingly, Monte Carlo methods have been used to determine the system transfer function. Every location at which the detectors have been, represents a set of possible events, and among these possible events a Monte Carlo simulation is used to generate the system transfer function. Alternatively, the set of all possible events can be tabulated at each position in order to generate a system transfer function. Since with a handheld detector the user can vary the dwell time that is spent at each detector location, a dwell time factor also can be incorporated in the system transfer function incrementing the set of possible events at each detector location for the dwell time at that location.
[0077 Once a system transfer function is determined, iterative reconstruction methods listed above, or other iterative reconstruction methods, are used to generate estimates for the actual source distribution. For example, in one reconstruction method, a user prescribes a volume of reconstruction, so that events outside this volume are not allowed in the transition matrix. In an alternative reconstruction method, zero elements of the transition matrix are discarded using a "sparse matrix" approach.
[0078] Another application for position-integrating technology is in the surveillance of radioactive waste. There is a need for a method to describe the distribution (in three dimensions) of radioactive materials generated by, for example, a nuclear power station. These materials are occasionally buried underground as nuclear waste. High-level radioactive waste has been generated in the form of spent fuel from civilian nuclear power plants in the U.S. The task of cleaning up high-level and low-level radioactive waste has been investigated by the Environment Management program at the Department of Energy, with a Base Case estimate of $227 billion required for cleanup.
[0079] Existing tools for surveillance of decommissioned sites include Geiger counters and portable gamma cameras. Robotic devices have been made and rapid response vehicles for radioactive surveillance are also known. It is believed that no existing product provides depth information about radioactive sources.
This depth information is important not only in locating the radioactive source for removal, but also in estimating the amount of buried radioactive material.
[00807 It is contemplated that methods developed for medical applications can be used for radioactive waste surveillance, provided that the one-inch diameter gamma camera for surgical applications is replaced with a larger (e. g., 10-inch diameter) gamma camera, and the sensor (e. g., a 1-cm Polhemus electromagnetic sensor). is replaced with one or more Global Position Sensing detectors. ' [0081] It is contemplated that detectors that are sensitive to radioactive emissions other than gamma-rays also can be used. For example a gamma ray detector array can be replaced with a neutron detector such as a lutetium orthosilicate, which has a high stopping power for neutrons, and the tungsten pinhole face replaced or augmented with cadmium. Through the above extension to neutron radiation, it can be seen that the detector system of this invention can detect many forms of radiation, including gamma radiation, particulate radiation (e. g., neutrons, alpha particles) or electromagnetic radiation.
[0082] In the present disclosure, the words "a" or " an" are to be taken to include both the singular and the plural. Conversely, any reference to plural items shall, where appropriate, include the singular.
[0083] In the present disclosure, the word "handheld" is to be taken to include any method of moving a part without a fixed predetermined orbit. Thus the use of a human hand to move the part is not required, and computer-directed or servo-controlled positioning is also intended to be covered by the invention disclosure, so long as the orbit is determined according to (or influenced by) the specific requirements determined at the time of examination.
[0084] From the foregoing it will be observed that numerous modifications and variations can °be effectuated without departing from the true spirit and scope of the novel concepts of the present invention. It is to be understood that no limitation with respect to the specific embodiments illustrated is intended or should be inferred. The disclosure is intended to cover by the appended claims all such modifications as fall within the scope of the claims.
[0013] The present invention also contemplates a method of imaging a radiation emitting lesion located in a volume of,interest comprising the steps of positioning a detector that is sensitive to radiation at multiple locations in relation to the lesion to obtain information as to the energy and distribution of the radiation detected by the detector;
recording the positions and angulations of the detector in relation to the lesion; integrating the positions and angulations of the detector with the energy and distribution information; and deriving a three dimensional representation of the lesion based on the integration.
[0014] The present invention also contemplates a method of imaging a gamma ray emitting lesion in a patient comprising the steps of positioning a gamma camera that is sensitive to gamma radiation at multiple locations in relation to the lesion to obtain information as to the energy and distribution of gamma rays detected by the detector; recording the positions and angulations of the camera in relation to the lesion;
integrating the positions and angulations of the camera with the energy and distribution information; and deriving a three dimensional representation of the lesion based on the integration.
[0015] The present invention also contemplates a method of radioactive waste surveillance comprising the steps of positioning a detector at multiple locations in relation to the waste, the detector being sensitive to radiation emitted by the waste and providing information as to the energy and distribution of the radiation detected by the detector;
recording the positions and angulations of the detector in relation to the waste; integrating the positions and angulations of the detector with the source and distribution information; and deriving a three dimensional representation of the waste based on the integration.
BRTEF DESCRIPTION OF THE DRAWINGS
[0016] The benefits and advantages of the present invention will become more readily apparent to those of ordinary skill in the relevant art after reviewing the following detailed description and accompanying drawings, whereine [0017] FIG. 1 shows a schematic of the inventive imaging system using a single hand-held pinhole gamma camera.
[0018] FIG. 2 shows a comparison of the geometry for positron emission mammography backprojection and hand-held pinhole gamma camera backprojection.
[0019] FIG. 3 shows an example of a graphic user interface supplying a user with feedback about the location of a hand-held pinhole gamma camera.
[0020] FIG. 4 illustrates a comparison of simulated backprojection images derived from a hand-held pinhole gamma camera. Figure 4B is inverted with respect to Figure 4C.
[0021] FIG. 5 illustrates an efficiency correction calculation for a hand-held pinhole gamma camera.
[0022] FIG. 6 compares positron emission mammography reconstruction geometry and hand-held pinhole gamma camera backprojection geometry.
[0023] FIG. 7 illustrates an alternative embodiment of the imaging system using a multiple gamma-ray sensitive cameras to define a line of backprojection between two of the cameras.
[0024] FIG. 8 illustrates an alternative embodiment utilizing a Compton scattering approach to define a line of backprojection.
[0025] FIG. 9 illustrates an alternative embodiment of the imaging system using a multiple cameras to define a line of backprojection, in which the cameras are of different sizes.
[0026]' FIG. 10 illustrates an alternative embodiment in which two gamma-ray sensitive cameras can obtain an image of the prostate.
DETAILED DESCRIPTION OF THE INVENTION
[0027] Referring now to Figure 1, it is seen that in one embodiment of the inventive imaging system, a moveable detector 1 that is sensitive to radiation 3 emitted by a source 4 in a volume of interest 5 is provided.
[0028] In a preferred embodiment, the detector can be a hand held gamma camera that provides a two dimensional image of radiation that enters the camera through an aperture 7 and strikes material on a backplane 6, which material is sensitive to the deposition of energy from incident gamma rays.
[0029) Affixed rigidly to the camera body is a position sensor 2, or some other device for recording the position and angulation of the camera in time with respect to the source 4 and volume of interest 5. Information regarding the camera's position and angulation, and the detected radiation are transmitted to a data acquisition device 8 and sent digitally to a computer 9 or other computational device with a display also sometimes referred to as a graphical user interface (an example of which is shown in Figure 3). ' [0030] The camera 1 may contain shielding material to reduce the number of events detected on the backplane that do not traverse the aperture. The aperture may be a single hole (i.e, "pinhole") or multiple pinholes ("i.e., coded aperture"), or many pinholes in a grid ("parallel hole collimator"). The pinhole grid pattern may converge ("converging hole collimator"), diverge ("diverging hole collimator"), or slant ("slant hole collimator").
[0031] In one embodiment, a one-inch square field-of-view portable gamma camera was built, using a Hamamatsu 85900 C8 position-sensitive photomultiplier.
The camera head weighed less than two pounds, and included a tungsten housing and pinhole assembly. A
pinhole size of 2.5mm diameter was selected to provide good spatial resolution. A GSO (Gd~Si05) crystal array (2mm pitch) was affixed to the photomultiplier face.
[0032] The camera head is placed on a plastic handle in which a Polhemus electromagnetic position sensor is disposed. The Polhemus system includes an electromagnetic transmitter and one or more receivers. The receiver is about one cubic centimeter in size. The Polhemus transmitter and the camera are attached via cables to an electronics controller assembly, which includes analog-to-digital conversion and power circuitry. The Polhemus circuitry interfaces with the computer via a serial port, using ASCII
commands that are interpreted by the controller circuitry. The Polhemus Euler angle specification system is similar to the standard "Goldstein" type angles.
[0033] The inventive imaging system and method integrates information as to the energy and distribution of gamma rays detected by the detector backplane 6 as well as the position and angulation of the detector backplane 6 in time. Integrating positions and angulations of the detector backplane with source information solves several problems. For example, the angulation information provided by the imaging system is needed to perform a tomographic backprojection or reconstruction with the data generated from the detector. The angulation information also is needed to determine the depth of a source or lesion. ;Further, data is not discarded when the detector is moved to different positions and angulations and the field of view of the hand-held camera is determined by the position of the camera in relation to the lesion of interest, so that the size of the camera is not a critical determinant of the field-of-view. A small camera that can fit in close quarters can be moved by the user to provide information relating to a much large field-of-view.
[0034] Referring now to Figures 2A and 2B, it is seen that a "backprojection method" can be used to determine the position of the source and to provide a three dimensional mapping or representation of the source, which for example can be a lesion in a patient that has been injected with a gamma ray emitting radiotracer. In a backprojection, a line is drawn (or calculated) from the detector in which energy was deposited, back towards the location where the gamma ray could have originated.
[0035] Referring now to Figure 2A, it is seen that in a typical positron emission mammography geometry, a positron annihilation 14 results in rays 16 and 18 from coincident events detected on two parallel detector planes 10 and 12, which are projected onto r virtual planes (plane 1-plane 6) between the two detectors. Pixel values are incremented at the locations on each virtual plane intersected by the rays.
A coincident angle a is defined by the angle of incidence of the ray 18 on the detector plane 22.
[0036] An example of a hand-held camera backprojection is shown in Figure 2B. In a typical sequence of data and image processing as shown in Figure 2B, each event recorded by the detector 20, which includes a scintillator in an array mounted on a position-sensitive photomultiplier, or 'PMT', is assigned to a particular crystal in the detector array, ' hereafter sometimes referred to as the -'event origin".
Both the location of the event origin 28 and the location of the center of the pinhole 24 are dependent on the position and angulation of the receiver. Because the pinhole, detector array, and position sensor are all on a single rigid body, the Euler angle formalism allows transformation from the rotating and translating camera frame of reference to the stationary frame of reference as for example is described in Goldstein, Classical Mechanics, Addison-Wesley Publishing Co., Copyright 1950, pages 107-109. An algorithm implemented in a computer program utilizes the angular and position information for each detected event to calculate the location of the event origin and center of the pinhole in the stationary frame of reference.
[0037] As shown in Figure 2B, ray 26 originating at event origin 28 is projected from the detector plane 20 through the center of a pinhole 24 disposed on plane 22 and onto a set of stationary virtual planes (plane 1-plane 6). Again, pixel values are incremented at the locations on each stationary virtual plane intersected by the rays. In the'hand-held case, unlike in typical positron emission mammography as for example is shown in Figure 2A, the detector plane 20 is free to move with respect to the stationary virtual planes. The output of the Polhemus receiver that is attached to the hand-held camera can be used to determine the ray angle I3.
[0038 Referring now to Figure 3 it is seen that a graphical user interface supplies a user with a virtual image of the hand-held camera 31, Polhemus transmitter 32 and a plane 40 in the backprojection volume. A line 33 is drawn between the Polhemus transmitter 32 to the Polhemus receiver 30, which is located at the base of the hand-held camera 31.
Although the origin of the stationary frame of reference can be placed in many different locations, for clarity we are placing it at the location of the Polhemus transmitter in Figure 3. A virtual ray 38 is generated arising from the event origin 35, where the energy of the gamma ray was deposited in the backplane 34 of the camera 31. A second line 37 shows the location of the rigid body hand held camera with respect to the Polhemus receiver 30. The direction of the virtual ray is determined using the two point formula for a line.
These two points (in the stationary reference frame) include the event location 35, and the center of the pinhole 36.
[0039) The virtual ray 38 is allowed to project onto one or more virtual planes, for example plane 40. At every point on each virtual plane that the virtual ray intersects (e.g., 42 in Figure 3), the virtual plane s pixel value is incremented by one. In an alternative embodiment the pixel value may be incremented by a number different from one (e.g., in excess of one) and may factor in such weighting factors as efficiency and dwell time.
[0040] In Figures 2B and 3, the virtual planes (Planes 1-6 in Figure 2B) or plane (40 in Figure 3) are set at varying distances X from the origin of the stationary frame of reference. The origin of the state of reference is shown as 32 in Figure 3, which in this figure is also the location of the Polhemus transmitter.
Typically a physician performing a biopsy will select one coordinate to be referred to as depth, and in this case "X" is the coordinate associated with the depth.
In an alternative embodiment, Y or Z can be the coordinates designated for depth. Position sensing can be implemented by utilizing a transmitter (set at the origin of the stationary frame of reference 32 in Figure 3) that creates an electromagnetic field, and one or more receivers (31 in Figure 3) that send information about the strength and direction of the electromagnetic field in order to tell the computer where the receiver is located. Alternatively, position sensing can be implemented with optical detectors or with potentiometers.
[0041] The depth information presented by this backprojection method is similar to a tomographic image in which a structure appears sharpest in the plane in which the structure is located. The positron emission mammography backprojection method reproduces this behavior, since if a source is actually present at the location along X of a particular virtual plane, rays passing through the source will converge in the same source virtual plane and form a sharp, bright point there. For positron emission mammography, the intersection of source rays form a blur in virtual planes far from the point source's plane, and are sharpest and brightest in the source's actual plane.
[0042] The behavior of a moving handheld camera in the present invention is different from the typical behavior of a stationary gamma camera. If the handheld camera were stationary during the entire data acquisition, rays from a source would form a disk on each virtual plane, with the size of the disk depending on the distance to the detector plane of the pinhole camera. This stationary handheld camera is equivalent to a simple stationary pinhole camera projection image, and would not provide any depth information.
[0043] However, when the handheld camera is moved, the motion of the handheld camera causes a blurring in each virtual plane that does not actually contain the source. In the actual source plane, rays still converge to a point. Backprojecting from a moving pinhole camera thus has the salutary effect of distinguishing overlapping sources from one another as is shown in Figures 4A-4C.
[0044] Figure 4A depicts simulation conditions in which three coplanar sources (a,b,c) are placed in plane X=100, and one source (d) is placed in plane X=0.
[0045] Figure 4B depicts an image as seen by a pinhole camera when no camera motion occurs. The sources appear inverted because of pinhole geometry.
All four sources (a,b,c,d) contribute to the image.
Because of overlap from source d, images of sources a and b cannot be individually resolved.
[0046] Figure 4C shows a backprojected image at plane X=200 obtained with a simulation of a moving pinhole camera. Figure 4C is not inverted because the software has corrected for inversion in order to present more intuitive information to the surgeon using the moving pinhole camera. In Figure 4C, backprojection removes the image of source d (at X=0) from the images of coplanar sources a, b, c.
[0047] Because the hand-held camera by definition has no predetermined orbit, the efficiency of data collection must be calculated for each actual position and angle in the orbit. The implemented efficiency correction involves backprojecting onto virtual planes the influence of all possible events that could have been detected on the detector plane for the time that the handheld detector was in each particular position. The actual event backprojection is then divided by this efficiency backprojection to get an efficiency-corrected backprojection data set. An approximate efficiency correction for the handheld detector was developed by calculating a simulated efficiency backprojection that was dependent on the positions and angulations encountered along the hand-held camera's orbit. This efficiency correction is shown schematically in Figures 5A and 5B.
[0048] Referring now to Figures 5A and 5B, it is seen that a first-order correction for inhomogeneity in efficiency has been implemented, by calculating the efficiency of the acquisition orbit in collecting counts from a simulated uniformly distributed source filling the specified volume of interest. For each simulated event with position and angulation encountered per unit time by the hand-held camera's orbit, each element of the detector backplane projects a ray through a pinhole.
[0049] Figure 5A shows an example of three detector elements 51, 53 and 55 disposed on moving detector backplane 50 contributing projections 58, 60 and 62, which pass through moving pinhole 52. A
mathematical artifice is employed to ensure that only the set of virtual rays traversing the volume of interest 59 will enter into the image formation process.
The artifice includes drawing the boundaries 54 and 56 of the volume of interest 59, backprojecting the virtual rays 58, 60 and 62 through the volume of interest 59, and determining mathematically where the virtual rays cross the boundaries of the volume of interest. The two locations on the boundaries 54 and 56 of the volume of interest 59 are then considered to have generated new lines. These new lines can be considered to be like the lines of response generated by coincident events in a positron emission tomography ("PET") scanner. These are called "pseudo-coincident" events.
[0050 Figure 5B shows an illustration of multiple simulated projections on boundary planes and reconstruction volume. It is seen that the accumulated pseudo-coincident events from all rays that intersect the boundary planes on both sides of the reconstruction volume 64 are binned into a sinogram, just as would be the case in a PET scanner. Because the orbit of the hand-held camera is not fixed, it is necessary to correct the sinogram for the effect of variations in collection efficiency along the orbit. This can be done as follows: the efficiency of the collection of radiation by the hand-held camera is estimated by simulating events in the volume of interest numerically (e.g., through Monte-Carlo statistical methods, or through analytic approaches), and using a computer to calculate the events that would have been detected by the hand-held camera as the camera ran through the actual orbit of the hand-held camera around the volume of interest. This numerical result yields an efficiency sinogram. The sinogram from actual gamma-ray events can be divided by the efficiency sinogram before proceeding to backprojection or reconstruction.' The same correction strategy follows for multiple camera heads as well, in which simulations of events arising in the volume of interest are combined with the known orbit of the camera in order to obtain an efficiency correction sinogram or other correction factor or matrix. The simulations can be updated as the camera moves to different positions.
[0051] As shown in Figures 6A and 6B, it is possible to cast the handheld pinhole camera configuration (Figure 6B) into a geometry that is similar to positron emission mammography (Figure 6A).
This is accomplished by representing the boundaries of the volume of reconstruction prescribed by the user as if there were detectors at these boundary locations.
[0052] Figure 6A shows geometry far reconstruction of positron emission mammography data sets. A positron annihilation 74 results in coincident gamma rays 71 and 73 that are detected at coordinates (Y1;Z1) on detector 70 and at coordinates (Y2,Z2) on detector 72. A reconstruction volume 76 is defined between the two detector planes 70 and 72. A sinogram of pixels on the two detector planes is submitted for reconstruction by an iterative Monte Carlo algorithm.
[0053] Figure 6B shows a formalism that converts rays from a moving pinhole into the equivalent of the dual detector geometry depicted in Figure 6A.
Boundary planes 88 and 96 are defined, which are traversed by rays passing through the pinhole 84. The detector backplane 82 is beyond the volume of reconstruction. The event origin on the detector backplane 82 is backprojected onto the boundary 88 of the reconstruction volume 90. The location of intersection 92 with this boundary 88 is then treated as if the ray actually arose from this location 92 on the boundary volume. Thus even though the physical locations of the detectors are beyond the boundaries of the volume of reconstruction, the event data is treated as if it was from lines which originated and terminated in "virtual" detectors 88 and 96 placed at the boundary locations 92 and 94.
[0054] The pixels on the boundary planes 88 and 96 that are struck by backprojected rays are considered to arise from two "pseudo-coincident " events.
The event can be considered to be equivalent to coincident gamma rays detected at coordinates (Y1,Z1) on planar detector 88 and coordinates(Y2,Z2) on planar detector 96. A sinogram is then formed from coincident events in these planar or "pseudo" detectors 88 and 96, and the sinogram is reconstructed as in the positron emission mammography case.
[0055] Creating virtual detectors on the boundary planes, which fit the handheld camera data into a positron emission mammography-like geometry (as in Figure 2), results in a loss of some angular information, because rays that do not hit both of the virtual boundary planes are discarded. In an alternative embodiment, it is possible to use boundary surfaces of other shapes and locations than parallel planes to reduce the amount of angular information that is lost.
For example, the boundary surfaces could form a six-sided cube enclosing the reconstruction volume.
[0056] Figure 7 shows an alternative embodiment where multiple detectors 100 and 116 are used. For example, both detectors can be handheld gamma cameras sensitive to radiation of 511 keV energy emitted coincidently by positron annihilation. Alternatively, each gamma camera could be sensitive to prompt radiation emitted in two directions by Indium-111. Alternatively, one of the detectors can be a handheld gamma camera and the other detector can be a stationary detector.
[0057] A positron annihilation occurring in source 104 within a volume of interest 106 emits coincident radiation 108 that is deposited simultaneously in two gamma-ray detectors at positions 112 and 118. Volume of interest 106 can be prescribed by the user of the system, or can be defined on the basis of other mathematical reasons. Position sensors 102 and 114 are affixed to the respective detectors~100 and 116.
The detectors 100 and 116 are connected to a data acquisition device 110 and a computer 120 with display.
The data acquisition device 110 gets input from both cameras as well as from position sensors 102 and 114 mounted on both cameras, and can determine whether coincidence occurred in hardware (via an AND gate) or software (by inspection of a time-stamped list mode file) or both.
[0058] Where one gamma camera and one stationary detector are used, the position-sensing receiver can be affixed to the hand-held component alone. In that case, since the other detector components) would not move, the position of these detector components can be determined through a rigid-body relationship to the Polhemus transmitter or other position-sensing apparatus.
[0059] When the lines from many energy depositions are accumulated, they generally cross at the location of the actual source (backprojection). This process can be modeled numerically as in Figure 2B, in which a set of virtual planes are imagined, and every time a plane is struck by a backprojected line, the 1s pixel on the plane which is struck is incremented by a value. The net effect is that in the planes that contain the actual source object (i.e., source plane), the pixel at the source's location in that plane has a high value.
[0060] Tt also is possible to simultaneously detect "prompt" coincidences in the two detector heads rather than annihilation coincidences, to reduce the effect of background radiation. For annihilation coincidences, gamma rays are formed from the annihilation of positions and electrons, while for prompt coincidences gamma rays are emitted essentially simultaneously but are not approximately 180 degrees apart as for annihilation radiation.
[0061] Referring now to Figure 8, it is seen that a "Compton camera" alternatively can be used to detect radiation from a source 124. A gamma ray 126 is emitted by a source 124 located in a volume of interest 125. Energy from the gamma ray 126 is deposited into a first detector plane 128 and then into a second detector plane 134 of a Compton camera 136 with position sensor 130, taking into account angular deviation 132 due to Compton interaction. Also included in the system are a data acquisition device 138 and a computer 140 with display. As in the case of the pinhole camera (Figure 1), a set of virtual rays can be generated for each event which deposits energy in both detector planes and each of these virtual rays can be backprojected to derive an image of the source 124.
[0062] Alternatively, the camera can be a mobile nonimaging detector or probe that records the amount of radiation detected at each point in time. The term " camera" is used loosely in reference to a nonimaging detector, because no image is formed at the backplane. In the case of a nonimaging detector, the number of counts received by the backplane can be integrated with the position and angulation data to backproject rays into a volume, thereby creating an image of the source. The backprojection method is similar to the method shown in Figure 2B and subsequent figures, except that the "event location" 28 of Figure 2B is the location of the detector element of the nonimaging detector, rather than a single pixel on a detector plane as represented in Figure 2B.
[0063] Alternatively, one component of the camera may emit radiation that passes through an area of interest and is then recorded by another component of the detector.
[0064] Alternatively it is possible to register data from anatomic imaging modalities such as ultrasound or magnetic resonance imaging to demonstrate needle tips or other landmarks of interest. The registration can be accomplished using position sensors that are touched by the user to fiduciary markers visible in the anatomic imaging modality or to the naked eye, and having these position sensors relay information to the computer that is displaying the images derived from data generated by the handheld scanner.
[0065] Using the principles of the invention as discussed herein, it is possible for the single camera head embodiments of Figures 1 and 8 to be placed within a body cavity or to be used intraoperatively.
[0066] It also is possible for one of the two gamma camera heads shown in Figure 7 to be miniaturized so that one or many of the camera heads can be in a body cavity and other camera heads be outside the body cavity, or for all camera heads to be inside the body cavity. For example, one of the camera heads of Figure 7 can be small, and can be inserted via a rectal probe to examine a patient's prostate, while the second detector head of Figure 7 can be larger and can be placed anterior to the patient's pelvis. An example of such an asymmetric configuration is illustrated in Figure 9. Referring now to Figure 9, it is seen that a first camera 156, which can inserted into a body cavity or orifice, includes a position sensor 158 and a backplane 160. A second camera 142 is disposed outside the patient°s body. The second camera includes a second position sensor 146 and a second backplane 144. The cameras are connected to a data acquisition device 162 and a computer 164 with a display. A gamma-ray emitting source 152 is located in a volume of interest 150 between the two cameras. The cameras detect emitted coincident or prompt gamma rays 148 and 154.
[0067] Figure 10 illustrates the use of the camera arrangement of Figure 9 to obtain an image of the prostate. Referring now to Figure 10, it is seen that a first camera 166 is disposed anterior to the urinary bladder 170 and pubic symphysis 168. A second smaller camera 172 is disposed inside the rectal vault near the prostate 180. In some cases it may be preferable to decrease the amount of background activity detected by the camera by, for example, draining the bladder with a catheter inserted through the urethra. The second camera 172 includes a handle 176 with an aperture 178 through which an optional biopsy or intervention needle 174 can be inserted. The needle, which is thus rigidly attached to camera 172, can be used to help the user localize objects of interest in space. By using rigid-body mechanics it is possible to display a "virtual needle tip" in the same graphical user interface that is showing the detector positions and the image of the object of interest.
[0066] Alternatively, the needle need not be inserted through the handle of one of the detectors and instead can be separate from the~detector. In such a case, a position sensor can be attached to the needle or to a holder for the needle.
[0069] Alternatively the needle 174 can be replaced by some other form of focal intervention such as for example, a focused ultrasound beam, radiotherapy beam, heating element, or a light guide.
[0070] Another way of deriving an image is through reconstruction rather than backprojection. In a backprojection, when the lines cross planes that do not contain the actual source (out-of-source planes), pixels on the planes are incremented. This has the effect of making it appear that background activity is present in the out-of-source planes (i.e., that do not actually contain a source). This false background activity in each plane from sources in other planes can obscure the effect of true sources in each plane. In a reconstruction, the contributions of backprojected lines from source planes to out-of-source planes are minimized or eliminated. Thus there is less background activity outside the source plane, thereby increasing detectability of sources.
[0071] ~ The classic method of performing reconstructions is called filtered backprojection, and is employed in x-ray computed tomography. In this method projection sets are decomposed into Fourier components, which are then modified with frequency-dependent filters. The filtered components are reverse transformed and backprojected.
[0072 Because a handheld gamma camera generally does not give a full set of angular ' information about the source distribution, filtered backprojections are not as useful as they would be for x-ray computed tomography where the detector goes around the source in a complete rotation. Therefore, for limited angle projection sets, which for example are provided by the inventive imaging system when a hand held gamma camera is used, it is preferable to use an iterative reconstruction, which is generally more flexible than a filtered backprojection approach.
[0073] Iterative reconstructions take an initial estimate as to the source distribution (i.e., "estimated source distribution"), and perform a projection of this estimate onto the known positions of the detectors. The projection of the estimated source distribution is then compared to the actual projection from the source, and the estimated source distribution is modified based on this comparison. In most cases, assuming that the estimate and the method of modifying the estimated source distribution were reasonable, an estimated source distribution is arrived at that is very close to the actual source distribution.
[0074) There are several methods of modifying the estimated source distribution. One such method is the maximum likelihood estimation method ("MZE"), as for example is described in Worstell et al., "Monte Carlo-based Implementation of the Mh-EM Algorithm for 3-D PET Reconstruction,°° Proceedings IEEE Nucl. Sci.
Symp. 1998. But other methods also are contemplated.
[0075, Another way of representing the process in which the estimated source distribution would cause a particular set of projections is with a matrix.
This matrix is typically referred to as a system transfer function, or as a transition matrix. The system transfer matrix can be thought of as a method of establishing which detectors would be excited by a source at which location. The system transfer function can be assembled analytically for cameras with detectors at known fixed positions. For rings of detectors, as in conventional PET scanners or SPECT scanner heads with.
fixed orbits, it is possible to solve for system transfer functions using analytical formulas, and to apply symmetries to reduce the size of the transfer functions (so-called "ordered-subset" methods).
[0076 In the case where one or more of the detectors are hand-held, the system transfer function is not readily calculated using analytical formulas, and there are few symmetries that can be used to reduce the matrix size. Accordingly, Monte Carlo methods have been used to determine the system transfer function. Every location at which the detectors have been, represents a set of possible events, and among these possible events a Monte Carlo simulation is used to generate the system transfer function. Alternatively, the set of all possible events can be tabulated at each position in order to generate a system transfer function. Since with a handheld detector the user can vary the dwell time that is spent at each detector location, a dwell time factor also can be incorporated in the system transfer function incrementing the set of possible events at each detector location for the dwell time at that location.
[0077 Once a system transfer function is determined, iterative reconstruction methods listed above, or other iterative reconstruction methods, are used to generate estimates for the actual source distribution. For example, in one reconstruction method, a user prescribes a volume of reconstruction, so that events outside this volume are not allowed in the transition matrix. In an alternative reconstruction method, zero elements of the transition matrix are discarded using a "sparse matrix" approach.
[0078] Another application for position-integrating technology is in the surveillance of radioactive waste. There is a need for a method to describe the distribution (in three dimensions) of radioactive materials generated by, for example, a nuclear power station. These materials are occasionally buried underground as nuclear waste. High-level radioactive waste has been generated in the form of spent fuel from civilian nuclear power plants in the U.S. The task of cleaning up high-level and low-level radioactive waste has been investigated by the Environment Management program at the Department of Energy, with a Base Case estimate of $227 billion required for cleanup.
[0079] Existing tools for surveillance of decommissioned sites include Geiger counters and portable gamma cameras. Robotic devices have been made and rapid response vehicles for radioactive surveillance are also known. It is believed that no existing product provides depth information about radioactive sources.
This depth information is important not only in locating the radioactive source for removal, but also in estimating the amount of buried radioactive material.
[00807 It is contemplated that methods developed for medical applications can be used for radioactive waste surveillance, provided that the one-inch diameter gamma camera for surgical applications is replaced with a larger (e. g., 10-inch diameter) gamma camera, and the sensor (e. g., a 1-cm Polhemus electromagnetic sensor). is replaced with one or more Global Position Sensing detectors. ' [0081] It is contemplated that detectors that are sensitive to radioactive emissions other than gamma-rays also can be used. For example a gamma ray detector array can be replaced with a neutron detector such as a lutetium orthosilicate, which has a high stopping power for neutrons, and the tungsten pinhole face replaced or augmented with cadmium. Through the above extension to neutron radiation, it can be seen that the detector system of this invention can detect many forms of radiation, including gamma radiation, particulate radiation (e. g., neutrons, alpha particles) or electromagnetic radiation.
[0082] In the present disclosure, the words "a" or " an" are to be taken to include both the singular and the plural. Conversely, any reference to plural items shall, where appropriate, include the singular.
[0083] In the present disclosure, the word "handheld" is to be taken to include any method of moving a part without a fixed predetermined orbit. Thus the use of a human hand to move the part is not required, and computer-directed or servo-controlled positioning is also intended to be covered by the invention disclosure, so long as the orbit is determined according to (or influenced by) the specific requirements determined at the time of examination.
[0084] From the foregoing it will be observed that numerous modifications and variations can °be effectuated without departing from the true spirit and scope of the novel concepts of the present invention. It is to be understood that no limitation with respect to the specific embodiments illustrated is intended or should be inferred. The disclosure is intended to cover by the appended claims all such modifications as fall within the scope of the claims.
Claims (52)
1. A tomographic imaging system comprising:
a moveable detector that is capable of detecting gamma radiation;
a first position sensor for determining the position and angulation of the detector in relation to a gamma ray emitting source; and a computational device for integrating the position and angulation of the detector with information as to the energy and distribution of gamma rays detected by the detector and deriving a three dimensional representation of the source based on the integration.
a moveable detector that is capable of detecting gamma radiation;
a first position sensor for determining the position and angulation of the detector in relation to a gamma ray emitting source; and a computational device for integrating the position and angulation of the detector with information as to the energy and distribution of gamma rays detected by the detector and deriving a three dimensional representation of the source based on the integration.
2. The system of claim 1 further comprising a graphical user interface for displaying the three dimensional representation.
3. The system of claim 1, wherein the detector is a gamma camera.
4. The system of Claim 1, wherein the detector includes shielding material.
5. The system of Claim 4, wherein the shielding material is interrupted by a pinhole, to form a pinhole collimator.
6. The system of Claim 4, wherein the shielding material is interrupted by a set of pinholes, to form a coded-aperture collimator.
7. The system of Claim 4, wherein the shielding material is interrupted by a set of pinholes configured to form a parallel hole collimator.
8. The system of Claim 4, wherein the shielding material is interrupted by a set of pinholes configured to form a converging hole collimator.
9. The system of Claim 4, wherein the shielding material is interrupted by a set of pinholes configured to form a diverging hole collimator.
10. The system of Claim 4, wherein the shielding material is interrupted by a set of pinholes configured to form a slant hole collimator.
11. The system of claim 1, further comprising a second detector that is capable of detecting gamma radiation.
12. The system of claim 11, wherein the computational device integrates information as to the energy and distribution of gamma rays detected by each detector with the position and angulation of each detector.
13. The system of Claim 11, further comprising a data acquisition apparatus for determining whether the detectors are struck by the gamma rays coincidently.
14. The system of claim 11, wherein one or more of the detectors is disposed in a body cavity.
15. The system of claim 11 further comprising a needle inserted through an aperture in one of the detectors, the needle providing a virtual needle tip that can be displayed on a graphic user interface.
16. The system of claim 2 further comprising a needle with a second position sensor affixed thereto, the second position sensor capable of generating a virtual needle tip that can be displayed on the graphic user interface.
17. The system of claim 1, wherein the detector is a Compton camera.
18. A tomographic imaging system comprising:
a moveable detector that is capable of detecting radiation;
a position sensor for determining the position and angulation of the detector in relation to a radiation emitting source; and a computational device for integrating the position and angulation of the detector with information as to the energy and distribution of radiation detected by the detector and deriving a three dimensional representation of the source based on the integration.
a moveable detector that is capable of detecting radiation;
a position sensor for determining the position and angulation of the detector in relation to a radiation emitting source; and a computational device for integrating the position and angulation of the detector with information as to the energy and distribution of radiation detected by the detector and deriving a three dimensional representation of the source based on the integration.
19. The system of claim 18 further comprising a graphical user interface for displaying the three dimensional representation.
20. The system of claim 18 further comprising multiple radiation detectors.
21. The system of claim 20, wherein the position sensor records the position and angulation of one or more of the detectors in relation to the source.
22. The system of claim 20, wherein the computational device integrates the position and angulation of one or more of the detectors with information as to the source and distribution of radiation detected by one or more of the detectors and derives a three dimensional representation of the source based on the integration.
23. The system of claim 18, wherein the detector is sensitive to gamma radiation.
24. The system of claim 18, wherein the detector is sensitive to electromagnetic radiation.
25. The system of claim 18, wherein the detector is sensitive to particulate radiation.
26. The system of claim 18, wherein the detector is a handheld gamma camera
27. The system of claim 18, wherein the detector is a Compton camera.
28. A tomographic imaging system comprising:
a moveable means for detecting radiation;
a means for determining the position and angulation of the detector means in relation to a radiation emitting source; and a means for integrating the position and angulation of the detector means with information as to the energy and distribution of radiation detected by the detector means and deriving a three dimensional representation of the source based on the integration.
a moveable means for detecting radiation;
a means for determining the position and angulation of the detector means in relation to a radiation emitting source; and a means for integrating the position and angulation of the detector means with information as to the energy and distribution of radiation detected by the detector means and deriving a three dimensional representation of the source based on the integration.
29. The system of claim 28 further comprising a means for displaying the three dimensional representation as an image or set of images.
30. The system of claim 28 wherein the detector means is sensitive to gamma radiation.
31. The system of claim 28 wherein the detector means is sensitive to electromagnetic radiation.
32. The system of claim 28 wherein the detector means is sensitive to particulate radiation.
33. A method of imaging a radiation emitting lesion located in a volume of interest comprising the steps of:
positioning a detector that is sensitive to radiation at multiple locations in relation to the lesion to obtain information as to the energy and distribution of the radiation detected by the detector;
recording the positions and angulations of the detector in relation to the lesion;
integrating the positions and angulations of the detector with the energy and distribution information;
and deriving a three dimensional representation of the lesion based on the integration.
positioning a detector that is sensitive to radiation at multiple locations in relation to the lesion to obtain information as to the energy and distribution of the radiation detected by the detector;
recording the positions and angulations of the detector in relation to the lesion;
integrating the positions and angulations of the detector with the energy and distribution information;
and deriving a three dimensional representation of the lesion based on the integration.
34. The method of Claim 33, wherein the representation of the lesion is derived using backprojection.
35. The method of Claim 33, wherein the representation of the lesion is derived using reconstruction.
36. The method of Claim 35, wherein the reconstruction is performed using an iterative method.
37. The method of claim 36 wherein the iterative method includes performing a Monte Carlo simulation.
38. The method of claim 33, wherein a virtual ray within a volume of interest is defined by two points upon the boundary of the volume of interest traversed by a backprojected ray, and where the virtual ray can be reconstructed or backprojected.
39. The method of claim 33 wherein the detector is sensitive to gamma radiation and the lesion emits gamma radiation.
40. The system of claim 33 wherein the detector is sensitive to electromagnetic radiation and the lesion emits electromagnetic radiation.
41. The system of claim 33 wherein the detector is sensitive to particulate radiation and the lesion emits particulate radiation.
42. The method of claim 39 wherein the detector comprises first and second detector components.
43. The method of claim 42 including the step of generating a set of virtual gamma rays using an event location on the first detector component as one point along the ray, and an event location on the second detector component as a second point along the ray.
44. The method of claim 43, further comprising the step of determining the amount of energy deposited in each of the detector components when energy from a gamma ray is deposited in both the detection components.
45. The method of claim 42 further comprising the step of using the difference in energy deposited in each detection component to determine a possible angular deviation direction for the gamma ray.
46. The method of claim 33, further comprising the step of inserting the detector into a body cavity.
47. The method of claim 33, wherein a virtual ray within a volume of interest is defined by two points upon a boundary of the volume of interest traversed by a backprojected ray.
48. The method of claim 47 where the virtual ray can be reconstructed or backprojected.
49. A method of imaging a gamma ray emitting lesion in a patient comprising the steps of:
positioning a gamma camera that is sensitive to gamma radiation at multiple locations in relation to the lesion to obtain information as to the energy and distribution of gamma rays detected by the detector;
recording the positions and angulations of the camera in relation to the lesion;
integrating the positions and angulations of the camera with the energy and distribution information; and deriving a three dimensional representation of the lesion based on the integration.
positioning a gamma camera that is sensitive to gamma radiation at multiple locations in relation to the lesion to obtain information as to the energy and distribution of gamma rays detected by the detector;
recording the positions and angulations of the camera in relation to the lesion;
integrating the positions and angulations of the camera with the energy and distribution information; and deriving a three dimensional representation of the lesion based on the integration.
50. The method of claim 49 wherein the representation of the lesion is derived using backprojection.
51. The method of claim 49 wherein the representation of the lesion is derived using reconstruction.
52. A method of radioactive waste surveillance comprising the steps of:
positioning a detector at multiple locations in relation to the waste, the detector being sensitive to radiation emitted by the waste and providing information as to the energy and distribution of the radiation detected by the detector;
recording the positions and angulations of the detector in relation to the waste;
integrating the positions and angulations of the detector with the source and distribution information;
and deriving a three dimensional representation of the waste based on the integration.
positioning a detector at multiple locations in relation to the waste, the detector being sensitive to radiation emitted by the waste and providing information as to the energy and distribution of the radiation detected by the detector;
recording the positions and angulations of the detector in relation to the waste;
integrating the positions and angulations of the detector with the source and distribution information;
and deriving a three dimensional representation of the waste based on the integration.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US19665400P | 2000-04-12 | 2000-04-12 | |
US60/196,654 | 2000-04-12 | ||
PCT/US2001/011708 WO2001079884A2 (en) | 2000-04-12 | 2001-04-11 | Hand held camera with tomograhic capability |
Publications (1)
Publication Number | Publication Date |
---|---|
CA2405592A1 true CA2405592A1 (en) | 2001-10-25 |
Family
ID=22726279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002405592A Abandoned CA2405592A1 (en) | 2000-04-12 | 2001-04-11 | Hand held camera with tomograhic capability |
Country Status (5)
Country | Link |
---|---|
EP (1) | EP1272108A4 (en) |
JP (1) | JP2003532870A (en) |
AU (1) | AU2001255305A1 (en) |
CA (1) | CA2405592A1 (en) |
WO (1) | WO2001079884A2 (en) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ATE533071T1 (en) * | 2002-08-12 | 2011-11-15 | Naviscan Pet Systems Inc | MISSION-SPECIFIC POSITRONEMISSION TOMOGRAPHY DEVICE |
ES2204322B1 (en) | 2002-10-01 | 2005-07-16 | Consejo Sup. De Invest. Cientificas | FUNCTIONAL BROWSER. |
GB0310998D0 (en) * | 2003-05-14 | 2003-06-18 | British Nuclear Fuels Plc | Detection of solid deposits |
EP1795142B1 (en) | 2005-11-24 | 2008-06-11 | BrainLAB AG | Medical tracking system using a gamma camera |
EP2140913A1 (en) * | 2008-07-03 | 2010-01-06 | Ion Beam Applications S.A. | Device and method for particle therapy verification |
CN108095761B (en) * | 2012-03-07 | 2021-10-15 | 齐特奥股份有限公司 | Spatial alignment apparatus, spatial alignment system and method for guiding a medical procedure |
JP5918093B2 (en) | 2012-09-21 | 2016-05-18 | 日立アロカメディカル株式会社 | Radiation measurement apparatus and radiation measurement method |
US9057684B2 (en) | 2013-04-05 | 2015-06-16 | The Arizona Board Of Regents On Behalf Of The University Of Arizona | Gamma ray imaging systems and methods |
JP6190302B2 (en) * | 2014-03-28 | 2017-08-30 | 国立研究開発法人国立がん研究センター | Biological function observation apparatus and radiation therapy system |
JP6456547B1 (en) * | 2018-10-08 | 2019-01-23 | 教裕 南郷 | Imaging equipment and image display equipment that are not easily affected by radiation |
CN110264565B (en) * | 2019-05-27 | 2021-07-30 | 浙江大学 | Three-dimensional reconstruction method based on half-peak probability density distribution |
US11191515B1 (en) * | 2020-10-23 | 2021-12-07 | Siemens Medical Solutions Usa, Inc | Internal dose assessment with portable single photon emission computed tomography |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3714429A (en) * | 1970-09-28 | 1973-01-30 | Afee J Mc | Tomographic radioisotopic imaging with a scintillation camera |
US4243884A (en) * | 1978-11-09 | 1981-01-06 | Actus, Inc. | Probe assembly |
US4899054A (en) * | 1988-01-19 | 1990-02-06 | General Electric Company | Gamma camera with image uniformity by energy correction offsets |
US5170055A (en) * | 1990-07-25 | 1992-12-08 | Care Wise Medical Products Corporation | Radiation detecting biopsy probe |
US5519221A (en) * | 1992-01-22 | 1996-05-21 | Ansel M. Schwartz | Dedicated apparatus and method for emission mammography |
US5606165A (en) * | 1993-11-19 | 1997-02-25 | Ail Systems Inc. | Square anti-symmetric uniformly redundant array coded aperture imaging system |
GB9512935D0 (en) * | 1995-06-24 | 1995-08-30 | British Nuclear Fuels Plc | Arrangements for detecting gamma radiation |
US5813985A (en) * | 1995-07-31 | 1998-09-29 | Care Wise Medical Products Corporation | Apparatus and methods for providing attenuation guidance and tumor targeting for external beam radiation therapy administration |
US5690113A (en) * | 1996-06-14 | 1997-11-25 | Acuson Corporation | Method and apparatus for two dimensional ultrasonic imaging |
US6064904A (en) * | 1997-11-28 | 2000-05-16 | Picker International, Inc. | Frameless stereotactic CT scanner with virtual needle display for planning image guided interventional procedures |
JPH11352233A (en) * | 1998-06-08 | 1999-12-24 | Toshiba Corp | Nuclear medicine diagnostic equipment |
-
2001
- 2001-04-11 CA CA002405592A patent/CA2405592A1/en not_active Abandoned
- 2001-04-11 JP JP2001576494A patent/JP2003532870A/en active Pending
- 2001-04-11 EP EP01928448A patent/EP1272108A4/en not_active Withdrawn
- 2001-04-11 AU AU2001255305A patent/AU2001255305A1/en not_active Abandoned
- 2001-04-11 WO PCT/US2001/011708 patent/WO2001079884A2/en active Application Filing
Also Published As
Publication number | Publication date |
---|---|
EP1272108A4 (en) | 2011-07-13 |
WO2001079884A2 (en) | 2001-10-25 |
EP1272108A2 (en) | 2003-01-08 |
AU2001255305A1 (en) | 2001-10-30 |
JP2003532870A (en) | 2003-11-05 |
WO2001079884A3 (en) | 2002-08-29 |
WO2001079884A9 (en) | 2002-12-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6628984B2 (en) | Hand held camera with tomographic capability | |
JP5554498B2 (en) | A stand-alone mini-gamma camera with a localization system for intraoperative use | |
US5961457A (en) | Method and apparatus for radiopharmaceutical-guided biopsy | |
US20070237290A1 (en) | Patient setup using tomosynthesis techniques | |
RU2401441C2 (en) | Positron emission tomography time-of-flight list mode reconstruction using detector response function | |
CN101842806B (en) | Dirty isotope PET reconstruction | |
US20120068076A1 (en) | Portable pet scanner for imaging of a portion of the body | |
US7840052B2 (en) | Restoration of the nuclear medicine 2D planar image by iterative constrained deconvolution | |
JP2000180550A (en) | Ml-em image reconfiguration method and medical image- forming device | |
JP2016510410A (en) | Scintigraphic imaging of automated 3D patient shapes with synthetic radiation-free. | |
CA2405592A1 (en) | Hand held camera with tomograhic capability | |
Matthies et al. | Mini gamma cameras for intra-operative nuclear tomographic reconstruction | |
CN106415317B (en) | Multiple emitted energies in Single Photon Emission Computed tomography | |
EP0844498B1 (en) | Radiation imaging apparatus and method | |
US7756244B2 (en) | Systems and methods for determining object position | |
US20080069414A1 (en) | System and method for segmentation | |
US9592408B2 (en) | Dose-volume kernel generation | |
WO2003086170A2 (en) | Flexible geometries for hand-held pet and spect cameras | |
US20130096421A1 (en) | Marker identification during gamma or positron imaging with application to interventional procedures | |
Yu et al. | Compton-camera-based SPECT for thyroid cancer imaging | |
Weinberg et al. | Implementing reconstruction with hand-held gamma cameras | |
CN102525521A (en) | Flashing layered phototropic vision | |
Weinberg et al. | Flexible geometries for hand-held PET and SPECT cameras | |
Zhang | Compton gamma-ray imaging probes for prostate and breast | |
Cebula | Design and Evaluation of Novel SPECT Systems for Organ Specific Imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
FZDE | Discontinued |