Reconstruction stabilizer and active vision

Information

  • Patent Grant
  • 8837793
  • Patent Number
    8,837,793
  • Date Filed
    Monday, January 9, 2012
    13 years ago
  • Date Issued
    Tuesday, September 16, 2014
    10 years ago
Abstract
A method for stabilizing the reconstruction of an imaged volume is presented. The method includes the steps of performing an analysis of the reliability of reconstruction of a radioactive-emission density distribution of the volume from radiation detected over a specified set of views, and defining modifications to the reconstruction process and/or data collection process to improve the reliability of reconstruction, in accordance with the analysis.
Description
FIELD AND BACKGROUND OF THE INVENTION

The present invention relates to radioactive-emission measurements of a volume. More particularly, the present invention relates to the accurate reconstruction of the volume, based on measurements from non-uniform views of a volume, and on the dynamic selection of views during the data acquisition process. Of particular interest is view selection for medical imaging and/or in conjunction with medical instruments, such as guided minimally invasive surgical instruments.


Radionuclide imaging is one of the most important applications of radioactivity in medicine. Its purpose is to obtain a distribution image of a radioactively labeled substance, e.g., a radiopharmaceutical, within the body following administration thereof to a patient. Radioactive-emission imaging relies on the fact that in general, pathologies, such as malignant tumors, malfunctioning organs, and inflammations, display a level of activity different from that of healthy tissue. Thus, radiopharmaceuticals, which circulate in the blood stream, are picked up by the active pathologies to a different extent than by the surrounding healthy tissue; in consequence, the pathologies are operative as radioactive-emission sources and may be detected by radioactive-emission imaging. It will be appreciated that the pathology may appear as a concentrated source of high radiation, a hot region, as may be associated with a tumor, or as a region of low-level radiation, which is nonetheless above the background level, as may be associated with carcinoma.


A reversed situation is similarly possible. Dead tissue has practically no pick up of radiopharmaceuticals, and is thus operative as a cold region.


Thus radiopharmaceuticals may be used for identifying active pathologies as well as dead tissue.


In the discussion that follows, the term body structure is intended to include organs, portions of organs, a part of the body, and a whole body. The term organ target is intended to include pathological features within organs. These pathological features may be expressed, by radioactive-emission imaging, as any one of the following:


i. hot regions, of a radioactive emission intensity higher than the background level;


ii. regions of low-level radioactive emission intensity, which is nonetheless above the background level; and


iii cold regions, of a radioactive emission intensity, lower than the background level.


Examples of radiopharmaceuticals include monoclonal antibodies or other agents, e.g., fibrinogen or fluorodeoxyglucose, tagged with a radioactive isotope, e.g., 99Mtechnetium, 67gallium, 201thallium, 111indium, 123iodine, 125iodine and 18fluorine, which may be administered orally or intravenously. The radiopharmaceuticals are designed to concentrate in the area of a tumor, and the uptake of such radiopharmaceuticals in the active part of a tumor, or other pathologies such as an inflammation, is higher and more rapid than in the tissue that neighbors the tumor. Thereafter, a radiation-emission-measuring-probe, which may be configured for extracorporeal or intracorporeal use, is employed for locating the position of the active area. Another application is the detection of blood clots with radiopharmaceuticals such as ACUTECT from Nycomed Amersham for the detection of newly formed thrombosis in veins, or clots in arteries of the heart or brain, in an emergency or operating room. Yet other applications include radioimaging of myocardial infarct using agents such as radioactive anti-myosin antibodies, radioimaging specific cell types using radioactively tagged molecules (also known as molecular imaging), etc.


The usual preferred emission for such applications is that of gamma rays, which emission is in the energy range of approximately 11-511 KeV. Beta radiation and positrons may also be detected.


Radioactive-emission imaging is performed with a radioactive-emission-measuring detector, such as a room temperature, solid-state CdZnTe (CZT) detector, which is among the more promising that are currently available. It may be configured as a single-pixel or a multi-pixel detector. Alternatively, another solid-state detector such as CdTe, HgI, Si, Ge, or the like, or a combination of a scintillation detector (such as NaT(Tl), LSO, GSO, CsI, CaF, or the like) and a photomultiplier, or another detector as known, may be used.



1
a-1i schematically illustrate detecting units 102 and detecting blocks 101 of various geometries and constructions, and radioactive-emission-measuring probes associated with them.



FIG. 1
a schematically illustrates a detecting unit 102, formed as a single-pixel detector 104, for example, a room-temperature solid-state CdZnTe (CZT) detector, having a diameter D and a thickness τd. Both the detector diameter D, or a diameter equivalent in the case of a non-circular detector, and the detector thickness τd affect the detecting efficiency. The detector diameter determines the surface area on which radioactive emission impinges; the greater the surface area, the greater the efficiency. The detector thickness affects the stopping power of the detector. High energy gamma rays may go through a thin detector, and the probability of their detection increases with detector thickness. By itself, a single-pixel detector cannot generate an image; rather, all counts are distributed over the surface area of the detector.



FIG. 1
b schematically illustrates the detecting unit 102 with a collimator 108, formed as a single cell of a diameter D, a length L, and a septa thickness τ, attached to the detector 104. The collimator 108 may be, for example, of lead, tungsten or another material which substantially blocks gamma and beta rays.


The collimator's geometry, and specifically, the ratio of D/L, provides the detecting unit 102 with a collection angle δ analogous to a viewing angle of an optical camera. The collection angle δ limits the radioactive-emission detection to substantially only that radioactive emission, which impinges on the detector 104 after passing through a “corridor” of the collimator 108 (although in practice, some high-energy gamma rays may penetrate the collimator's walls).



FIG. 1
c schematically illustrates a block 101 of the detecting units 102, with the collimator 108, formed as a multi-cell collimator, of a cell diameter D. The collection angle δ is defined for each of the detecting units 102 in the block, and each of the detecting units 102 forms a pixel in the block 101.



FIG. 1
d schematically illustrates a radioactive-emission-measuring probe 100 which comprises several detecting units 102, of different geometries and different collection angles δ, within a housing 107.



FIGS. 1
e-1i schematically illustrate the block 101, formed as a combination of a scintillation detector (such as NaI(Tl), LSO, GSO, CsI, CaF, or the like), a collimator grid, and photomultipliers.


As seen in FIG. 1e, the block 101, having proximal and distal ends 109 and 111, respectively, vis a vis an operator (not shown), is formed of the scintillation detector 104, of a single pixel, and the collimators 108, to create the detecting units 102. A plurality of photomultipliers 103 is associated with the single pixel scintillation detector 104, and with proper algorithms, as known, their output can provide a two dimensional image of the scintillations in the single pixel scintillation detector 104. In essence, this is an Anger camera, as known.


The distal view 111 of the collimator grid is seen in FIG. 1f.


Two optional proximal views 109 of the photomultipliers 103 are seen in FIGS. 1g and 1h, as a square grid arrangement, and as an arrangement of tubes.


An Anger camera 117, of the block 101 in the housing 107 is seen in FIG. 1i.


In each of the cases of FIGS. 1a-1i, the geometry of the collimator 108 determines the collection angle δ, wherein with no collimator, the collection angle δ, is essentially a solid angle of 4π steradians. Thus, the collimator's geometry affects both the detection efficiency and the image resolution, which are defined as follows:


i. The detection efficiency is the ratio of measured radiation to emitted radiation; and


ii. The image resolution is the capability of making distinguishable closely adjacent radioactive-emission organ targets, or the capability to accurately determine the size and shape of individual radioactive-emission organ targets.


Naturally, it is desired to optimize both the detection efficiency and the image resolution. Yet, they are inversely related to each other. The detection efficiency increases with increasing collimator's collection angle, and the image resolution decreases with increasing collimator's collection angle. For example, when the ratio of D/L is 1/2, the collection angle δ is substantially 2.5 steradians, so the cell views incident radiation within the confinement of about a 2.5-steradian sector. However, when the ratio of D/L is 1/12, the collection angle δ is substantially 0.31 steradians, so the cell views incident radiation within the confinement of about a 0.31-steradian sector.


Once the emission data is obtained, the data is processed to reconstruct the intensity distribution within the measured volume. The reconstruction process is generally complex, due to the large quantity of data which must be processed in order to obtain an accurate reconstruction. The following statistical model may be used to perform reconstruction.


We assume an intensity distribution, I, defined over an input space U, where U comprises a set of basic elements (e.g., pixels in two dimensional spaces, voxels in three dimensional spaces), and I(u) is the intensity of a given basic element uεU. A detecting unit positioned on a radiation-emission-measuring-probe takes a series of measurements y=(yi)i=1T from different positions and orientations around the volume U. The geometrical and physical properties of the detecting unit, together with its position and orientation in a given measurement i, determine the detection probability φi(u) of a photon emitted from location u. Thus the effective intensity of location u as viewed by the detecting unit during measurement i is φi(u)I(u).


The random count Xi(u) of photons that are emitted from location u and detected in measurement i is modeled by a Poisson process with mean φi(i)I(u). The total count of photons detected in measurement i is thus:

yi˜Poisson(ΣuεUφi(u)I(u))  (1a)

or in matrix notation:

y=Poisson(ΦI)  (1b)

where y is the vector of measurements yi, Φ is a matrix of detection probabilities over measurements i and voxels u, and I is a vector of intensity per voxel u. The reconstruction problem is to reconstruct the intensities I from the measurements y.


The 2-D Radon transform is a mathematical relationship which may be used to reconstruct the emission intensities of volume U when the set of measurements (yt)t=1T is unconstrained. The Radon transform is not statistical and does not take into account the Poissonian nature of the counts. In addition, it models the views as line projections. The Radon transform maps the spatial domain (x,y) to the Radon domain (p,φ). For a fixed projection angle, the Radon transform is simply a projection of the object. A technique known in the art as filtered back-projection (FBP) uses a back-projection operator and the inverse of the Radon transform to reconstruct the intensity distribution in volume U from measurements (yt)t=1T.


The basic, idealized problem solved by the FBP approach is to reconstruct an image from its Radon transform. The Radon transform, when properly defined, has a well-defined inverse. However, in order to invert the transform one needs measured data spanning 180°. In many medical imaging situations, the positioning of the detecting unit relative to the emitting object is constrained, so that complete measured data is not available. Reconstruction based on filtered back-projection is therefore of limited use for medical imaging. Maximum likelihood (ML) and Maximum A Posteriori (MAP) estimation methods, which address the statistical nature of the counts, have been found to provide better image reconstructions than FBP.


Limited-angle tomography is a reconstruction technique in the related art which reconstructs an image from projections acquired over a limited range of angular directions. The success of the reconstruction process depends upon the extent of the angular range acquired compared with the angular range of the missing projections. Any reconstruction from a limited range of projections potentially results in spatial distortions (artifacts) in the image. Limited angle techniques can be applied for both the Radon transform and the statistical models, but better results are generally achieved within the statistical framework. While it is known that the severity of the artifacts increases with the increasing angular range of the missing projections, limited-angle tomography does not provide information on which projections should be used in order to most effectively reconstruct the image.


Maximum likelihood (ML) estimation is a widely used method in the related art for reconstructing an image from a constrained set of measurements. A parameterization of the generative model described above is obtained by assigning an intensity I(u) to every voxel in U. The likelihood of the observed data y=(yt)t, given the set of parameters I={I(u):uεU} is:













L


(

y

I

)


=



ln






P


(

y

I

)









=



ln




t



P


(


y
t


I

)










=





t



ln






P
(




u




x
t



(
u
)




I

)









=





t



ln






Poisson
(


y
t





u





ϕ
t



(
u
)




I


(
u
)





)









=





t



{


-



u





ϕ
y



(
u
)




I


(
u
)





+


y
t


ln




u





ϕ
t



(
u
)




I


(
u
)





-

ln


(


y
t


1

)



}









(
2
)








Note that the lower and upper bound of an indexing variable (such as voxels u and time index t) are omitted in the following description, when they are clear from the context.


There is currently no analytic way to solve Eqn. 2 for the maximum of the likelihood function. However, optimization methods that find local maxima of the likelihood are known. One such method is the Expectation-Maximization (EM) process. In EM estimation, there is no guarantee that the sequence converges to a maximum likelihood estimator. For multimodal distributions, this means that an EM algorithm will converge to a local maximum (or saddle point) of the observed data likelihood function, depending on starting values.


Since the data generated by the model is only partially observable by our measurements, a basic ingredient of the Expectation-Maximization formalism is to define a set of random variables that completely define the data generated by the model. In the current case, since YtuXt(u), the set of variables {Xu(t):uεU; t=1, . . . , T} is such a set; the generated data is x=(xt)t where xt=(xt(u))u, and the observed data y is completely determined by x. The main tool in the EM formalism is the complete data likelihood:













ln






P


(

x

I

)



=



ln




t



P


(


x
t


I

)










=





t



ln








u



Poisson


(



x
t



(
u
)






ϕ
t



(
u
)




I


(
u
)




)











=





t





u



{



-


ϕ
t



(
u
)





I


(
u
)



+



x
t



(
u
)




ln


(



ϕ
t



(
u
)




I


(
u
)



)



+

ln


(



x
t



(
u
)



1

)



}










(
3
)







Since the likelihood depends on the complete data, which is only partially observable, we take its expectation with respect to the space of the unobserved data, given the current set of hypothesized parameters (i.e. the current estimator). The result is a function Q(I|I′) which assigns likelihood to sets I of model parameters, given the current set I′, and given the observed data y:













Q


(

I


I



)


=



E


[



ln






P


(

x

I

)




y

;

I



]








=





t





u



{



-


ϕ
t



(
u
)





I


(
u
)



+


E


[




x
t



(
u
)




y
t


;

I



]




ln


(



ϕ
t



(
u
)




I


(
u
)



)



+
C

}










(
4
)








where C is a term which is independent of the intensities I. The function Q(I|I′) is maximized by the following new estimates:











I


(
u
)


=


1



t




ϕ
t



(
u
)








t



E


[




x
t



(
u
)




y
t


;

I



]





;



u


U
.







(
5
)







The expectation in Eqn. 5 is obtained as follows:













P


X
t



(
u
)





(





(
u
)




;

I



)



=





(






(
u
)



;

I



)




(




(
u
)




I



)





(



I



)



=






Poisson
(



y
t

-


x
t



(
u
)








v

u










ϕ
t



(
v
)





I




(
v
)





)






Poisson


(




(
u
)





ϕ
t




(
u
)




I




(
u
)




)






Poisson


(


y
t






ϕ
t




(
v
)




I




(
v
)





)



=

Binomial


(





(
u
)







(
u
)




I




(
u
)






ϕ
t




(
v
)




I




(
v
)





;

y
t


)








(
6
)







It follows that E[xt(u)|yt;









I


]

=


y
t






ϕ
t



(
u
)





I




(
u
)






v









ϕ
t



(
v
)





I




(
v
)







,





and hence the EM iteration is:










I


(
u
)


=


1



t








ϕ
t



(
u
)








t




y
t






ϕ
t



(
u
)





I




(
u
)






v









ϕ
t



(
v
)





I




(
v
)












(
7
)







It is provable that each EM iteration improves the likelihood. Thus, given a random starting estimator, the EM algorithm iterates the above improvement step until it converges to a local maximum of the likelihood. Several random starts increase the chance of finding a globally good estimator.


It is usually desired to maximize the expected posterior probability (given a proper prior) rather than the expected likelihood. In that case we assume a prior probability on the intensities P(I)=ΠuP(I(u)). A proper conjugate prior for the Poisson distribution is the Gamma distribution:










P


(

I


(
u
)


)


=


Gamma


(


I


(
u
)




)


=



Γ
(




I


(
u
)









(
8
)







Now the maximization is done on Q(I|I′)=E[ ln P(x|I)p(I)|y; I′]. Plugging the Gamma prior into Q, and solving for I(u), we get the following EM iteration for the maximum posterior estimation:













I


(
u
)


=




+



t







E


[




x
t



(
u
)




y
t


;

I



]





+








ϕ
t



(
u
)











=





1




+



t





ϕ
t



(
u
)







[

+



t








y
t




ϕ



(
u
)




I




(
u
)






ϕ
t




I








]



(
10
)









(
9
)







The EM update step can be formulated in matrix notation as follows. Let Φ be the matrix of the projections [φl(u)]t,u, and let I, I′, y, α and β be represented as column vectors. Eqn. 10 can be written in vector and matrix notations as:









I
=


α
+


I


·

(


Φ
T



y

Φ






I





)




β
+


Φ
T


1







(
11
)








where the explicit multiplication and division denote element-wise operations, and where 1 is a vector (of the appropriate length) consisting solely of 1's.


Limited computational resources (i.e., when the entire projection matrix Φ cannot be kept in memory) may require breaking the update computation according to a partition of Φ into a set of sub-matrices (Φi). In that case the intensities can be updated gradually (using only one sub-matrix at each step) according to the following computation:









I
=


α
+


I


·



i








Φ
i
T




y
i



Φ
i



I









β
+



i








Φ
i
T


1








(
12
)








where yi is the vector of observations that are obtained using the views of Φi.


In the context of image reconstruction from a constrained set of views, the utility of the reconstruction algorithms described above is limited. Many estimation algorithms, including EM and the Radon transform, require measurements from a complete set of views surrounding the imaged object. Although algorithms for EM with missing views have been developed, these algorithms are based on equally spaced views surrounding the imaged object, of which a number of views are not available. These algorithms do not provide a generalized solution for an unconstrained set of views, which may not be equally spaced or available for all directions surrounding the object.


Singular value decomposition (SVD) is a known technique for factorizing a rectangular real or complex matrix, with applications in signal processing and statistics. SVD may be considered a generalization of Eigenvalue decomposition to m*n matrices, whereas Eigenvalue decomposition is applicable only to square matrices.


SVD states that given the m-by-n matrix M whose entries are either from the field of real numbers or the field of complex numbers, there exists a factorization of the form:

M=UDVT  (13)

where U is an m-by-m unitary matrix, D is m-by-n with nonnegative numbers on the diagonal and zeros off the diagonal, and VT, the conjugate transpose of V, is an n-by-n unitary matrix. The elements along the diagonal of D are denoted the singular values. Such a factorization is called a singular-value decomposition of M. A common convention is to order the singular values Di,i in non-increasing fashion, so that the diagonal matrix D is uniquely determined by M.


The condition number of a matrix is defined as the ratio of the matrix's largest singular value to its smallest singular value. In numerical analysis, the condition number associated with a problem is a measure of that problem's amenability to digital computation. A problem with a low condition number is said to be well-conditioned, while a problem with a high condition number is said to be ill-conditioned. For example, the condition number associated with the linear equation x=My gives a bound on how accurate the solution y will be after approximate solution.


SVD may be employed for the solution of linear inverse problems. The inverse problem for a set of linear equations, x=My, is to calculate vector y from a known x and M. Using SVD, with M=UDV*, the problem may be restated as x=UDV*y. If D is an invertible matrix, y is obtained as:

y=VD−1UTx  (14)

where U and Vt are unitary matrices, and D is a diagonal matrix containing the singular values of M. Since U and Vt are easily transposable, a solution to Eqn. 13 is obtainable only if D is invertible. As D is a diagonal matrix of singular values e1 to en (where ek denotes the k-th singular value Dk,k), D−1 is a diagonal matrix of the reciprocals of the singular values. That is:










D

-
1


=

[




e
1

-
1



























e
2

-
1




















































e
n

-
1





]





(
15
)







If the dimension of y is larger than the dimension of x, y cannot be determined directly but rather can only be estimated, which is often performed by an iterative procedure. If the condition number of D is very large, the multiplicative factors ei−1 will vary greatly, thus multiplying any measurement or calculation errors based on individual or linear combinations of elements of y, and decreasing the likelihood of convergence of the iterative estimation procedure.


Truncated SVD is a known technique for reducing sensitivity to inaccuracies or noise when solving a set of linear equations. In truncated SVD, the constraints associated with the smaller singular values are eliminated from the estimation. In terms of the inverse problem of Eqn. 14, this is accomplished for the i-th singular value by setting equal to zero in matrix D−1. Thus the lower valued singular values no longer affect the estimation process.


In order to achieve a reconstructed image which is adequate for medical diagnostic and treatment purposes, a reliable, high-resolution image of the tested object (i.e. body structure) must be obtained. Currently, reliable reconstruction algorithms are available only for complete data sets, which provide coverage of the entire volume. Such data is generally not available during medical imaging. Additionally, when high-resolution detecting units are used, their efficiency is relatively low, and the detecting units must remain at each position for a relatively long time in order to achieve a high probability of detection. Since during medical testing, measurements are generally performed at many locations as the detecting unit is moved relative to the observed body structure, the testing procedure generally requires a long time and is physically and emotionally difficult for the patient. Additionally, reconstruction is based upon a large quantity of data, and is a lengthy and computationally complex process.


There is thus a widely recognized need for, and it would be highly advantageous to have, an apparatus, system and method devoid of the above limitations.


SUMMARY OF THE INVENTION

According to a first aspect of the present invention there is provided a method for stabilizing the reconstruction of an imaged volume. The method includes the steps of performing an analysis of the reliability of reconstruction of a radioactive-emission density distribution of the volume from radiation detected over a specified set of views, and defining modifications to the reconstruction process and/or data collection process to improve the reliability of reconstruction, in accordance with the analysis.


According to a second aspect of the present invention there is provided a reconstruction stabilizer, for improving the reliability of reconstruction of an imaged volume, according to a preferred embodiment of the present invention. The reconstruction stabilizer includes a reliability analyzer for performing an analysis of the reliability of reconstruction of a radioactive-emission density distribution of the volume from radiation detected over a specified set of views, and a modifier associated with the reliability analyzer, for defining modifications to at least one of a reconstruction process and a data collection process to improve the reliability of reconstruction, in accordance with the analysis.


According to a third aspect of the present invention there is provided a system for generating a three-dimensional image of volume, according to a preferred embodiment of the present invention. The system includes a radiological imaging camera comprising a plurality of detectors configured for independent movement during data acquisition, and configured for detecting radiation emitted from the volume thereby to provide radiation data, a reconstructor, configured for performing an analysis of the radiation data so as to reconstruct a three-dimensional image of the volume, and a reconstruction stabilizer associated with the camera and the reconstructor. The reconstruction stabilizer includes a reliability analyzer for performing an analysis of the reliability of reconstruction of a radioactive-emission density distribution of the volume from radiation detected over a specified set of views, and a modifier associated with the reliability analyzer, for defining modifications to at least one of a reconstruction process and a data collection process to improve the reliability of reconstruction, in accordance with the analysis, and for providing the modifications to at least one of the camera and the reconstructor.


According to a fourth aspect of the present invention there is provided a method of radioactive-emission measurements of a body structure, according to a preferred embodiment of the present invention. The method includes the steps of performing radioactive-emission measurements of the body structure, at a predetermined set of views, analyzing the radioactive-emission measurements, and dynamically defining further views for measurements, based on the analyzing.


According to a fifth aspect of the present invention there is provided a measurement unit for performing radioactive-emission measurements of a body structure, according to a preferred embodiment of the present invention. The measurement unit includes a probe for performing the radioactive-emission measurements of the body structure, where the probe is controllable to perform the measurements at a predetermined set of views, an analysis unit for analyzing the radioactive-emission measurements, and a view definer for dynamically defining further views for measurements, based on the analyzing.


The present invention successfully addresses the shortcomings of the presently known configurations by providing methods for stabilizing the reconstruction of an imaged volume and of performing radioactive-emission measurements of a body structure.


Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, suitable methods and materials are described below. In case of conflict, the patent specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and not intended to be limiting.


Implementation of the method and system of the present invention involves performing or completing selected tasks or steps manually, automatically, or a combination thereof. Moreover, according to actual instrumentation and equipment of preferred embodiments of the method and system of the present invention, several selected steps could be implemented by hardware or by software on any operating system of any firmware or a combination thereof. For example, as hardware, selected steps of the invention could be implemented as a chip or a circuit. As software, selected steps of the invention could be implemented as a plurality of software instructions being executed by a computer using any suitable operating system. In any case, selected steps of the method and system of the invention could be described as being performed by a data processor, such as a computing platform for executing a plurality of instructions.





BRIEF DESCRIPTION OF THE DRAWINGS

The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.


The invention is herein described, by way of example only, with reference to the accompanying drawings. With specific reference now to the drawings in detail, it is stressed that the particulars shown are by way of example and for purposes of illustrative discussion of the preferred embodiments of the present invention only, and are presented in the cause of providing what is believed to be the most useful and readily understood description of the principles and conceptual aspects of the invention. In this regard, no attempt is made to show structural details of the invention in more detail than is necessary for a fundamental understanding of the invention, the description taken with the drawings making apparent to those skilled in the art how the several forms of the invention may be embodied in practice.


In the drawings:



FIGS. 1
a-1i show detecting units and blocks of various geometries and constructions and radioactive-emission-measuring probes, associated with them.



FIGS. 1
j and 1k pictorially illustrate a view and viewing parameters associated with it, in accordance with definitions of the present invention.



FIG. 2
a is a simplified flowchart of a method for stabilizing the reconstruction of an imaged volume, according to a preferred embodiment of the present invention.



FIG. 2
b is a simplified flowchart of a method for analyzing a detection probability matrix, according to a preferred embodiment of the present invention.



FIG. 2
c is a simplified block diagram of a reconstruction stabilizer, according to a preferred embodiment of the present invention.



FIG. 2
d is a simplified flowchart of a method of performing radioactive-emission measurements of a body structure, according to a preferred embodiment of the present invention.



FIG. 3 shows an object shaped as a cylinder with a front protrusion, and having a high-emittance portion (hotspot).



FIG. 4
a illustrates an object having two high-emission regions of interest.



FIG. 4
b illustrates the added information provided by each of views VA to VF.



FIGS. 5
a and 5b are simplified flowcharts of iterative methods of performing radioactive-emission measurements of a body structure, according to a first and a second preferred embodiment of the present invention.



FIGS. 6
a and 6b are simplified flowcharts of methods for dynamically defining further views, according to a first and a second preferred embodiment of the present invention.



FIGS. 7
a and 7b present the principles of modeling, for obtaining an optimal set of views, in accordance with the present invention;



FIG. 8 shows four models of a prostate emittance model.



FIGS. 9
a-9f show emittance models of a given volume to illustrate view selection using the separability criterion.



FIGS. 9
g-9i show emittance models of a given volume to illustrate view selection using a weighted-combination criterion.



FIG. 10 is a simplified flowchart of an iterative method for selecting further views, according to a preferred embodiment of the present invention.



FIG. 11 is a simplified flowchart of a single iteration of a view selection method, according to a preferred embodiment of the present invention.



FIG. 12 is a simplified flowchart of a method for dynamically defining further views, according to a third preferred embodiment of the present invention.



FIG. 13 is a simplified block diagram of measurement unit for performing radioactive-emission measurements of a body structure, according to a preferred embodiment of the present invention.



FIG. 14 is a simplified flowchart of a method for selecting an optimal initial predetermined set of views of a volume to be imaged, according to a preferred embodiment of the present invention.



FIG. 15 shows an example of a volume with two views.



FIG. 16 illustrates the concept of uniform coverage of a volume.



FIG. 17 is a simplified flowchart of method for stabilizing the reconstruction of an imaged volume with dynamic view definition, according to a preferred embodiment of the present invention.





DESCRIPTION OF THE PREFERRED EMBODIMENTS

The present embodiments teach providing modifications of the reconstruction and/or imaging processes in order to obtain a reliable reconstruction of the imaged volume. Specifically, the methods teach analyzing the reliability of the reconstruction obtainable from collected emission data over a set of views to determine modifications which are expected to improve reliability. The present embodiments further teach using radioactive-emission measurements to define views for further radioactive-emission measurements of a body structure, to be performed during the current measurement process.


With non-uniform scanning, the amount of information available for different voxels is not uniform, thereby constraining the ability to obtain an accurate reconstruction of the volume. In addition, there may be a need to focus on a region or regions of interest, and to control reconstruction resolution and accuracy where and when necessary. In such cases it is important to provide features and algorithmic components which handle or compensate for the lack of uniformity of information, and are capable of focusing on features or regions of interest, while maintaining short and efficient data acquisition.


In the context of medical imaging the imaged volume corresponds to a body structure, which may include a whole body, portion of a body, target organ and so forth. Those non-limiting embodiments presented below which are directed at imaging a body structure are to be understood as applying to any imaged volume.


Before explaining at least one embodiment of the invention in detail, it is to be understood that the invention is not limited in its application to the details of construction and the arrangement of the components set forth in the following description or illustrated in the drawings. The invention is capable of other embodiments or of being practiced or carried out in various ways. Also, it is to be understood that the phraseology and terminology employed herein is for the purpose of description and should not be regarded as limiting.


Referring now to the drawings, FIGS. 1j and 1k pictorially illustrate a view and viewing parameters associated with it, in accordance with definitions of the present invention.


Seen in FIG. 1j is a volume U, subdivided into voxels u. The volume U is defined in a six-degree coordinate system (x;y;z;ω;θ;σ) and has a point of origin P0(x0; y0; z0; ω0; θ0; σ0). A detecting unit 102 is positioned at a location and orientation P1(x1; y1; z1; ω1; θ1; σ1). The detecting unit 102 has a detector 104 of a specific detector material of a thickness τd, and a collimator 108 of a diameter D and a length L, so as to define a collection angle δ. The location, orientation, and collection angle parameters determine a three-dimensional sector S, which is the portion of the volume U that is within the detector unit's field of view. The probability of detection of a given voxel outside of sector S is negligible, since a photon emitted from the given voxel outside of sector S will have a very low probability of reaching detector 104 through collimator 108.



FIG. 1
k schematically illustrates the emission rate of the volume U, as a function of time, given that a radioactive material of a specific half-life has been administered at a time T0.


A view may thus be defined as a group of nonzero probabilities of detecting a radioactive emission associated with all the voxels that form sector S (FIG. 1j).


A view is sometimes referred to as a projection, and the two terms are synonymous. Furthermore, a view defined over the sector S can be naturally extended to be defined over the group of all voxels in the volume U, by simply associating a zero probability with every voxel outside the sector S. This enables applying mathematical operations over the entire volume U.


The viewing parameters, which are the factors affecting the detection of radioactive emissions, are as follows:


i. Location and Orientation Parameters:






    • a location and an orientation in a six-dimensional space, P1(x1; y1; z1; ω1; θ1; σ1), with respect to the origin P0(x0; y0; z0; ω0; θ0; σ0) of the volume U, in which the detecting unit 102 is positioned;


      ii. Detector-Unit Parameters:

    • the detector unit geometry (e.g. length L, diameter D, and/or collection angle δ), which together with the location and orientation parameters, P1(x1; y1; z1; ω1; θ1; σ1) with respect to the origin P0(x0; y0; z0; ω0; θ0; σ0) define the sector S;

    • the septa thickness τ, which affects the probability that a photon that enters the collimator will reach the detector as well as crosstalk effects (which occur when a photon which entered a neighboring cell penetrates the collimator and reaches the detector), hence, the detector efficiency

    • the detector material, which affects the detector efficiency; and

    • the detector thickness τd, which affects the detector's stopping power, hence, its efficiency;


      iii. Attenuation Parameters:

    • attenuation properties of all the voxels within the sector S, as they affect the probabilities that radioactive emissions from a specific voxel within the sector S will reach the detector, wherein different voxels within the sector S may have different attenuation properties, since several types of tissue may be involved;


      iv. Time Parameters:

    • since the radiopharmaceutical decays with a specific half-life, the time t1 since administration, and the duration of the measurement Δt1, affect the number of emissions that occur during the radioactive-emission measurement.


      v. Radiopharmaceutical Parameters:





The half-life t1/2, of the radiopharmaceutical, the types of radioactive emission, whether gamma or beta, and the energies of the radioactive emission affect the probability of detection.


Some of these viewing parameters are fixed for a particular situation. Specifically, the tissue attenuation parameters are given. Additionally, the time t1 since administration of the radiopharmaceutical is generally governed by the blood pool radioactivity, since it is generally necessary to wait until the blood pool radioactivity dies out for low-level detection to be possible. For the remaining viewing parameters, optimization may be carried out.


To recapitulate the problem described above, an intensity distribution/, in terms of radioactive emissions per seconds, is defined over the volume U, forming our input space. Volume U comprises a set of basic elements u (e.g., pixels in two dimensional spaces, voxels in three dimensional spaces), and I(u) is the intensity in a given basic element uεU. A view (also denoted a projection) φεΦ is defined by the set of probabilities {φ(u):uεU}, where φ(u) is the probability of detecting a radioactive emission from a voxel u, as defined by viewing parameters, such as the physical and geometrical properties of the detecting unit, as well as the attenuation parameters of the viewed volume U, and the time parameters. A measurement is obtained by choosing a view φεΦ, and then sampling according to the viewing parameters.


As shown in Eqn. 1b, the relationship between the measurements, y, and the emission intensities over the body structure may be represented as y=Poisson(ΦI). The reconstruction problem we are faced with is to calculate the intensity vector I from the measurement vector y, given a known probability matrix Φ.


Applying SVD analysis to the probability matrix Φ, we obtain:

Φ=UDVt=>y=Poisson(UDVtI)  (16)

where U and Vt are unitary matrices, and D is a diagonal matrix containing the singular values of Φ. Any destabilizing singular values in D are likewise destabilizing to the Poisson process. Specifically, those elements of I (or linear combinations of elements) for which insufficient data is available are equivalent, regardless of whether or not the photon emissions are modeled as Poissonian.


The following description is directed at a non-limiting preferred embodiment wherein the relationship between the measurements y and intensities I are modeled as a set of linear equations, y=ΦI, for the purposes of reconstruction. This model is suitable for least squares optimization when photon emissions have a Gaussian distribution, and thus may be reasonably applied when the measured counts (i.e. y) are relatively high. However the embodiments are extendable to modeling the emissions as Poissonian. In a preferred embodiment, when a given voxel or voxels are not within Gaussian range, the y values of several low-intensity views are combined and set equal to the sum of all the related linear equations, in order to bring the total count into Gaussian range.


Assuming intensities with a Gaussian distribution, then by applying SVD to the present inverse problem, I is obtained as:

I=VD−1Uty  (17)

As discussed above, the condition number of D may serve as an indicator of the stability of the reconstruction, particularly during an iterative reconstruction process. Furthermore, known techniques such truncated SVD may be used for optimization and reconstruction purposes.


The present embodiments are of a method, apparatus, and system which analyze the reliability of the reconstruction possible by a given imaging constellation, and define modifications in order to improve the quality of reconstruction. The modifications may be to the reconstruction and/or data collection aspects of the imaging process. As described below, the modifications may use active view selection during imaging and/or to guide non-uniform scanning of the imaged volume.


Reference is now made to FIG. 2a, which is a simplified flowchart of a method for stabilizing the reconstruction of an imaged volume, according to a preferred embodiment of the present invention. In step 115, an analysis of the reliability of reconstruction of a radioactive-emission density distribution of the volume is performed. The analysis is made to determine the reliability of the data which may be collected by radioactive-emission measurements taken over a specified set of views. In step 120, modifications to improve reliability of reconstruction are defined in accordance with the analysis. The modifications may involve a change or adaptation of the reconstruction process and/or a change or adaptation of the data collection process. Both the analysis and modification aspects of the present embodiment are described in detail below.


The present embodiment tailors the reconstruction and/or data collections processes in order to obtain a more accurate reconstruction of the imaged volume. Low reliability data may cause large errors in the reconstruction process. As a result the reconstruction process may be noisy and unstable, and fail to converge properly. The reconstructed image may contain artifacts, unsupported frequencies and other errors.


The present embodiments are applicable to all stages of the imaging and/or reconstruction processes. The analysis may be performed prior to data collection, in order to define a scanning procedure which provides high-reliability data. During data collection, the scan pattern may be adapted or views may be added (e.g. active vision/adaptive scanning), utilizing the new information provided by the collected data. After data collection, modifications may be made to the reconstruction process in order to counteract the effects of unreliable data and non-uniformities, and to improve reconstruction accuracy and stability.


In a first preferred embodiment, the analysis of the reliability of reconstruction is based upon on analysis of the singular values of the probability matrix Φ. A non-limiting preferred embodiment of such an analysis is presented below.


Reference is now made to FIG. 2b, which is a simplified flowchart of a method for analyzing a detection probability matrix, according to a preferred embodiment of the present invention. In step 125, a detection probability matrix is provided. As described above, the detection probability matrix gives the detection probability for each of the views over the voxels of the volume. More explicitly, in the preferred embodiment each view has an associated row in the detection probability matrix. The values in a given row are derived from the probability of detection of a photon emerging from each of the voxels of the volume, by a detector at the associated view. In step 130, the singular values of the detection probability matrix are calculated. In the preferred embodiment, the singular values are calculated by employing SVD on the detection probability matrix. Performing SVD on the detection probability matrix enables easy implementation of stabilization methods such as truncated SVD for stabilizing the reconstruction process. In step 135, those singular values which may have a destabilizing effect on the reconstruction process are identified as destabilizing singular values.


In the preferred embodiment, a measure of the reconstruction reliability (denoted herein the reliability measure) is determined. The reliability measure is used for determining if any modification of the imaging and/or reconstruction processes is required, or if a stable reconstruction is possible without such modification.


Preferably, the condition number of the probability matrix is utilized as a reliability measure. If the condition number is satisfactory the reconstruction may be considered stable, with no need for modification. Preferably, the identifying of destabilizing singular values is performed only if the condition number is below a specified magnitude.


In an alternate preferred embodiment, the information-theoretic Fisher information is calculated as a reliability measure from intensity distributions provided of the imaged volume. An intensity distribution may be obtained either as an emittance model created prior to the data measurement process, or as an intermediate result of the reconstruction process. The Fisher information is described in detail below, in the context of the reliability criterion for active vision.


In the preferred embodiment, the identification of destabilizing singular values is performed as follows. The ratio of the largest singular value to each of the remaining singular values is calculated. A singular value for which the ratio is above a specified threshold is considered destabilizing. One may consider that if the effects of the destabilizing singular values upon the reconstruction are eliminated, the condition number of the detection probability matrix is effectively below the threshold.


Preferably, unreliable voxels or combinations of voxels associated with a given destabilizing singular value are identified. The intensity levels of such unreliable voxels or combination of voxels are unsupported in the collected emission data, due to the destabilizing effects of the associated singular value. Therefore the accurate reconstruction of such voxels and/or combinations is unlikely, and may destabilize the reconstruction process. The destabilization of the reconstruction by unsupported voxels may lead to artifacts and unsupported frequencies in the reconstructed intensities. These effects are not necessarily limited to the unsupported regions, but may extend into other portions of the volume and create artifacts there as well.


The voxels and linear combinations of voxels associated with a singular may be determined via SVD decomposition of the intensity distribution matrix Φ. The matrix V obtained by SVD decomposition of Φ indicates voxels or linear combinations of voxels whose intensity after reconstruction is most highly affected by each of the singular values in D. Each row of V indicates a weighting associated with a given singular value for each voxel in I. Consider the following example of a volume with four voxels i1-i4:









y
=



UDV
T


I

=




U


[




e
1


























e
2


























e
3


























e
4




]




[




v
11




v
12




v
13




v
14










































v
41




v
42




v
43




v
44




]




[




i





1






i





2






i





3






i





4




]


=

U








[




















e
4



(



v
41


i





1

+


v
42


i





4

+


v
43


i





3

+


v
44


i





4


)





]







(
18
)







In a first example, assume that e4 is a destabilizing singular value, and that the fourth row of V, v4, equals [0.001 0.01 0.97 0.015]. Examining v4 it is seen that the coefficient for the single voxel i3 is significantly larger the coefficients of the other voxels. It is therefore quite apparent that the reliability of data regarding voxel i3 is particularly low, and it is unlikely that i3 can be accurately reconstructed. In this example, the preferred embodiment is to perform smoothing of i3 relative to the neighboring voxels. Smoothing may be a practical approach where an inaccuracy occurs in spatially localized regions of the volume. An exemplary criterion to determine which voxels should be smoothed, is to select voxels whose corresponding coefficient of vi is above a certain threshold.


In a second example, assume that e3 is a destabilizing singular value, and that the third row of V, v3, equals [0.5 −0.5 0.5 −0.5]. In the present example, the unreliability of the data relates to a combination of voxels which may be distributed over the volume, rather than to an independent, spatially located voxel or group of voxels. In the present example, the reliability of data regarding the linear combination 0.5i1−0.5i2+0.5i3−0.5i4, as reflected by measurements y, is low, which may lead to instability and inaccuracies during reconstruction. Since there are no dominant voxels it may be less effective to perform smoothing. Instead, the preferred embodiment is to constrain intensities of the combination of voxel intensities between iterations.


A proposed approach for applying such a constraint is as follows. Given a reconstructed intensity vector I′, modify the intensities prior to performing the next iteration as follows:

I′=>I′−vkα(vkTI)  (18)

where vk is the row of V associated with a destabilizing singular number, ek, and α is a proportionality factor which may be related to the magnitude of ek. The effect of the voxel combination on the following iteration is therefore reduced as desired, where the reduction may be correlated with the dominance of the associated singular value.


Preferably, constraints are added during the reconstruction process to reduce or eliminate the effects of the unreliable data associated with destabilizing singular values.


In the preferred embodiment, the modifications include defining constraints on the reconstruction process. The reconstruction may be considered unreliable if one or more destabilizing singular values have been found, and the modifications may be made in order to reduce the effect of the destabilizing singular values upon the reconstruction. Examples of such constraints may include:

    • 1) Smoothing or piecewise smoothing for all or portions of the volume
    • 2) Uniting voxels
    • 3) Subdividing voxels
    • 4) Adjusting the levels of one or more voxels after a reconstruction iteration
    • 5) Constraining the input data
    • 6) Performing reconstruction with differing resolutions for different portions of the volume


      The above exemplary constraints are now discussed in turn.


Smoothing a reconstructed image is performed by calculating the intensity level of a voxel based on the intensities of surrounding voxels, in order to control the magnitude and rate of the fluctuations in intensities between voxels. For example, the intensity of a given voxel may be corrected to better reflect its value as an average of the surrounding voxels. Smoothing may be performed over the entire volume, or only on portions of the volume deemed unreliable.


When a non-uniform scan is performed, it is possible that not all directions surrounding a voxel will have equal support. It may therefore be easier to distinguish voxel from some neighbors than from others. In the preferred embodiment, smoothing is performed directionally. In directional smoothing, the smoothing is weighted more strongly in certain directions than in others.


In a first example, directional smoothing is used to overcome the limitations of unreliable data by smoothing the voxel in the weaker directions. In a second example, a higher resolution is available for some portions of the volume than is available in other portions. If it is important to obtain a reconstruction having a uniform resolution, directional smoothing may be applied reduce the resolution in highly reliable portions of the volume to the resolution obtainable elsewhere.


A directional smoothing policy may be based on the Fisher information measure, which provides directionality information. The Fisher information matrix indicates directionality and cross-relationships between voxels, but may be difficult to calculate. In the preferred embodiment, a scalar Fisher information is calculated for groups of views, where each group is localized in a given direction. The scalar Fisher information measures may then be analyzed together, in order to determine the relationships between the different directions.


Low reliability voxels may be united or combined with surrounding voxels in order to increase their reliability. Voxels may combined by assigning the intensity level of a neighboring high-reliability voxel to the low-reliability voxel. Voxel merging may be carried out more extensively, by repeatedly merging the lowest reliability voxel with one of the neighboring voxels to form an aggregate of voxels, until no voxel is left with a reliability below a threshold. It will be appreciated that the aggregate voxels have a lower resolution locally, but all other regions with good coverage remain of high resolution, according to their coverage.


An opposite approach for high-reliability portions of the volume is to subdivide high-reliability voxels into smaller voxels in order to improve resolution while still achieving a required reliability. For example, if there are 1000 different views independently covering 1 cubic cm and almost not affected by the surrounding voxels, then theoretically that volume can be divided up to about 1000 voxels, if the views create a linearly independent set of equations with a good condition number.


A further possible modification of the reconstruction process may be to adjust the effect of certain unreliable voxels or linear combinations of voxels after each of one or more iterations of the reconstruction process. For example, the magnitude of unreliable voxels or combinations of voxels associated with a destabilizing singular value may be reduced towards zero or towards the value of neighboring voxels by a specified factor, possibly in proportion to their respective weightings in the associated row of matrix V and/or to the associated singular value.


An additional possible modification is to constrain the value of the input. For example, a realistic range of counts may be defined overall or per view based on view parameters, typical spatial structures and the like. Examples include adding a prior constraint such as a gamma distribution of intensities, a constant or piecewise-constant progression of intensities, a linear or piecewise-linear progression of intensities, smoothing or piecewise smoothing constraints, an intensity distribution based on the shape of or magnitude of the object being scanned (possibly determined during a pre-scan), or a maximal or known range of expected intensities,


The reconstruction process may be performed in a manner that obtains varying resolutions over the volume. In one approach, the entire volume is defined as a single huge voxel, and split over the reconstruction iterations to form smaller voxels. The process is performed repeatedly over all or some of the voxels of the volume, as long as the result of the split maintains a reconstruction reliability high enough for stable results.


In a further preferred embodiment, the information-theoretic Fisher information is used as a reliability measure. Modifications to the reconstruction and/or data collection processes may be performed when the Fisher information is deemed to be outside a specified range. The calculation of the Fisher information is described in detail below, in the context of the reliability criterion for active vision. The Fisher information may be calculated from the results of a previous reconstruction iteration, or, initially, from an emittance model provided of the imaged volume.


The following addresses a preferred embodiment in which the modification is implemented by defining views for imaging the volume. The defined views serve to guide the data collection process in order to obtain measurements which enable performing a stable and accurate reconstruction of the intensity distribution of the volume. In this manner, detecting resources may be invested effectively in order to improve reconstruction reliability. Scanning resources include, for example, detector dwell time, number of detectors, angular and translational increments, and the like—features that increase the amount of data collection.


Such modification may be performed when a region is interesting but is still left with coarse resolution, so as to allocate more scanning resources, such as dwell time, number of detectors, angular and translational increments, to cover that region and to form more independent views such as to increase the reliability of the reconstruction of that area. Preferably, the views are defined so as to obtained a desired resolution over the region of interest.


The defined views may yield a non-uniform scanning procedure. A scanning density may be specified by the angular and translational increment size, or steps, the smaller the increment, the higher the density. Additionally, the scanning density may be defined by the acquisition time at each position—the longer the time, the higher the density.


Non-uniform scanning may be defined by specifying varying scan densities for imaging the volume. For example, the modification may entail adjusting a local scanning density to scan a region of interest with high density, and to scan other regions with low density.


The non-uniform scans may relate to non-uniform angular steps of the detector along a sweep, non-uniform detector translational steps, or different steps by different detectors. Some detectors may employ dense steps and others may employ sparse steps, for example, based on active vision as taught hereinbelow.


The scan density may be adapted to the distance to the object of interest. Since resolution decreases with distance, the higher density may compensate for increased distance.


Additionally or alternatively, the angular steps may increase in density when scanning a region of interest, and decrease in density, when scanning other portions of the volume.


Furthermore, more than one region of interest may be scanned with dense steps, simultaneously. The two regions of interest may be, for example, a tissue region and a blood pool region. This has applicability, for example, to dynamic studies of blood perfusion, by providing even scanning resources both to the blood and to the tissue.


Additionally, convex scans may be employed.


Variable scans, where a same region is scanned first with a first density and then with another density, may be employed. Alternatively, a same region may be scanned by a first group of detectors with a first density and then by a second group of detectors with another density, concurrently, or at different times. Thus the same region is scanned with at some density by a given detector and at a different density by another detector.


In a further preferred embodiment, view definition is performed dynamically during radioactive-emission measurements of the volume. The definition of further views for measurement during the data collection process is denoted active vision herein and is described in detail below. Active vision may be performed independently or in conjunction with the present preferred embodiment of stabilizing the reconstruction of an imaged volume.


In the preferred embodiment, the method includes the further step of iteratively reconstructing the radioactive-emission density distribution of the volume, preferably by EM estimation. Reconstruction reliability is preferably evaluated after every iteration.


Reference is now made to FIG. 2c, which is a simplified block diagram of a reconstruction stabilizer, according to a preferred embodiment of the present invention. Reconstruction stabilizer 140 includes reliability analyzer 145 which analyses the reconstruction reliability of a radioactive-emission density distribution of a volume. The reconstruction reliability may be determined from measures such as the condition number of the probability density matrix or the Fisher information, as described above. Modifier 150 defines modifications to improve the reliability of the reconstruction, based on the analysis performed by reliability analyzer 145. Modifier 150 preferably applies at least one constraint to the reconstruction process and/or data collection. The constraints may include applying smoothing, uniting voxels, adjusting the levels of voxels, defining new views or non-uniform scanning, and other constraints described above. The views may be defined before or during emission data collection.


In the preferred embodiment, reconstruction stabilizer 140 is included in a system for generating a three-dimensional image of volume. In the preferred embodiment, the system further includes radiological imaging camera 155 and reconstructor 160. Camera 155 includes detectors capable of independent movement during data acquisition, and which are capable of detecting radiation emitted from the volume thereby to provide radiation data. In accordance with embodiments of the present invention, each block of the imaging camera construction or each detecting unit, where single-pixel detecting units are used, may be provided with at least one, and preferably, two, three, or as many as six degrees of motion such as, for example, rotational motion around the x, y, and z, axis, oscillatory motion about these axes, or translational motion along these axes. Camera 155 thus enables flexible view definition by modifier 150, increasing the likelihood that a reliable reconstruction will be achieved. Reconstructor 160 analyzes the radiation data provided by camera 155 so as to reconstruct a three-dimensional image of the volume.


Reconstructor stabilizer 140 guides camera 155 and reconstructor 160, by providing one or both with the modifications determined to improve reconstruction reliability.


In practice, it is generally necessary to reconstruct volumes at a high resolution (i.e. a large number of voxels) from many measurements. Implementing the techniques described above then requires the manipulation of extremely large matrices. The large size of these matrices poses difficulties both for performing the calculations and for memory management. These difficulties are particularly problematic when active view selection is performed.


In the preferred embodiment, calculations are performed in a localized manner. That is, calculations are performed for different sections of the volume in turn, progressing through the volume until it is covered in its entirety. When calculations are being made of a particular section of the volume, the variables relating to the other sections are “frozen”. Thus first a selected section of the volume may be stabilized, so that the reconstruction of following sections is based on data with improved reliability. Localized calculations, such as SVD decomposition, are performed on smaller sub-matrices rather than on a single large matrix, alleviating the difficulties of manipulating very large matrices.


The division of the volume into sections may be devised in any practical way, for example as sequential slices or by another spatial connection. The division may be based on the scan pattern, so that well-supported sections may be stabilized first, followed by less supported section. The division may also or alternately be based on knowledge of the body being imaged, such as its shape or composition. When imaging a body structure, the sections may be based on the known axes or structures of the organ


The following embodiments are of a method for determining further views for the imaging of a body structure, denoted active vision herein. Active vision addresses the problem of ensuring that the quality of data gathered during the measurement process is adequate to provide a high quality image. The collected data and/or the image reconstructed from the collected data is analyzed the while the measurement process is taking place. Based on the analysis, further views are defined. Since each view is associated with known values of the viewing parameter(s), selecting a view effectively specifies known viewing parameter values. The defined further views thus define a set of viewing parameter values, which are used during the current measurement process in order to collect data which yields a high-quality reconstruction of the body structure.


Active vision may be performed independently, as described in the embodiments presented below. Additionally or alternatively, active vision may be performed in conjunction with the above-described method for stabilizing reconstruction, by modifying data collection by dynamically providing additional views to improve reconstruction reliability (see FIG. 17).


The following embodiments are not confined to a specific reconstruction algorithm. Further views are preferably defined based on one or more of the following:

    • 1) Detector photon count
    • 2) Geometric properties of the reconstructed body structure
    • 3) Information theoretic measures that quantify the quality of the data fed to the reconstruction algorithm


      Each of these criteria is discussed in detail below.


Reference is now made to FIG. 2d, which is a simplified flowchart of a method of performing radioactive-emission measurements of a body structure, according to a preferred embodiment of the present invention. In step 200, radioactive-emission measurements of the body structure are performed at predetermined views, preferably in vivo. Preferably the measurements are performed for diagnostic purposes. These predetermined views are selected prior to the measurement process, based on a model of the body structure being imaged. In the model more and less informative viewing directions have been identified. The predetermined views of step 200 preferably include those views expected to be informative, based on an analysis of the model.


Preferably the body structure is all or a portion of: a prostate, a heart, a brain, a breast, a uterus, an ovary, a liver, a kidney, a stomach, a colon, a small intestine, an oral cavity, a throat, a gland, a lymph node, the skin, another body organ, a limb, a bone, another part of the body, and a whole body.


In step 210 the radioactive-emission measurements are analyzed. Preferably the analysis includes one or more of:

    • 1) Analyzing detector photon count(s)
    • 2) Analyzing detector photon count rate(s)
    • 3) Identifying detector saturation
    • 4) Reconstructing a body structure image from emission measurements
    • 5) Identifying geometric properties of the reconstructed image
    • 6) Applying information-theoretic measures to the reconstructed image


In step 220, further views for measurements are dynamically defined, based on the analysis performed in step 210. Preferably, each of the views is associated with viewing parameters selected from the group consisting of: detector unit location, detector unit orientation, collection angle, and measurement duration. Defining a view consists of providing a value for each of the parameters associated with the given view. The analysis (step 210) and/or dynamic view definition (step 220) may take into account external parameters including: measurement duration, time elapsed from the administration of the pharmaceutical to the measurement, radiopharmaceutical half life, radioactive emission type, and radioactive emission energy.


Each of these analysis techniques, and their application to view definition, is now discussed in turn. While each of the analysis/view determination techniques is discussed as a separate embodiment, multiple techniques may be used together to obtain the desired image quality.


In a first preferred embodiment, a photon count analysis ensures that the photon count at a given view yields an acceptable measurement error. As discussed above, the radiative emissions of the body structure being imaged is a Poisson process. In a Poisson process the Poisson noise grows inversely to the square root of the number of photons detected. In other words, if N photons are collected from a given view, the resulting signal to noise ratio (SNR) equals:

SNR=n/√{square root over (N)}=√{square root over (N)}  (19)


The unprocessed detector photon count at a given view thus provides significant information regarding the quality of the information obtained at a given view. If the photon count is too low, it may be desired to continue to collect photons at the current location/orientation in order to obtain a satisfactory SNR. Alternatively, it may be determined that enough photons have already been collected, and to terminate the current view and move on to the next view.


The analysis is preferably performed by defining a global or local required measurement error, and comparing the square root of the obtained photon count to the required measurement error. Photon count analysis can be applied to the current and/or previous views. When a photon count of a current view is found to be too low, the duration of the current view is preferably extended in order to obtain the required error value. When a photon count of a past view is found to be too low, an emission measurement at substantially the same location and orientation but having a longer duration than previously is preferably performed. Alternately or additionally, the collection angle at the given location/orientation is preferably increased.


In an additional preferred embodiment, a detector photon count is analyzed to identify detector saturation at a given view. Preferably, when a detector is determined to have saturated, a new view or views are selected to reinforce those views that have saturated. In an alternate preferred embodiment, further views are defined to avoid highly-radiating portions of the body structure.


In a second preferred embodiment, a photon collection rate at a given view is analyzed to determine if it is within a specified range. In the preferred embodiment, the photon count rate is used to identify regions of high or low interest. In prostate imaging, for example, a region of high interest may be identified by a high photon rate, indicative of a tumor. In a second example, a region of high interest may be identified in heart imaging by a low photon rate, indicative of non-functional tissues. After one or more areas of high and/or low interest are found, further views are preferably defined by selecting views to concentrate on regions of high interest and/or to avoid regions of low interest. It is thus possible to zoom in on a suspected pathology without repeating the emission measurement process.


In a further preferred embodiment, the analyzing of step 210 includes reconstructing a radioactive-emission density distribution of the body structure. Reconstruction may be performed according to any applicable technique known in the art. The reconstruction is then used as the basis for further analysis.


Reconstruction based on the data collected from the predetermined views provides information regarding the quality of information obtained from the preceding measurements, and which further views are likely to be most informative. Selecting new views based on reconstruction is intended to bring us into viewing from the more informative views or combinations of views.


Reference is now made to FIGS. 3 and 4a-4b, which pictorially illustrate how different views provide differing types and quality of information. FIG. 3 shows an object 300 shaped as a cylinder with a front protrusion, and having a high-emittance portion (hotspot) 310. Four views of object 300 are shown, which can be seen to provide different levels of information. Front views, such as VI, provide little information regarding the shape of object 300 and have relatively little attenuation between the detector and hotspot 310. Side views, such as V2, provide edge information regarding the object shape or profile, and when correlated with front views help locate hotspot 310 spatially within object 300. Top views, such as V3, provide information regarding the cylinder edge region 320. Finally, rear views, such as V4, are uninformative about the shape of object 300 and have high attenuation relative to hot region 310.



FIGS. 4
a and 4b demonstrate how the proper selection of views may improve the quality of information obtained for the body structure, for example in distinguishing between two regions of interest within a given volume.



FIG. 4
a illustrates an object 400 having two high-emission regions of interest (ROI), 410 and 420. For clarity the views VA to VF are shown as lines in FIG. 4a, however in practice they will each have a finite collection angle δ. The position of ROIs 410 and 420 are assumed to have been estimated based on a model of object 400 and/or a previously performed prescan. The goal of the present invention is to select an additional new view or views which increase the information we have regarding the separation of ROIs 410 and 420 within object 400.


In simple terms, consider the object as having three regions: ROI 410 with intensity I1, ROI 420 with intensity I2, and a low-emission region 430 between the two ROIs with intensity I3. The detected intensity at a given detector is proportional to








I
n


r
ni
2


,





where In is the emission intensity of region n and ri is the distance of region n from detector Vi.



FIG. 4
b illustrates the added information provided by each of the shown views, VA to VF. Views VB and VC collect emissions from all three regions, and are therefore least informative. Views VD and VE collect emissions from only low emittance region 430, and therefore provide most information regarding the location of each ROI within the volume and the separation between ROIs 410 and 420. Views VA and VF pass only through a single ROI, and therefore provide an intermediate level of information. It is a goal of the present invention to determine, while the emission measurements of the body structure are taking place, that views in the vicinity of VD and VE are highly informative, and to add these further views to the measurement process.


A body structure reconstruction can be utilized in several ways to define further views. A first way is to identify interesting portions of the contour and structure of the reconstruction. For example, it is seen in FIG. 3 that top views are informative about edge region 320. Further top view measurements will therefore be informative re edge region 320, and may enable defining the edge more accurately.


In a preferred embodiment, the reconstruction is analyzed to identify textural edges within the reconstruction, and view definition preferably includes selecting views at an angle to the textural edges. In the preferred embodiment, the angle is a substantially sharp angle in order to provide information regarding the edge.


In another preferred embodiment, the reconstruction is analyzed to identify volumetric boundaries within the reconstruction, and view definition preferably includes selecting views at an angle to the volumetric boundaries. It is expected that the defined views will provide information regarding the boundary and differences in surrounding tissues on either side of the boundary. In the preferred embodiment, the angle is a substantially sharp angle.


Another way to utilize the reconstruction to define further views is to identify suspected organ targets within the reconstruction, and to select further view(s) in close proximity to the suspected organ targets. A suspected organ target is typically detected by identifying portions of the reconstruction whose emission intensity distribution and spatial characteristics are typical of a suspect region.


In a first preferred embodiment, a suspected organ target is defined as a high-emittance portion of the reconstruction. In prostate imaging, for example, a suspected organ target may be a high-emittance portion of the structure, indicating a tumor. A high-emittance portion is characterized by an intensity that is greater than the background intensity by a factor of at least (1+α), where α is a parameter specified by the user. In practice, a hotspot is usually detectable only if the radiation levels within the hotspot are higher than the background level by a factor of 1.5-2. α is therefore typically defined between 0.5-1. However, the detectability of a hotspot rises as the radioactive intensity of the body rises, raising the photon count. Thus, a lower value of α may be used when the measured body structure is in a state of high-intensity emittance. For example, a body structure may be characterized by relatively high emittance immediately following the administration of the radiopharmaceutical, and be characterized by lower emittance at a later time.


In a second preferred embodiment, a suspected organ target is defined as a low-emittance portion of the reconstruction. In heart imaging, for example, a suspected organ target is defined as a low-emittance portion of the reconstruction, indicating non-functional tissues. A low-emittance portion is characterized by an intensity that is lower than the background intensity by a factor of at least (1+β), where β is a parameter specified by the user.


In the preferred embodiment the further views are used immediately for radioactive-emission measurements. The results of the new measurements are then used in another analysis to define new further views for additional measurements. The radioactive-emission measurements may then be said to be performed iteratively.


Reference is now made to FIG. 5a, which is a simplified flowchart of a iterative method of performing radioactive-emission measurements of a body structure, according to a first preferred embodiment of the present invention. In step 500, radioactive-emission measurements of the body structure are performed at predetermined views. In step 510, an analysis is performed of the previously performed emission measurements. In step 520 a decision is made whether to continue with further measurements. If yes, in step 530 further views are defined based on the analysis. Subsequent iterations continue until the decision to end the emission measurement process. After the first iteration, the analysis performed at a given stage may include consideration of all or on part of the measurements performed during one or more previous iterations, in addition to the new measurements.


Reference is now made to FIG. 5b, which is a simplified flowchart of a iterative method of performing radioactive-emission measurements of a body structure, according to a second preferred embodiment of the present invention. In the present preferred embodiment, a reconstruction of the body structure is formed in step 505. The analysis step 510 is then performed utilizing data provided by the reconstruction(s).


Referring again to FIG. 2, preferably, analysis step 210 includes determining an accuracy of the reconstruction. Accuracy is preferably determined by analyzing the variance of the reconstructions formed over multiple iterations. Preferably, further views are defined in step 220 to concentrate on the region for which higher accuracy is required. Regions of the reconstruction having low variance provide a high degree of confidence regarding the accuracy of the reconstruction in the given region (where a portion may include the entirety of the body structure being imaged). Further views may be added to the current measurements until the variance is reduced to a required level.


Preferably, analysis step 210 includes determining a resolution of the reconstruction. Resolution is preferably determined by analyzing the full width at half maximum (FWHM) of peak values of the reconstruction. The FWHM is given by the distance between points at which the reconstructions reaches half of a peak value. Preferably, further views are defined in step 220 to concentrate on the region for which higher resolution is required.


An additional way to define future views using the reconstruction is on an information-theoretic basis. A quality function expressing an information theoretic measure is defined. The quality function rates the information that is obtainable from the body structure when one or more permissible views are added to current measurement process. Several examples of quality functions based on information-theoretic measures are discussed in detail below. The quality function is used to rate potential further views. The measurement process may then continue at those further views whose addition to the previous views yields a high rating.


Reference is now made to FIG. 6a, which is a simplified flowchart of a method for dynamically defining further views, according to a first preferred embodiment of the present invention. In step 610 a quality function is provided. The quality function expresses an information-theoretic measure which rates the quality of information obtainable from potential further views. In step 620 a set of further views is selected to maximize the quality function. Preferably the selected further views fulfill certain constraints; for example the further views may be selected from a predefined set or may be located in the vicinity of a region of interest within the body structure.


In the abovedescribed reconstruction-based analyses, the quality function is evaluated independently for a single reconstruction of the emission intensity of the body structure. However, quality functions may be defined which calculate the score for a given set in relation to one or more reconstructions and/or emittance models. An emittance model is a representation of a specific radioactive-emission intensity distribution within the volume U, so as to model organ targets, such as hot regions, of a radioactive emission intensity, higher than the background level, regions of low-level radioactive emission intensity, which is nonetheless above the background level, and cold regions, of a radioactive emission intensity, lower than the background level. Given an object or class of objects, emittance models may be devised to reflect expected or typical emission patterns for the given object.


Developing an emittance model for a particular body structure involves analyzing known information about the body structure to determine expected emission patterns of the body structure. In order to develop a model of a particular body structure, for example a prostate, many factors may be considered. Physical aspects, such as the size and shape of the prostate and the position of the prostate within the torso may be considered, as well as medical knowledge regarding typical emissions from healthy and diseased prostates. Additional information may concern variations between individuals, such as age, weight, percentage of body fat, and the like.


For simplicity, the following discussion describes the evaluation of information-theoretic quality functions based on emittance models only. It is to be understood that at least one of the emittance models is a reconstruction of the body structure based on past measurements. Any remaining emittance models are provided externally, and may be based on general medical knowledge or on information gathered during a previous round of emission measurements of the body structure.


Reference is now made to FIG. 6b, which is a simplified flowchart of a method for dynamically defining further views, according to a second preferred embodiment of the present invention. The current method differs from the method of FIG. 6a by the addition of steps 605-606. In step 605 a set of one or more emittance models is provided (where the set includes one or more reconstructions of the body structure). An emittance model specifies the radiative intensity of each voxel in the body structure. As discussed above, some of the viewing parameters affect the radiative intensity of the voxels in the volume, for example the type of radiopharmaceutical and the time since administration of the radiopharmaceutical. Therefore, the emittance models provided in step 605 preferably correspond to the relevant viewing parameters. In step 606 a collection of possible further views of the body structure is provided. The collection of views includes possible further views for future measurements, preferably based on anatomical and other constraints. Furthermore, the quality function provided in step 610 may utilize multiple emission models.


Reference is now made to FIGS. 7a and 7b which present the principles of modeling of a body structure, in accordance with the present invention. FIG. 7a schematically illustrates a body section 720, having a region of interest (ROI) 730. The region of interest 730 is associated with a body structure 740, with a specific radioactive-emission-density distribution, possibly suggestive of a pathological feature 750, termed herein an organ target 750. Additionally, there may be certain physical viewing constraints, associated with the region of interest 730. Model 760, as illustrated in FIG. 7b, represents body section 720 as a collection of voxels, having a specified radioactive-emission intensity distribution.


It will be appreciated that the model 760 of the region of interest 730 may be based on general medical information of the body structure 740 and common pathological features associated with it. Additionally, the model may be based on information related to a specific patient, such as age, sex, weight, and body type. Furthermore, a structural image, such as by ultrasound or MRI, may be used for providing information about the size and location of the body structure 740 in relation to the body section 720, for generating the model 760.


Reference is now made to FIG. 8, which shows four models of a prostate emittance model. Three views are given for each of the emittance models shown. Each of the views is a maximum intensity projection (MIP) along a different one of the three axes. For example, in an X-Y projection the intensity of a given point is taken as the maximum intensity for that point along a line parallel to the z-axis. Thus the volume in effect becomes “transparent”, and the maximum intensity regions are shown clearly. Variations in emission levels are indicated by differences in shading. The emittance model is shaped as an ellipsoid, the typical shape of a prostate. Since a diseased prostate is generally characterized by one or more hot regions, each of the emittance models shown has a number of high-emittance portions (two hot regions in Models A-C and three hot regions in Model D). Each high-emittance area is an ellipsoid, with a size of approximately 1 cubic centimeter. The intensity of the modeled organ targets varies randomly from two to eight times brighter than the surrounding volume.


In the preferred embodiment, one or more of the emittance models contains at least one high-emittance portion (i.e. hot region). A prostate containing a tumor, for example, may be modeled as an ellipsoid volume with one or more high-emittance portions.


In the preferred embodiment, one or more of the emittance models contains at least one low-emittance portion. A diseased heart may therefore be modeled as a heart-shaped volume with low-emittance portions.


Note that an emittance model need not contain high- or low-emittance portions, but may have a uniform intensity or a slowly varying intensity.


Following are two embodiments of quality functions, both of which include considerations of the emission distribution of the volume, in light of the emittance models provided (step 605 of FIG. 6b). The information quality given by one or more further views is evaluated in relation to multiple reconstructions and/or emittance models of the body structure. In a first preferred embodiment the quality function is based on an information theoretic measure of separability. In a second preferred embodiment, the quality function is based on an information-theoretic Fisher information measure, for ensuring reliable reconstruction of a radioactive-emission density distribution.


In a first preferred embodiment, the quality function implements a separability criterion. Separability is a measure of the extent to which the measurements that are obtained from each pair of models can be distinguished from one another.


The concept of separability is illustrated in FIGS. 9a-9c, each of which shows an emittance model of a given volume 900. Each of the emittance models contains two high-emittance portions (hot regions), 910 and 920, which are located in respectively different locations in each of the models. It can be seen that the hot regions in emittance models 9b and 9c are in similar portions of the volume, as opposed to the hot regions in model 9a. It may therefore be difficult to distinguish between reconstructions of emittance models 9b and 9c. The separability criterion ensures that the selected set includes views which provide reconstructions which distinguish between emittance models 9b and 9c.


Letting I be the emittance model set, a measure for the dissimilarity of any two given densities in I is defined. Since most state-of-the-art estimating algorithm are aimed at finding ML estimators, in the current example the quality function is based on the likelihood function. The likelihood of an estimator of I, given a set of Poissonian measurements y is:













(
I
)


=



t







{


-



u





ϕ
t



(
u
)




I


(
u
)





+


y
t


ln




u





ϕ
t



(
u
)




I


(
u
)





-

ln


(


y
t

!

)



}






(
19
)







For separability, it is desired that this measure be different for each IεI. Since the measure is a random variable that depends on the actual measurements, all possible pairings of emittance models should be examined to ensure that the resulting distributions are separable. A quality function that captures this separability is given by the square of the difference between the means of the distributions normalized by the sum of their variances:











SEPARABILITY
Φ



(


I
1

,

I
2


)


=






E









(

I
1

)



-

E









(

I
2

)






2



Var









(

I
1

)



+

Var









(

I
2

)









(
20
)







The expectations and variances in Equation 20 are taken over random measurements y, sampled from the true intensity I1 (note that the measure is not symmetric).


Since the true (unknown) intensity can be any IεI, a projection set Φ* that maximizes the worst-case separability is desired. That is:

Φ*=arg maxΦminI1,I2εISEPARABILITYΦ(I1,I2)  (21)


Scoring for separability is based on the minimum separability obtained with a given set of views for all of the possible pairings of emittance models from the set of emittance models, thereby enabling defining a desired resolution in more than one direction, or in more than one portion of the volume. All of the emittance models are modeled on a substantially identical volume. The emittance models preferably differ from one another in the modeled organ targets, where the modeled organ targets are separated by a difference of at least the required resolution (where the displacement which produces the required resolution is denoted delta herein). Substantially identical sets of views are formed from the collection of views, and each of the formed sets is scored with respect to each of the pairs. One of the sets of views is selected, based on the minimum or average score for the plurality of pairs.


For example, assume the set of emittance models contains the three models 9a-9c. A separability score is calculated for a given formed set of views by applying Equation 20 to all three pairs 9a/9b, 9a/9c, and 9b/9c. The lowest of the three calculated values is taken as the separability score for the formed set. Once a separability score has been calculated in such manner for each of the formed sets of views, the view set having the highest separability is selected.


The separability criterion may be used to ensure that a required resolution is obtained in all or a portion of the body. In a preferred embodiment, view set selection for separability is performed utilizing a set of emittance models consisting of one pair of emittance models having substantially identical volumes but with different modeled organ targets. The modeled organ targets are separated by a delta in a given direction so as to define a required resolution in that direction and portion of the volume U. Substantially identical sets of views are formed from the collection of views, and scored with respect to the pair of emittance models, using a quality function based on the separability criterion, and one of the sets of views is selected based on the separability scores. The selected set is thus the set which provides the optimum resolution in the given direction and in the vicinity of the modeled organ targets.


A similar approach may be used to ensure resolution in more than one direction and/or portion of the volume U. Consider for example, a pair of models of substantially identical volumes, as follows: The model of FIG. 9d, which schematically illustrates the volume U, having the modeled organ target HS, whose center is at a location (x;y;z)HS, and the model of FIG. 9e, which schematically illustrates the volume U, having the modeled organ target HS′, whose center is at a location (x;y;z)HS′. FIG. 9e also shows modeled organ target HS of FIG. 9d superimposed over the present emittance model to illustrate that HS′ is somewhat displaced from HS, along the x-axis, and the displacement, is denoted as delta1 in the present example. The displacement delta may be measured, for example, in mm.


An optimal set of views, from the standpoint of separability, is that which will best distinguish between the model of FIG. 9d and that of FIG. 9e. Thus, a score, in terms of separability is given for the pair of models, and relates to a resolution as defined by the difference between the models of the pair. In the present example, the difference is delta1 along the x-axis, around the locations of HS and HS′, so the score given by the information theoretic measure of separability will relate specifically to a resolution as defined by delta1 along the x-axis, around the locations of HS and HS′. Other portions of the volume U and other directions may have different resolutions.


For example, consider the model of FIG. 9f, which schematically illustrates the volume U, having the modeled organ target HS″, whose center is at a location (x;y;z)HS″, and shows HS of FIG. 9d and HS′ of FIG. 9e superimposed over the present emittance model. HS″ is displaced from HS, along the z-axis, a displacement delta2, and additionally, HS″ is displaced from HS′, along the x- and z-axes, a displacement delta3.


Scores, in terms of separability, may be given to all the paring combinations, that is the models of FIGS. 9d-9e, relating to delta 1; the models of FIGS. 9d-9f, relating to delta2, and the models of FIGS. 9e-9f, relating to delta3. An optimal set of views may be selected based on its minimum or average scores for all the pairing combinations; for example, the optimal set of views may be that whose average score for all the pairing combinations is the highest. Alternatively, a weighted average may be applied.


In a second preferred embodiment, the quality function implements a reliability criterion, which is a measure of how reliably the intensity distribution of a given object may be reconstructed from the sampled views. Since the input to the reconstructed algorithm is a random sample, the output estimator is also random. A desired property of this output is that it be reliable in the sense that similar estimators for different projected samples (i.e. different sets of measurements) of the same input intensity are obtained with high confidence.


The Fisher Information, Fφ(I) is a measure, known in the art, which is used to evaluate the expected curvature of the likelihood of a model/(taken over measurements sampled from the model I). The Fisher Information is defined as:

Fφ(I)=−E∇2L(I)  (22)


The derivatives are taken with respect to the parameters of the intensity I, and the expectation is taken with respect to the random counts. Intuitively, a sharper curvature means that a maximum-likelihood estimation algorithm is more likely to produce a low-variance estimator. Indeed, this property is captured by the Cramer-Rao lower bound, which states that the inverse of the Fisher Information is a lower bound for the variance of any unbiased estimator, where an estimator f of the intensity I is unbiased if Ef(y)=I.


The Fisher information provides a value, FΦ(I)u,v, for each pair of voxels u and v. To provide a single measure for the reliability of the estimator, the average level, worst case, or other reasonable measure may be taken over the voxels. In the current example, the quality function is based on the average Fisher information. Starting with a set of emittance models to be used by the estimating algorithm to provide a reliable estimator, the average of the above measure is maximized over the entire set, defining the selected set Φ* as:










Φ
*

=


argmax
ϕ






I

I











u

U





[


F
Φ



(
I
)


]


u
,
u









(
23
)







Alternately, a set may be chosen to minimize the inverse of the Fisher information, by selecting for:










Φ
*

=


argmin
ϕ






I

I











u

U





[



F
Φ



(
I
)



-
1


]


u
,
u









(
24
)







Since inverting FΦ(I) may be computationally expensive, the [FΦ(I)−1]u,u term in Equation 24 may be replaced with 1/[FΦ(I)]u,u (thus neglecting the off-diagonal elements of the Fisher Information matrix FΦ(I)). Note that Equations 23 and 24 are not mathematically equivalent, and may therefore yield different selected sets.


In the preferred embodiment, evaluating the quality function is performed using the reliability criterion, and two or more emittance models are provided, having substantially identical volumes, but different modeled organ targets. Substantially identical sets of views are formed for all the emittance models, and each set is scored for reliability. One of the sets of views is then selected based on the average score for all the of the emittance models.


In a further preferred embodiment, a weighted combination of several information theoretic measures is used. For example, a plurality of models may be provided, all having substantially identical dimensions and volumes, as follows:


i. a pair of models with slightly different distributions of radioactive emission sources, as seen in FIGS. 9g and 9h, for scoring sets of views in terms of separability;


ii. a model with a given distribution of radioactive emission sources, as seen in any one of FIG. 9g, 9h or 9i, for scoring sets of views in terms of reliability.


Identical sets of views may be applied to all the models, and each view may be scored in terms of separability and reliability. An optimal set of views may be selected based on a summation of the two scores, or based on a weighted average of the two scores, where the relative weight given to each criterion reflects the relative importance of each measure per the given application.


The quality function is preferably defined in accordance with one of the following: worst case effectiveness for the given further view over the volume, average effectiveness for the given further view over the volume, worst case effectiveness for the given further view over the set of emittance models and average effectiveness for the given further view over the set of emittance models.


Maximization of the quality function may be performed utilizing any method known in the art such as simulated annealing and gradient ascent. In the simulated annealing (SA) method, each point of the search space is compared to a state of some physical system. The quality function to be maximized is interpreted as the internal energy of the system in that state. Therefore the goal is to bring the system from an arbitrary initial state to a state with the minimum possible energy.


The neighbors of each state and the probabilities of making a transition from each step to its neighboring states are specified. At each step, the SA heuristic probabilistically decides between moving the system to a neighboring state s′ or staying put in state s. The probabilities are chosen so that the system ultimately tends to move to states of lower energy. Typically this step is repeated until the system reaches and acceptable energy level.


Gradient ascent, on the other hand, is based on the observation that if a real-valued function F(x), such as the quality function of the present embodiments, is defined and differentiable in a neighborhood of a point a, then F(x) increases fastest if one goes from a in the direction of the gradient of F at a, ∇F(a). It follows that if:

b=a+γF(a)  (25)

for γ>0 a small enough number, then F(a)≦F(b). Gradient ascent starts with a guess x0 for a local maximum of F, and considers the sequence x0, x1, x2, . . . such that:

xn+1=xn+γ∇F(xn), n≧0.  (26)


Since F(x0)≦F(x1)≦F(x2)≦ . . . , the sequence (xn) is expected converges to a local maximum.


Preferably, the set of views selected with the quality function is increased by at least one randomly selected view. The randomly selected view(s) increase the probability that the quality of information obtained with the further views is maximized globally rather than locally.


As discussed above, selecting the best set of size N from amongst a large set of candidate projections is computationally complex. Since the size of the collection of views and of the required set may be large, a brute force scheme might not be computationally feasible.


In an additional preferred embodiment, a so-called “greedy algorithm” is used to incrementally construct larger and larger sets, until the required number of further views is defined. When multiple further views are required, it is computationally complex to maximize the quality function over all possible combinations of further views. The greedy algorithm reduces the computational burden by selecting the further views one at a time. The algorithm starts with a current set of views, and for each iteration determines a single view that yields the maximum improvement of the set score (hence the name “greedy”).


In theoretical terms, assume ρ(•) is the quality measure we are using for the view selection, and assume without loss of generality that we are trying to maximize this measure. We gradually build a set W of projections as follows. We start with an empty set W=Ø, and at every stage choose the projection that maximizes the quality measure when added to the current set:

W←arg maxW′{ρ(W′)|W′=W∪{φ},φεΦ}  (27)


In other words, during a given iteration, a respective score is calculated for a combination of the previous set with each of the views which is not a member of the current set. The current set is then expanded by adding the view which yielded the highest respective score, and the expanded current set serves as the input to the following iteration. Thus the number of times the scoring function is calculated per iteration drops from iteration to iteration. For a large collection of possible views, the greedy algorithm reduces the total number of computations required for set selection.


Reference is now made to FIG. 10, which is a simplified flowchart of an iterative “greedy” method for defining further views, according to a preferred embodiment of the present invention. In step 1000 a collection of views of the body structure is provided. The collection of views includes possible further views for future measurements, preferably based on anatomical and other constraints. In step 1010, the set of views used for the previous emission measurements is established as a current set of views. In step 1020 the view set is incrementally increased by a single further view during each iteration, until the required number of further views has been selected.


Reference is now made to FIG. 11, which is a simplified flowchart of a single iteration of the view selection method of FIG. 10, according to a preferred embodiment of the present invention. The method of FIG. 11 expands the current set of views by a single view. The method begins with a current set of views, which is the predetermined set (step 1010 above) for the first iteration of the greedy algorithm, or the set formed at the end of the previous iteration (step 1120 below) for all subsequent iterations. In step 1100, a respective expanded set is formed for each view not yet in the current set of views. A given view's expanded set contains all the views of the current set of views as well as the given view. In step 1110, a respective score is calculated for each of the expanded sets using the quality function. In step 1120, the view which yielded the highest-scoring expanded set is selected as a further view, to be used for further radioactive emission measurements. Finally, in step 1130, the current set is equated to the highest-scoring expanded set by adding the selected view to the current set. The newly formed current set serves as an input to the subsequent iteration, until the desired number of views is attained.


Reference is now made to FIG. 12, which is a simplified flowchart of a method for dynamically defining further views, according to a third preferred embodiment of the present invention. In step 1210, a collection of possible further views for performing radioactive-emission measurements of the body structure is provided. Each of the views is associated with at least one viewing parameter. Preferably the viewing parameters consist of at least one the following: detector unit location, detector unit orientation, collection angle, and measurement duration.


In step 1220 at least one quality function is provided. Each quality function is for evaluating sets of views, essentially as described above. A single quality function may be used to select several sets of views, where each set of views contains a different number of views.


In step 1230, multiple sets of further views (where a set may include a single further view) are formed from the collection of views, using the quality function(s) provided in step 1220. In a first preferred embodiment, each of the sets is formed using a different one of the quality functions. In an alternate preferred embodiment, one or more of the quality functions are used to form more than one set of views, where sets formed with the same quality function have differing numbers of views.


In step 1240, a selected set of views is obtained from the sets formed in step 1230.


In a first preferred embodiment, the final set of views is obtained by choosing one of the sets formed in step 1230 using a set selection criterion. For example, a respective set is formed in step 1230 for the separability and reliability criteria independently. A set selection criterion which calculates an overall performance rating for a given set taking both criteria into account is defined, and the formed set with the highest overall rating is selected as the final set.


In another preferred embodiment, the selected set of views is obtained by merging the sets formed in step 1230 according to the relative importance of the respective quality function used to form each set.


In the preferred embodiment, the method further consists of providing at least one emittance model and/or reconstruction representing the radioactive-emission density distribution of the volume, and of evaluating with at least one of the quality functions of step 1220 is performed in relation to the emittance models.


As discussed above, since each view is associated with one or more parameters, the selected set yields a group of parameter values for performing effective detection of the intensity distribution of the body structure. For example, if each view is associated with a view location parameter the selected set defines a set of locations for collecting emission data from an object, in order to provide a high-quality reconstruction of the intensity distribution of the body structure.


Reference is now made to FIG. 13, which is a simplified block diagram of measurement unit for performing radioactive-emission measurements of a body structure, according to a preferred embodiment of the present invention. Measurement unit 1300 includes probe 1310, analyzer 1320 and view definer 1330. Probe 1310 performs the radioactive-emission measurements of the body structure. Radioactive-emission-measuring probe 1310 preferably comprises several detecting units, which may be of different geometries and different collection angles δ, within a housing (see for example FIG. 1d). Preferably, the orientation and/or collection angle of the individual collimators is controllable. Analyzer 1320 analyzes the radioactive-emission measurements obtained from probe 1310. View definer 1330 dynamically defines further views for measurements, based on the analysis provided by analysis unit 1320. The analysis and view definition are performed substantially as described above.


The abovedescribed methods for radioactive-emission measurements of a body structure begin by performing measurements at a predetermined set of views. The results of the initial measurements are then analyzed and further views are defined.


The initial set of views is preferably selected based on information theoretic measures that quantify the quality of the data fed to the reconstruction algorithm, in order to obtain the best data for reconstructing a three-dimensional image of the body structure.


In accordance with the present invention, our approach is delineated by the following process:

  • i. modeling the body structure or region of interest as a model of a volume U, possibly with one or several modeled organ targets HS, within anatomical constraints AC (see FIGS. 7a-7b);
  • ii. obtaining an optimal and permissible initial set of views for the modeled volume U; and
  • iii. initiating the radioactive-emission measurements at the optimal set of views.


The following section concentrates on the second step of the process, namely, obtaining the optimal and permissible set of initial views for performing the radioactive-emission measurements of the body structure. The initial predetermined set of views is denoted herein the optimal set of views.


We consider here the following problem: Assume that there is a large pool of candidate views to choose from, but due to time restrictions or other restrictions we are limited to a specific number of views N. Which are the best N projections in terms of the quality of the reconstruction? It is further assumed that the pool of projections may be constrained, and hence general sampling theorems (e.g., Radon Transform) cannot be applied. For instance, we consider a scenario in emission tomography where the detecting unit can be located on top of one face of a given volume but not on the others. In such cases, methods such as EM estimation do not clearly establish what is the best scanning scheme, in order to provide the best reconstruction of the radioactive intensity distribution of the volume.


Reference is now made to FIG. 14, which is a simplified flowchart of a method for selecting an optimal initial predetermined set of views of a volume to be imaged, according to a preferred embodiment of the present invention. In step 1400, an emittance model of the body structure being imaged is provided. The body structure is modeled as a volume, and is preferably represented as a collection of voxels. In contrast to the use of emittance models for evaluation of the quality function, in the following all of the emittance models are provided prior to performing the radioactive-emission measurements, and thus do not include reconstructions of the body structure based on a current round of measurements. In step 1410, a collection of views of the body structure is provided. Each view reflects a possible data collection option, from which the goal is to select a set of views which together provide the desired image reconstruction quality. Each view is associated with one or more viewing parameters, such as the location and orientation of the detecting unit with respect to the volume, as described in detail above.


The view set is preferably tailored to correspond to the actual limitations and constraints likely to be encountered when collecting data for a given object, that is a given body structure. For example, in medical applications the detecting unit can generally view the body structure only from certain distances and orientations, which are determined by the patient's anatomy and the structure of the probe. The view set may therefore be designed to contain views having only those viewing parameter values consistent with attainable distances and orientations. The view set may also be designed to contain views having only those viewing parameter values suitable for a given measurement scenario, for example having identical values for the type of radiopharmaceutical and time since administration.


Preferably, the collection of views represents a quantized continuum of views. The view collection thus reflects the detection probability distribution for a detecting unit making periodic measurements while moving along a trajectory.


Reference is now made to FIG. 15, which shows a non-limiting example of a volume 1500 with two views, V1 and V2. Volume 1500 divided into twenty-seven cube-shaped voxels (labeled 1-27). Each view reflects a number of viewing parameters, in the case shown a location (in the XYZ coordinate system), orientation (ω, θ and σ) and collection angle δ. The viewing parameters determine which voxels in volume 1500 are within a detecting unit's collection angle, and generally affect the probability of detection for each voxel. The probability of detection may be affected by additional factors, including the attenuation coefficients within the volume.


During data collection, the probability of detection for each voxel is dependent on the parameters outlined above. In the preferred embodiment, a respective projection is calculated for each view, giving the view's detection probability distribution (i.e. the detection probability for each voxel of the volume). For a given view, the associated projection will have significant detection probabilities only for those voxels within the sector defined by the detecting unit's collection angle and location and orientation, as illustrated in conjunction with FIG. 1j.


The detection probability distribution for each view, that is the group of probabilities for each voxel, for a given view, is calculated according to techniques known in the art for determining the detection probability of a radioactive emission from a given voxel, under the constraints specified by the viewing parameters, for example, based on computer simulations of the geometry and the other viewing parameters, delineated hereinabove.


Generally, more distant voxels will have a lower probability of detection than closer voxels along the same line of sight. The volume attenuation coefficient may be constant or may vary over the volume. Thus, different sets of views may be produced for a single sector, by defining volumes with differing attenuation coefficients. For example, bone tissue and muscle tissue have different attenuations. Anatomical knowledge of the body structure being imaged may be used to develop a model of the volume U with a non-uniform attenuation that reflects the expected attenuation of the given body structure.


Referring again to FIG. 14, in step 1420 a scoring function is provided. The scoring function rates the information that is obtainable from the volume using any set of views (containing at least one view) taken from the collection of views. Evaluation of the scoring function does not require data from previously performed measurements. This is in contrast to the quality function, which rates the information quality on the basis of previous emission measurements as well as by the characteristics of possible further views. Preferably, the scoring function is a function of the detection probability distribution, as given by the projections. Several preferred embodiments of the scoring function, are discussed in detail below.


In step 1425, sets of views are formed from the collection of views. A score is then calculated for each of the sets.


In step 1430, the scores calculated in step 1425 are used to select one of the formed sets of views as the optimal set of views for the abovedescribed methods of performing radioactive-emission measurements of a body structure. A given scoring function may be used to select a set in a number of ways. In a first preferred embodiment, a required number of views is specified, and the highest scoring set with the specified number of views is selected. In a second preferred embodiment, the user may specify a minimal score which is known to provide satisfactory information quality, and select the smallest set which provides the specified score. However given a large collection of views the required number of calculations may be prohibitive. A third preferred embodiment used to reduce the computational burden is the greedy algorithm embodiment, similar to the abovedescribed method of FIGS. 10-11.


As discussed above, the scoring function is a measure of the quality of information which may be gathered for the volume using the given set of views. Both the separability and reliability criteria discussed above may be utilized as scoring functions, substantially as described above.


In an additional preferred embodiment, the scoring function implements a uniformity criterion, to ensure uniform coverage of the volume. It is often desired to obtain a uniform reconstruction quality among a given set of voxels WU, where W is the set of voxels for which it is desired to obtain uniform detection. Note that by selecting the appropriate W, the uniformity criterion is applied to all or a portion of the body. The uniformity criterion ensures that the spread of the total influence of each element on the set of measurements is as uniform as possible. The uniformity criterion depends only on the collection of views Φ and requires no assumptions on the distribution I.


For a set Φ of projections, the total influence of an element, u, is given by ΣφεΦφ(u). Normalizing these values to PΦ(u), such that











u

W





P
Φ



(
u
)



=
1

,





a probability measure is obtained for which the entropy H(Φ) can serve as a uniformity measure:










H


(
Φ
)


=

-




u

W





P
Φ


log







P
Φ



(
u
)









(
28
)







The selected set Φ* is the set (containing the required number of views) that satisfies:

Φ*=arg maxΦH(Φ)  (29)


Reference is now made to FIG. 16, which illustrates the concept of uniform coverage of a volume. Volume 1600 consists of twenty-seven cube-shaped voxels, including voxels 1610 and 1620. A, B, C and D are three views of volume 1600, showing the detector position and orientation for each view. For the uniformity criterion, it is desired that the overall influence of voxels 1610 and 1620 be approximately equal.


Assume that the probabilities of detection are as follows:














View
Voxel 1610
Voxel 1620

















A
0.6
0


B
0.2
0.5


C
0
0.3


D
0
0.1









Consider two possible sets of views: set {A, B, C} and set {B, C, D}. For set {A, B, C}, the total contribution of voxel 1610 is 0.8 (0.6+0.2+0) and of voxel 1620 is 0.8 (0+0.5+0.3). Normalizing these values for set {A, B, C} gives a probability set of [0.5,0.5]: For set {B, C, D}, the total contribution of voxel 1610 is 0.2 (0.2+0+0) and of voxel 1620 is 0.9 (0.5+0.3+0.1). Normalizing these values for set for set {B, C, D} gives a probability set of [0.18,0.82]. Thus:

H({A,B,C})=−(0.5*log2 0.5+0.5*log2 0.5)=−(−0.5−0.5)=1
H({B,C,D})=−(0.18*log2 0.18+0.82*log2 0.82)=−(−0.44−0.07)=0.51

Set {A, B, C} is thus seen to provide a more uniform coverage of volume 1600 than set {B, C, D}.


In an additional preferred embodiment, the greedy algorithm is used to incrementally construct larger and larger sets, until a set containing the required number of views is obtained. As described for FIGS. 10-11 above, the algorithm starts with an initial set of views (which may be the empty set), and in each iteration adds the view that yields the maximum improvement of the set score, until the required optimal set of views is obtained.


In a further preferred embodiment of a method for the selection of an optimal set, multiple view sets are first formed from one or more scoring functions, and then a final selection is made of one of the resulting sets, substantially as described for FIG. 12 above. Preferably the viewing parameters consist of at least one the following: detector location, detector orientation, viewing angle, material, thickness, collimator length, septa thickness, cell size, detection duration, time of detection, and a type of radiopharmaceutical. Knowledge about the conditions under which a particular set of measurements will be taken may be used to determine the optimal final set of views for a given data collection scenario.


The abovedescribed methods may each be embodied as a computer program stored on a computer-readable storage medium. In the preferred embodiment, computer-readable storage medium contains a set of instructions for defining views for radioactive-emission measurements of the body structure. An analysis routine analyzes the radioactive-emission measurements obtained from a radioactive-emission-measuring probe, and a view definition routine dynamically defines further views for measurements, based on the analyzing.


The active vision embodiments discussed above may be provided as modifications to improve reconstruction reliability. Reference is now made to FIG. 17, which is a simplified flowchart of a method for stabilizing the reconstruction of an imaged volume with dynamic view definition, according to a preferred embodiment of the present invention. FIG. 17 presents an embodiment in which the modifications to improve reconstruction reliability are performed by the active vision techniques described herein. In step 1700 radioactive emission measurements are performed. In step 1710, the reliability of reconstruction is analyzed, for example utilizing the SVD analysis described above. In step 1720, a determination is made whether measurements are finished, preferably based on the reliability analysis performed in step 1710. If it is determined that modifications are required, further views are defined without ceasing the measurement process, and preferably utilizing active vision techniques. Further measurements are then performed at the newly-defined views, until an adequate reconstruction reliability is achieved or it is decided to terminate the measurement process for other reasons.


The above-described improved reconstruction reliability and active vision techniques enable high-quality reconstruction based on data collected from a limited collection of views. The ability to improve reconstruction reliability, both before and during imaging, reduces the need for repeated imaging to obtain additional measurements. Reconstructing the intensity distribution from a smaller quantity of collected data utilizing a stable reconstruction process simplifies the computational process. Furthermore, these methods enable resolving the current conflict between the relatively large-pixel detectors needed for measurement speed and data processing considerations, with the small-pixel detectors needed until now to obtain a high-resolution reconstruction. The abovedescribed embodiments are particularly suitable for medical imaging purposes, where a high-resolution image is needed and it is desired to minimize the difficulties of the patient undergoing the diagnostic testing or treatment.


This application claims the benefit of:

  • International Application PCT/IL2005/001215, filed Nov. 16, 2005;
  • International Application PCT/IL2005/001173, filed Nov. 9, 2005,
  • U.S. Application 60/700,318, filed Jul. 19, 2005;
  • U.S. Application 60/700,299, filed Jul. 19, 2005;
  • U.S. Application 60/700,317, filed Jul. 19, 2005;
  • U.S. Application 60/700,753, filed Jul. 20, 2005;
  • U.S. Application 60/700,752, filed Jul. 20, 2005;
  • U.S. Application 60/702,979, filed Jul. 28, 2005;
  • U.S. Application 60/720,034, filed Sep. 26, 2005;
  • U.S. Application 60/720,652, filed Sep. 27, 2005;
  • U.S. Application 60/720,541, filed Sep. 27, 2005,
  • U.S. Application 60/750,287, filed Dec. 13, 2005;
  • U.S. Application 60/750,334, filed Dec. 15, 2005;
  • U.S. Application 60/750,597, filed Dec. 15, 2005;
  • U.S. Application 60/800,845, filed May 17, 2006;
  • U.S. Application 60/800,846, filed May 17, 2006;
  • Israel Application 171346, filed Oct. 10, 2005;
  • Israel Application 172349, filed Nov. 27, 2005;
  • U.S. Application 60/741,440, filed Dec. 2, 2005;
  • International Application PCT/IL2006/000059, filed Jan. 15, 2006;
  • U.S. Application 60/763,458, filed Jan. 31, 2006;
  • International Application PCT/IL2006/000562, filed May 11, 2006;
  • U.S. Application 60/799,688, filed May 11, 2006; and
  • U.S. Application 60/816,970, filed Jun. 28, 2006.


    Information of all of which is herein incorporated in entirety by reference.


This application further incorporates by reference all the information of the International Application entitled “Imaging Protocols” which is being co-filed by the same assignee of the present invention on Jul. 19, 2006.


It is expected that during the life of this patent many relevant detection probes, detector types, radiation-based detection systems, algorithms for reconstruction and analyzing the reliability of reconstruction, and algorithms will be developed and the scope of the corresponding terms are intended to include all such new technologies a priori.


It is appreciated that certain features of the invention, which are, for clarity, described in the context of separate embodiments, may also be provided in combination in a single embodiment. Conversely, various features of the invention, which are, for brevity, described in the context of a single embodiment, may also be provided separately or in any suitable subcombination.


Although the invention has been described in conjunction with specific embodiments thereof, it is evident that many alternatives, modifications and variations will be apparent to those skilled in the art. Accordingly, it is intended to embrace all such alternatives, modifications and variations that fall within the spirit and broad scope of the appended claims. All publications, patents and patent applications mentioned in this specification are herein incorporated in their entirety by reference into the specification, to the same extent as if each individual publication, patent or patent application was specifically and individually indicated to be incorporated herein by reference. In addition, citation or identification of any reference in this application shall not be construed as an admission that such reference is available as prior art to the present invention.

Claims
  • 1. A method of radioactive-emission measurements of a body structure, comprising: performing radioactive-emission measurements of the body structure, at a predetermined set of views taken from different directions;analyzing the radioactive-emission measurements from the different view directions including extracting a detected photon count; anddynamically defining at least one view parameter for further radioactive measurements based on said photon count detected from the different view directions from at least a target region, andmatching said dynamically defined at least one view parameter to one or more particular views and view directions.
  • 2. The method of claim 1, wherein said performing is in vivo.
  • 3. The method of claim 1, comprising reconstructing a 3D non-coincidence image from said measurements.
  • 4. The method of claim 1, wherein said performing is for diagnostic purposes.
  • 5. The method of claim 1, wherein said performing is as part of a pre-scan.
  • 6. The method of claim 1, wherein said at least one acquisition parameters comprises a view selection of a plurality of views.
  • 7. The method of claim 6, further comprising performing radioactive-emission measurements of the body structure, at said plurality of views.
  • 8. The method of claim 6, wherein each of said views is associated with viewing parameters selected from the group consisting of: detector unit location, detector unit orientation and collection angle.
  • 9. The method of claim 6, wherein each of said views is associated with a measurement duration.
  • 10. The method of claim 1, wherein dynamically defining comprises defining at least one view measurement duration to achieve a desired photon count.
  • 11. The method of claim 1, wherein dynamically defining comprises defining at least one view measurement duration to achieve a desired total photon count from at least a region of interest.
  • 12. A method according to claim 1, wherein said analyzing comprises estimating a measurement error provided by a photon count at a given view.
  • 13. A method according to claim 12, wherein said defining comprises extending a duration of a view to obtain a required error value or better.
  • 14. A method according to claim 12, wherein said defining comprises adding views to obtain a required error value or better.
  • 15. A method according to claim 1, wherein said analyzing comprises determining that a photon collection rate at a given view is within a specified range.
  • 16. A method according to claim 1, wherein said defining comprises selecting views to reduce acquisition time from regions of low interest.
  • 17. A method according to claim 1, wherein said defining comprises selecting views to increase acquisition time viewing regions of high interest.
  • 18. A method according to claim 1, wherein said analyzing comprises identifying detector saturation at a given view.
  • 19. A method according to claim 18, wherein said defining comprises selecting views to reinforce views that have saturated.
  • 20. A method according to claim 1, wherein said analyzing comprises forming a reconstruction of a radioactive-emission density distribution of said body structure.
  • 21. A method according to claim 20, wherein said analyzing comprises identifying edges within said a volume reconstructed by said reconstructing.
  • 22. A method according to claim 21, wherein said defining comprises selecting views at an angle to said edges.
  • 23. A method according to claim 22, wherein said angle is a substantially sharp angle.
  • 24. A method according to claim 20, wherein said analyzing comprises identifying volumetric boundaries within said reconstruction.
  • 25. A method according to claim 24, wherein said defining comprises selecting views at an angle to said volumetric boundaries.
  • 26. A method according to claim 25, wherein said angle is substantially sharp angle.
  • 27. A method according to claim 20, wherein said defining comprises selecting views covering regions in close proximity to a suspected target organ target.
  • 28. The method of claim 1, wherein said method is performed iteratively.
  • 29. The method of claim 20, wherein said method is performed iteratively.
  • 30. A method according to claim 29, wherein said analyzing comprises determining a variance of a plurality of reconstructions.
  • 31. A method according to claim 30, wherein said defining comprises selecting views to reduce the variance between iterations of reconstruction of a region of high interest to a required level.
  • 32. A method according to claim 20, wherein said analyzing comprises determining an accuracy of said reconstruction.
  • 33. A method according to claim 20, wherein said analyzing comprises determining a resolution of said reconstruction.
  • 34. A method according to claim 1, wherein said dynamically defining comprises: providing a quality function for rating a quality of information obtainable from said measurements by said quality function; andselecting a set of views to maximize said quality function.
  • 35. A method according to claim 34, wherein the quality function is based on an information theoretic measure of separability.
  • 36. A method according to claim 34, wherein the quality function is based on the reliability of reconstruction of a radioactive-emission density distribution.
  • 37. A method according to claim 34, wherein the quality function is based on an information-theoretic Fisher information measure.
  • 38. A method according to claim 37, wherein the Fisher information measure is a non-scalar measure.
  • 39. A method according to claim 34, further comprising increasing said selected set of views by at least one randomly selected view.
  • 40. The method of claim 34, wherein said dynamically defining comprises: providing a collection of views of said body structure;establishing said predetermined set of views as a current set of views; anditeratively expanding said current set until a predefined number of further views are obtained, wherein a set expansion iteration comprises:forming a respective expanded set for said collection of views and not in said current set, wherein an expanded set comprises said current set and a respective view;calculating a respective score for each of said expanded sets using said quality function;selecting said respective view of a highest-scoring expanded set as a further view; andequating said current set to said highest-scoring expanded set.
  • 41. The method of claim 1, wherein dynamically defining comprises selecting a set of views, selected by a method of: providing a model of a volume;providing a collection of views of said volume;providing a scoring function, by which any set of at least one view from said collection is scorable with a score that rates information obtained from said volume by said set;forming sets of views and scoring them, by said scoring function; andselecting a set of views from said sets, based on said scoring function.
  • 42. The method of claim 41, wherein each of said views is associated with viewing parameters comprising a detector location and orientation.
  • 43. The method of claim 41, wherein each of said views is associated with viewing parameters selected from the group consisting: of measurement duration, time elapsed from the administration of the pharmaceutical to the measurement, radiopharmaceutical half life, radioactive emission type, and radioactive emission energy.
  • 44. The method of claim 41, wherein said selecting is to obtain a minimal-size set which attains a predefined score.
  • 45. The method of claim 41, wherein providing a model comprises providing at least one emittance model, wherein an emittance model comprises a representation of a radioactive-emission density distribution of said volume.
  • 46. The method of claim 45, wherein said providing at least one emittance model includes modeling a radioactive-emission density distribution to form at least one modeled organ target.
  • 47. The method of claim 45, wherein: said providing at least one emittance model comprises providing a plurality of emittance models of substantially identical volumes, but different modeled organ targets;the forming sets of views from the collection of views and scoring them, comprises forming substantially identical sets of views for all the models and scoring the sets with respect to each model, based on the information theoretic measure of reliability; andthe selecting a sets of views, based on its score, comprises selecting based on the average score for the plurality of models.
  • 48. The method of claim 45, wherein: said providing at least one emittance model comprises providing a pair of models of substantially identical volumes, but different modeled organ targets, wherein the difference between the modeled organ targets is defined by at least one delta;the scoring function is based on an information theoretic measure of separability;the forming the sets of views from the collection of views and scoring them, comprises forming substantially identical sets of views for the pair of models and scoring the sets with respect to the pair; andthe selecting a sets of views, based on its score, comprises selecting based on the score for the pair.
  • 49. The method of claim 45, wherein: said providing at least one emittance model comprises providing a plurality of pairs of models of substantially identical volumes, but different modeled organ targets, wherein the difference between the modeled organ targets for each pair is defined by at least one delta;the scoring function is based on an information theoretic measure of separability;the forming the sets of views from the collection of views and scoring them, comprises forming substantially identical sets of views for the plurality of pairs of models and scoring the sets with respect to each of the pairs; andthe selecting a sets of views, based on its score, comprises selecting based on an average score for the plurality of pairs.
  • 50. The method of claim 45, wherein: the modeling the body-structure comprises providing at least one model of no radiation emission density distribution and at least two models of different radiation emission density distributions, so as to form different modeled organ targets for those two models; andthe providing the scoring function comprises providing a scoring function as a combination of uniformity, separability and reliability.
  • 51. The method of claim 45, wherein said scoring function comprises a worst-case effectiveness for said given view over said set of emittance models.
  • 52. The method of claim 45, wherein said scoring function comprises an average effectiveness for said given view over said set of emittance models.
  • 53. The method of claim 45, wherein said set of emittance models comprises a single emittance model.
  • 54. The method of claim 45, wherein said at least one emittance model corresponds to said body structure.
  • 55. The method of claim 54, wherein said body structure comprises one of a group of body structures comprising: a prostate, a heart, a brain, a breast, a uterus, an ovary, a liver, a kidney, a stomach, a colon, a small intestine, an oral cavity, a throat, a gland, a lymph node, the skin, another body organ, a limb, a bone, another part of the body, and a whole body.
  • 56. A method of radioactive-emission measurements of a body structure, comprising: performing radioactive-emission measurements of the body structure, at a predetermined set of views;analyzing the radioactive-emission measurements including extracting a detected photon count; anddynamically defining at least one view parameter for further radioactive measurements based on said photon count detected from at least a target region,wherein said at least one view parameter is one or more of a detector position, a detector angular orientation, and the duration of particular views.
Priority Claims (2)
Number Date Country Kind
171346 Oct 2005 IL national
172349 Nov 2005 IL national
RELATED APPLICATIONS

This application is a continuation of U.S. patent application Ser. No. 11/988,926 filed on Feb. 5, 2009, which is a National Phase of PCT Patent Application No. PCT/IL2006/000840 having International Filing Date of Jul. 19, 2006, which claims the benefit of priority of U.S. Provisional Patent Application Nos. 60/816,970 filed on Jun. 28, 2006, 60/800,846 filed on May 17, 2006, 60/800,845 filed on May 17, 2006, 60/799,688 filed on May 11, 2006, 60/763,458 filed on Jan. 31, 2006, 60/750,597 filed on Dec. 15, 2005, 60/750,334 filed on Dec. 15, 2005, 60/750,287 filed on Dec. 13, 2005, 60/741,440 filed on Dec. 2, 2005,60/720,652 filed on Sep. 27, 2005, 60/720,541 filed on Sep. 27, 2005, 60/720,034 filed on Sep. 26, 2005, 60/702,979 filed on Jul. 28, 2005, 60/700,753 filed on July 20; 2005, 60/700,752 filed on Jul. 20, 2005, 60/700,318 filed on Jul. 19, 2005, 60/700,299 filed on Jul. 19, 2005 and 60/700,317 filed on Jul. 19, 2005. PCT/IL2006/000840 is also a Continuation-in-Part of PCT Patent Application Nos. PCT/IL2006/000562 filed on May 11, 2006, PCT/IL2006/000059 filed on Jan. 15, 2006, PCT/IL2005/001215 filed on Nov. 16, 2005 and PCT/IL2005/001173 filed on Nov. 9, 2005. PCT/IL2006/000840 claims the benefit of priority of Israel Patent Application Nos. 172349 filed on Nov. 27, 2005 and 171346 filed on Oct. 10, 2005. The contents of the above Applications are all incorporated herein by reference.

US Referenced Citations (749)
Number Name Date Kind
630611 Knapp et al. Aug 1899 A
2776377 Anger Jan 1957 A
3340866 Noeller Sep 1967 A
3446965 Ogier et al. May 1969 A
3535085 Shumate et al. Oct 1970 A
3684887 Hugonin Aug 1972 A
3690309 Pluzhnikov et al. Sep 1972 A
3719183 Schwartz Mar 1973 A
3739279 Hollis Jun 1973 A
3971362 Pope et al. Jul 1976 A
3978337 Nickles et al. Aug 1976 A
3988585 O'Neill et al. Oct 1976 A
4000502 Butler et al. Dec 1976 A
4015592 Bradley-Moore Apr 1977 A
4055765 Gerber et al. Oct 1977 A
4061919 Miller et al. Dec 1977 A
4095107 Genna et al. Jun 1978 A
4165462 Macovski et al. Aug 1979 A
4181856 Bone Jan 1980 A
4278077 Mizumoto Jul 1981 A
4289969 Cooperstein et al. Sep 1981 A
4291708 Frei et al. Sep 1981 A
4296785 Vitello et al. Oct 1981 A
4302675 Wake et al. Nov 1981 A
4364377 Smith Dec 1982 A
4383327 Kruger May 1983 A
4476381 Rubin Oct 1984 A
4503331 Kovacs, Jr. et al. Mar 1985 A
4521688 Yin Jun 1985 A
H12 Bennett et al. Jan 1986 H
4580054 Shimoni Apr 1986 A
4595014 Barrett et al. Jun 1986 A
4674107 Urban et al. Jun 1987 A
4679142 Lee Jul 1987 A
4689041 Corday et al. Aug 1987 A
4689621 Kleinberg Aug 1987 A
4709382 Sones Nov 1987 A
4710624 Alvarez et al. Dec 1987 A
4731536 Rische et al. Mar 1988 A
4773430 Porath Sep 1988 A
4782840 Martin, Jr. et al. Nov 1988 A
4791934 Brunnett Dec 1988 A
4801803 Denen et al. Jan 1989 A
4828841 Porter et al. May 1989 A
4834112 Machek et al. May 1989 A
4844067 Ikada et al. Jul 1989 A
4844076 Lesho et al. Jul 1989 A
4853546 Abe et al. Aug 1989 A
4854324 Hirschman et al. Aug 1989 A
4854330 Evans, III et al. Aug 1989 A
4867962 Abrams Sep 1989 A
4893013 Denen et al. Jan 1990 A
4893322 Hellmick et al. Jan 1990 A
4919146 Rhinehart et al. Apr 1990 A
4924486 Weber et al. May 1990 A
4928250 Greenberg et al. May 1990 A
4929832 Ledly May 1990 A
4938230 Machek et al. Jul 1990 A
4951653 Fry et al. Aug 1990 A
4959547 Carroll et al. Sep 1990 A
4970391 Uber, III Nov 1990 A
4995396 Inaba et al. Feb 1991 A
5014708 Hayashi et al. May 1991 A
5018182 Cowan et al. May 1991 A
5032729 Charpak Jul 1991 A
5033998 Corday et al. Jul 1991 A
5039863 Matsuno et al. Aug 1991 A
5042056 Hellmick et al. Aug 1991 A
5070877 Mohiuddin et al. Dec 1991 A
5070878 Denen Dec 1991 A
5088492 Takayama et al. Feb 1992 A
5115137 Andersson-Engels et al. May 1992 A
5119818 Carroll et al. Jun 1992 A
5132542 Bassalleck et al. Jul 1992 A
5145163 Cowan et al. Sep 1992 A
5151598 Denen Sep 1992 A
5170055 Carroll et al. Dec 1992 A
5170439 Zeng et al. Dec 1992 A
5170789 Narayan et al. Dec 1992 A
5196796 Misic et al. Mar 1993 A
5210421 Gullberg et al. May 1993 A
5243988 Sieben et al. Sep 1993 A
5246005 Carroll et al. Sep 1993 A
5249124 DeVito Sep 1993 A
5252830 Weinberg Oct 1993 A
5254101 Trombley, III Oct 1993 A
5258717 Misic et al. Nov 1993 A
5263077 Cowan et al. Nov 1993 A
5279607 Schentag et al. Jan 1994 A
5284147 Hanaoka et al. Feb 1994 A
5299253 Wessels Mar 1994 A
5304165 Haber et al. Apr 1994 A
5307808 Dumoulin et al. May 1994 A
5307814 Kressel et al. May 1994 A
5309959 Shaw et al. May 1994 A
5317506 Coutre et al. May 1994 A
5317619 Hellmick et al. May 1994 A
5323006 Thompson et al. Jun 1994 A
5329976 Haber et al. Jul 1994 A
5334141 Carr et al. Aug 1994 A
5349190 Hines et al. Sep 1994 A
5355087 Claiborne et al. Oct 1994 A
5365069 Eisen et al. Nov 1994 A
5365928 Rhinehart et al. Nov 1994 A
5367552 Peschmann Nov 1994 A
5377681 Drane Jan 1995 A
5381791 Qian Jan 1995 A
5383456 Arnold et al. Jan 1995 A
5383858 Reilly et al. Jan 1995 A
5386446 Fujimoto et al. Jan 1995 A
5387409 Nunn et al. Feb 1995 A
5391877 Marks Feb 1995 A
5395366 D'Andrea Mar 1995 A
5399868 Jones et al. Mar 1995 A
5404293 Weng et al. Apr 1995 A
5415181 Hofgrefe et al. May 1995 A
5431161 Ryals et al. Jul 1995 A
5435302 Lenkinski et al. Jul 1995 A
5436458 Tran et al. Jul 1995 A
5441050 Thurston et al. Aug 1995 A
5448073 Jeanguillaume Sep 1995 A
5451232 Rhinehart et al. Sep 1995 A
5472403 Cornacchia et al. Dec 1995 A
5475219 Olson Dec 1995 A
5475232 Powers et al. Dec 1995 A
5476095 Schnall et al. Dec 1995 A
5479969 Hardie et al. Jan 1996 A
5481115 Hsieh et al. Jan 1996 A
5484384 Fearnot Jan 1996 A
5489782 Wernikoff Feb 1996 A
5493595 Schoolman Feb 1996 A
5493805 Penuela et al. Feb 1996 A
5494036 Uber, III et al. Feb 1996 A
5501674 Trombley, III et al. Mar 1996 A
5517120 Misik et al. May 1996 A
5519221 Weinberg May 1996 A
5519222 Besett May 1996 A
5519931 Reich May 1996 A
5520182 Leighton et al. May 1996 A
5520653 Reilly et al. May 1996 A
5521506 Misic et al. May 1996 A
5536945 Reich Jul 1996 A
5545899 Tran et al. Aug 1996 A
5559335 Zeng et al. Sep 1996 A
5565684 Gullberg et al. Oct 1996 A
5569181 Heilman et al. Oct 1996 A
5572132 Pulyer et al. Nov 1996 A
5572999 Funda et al. Nov 1996 A
5579766 Gray Dec 1996 A
5580541 Wells et al. Dec 1996 A
5585637 Bertelsen et al. Dec 1996 A
5587585 Eisen et al. Dec 1996 A
5591143 Trombley, III et al. Jan 1997 A
5600145 Plummer Feb 1997 A
5604531 Iddan et al. Feb 1997 A
5610520 Misic Mar 1997 A
5617858 Taverna et al. Apr 1997 A
5629524 Stettner et al. May 1997 A
5635717 Popescu Jun 1997 A
5657759 Essen-Moller Aug 1997 A
5672877 Liebig et al. Sep 1997 A
5677539 Apotovsky et al. Oct 1997 A
5682888 Olson et al. Nov 1997 A
5687542 Lawecki et al. Nov 1997 A
5690691 Chen et al. Nov 1997 A
5692640 Caulfield et al. Dec 1997 A
5694933 Madden et al. Dec 1997 A
5695500 Taylor et al. Dec 1997 A
5716595 Goldenberg Feb 1998 A
5727554 Kalend et al. Mar 1998 A
5729129 Acker Mar 1998 A
5732704 Thurston et al. Mar 1998 A
5739508 Uber, III Apr 1998 A
5741232 Reilly et al. Apr 1998 A
5742060 Ashburn Apr 1998 A
5744805 Raylman et al. Apr 1998 A
5757006 De Vito et al. May 1998 A
5779675 Reilly et al. Jul 1998 A
5780855 Pare et al. Jul 1998 A
5781442 Engleson et al. Jul 1998 A
5784432 Kurtz et al. Jul 1998 A
5786597 Lingren et al. Jul 1998 A
5795333 Reilly et al. Aug 1998 A
5799111 Guissin Aug 1998 A
5800355 Hasegawa Sep 1998 A
5803914 Ryals et al. Sep 1998 A
5806519 Evans, III et al. Sep 1998 A
5808203 Nolan, Jr. et al. Sep 1998 A
5810742 Pearlman Sep 1998 A
5811814 Leone et al. Sep 1998 A
5813985 Carroll Sep 1998 A
5818050 Dilmanian et al. Oct 1998 A
5821541 Tuemer Oct 1998 A
5825031 Wong et al. Oct 1998 A
5827219 Uber, III et al. Oct 1998 A
5828073 Zhu et al. Oct 1998 A
5833603 Kovacs et al. Nov 1998 A
5838009 Plummer et al. Nov 1998 A
5840026 Uber, III et al. Nov 1998 A
5841141 Gullberg et al. Nov 1998 A
5842977 Lesho et al. Dec 1998 A
5843037 Uber, III Dec 1998 A
5846513 Carroll et al. Dec 1998 A
5847396 Lingren et al. Dec 1998 A
5857463 Thurston et al. Jan 1999 A
5871013 Wainer et al. Feb 1999 A
5873861 Hitchins et al. Feb 1999 A
5880475 Oka et al. Mar 1999 A
5882338 Gray Mar 1999 A
5884457 Ortiz et al. Mar 1999 A
5885216 Evans, III et al. Mar 1999 A
5891030 Johnson et al. Apr 1999 A
5893397 Peterson et al. Apr 1999 A
5899885 Reilly et al. May 1999 A
5900533 Chou May 1999 A
5903008 Li May 1999 A
5910112 Judd et al. Jun 1999 A
5911252 Cassel Jun 1999 A
5916167 Kramer et al. Jun 1999 A
5916197 Reilly et al. Jun 1999 A
5920054 Uber, III Jul 1999 A
5927351 Zhu et al. Jul 1999 A
5928150 Call Jul 1999 A
5932879 Raylman et al. Aug 1999 A
5938639 Reilly et al. Aug 1999 A
5939724 Eisen et al. Aug 1999 A
5944190 Edelen Aug 1999 A
5944694 Hitchins et al. Aug 1999 A
5947935 Rhinehart et al. Sep 1999 A
5953884 Lawecki et al. Sep 1999 A
5954668 Uber, III et al. Sep 1999 A
5961457 Raylman et al. Oct 1999 A
5967983 Ashburn Oct 1999 A
5973598 Beigel Oct 1999 A
5974165 Giger et al. Oct 1999 A
5984860 Shan Nov 1999 A
5987350 Thurston Nov 1999 A
5993378 Lemelson Nov 1999 A
5997502 Reilly et al. Dec 1999 A
6002134 Lingren Dec 1999 A
6002480 Izatt et al. Dec 1999 A
6017330 Hitchins et al. Jan 2000 A
6019745 Gray Feb 2000 A
6021341 Scibilia et al. Feb 2000 A
6026317 Verani Feb 2000 A
6037595 Lingren Mar 2000 A
6040697 Misic Mar 2000 A
6042565 Hirschman et al. Mar 2000 A
RE36648 Uber, III et al. Apr 2000 E
6046454 Lingren et al. Apr 2000 A
6048334 Hirschman et al. Apr 2000 A
6052618 Dahlke et al. Apr 2000 A
6055450 Ashburn Apr 2000 A
6055452 Pearlman Apr 2000 A
RE36693 Reich May 2000 E
6063052 Uber et al. May 2000 A
D426891 Beale et al. Jun 2000 S
D426892 Beale et al. Jun 2000 S
6072177 McCroskey et al. Jun 2000 A
6076009 Raylman et al. Jun 2000 A
6080984 Friesenhahn Jun 2000 A
D428491 Beale et al. Jul 2000 S
6082366 Andra et al. Jul 2000 A
6090064 Reilly et al. Jul 2000 A
6091070 Lingren et al. Jul 2000 A
6096011 Trombley, III et al. Aug 2000 A
6107102 Ferrari Aug 2000 A
6115635 Bourgeois Sep 2000 A
6129670 Burdette et al. Oct 2000 A
6132372 Essen-Moller Oct 2000 A
6135955 Madden et al. Oct 2000 A
6135968 Brounstein Oct 2000 A
6137109 Hayes Oct 2000 A
6145277 Lawecki et al. Nov 2000 A
6147352 Ashburn Nov 2000 A
6147353 Gagnon et al. Nov 2000 A
6148229 Morris, Sr. et al. Nov 2000 A
6149627 Uber, III Nov 2000 A
6155485 Coughlin et al. Dec 2000 A
6160398 Walsh Dec 2000 A
6162198 Coffey et al. Dec 2000 A
6172362 Lingren et al. Jan 2001 B1
6173201 Front Jan 2001 B1
6184530 Hines et al. Feb 2001 B1
6189195 Reilly et al. Feb 2001 B1
6194715 Lingren et al. Feb 2001 B1
6194725 Colsher et al. Feb 2001 B1
6194726 Pi et al. Feb 2001 B1
6197000 Reilly et al. Mar 2001 B1
6202923 Boyer et al. Mar 2001 B1
6203775 Torchilin et al. Mar 2001 B1
6205347 Morgan et al. Mar 2001 B1
6212423 Krakovitz Apr 2001 B1
6223065 Misic et al. Apr 2001 B1
6224577 Dedola et al. May 2001 B1
6226350 Hsieh May 2001 B1
6229145 Weinberg May 2001 B1
6232605 Soluri et al. May 2001 B1
6233304 Hu et al. May 2001 B1
6236050 Tumer May 2001 B1
6236878 Taylor et al. May 2001 B1
6236880 Raylman et al. May 2001 B1
6239438 Schubert May 2001 B1
6240312 Alfano et al. May 2001 B1
6241708 Reilly et al. Jun 2001 B1
6242743 DeVito Jun 2001 B1
6242744 Soluri et al. Jun 2001 B1
6242745 Berlad et al. Jun 2001 B1
6246901 Benaron Jun 2001 B1
6252924 Davantes et al. Jun 2001 B1
6258576 Richards-Kortum et al. Jul 2001 B1
6259095 Bouton et al. Jul 2001 B1
6261562 Xu et al. Jul 2001 B1
6263229 Atalar et al. Jul 2001 B1
6269340 Ford et al. Jul 2001 B1
6270463 Morris, Sr. et al. Aug 2001 B1
6271524 Wainer et al. Aug 2001 B1
6271525 Majewski et al. Aug 2001 B1
6280704 Schutt et al. Aug 2001 B1
6281505 Hines et al. Aug 2001 B1
6308097 Pearlman Oct 2001 B1
6310968 Hawkins et al. Oct 2001 B1
6315981 Unger Nov 2001 B1
6317623 Griffiths et al. Nov 2001 B1
6317648 Sleep et al. Nov 2001 B1
6318630 Coughlin et al. Nov 2001 B1
6322535 Hitchins et al. Nov 2001 B1
6323648 Belt et al. Nov 2001 B1
6324418 Crowley et al. Nov 2001 B1
RE37487 Reilly et al. Dec 2001 E
D452737 Nolan, Jr. et al. Jan 2002 S
6336913 Spohn et al. Jan 2002 B1
6339652 Hawkins et al. Jan 2002 B1
6339718 Zatezalo et al. Jan 2002 B1
6344745 Reisker et al. Feb 2002 B1
6346706 Rogers et al. Feb 2002 B1
6346886 de la Huerga Feb 2002 B1
RE37602 Uber, III et al. Mar 2002 E
6353227 Boxen Mar 2002 B1
6356081 Misic Mar 2002 B1
6368331 Front et al. Apr 2002 B1
6371938 Reilly et al. Apr 2002 B1
6375624 Uber, III et al. Apr 2002 B1
6377838 Iwanczyk et al. Apr 2002 B1
6381349 Zeng et al. Apr 2002 B1
6385483 Uber, III et al. May 2002 B1
6388244 Gagnon May 2002 B1
6388258 Berlad et al. May 2002 B1
6392235 Barrett et al. May 2002 B1
6396273 Misic May 2002 B2
6397098 Uber, III et al. May 2002 B1
6399951 Paulus et al. Jun 2002 B1
6402717 Reilly et al. Jun 2002 B1
6402718 Reilly et al. Jun 2002 B1
6407391 Mastrippolito et al. Jun 2002 B1
6408204 Hirschman Jun 2002 B1
6409987 Cardin et al. Jun 2002 B1
6415046 Kerut, Sr. Jul 2002 B1
6420711 Tuemer Jul 2002 B2
6425174 Riech Jul 2002 B1
6426917 Tabanou et al. Jul 2002 B1
6429431 Wilk Aug 2002 B1
6431175 Penner et al. Aug 2002 B1
6432089 Kakimi et al. Aug 2002 B1
6438401 Cheng et al. Aug 2002 B1
6439444 Shields, II Aug 2002 B1
6440107 Trombley, III et al. Aug 2002 B1
6442418 Evans, III et al. Aug 2002 B1
6448560 Tumer Sep 2002 B1
6453199 Kobozev Sep 2002 B1
6459925 Nields et al. Oct 2002 B1
6459931 Hirschman Oct 2002 B1
6468261 Small et al. Oct 2002 B1
6469306 Van Dulmen et al. Oct 2002 B1
6471674 Emig et al. Oct 2002 B1
6480732 Tanaka et al. Nov 2002 B1
6484051 Daniel Nov 2002 B1
6488661 Spohn et al. Dec 2002 B1
6490476 Townsend et al. Dec 2002 B1
6504157 Juhi Jan 2003 B2
6504178 Carlson et al. Jan 2003 B2
6504899 Pugachev et al. Jan 2003 B2
6506155 Sluis et al. Jan 2003 B2
6510336 Daghighian et al. Jan 2003 B1
6512374 Misic et al. Jan 2003 B1
6516213 Nevo Feb 2003 B1
6519569 White et al. Feb 2003 B1
6520930 Critchlow et al. Feb 2003 B2
6522945 Sleep et al. Feb 2003 B2
6525320 Juni Feb 2003 B1
6525321 Juni Feb 2003 B2
6541763 Lingren et al. Apr 2003 B2
6545280 Weinberg et al. Apr 2003 B2
6549646 Yeh et al. Apr 2003 B1
6560354 Maurer et al. May 2003 B1
6562008 Reilly et al. May 2003 B1
6563942 Takeo et al. May 2003 B2
6565502 Bede et al. May 2003 B1
6567687 Front et al. May 2003 B2
6575930 Trombley, III et al. Jun 2003 B1
6576918 Fu et al. Jun 2003 B1
6583420 Nelson et al. Jun 2003 B1
6584348 Glukhovsky Jun 2003 B2
6585700 Trocki et al. Jul 2003 B1
6587710 Wainer Jul 2003 B1
6589158 Winkler Jul 2003 B2
6591127 McKinnon Jul 2003 B1
6592520 Peszynski et al. Jul 2003 B1
6602488 Daghighian Aug 2003 B1
6607301 Glukhovsky et al. Aug 2003 B1
6611141 Schulz et al. Aug 2003 B1
6614453 Suri et al. Sep 2003 B1
6620134 Trombley, III et al. Sep 2003 B1
6627893 Zeng et al. Sep 2003 B1
6628983 Gagnon Sep 2003 B1
6628984 Weinberg Sep 2003 B2
6630735 Carlson et al. Oct 2003 B1
6631284 Nutt et al. Oct 2003 B2
6632216 Houzego et al. Oct 2003 B2
6633658 Dabney et al. Oct 2003 B1
6638752 Contag et al. Oct 2003 B2
6643537 Zatezalo et al. Nov 2003 B1
6643538 Majewski et al. Nov 2003 B1
6652489 Trocki et al. Nov 2003 B2
6657200 Nygard et al. Dec 2003 B2
6662036 Cosman Dec 2003 B2
6664542 Ye et al. Dec 2003 B2
6670258 Carlson et al. Dec 2003 B2
6671563 Engleson et al. Dec 2003 B1
6673033 Sciulli et al. Jan 2004 B1
6674834 Acharya et al. Jan 2004 B1
6676634 Spohn et al. Jan 2004 B1
6677182 Carlson et al. Jan 2004 B2
6677755 Belt et al. Jan 2004 B2
6680750 Tournier et al. Jan 2004 B1
6694172 Gagnon et al. Feb 2004 B1
6697660 Robinson Feb 2004 B1
6699219 Emig et al. Mar 2004 B2
6704592 Reynolds et al. Mar 2004 B1
6713766 Garrard et al. Mar 2004 B2
6714012 Belt et al. Mar 2004 B2
6714013 Misic Mar 2004 B2
6716195 Nolan, Jr. et al. Apr 2004 B2
6722499 Reich Apr 2004 B2
6723988 Wainer Apr 2004 B1
6726657 Dedig et al. Apr 2004 B1
6728583 Hallett Apr 2004 B2
6731971 Evans, III et al. May 2004 B2
6731989 Engleson et al. May 2004 B2
6733477 Cowan et al. May 2004 B2
6733478 Reilly et al. May 2004 B2
6734416 Carlson et al. May 2004 B2
6734430 Soluri et al. May 2004 B2
6737652 Lanza et al. May 2004 B2
6737866 Belt et al. May 2004 B2
6740882 Weinberg et al. May 2004 B2
6743202 Hirschman et al. Jun 2004 B2
6743205 Nolan, Jr. et al. Jun 2004 B2
6747454 Belt Jun 2004 B2
6748259 Benaron et al. Jun 2004 B1
6751500 Hirschman et al. Jun 2004 B2
6765981 Heumann Jul 2004 B2
6766048 Launay et al. Jul 2004 B1
6771802 Patt et al. Aug 2004 B1
6774358 Hamill et al. Aug 2004 B2
6776977 Liu Aug 2004 B2
6787777 Gagnon et al. Sep 2004 B1
6788758 De Villiers Sep 2004 B2
6798206 Misic Sep 2004 B2
6808513 Reilly et al. Oct 2004 B2
6809321 Rempel Oct 2004 B2
6813868 Baldwin et al. Nov 2004 B2
6821013 Reilly et al. Nov 2004 B2
6822237 Inoue et al. Nov 2004 B2
6833705 Misic Dec 2004 B2
6838672 Wagenaar et al. Jan 2005 B2
6841782 Balan et al. Jan 2005 B1
6843357 Bybee et al. Jan 2005 B2
6851615 Jones Feb 2005 B2
6866654 Callan et al. Mar 2005 B2
6870175 Dell et al. Mar 2005 B2
6881043 Barak Apr 2005 B2
6888351 Belt et al. May 2005 B2
6889074 Uber, III et al. May 2005 B2
6897658 Belt et al. May 2005 B2
6906330 Blevis et al. Jun 2005 B2
D507832 Yanniello et al. Jul 2005 S
6915170 Engleson et al. Jul 2005 B2
6915823 Osborne et al. Jul 2005 B2
6917828 Fukuda Jul 2005 B2
6921384 Reilly et al. Jul 2005 B2
6928142 Shao et al. Aug 2005 B2
6935560 Andreasson et al. Aug 2005 B2
6936030 Pavlik et al. Aug 2005 B1
6937750 Natanzon et al. Aug 2005 B2
6939302 Griffiths et al. Sep 2005 B2
6940070 Tumer Sep 2005 B2
6943355 Shwartz et al. Sep 2005 B2
6957522 Baldwin et al. Oct 2005 B2
6958053 Reilly Oct 2005 B1
6963770 Scarantino et al. Nov 2005 B2
6970735 Uber, III et al. Nov 2005 B2
6972001 Emig et al. Dec 2005 B2
6974443 Reilly et al. Dec 2005 B2
6976349 Baldwin et al. Dec 2005 B2
6984222 Hitchins et al. Jan 2006 B1
6985870 Martucci et al. Jan 2006 B2
6988981 Hamazaki Jan 2006 B2
6994249 Peterka et al. Feb 2006 B2
7009183 Wainer et al. Mar 2006 B2
7011814 Suddarth et al. Mar 2006 B2
7012430 Misic Mar 2006 B2
7017622 Osborne et al. Mar 2006 B2
7018363 Cowan et al. Mar 2006 B2
7019783 Kindem et al. Mar 2006 B2
7025757 Reilly et al. Apr 2006 B2
7026623 Oaknin et al. Apr 2006 B2
7043063 Noble et al. May 2006 B1
7102138 Belvis et al. Sep 2006 B2
7103204 Celler et al. Sep 2006 B1
7127026 Amemiya et al. Oct 2006 B2
7142634 Engler et al. Nov 2006 B2
7145986 Wear et al. Dec 2006 B2
7147372 Nelson et al. Dec 2006 B2
7164130 Welsh et al. Jan 2007 B2
7176466 Rousso et al. Feb 2007 B2
7187790 Sabol et al. Mar 2007 B2
7217953 Carlson May 2007 B2
7256386 Carlson et al. Aug 2007 B2
7291841 Nelson et al. Nov 2007 B2
7327822 Sauer et al. Feb 2008 B2
7359535 Salla et al. Apr 2008 B2
7373197 Daighighian et al. May 2008 B2
7394923 Zou et al. Jul 2008 B2
7444010 De Man Oct 2008 B2
7468513 Charron et al. Dec 2008 B2
7470896 Pawlak et al. Dec 2008 B2
7490085 Walker et al. Feb 2009 B2
7495225 Hefetz et al. Feb 2009 B2
7502499 Grady Mar 2009 B2
7570732 Stanton et al. Aug 2009 B2
7592597 Hefetz et al. Sep 2009 B2
7620444 Le et al. Nov 2009 B2
7627084 Jabri et al. Dec 2009 B2
7652259 Kimchy et al. Jan 2010 B2
7671331 Hefez Mar 2010 B2
7671340 Uribe et al. Mar 2010 B2
7672491 Krishnan et al. Mar 2010 B2
7680240 Manjeshwar et al. Mar 2010 B2
7705316 Rousso et al. Apr 2010 B2
7734331 Dhawale et al. Jun 2010 B2
7826889 David et al. Nov 2010 B2
7831024 Metzler et al. Nov 2010 B2
7835927 Schlotterbeck et al. Nov 2010 B2
7872235 Rousso et al. Jan 2011 B2
7894650 Weng et al. Feb 2011 B2
7968851 Rousso et al. Jun 2011 B2
8013308 Guerin et al. Sep 2011 B2
8055329 Kimchy et al. Nov 2011 B2
8111886 Rousso et al. Feb 2012 B2
8158951 Bal et al. Apr 2012 B2
8163661 Akiyoshi et al. Apr 2012 B2
8204500 Weintraub et al. Jun 2012 B2
8338788 Zilberstein et al. Dec 2012 B2
8440168 Yang et al. May 2013 B2
20010016029 Tumer Aug 2001 A1
20010020131 Kawagishi et al. Sep 2001 A1
20010035902 Iddan et al. Nov 2001 A1
20010049608 Hochman Dec 2001 A1
20020068864 Bishop et al. Jun 2002 A1
20020072784 Sheppard, Jr. et al. Jun 2002 A1
20020085748 Baumberg Jul 2002 A1
20020087101 Barrick et al. Jul 2002 A1
20020099295 Gil et al. Jul 2002 A1
20020099310 Kimchy et al. Jul 2002 A1
20020099334 Hanson et al. Jul 2002 A1
20020103429 DeCharms Aug 2002 A1
20020103431 Toker et al. Aug 2002 A1
20020145114 Inoue et al. Oct 2002 A1
20020148970 Wong et al. Oct 2002 A1
20020165491 Reilly Nov 2002 A1
20020168094 Kaushikkar et al. Nov 2002 A1
20020168317 Daighighian et al. Nov 2002 A1
20020172405 Schultz Nov 2002 A1
20020179843 Tanaka et al. Dec 2002 A1
20020183645 Nachaliel Dec 2002 A1
20020188197 Bishop et al. Dec 2002 A1
20020198738 Osborne Dec 2002 A1
20030001098 Stoddart et al. Jan 2003 A1
20030001837 Baumberg Jan 2003 A1
20030006376 Tumer Jan 2003 A1
20030013950 Rollo et al. Jan 2003 A1
20030013966 Barnes et al. Jan 2003 A1
20030038240 Weinberg Feb 2003 A1
20030055685 Cobb et al. Mar 2003 A1
20030063787 Natanzon et al. Apr 2003 A1
20030071219 Motomura et al. Apr 2003 A1
20030081716 Tumer May 2003 A1
20030135388 Martucci et al. Jul 2003 A1
20030136912 Juni Jul 2003 A1
20030144322 Kozikowski et al. Jul 2003 A1
20030147887 Wang et al. Aug 2003 A1
20030158481 Stotzka et al. Aug 2003 A1
20030183226 Brand et al. Oct 2003 A1
20030189174 Tanaka et al. Oct 2003 A1
20030191430 D'Andrea et al. Oct 2003 A1
20030202629 Dunham et al. Oct 2003 A1
20030208117 Shwartz et al. Nov 2003 A1
20030215122 Tanaka Nov 2003 A1
20030215124 Li Nov 2003 A1
20030216631 Bloch et al. Nov 2003 A1
20030219149 Vailaya et al. Nov 2003 A1
20040003001 Shimura Jan 2004 A1
20040010397 Barbour et al. Jan 2004 A1
20040015075 Kimchy et al. Jan 2004 A1
20040021065 Weber Feb 2004 A1
20040044282 Mixon et al. Mar 2004 A1
20040051368 Caputo et al. Mar 2004 A1
20040054248 Kimchy et al. Mar 2004 A1
20040054278 Kimchy et al. Mar 2004 A1
20040065838 Tumer Apr 2004 A1
20040075058 Blevis et al. Apr 2004 A1
20040081623 Eriksen et al. Apr 2004 A1
20040082918 Evans et al. Apr 2004 A1
20040084340 Morelle et al. May 2004 A1
20040086437 Jackson et al. May 2004 A1
20040101176 Mendonca et al. May 2004 A1
20040101177 Zahlmann et al. May 2004 A1
20040116807 Amrami et al. Jun 2004 A1
20040120557 Sabol Jun 2004 A1
20040122311 Cosman Jun 2004 A1
20040125918 Shanmugavel et al. Jul 2004 A1
20040138557 Le et al. Jul 2004 A1
20040143449 Behrenbruch et al. Jul 2004 A1
20040144925 Stoddart et al. Jul 2004 A1
20040153128 Suresh et al. Aug 2004 A1
20040162492 Kobayashi Aug 2004 A1
20040171924 Mire et al. Sep 2004 A1
20040183022 Weinberg Sep 2004 A1
20040184644 Leichter et al. Sep 2004 A1
20040193453 Butterfield et al. Sep 2004 A1
20040195512 Crosetto Oct 2004 A1
20040204646 Nagler et al. Oct 2004 A1
20040205343 Forth et al. Oct 2004 A1
20040210126 Hajaj et al. Oct 2004 A1
20040238743 Gravrand et al. Dec 2004 A1
20040251419 Nelson et al. Dec 2004 A1
20040253177 Elmaleh et al. Dec 2004 A1
20040263865 Pawlak et al. Dec 2004 A1
20050001170 Juni Jan 2005 A1
20050006589 Young et al. Jan 2005 A1
20050020898 Vosniak et al. Jan 2005 A1
20050020915 Bellardinelli et al. Jan 2005 A1
20050023474 Persyk et al. Feb 2005 A1
20050029277 Tachibana Feb 2005 A1
20050033157 Klein et al. Feb 2005 A1
20050049487 Johnson et al. Mar 2005 A1
20050055174 David et al. Mar 2005 A1
20050056788 Juni Mar 2005 A1
20050074402 Cagnolini et al. Apr 2005 A1
20050107698 Powers et al. May 2005 A1
20050107914 Engleson et al. May 2005 A1
20050108044 Koster May 2005 A1
20050113945 Engleson et al. May 2005 A1
20050121505 Metz et al. Jun 2005 A1
20050131270 Weil et al. Jun 2005 A1
20050145797 Oaknin et al. Jul 2005 A1
20050148869 Masuda Jul 2005 A1
20050149350 Kerr et al. Jul 2005 A1
20050156115 Kobayashi et al. Jul 2005 A1
20050173643 Tumer Aug 2005 A1
20050187465 Motomura et al. Aug 2005 A1
20050198800 Reich Sep 2005 A1
20050203389 Williams Sep 2005 A1
20050205792 Rousso et al. Sep 2005 A1
20050205796 Bryman Sep 2005 A1
20050207526 Altman Sep 2005 A1
20050211909 Smith Sep 2005 A1
20050215889 Patterson, II Sep 2005 A1
20050234424 Besing et al. Oct 2005 A1
20050247893 Fu et al. Nov 2005 A1
20050253073 Joram et al. Nov 2005 A1
20050261936 Silverbrook et al. Nov 2005 A1
20050261937 Silverbrook et al. Nov 2005 A1
20050261938 Silverbrook et al. Nov 2005 A1
20050266074 Zilberstein et al. Dec 2005 A1
20050277833 Williams, Jr. Dec 2005 A1
20050277911 Stewart et al. Dec 2005 A1
20050278066 Graves et al. Dec 2005 A1
20050288869 Kroll et al. Dec 2005 A1
20060000983 Charron et al. Jan 2006 A1
20060033028 Juni Feb 2006 A1
20060036157 Tumer Feb 2006 A1
20060072799 McLain Apr 2006 A1
20060074290 Chen et al. Apr 2006 A1
20060104519 Stoeckel et al. May 2006 A1
20060109950 Arenson et al. May 2006 A1
20060122503 Burbank et al. Jun 2006 A1
20060145081 Hawman Jul 2006 A1
20060160157 Zuckerman Jul 2006 A1
20060188136 Ritt et al. Aug 2006 A1
20060214097 Wang et al. Sep 2006 A1
20060237652 Kimchy et al. Oct 2006 A1
20060257012 Kaufman et al. Nov 2006 A1
20070081700 Blumenfeld et al. Apr 2007 A1
20070116170 De Man et al. May 2007 A1
20070133852 Collins et al. Jun 2007 A1
20070156047 Nagler et al. Jul 2007 A1
20070166227 Liu et al. Jul 2007 A1
20070189436 Goto et al. Aug 2007 A1
20070194241 Rousso et al. Aug 2007 A1
20070265230 Rousso et al. Nov 2007 A1
20080001090 Ben-Haim et al. Jan 2008 A1
20080029704 Hefetz et al. Feb 2008 A1
20080033291 Rousso et al. Feb 2008 A1
20080036882 Uemura et al. Feb 2008 A1
20080039721 Shai et al. Feb 2008 A1
20080042067 Rousso et al. Feb 2008 A1
20080128626 Rousso et al. Jun 2008 A1
20080137938 Zahniser Jun 2008 A1
20080230702 Rousso et al. Sep 2008 A1
20080230705 Rousso et al. Sep 2008 A1
20080237482 Shahar et al. Oct 2008 A1
20080260228 Dichterman et al. Oct 2008 A1
20080260580 Helle et al. Oct 2008 A1
20080260637 Dickman Oct 2008 A1
20080277591 Shahar et al. Nov 2008 A1
20090001273 Hawman Jan 2009 A1
20090018412 Schmitt Jan 2009 A1
20090078875 Rousso et al. Mar 2009 A1
20090112086 Melman Apr 2009 A1
20090152471 Rousso et al. Jun 2009 A1
20090190807 Rousso et al. Jul 2009 A1
20090201291 Ziv et al. Aug 2009 A1
20090236532 Frach et al. Sep 2009 A1
20090304582 Rousso et al. Dec 2009 A1
20100006770 Balakin Jan 2010 A1
20100021378 Rousso et al. Jan 2010 A1
20100102242 Burr et al. Apr 2010 A1
20100121184 Dhawale et al. May 2010 A1
20100140483 Rousso et al. Jun 2010 A1
20100202664 Busch et al. Aug 2010 A1
20100245354 Rousso et al. Sep 2010 A1
20120172699 Nagler et al. Jul 2012 A1
20120248320 Wangerin et al. Oct 2012 A1
20120326034 Sachs et al. Dec 2012 A1
20130114792 Zilberstein et al. May 2013 A1
20130308749 Zilberstein et al. Nov 2013 A1
20140151563 Rousso et al. Jun 2014 A1
Foreign Referenced Citations (63)
Number Date Country
1516429 Dec 1969 DE
19814199 Oct 1999 DE
19815362 Oct 1999 DE
0273257 Jul 1988 EP
0525954 Feb 1993 EP
0526970 Feb 1993 EP
0543626 May 1993 EP
0592093 Apr 1994 EP
0697193 Feb 1996 EP
0813692 Dec 1997 EP
0887661 Dec 1998 EP
1237013 Sep 2002 EP
2031142 Apr 1980 GB
59-141084 Aug 1984 JP
61-026879 Feb 1986 JP
01-324568 Jun 1986 JP
03-121549 May 1991 JP
04-151120 May 1992 JP
06-109848 Apr 1994 JP
07-059763 Mar 1995 JP
07-141523 Jun 1995 JP
08-292268 Nov 1996 JP
10-260258 Sep 1998 JP
11-072564 Mar 1999 JP
WO 9200402 Jan 1992 WO
WO 9816852 Apr 1998 WO
WO 9903003 Jan 1999 WO
WO 9930610 Jun 1999 WO
WO 9939650 Aug 1999 WO
WO 0010034 Feb 2000 WO
WO 0018294 Apr 2000 WO
WO 0022975 Apr 2000 WO
WO 0025268 May 2000 WO
WO 0031522 Jun 2000 WO
WO 0038197 Jun 2000 WO
WO 0189384 Nov 2001 WO
WO 0216965 Feb 2002 WO
WO 02058531 Aug 2002 WO
WO 02075357 Sep 2002 WO
WO 03073938 Sep 2003 WO
WO 03086170 Oct 2003 WO
WO 2004004787 Jan 2004 WO
WO 2004032151 Apr 2004 WO
WO 2004042546 May 2004 WO
WO 2004113951 Dec 2004 WO
WO 2005002971 Jan 2005 WO
WO 2005059592 Jun 2005 WO
WO 2005059840 Jun 2005 WO
WO 2005067383 Jul 2005 WO
WO 2005104939 Nov 2005 WO
WO 2005118659 Dec 2005 WO
WO 2005119025 Dec 2005 WO
WO 2006042077 Apr 2006 WO
WO 2006051531 May 2006 WO
WO 2006054296 May 2006 WO
WO 2006075333 Jul 2006 WO
WO 2006129301 Dec 2006 WO
WO 2007010534 Jan 2007 WO
WO 2007010537 Jan 2007 WO
WO 2007054935 May 2007 WO
WO 2007074467 Jul 2007 WO
WO 2008010227 Jan 2008 WO
WO 2008075362 Jun 2008 WO
Non-Patent Literature Citations (449)
Entry
Communication Pursuant to Article 94(3) EPC Dated Oct. 26, 2012 From the European Patent Office Re. Application No. 05803689.8.
Official Action Dated Aug. 5, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/309,479.
Official Action Dated Aug. 14, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/448,473.
Berman et al. “D-SPECT: A Novel Camera for High Speed Quantitative Molecular Imaging: Initial Clinical Results”, The Journal of Nuclear Medicine, 47(Suppl.1): 131P, 2006.
Berman et al. “Myocardial Perfusion Imaging With Technetium-99m-Sestamibi: Comparative Analysis of Available Imaging Protocols”, The Journal of Nuclear Medicine, 35: 681-688, 1994.
Borges-Neto et al. “Perfusion and Function at Rest and Treadmill Exercise Using Technetium-99m-Sestamibi: Comparison of One- and Two-Day Protocols in Normal Volunteers”, The Journal of Nuclear Medicine, 31(7): 1128-1132, Jul. 1990.
Kwok et al. “Feasability of Simultaneous Dual-Isotope Myocardial Perfusion Acquisition Using a Lower Dose of Sestamibi”, European Journal of Nuclear Medicine, 24(3): 281-285, Mar. 1997.
Patton et al. “D-SPECT: A New Solid State Camera for High Speed Molecular Imaging”, The Journal of Nuclear Medicine, 47(Suppl.1): 189P, 2006.
Restriction Official Action Dated Aug. 16, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/448,473.
Supplemental Notice of Allowability Dated Aug. 20, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,653.
Communication Pursuant to Article 94(3) EPC Dated Sep. 16, 2013 From the European Patent Office Re.: Application No. 06832278.3.
Communication Pursuant to Article 94(3) EPC Dated Sep. 23, 2011 From the European Patent Office Re.: Application No. 06832278.3.
Written Opinion Dated Nov. 1, 2007 from the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL06/00840.
Written Opinion Dated Jul. 11, 2008 from the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL06/01511.
Applicant-Initiated Interview Summary Dated Jan. 28, 2013 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Official Action Dated Feb. 7, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,653.
Official Action Dated Dec. 19, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/448,473.
Official Action Dated Apr. 16, 2012 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Advisory Action Before the Filing of an Appeal Brief Dated Feb. 26, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/989,223.
Official Action Dated Nov. 30, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/989,223.
Official Action Dated Jun. 12, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/087,150.
Bacharach et al. “Attenuation Correction in Cardiac Positron Emission Tomography and Single-Photon Emission Computed Tomography”, Journal of Nucelar Cardiology, 2(3): 246-255, 1995.
Uni Magdeburg “Attenuation Map”, University of Magdeburg, Germany, Retrieved From the Internet, Archived on Jul. 31, 2002.
Zaidi et al. “Determination of the Attenuation Map in Emission Tomography”, Journal of Nuclear Medicine, 44(2): 291-315, 2003.
Official Action Dated Jun. 21, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/309,479.
Notice of Allowance Dated Jul. 19, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/989,223.
Notice of Allowance Dated Jul. 25, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 13/345,719.
Communication Pursuant to Article 94(3) EPC Dated Sep. 17, 2012 From the European Patent Office Re. Application No. 06832278.3.
Ouyang et al. “Incorporation of Correlated Structural Images in PET Image Reconstruction”, IEEE Transactions of Medical Imaging, 13(4): 627-640, Dec. 1994.
Notice of Allowance Dated Sep. 28, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/792,856.
Summons to Attend Oral Proceedings Pursuant to Rule 115(1) EPC Dated Nov. 29, 2012 From the European Patent Office Re. Application No. 06756259.5.
Brzymialkiewicz et al. “Evaluation of Fully 3-D Emission Mammotomography With a Compact Cadmium Zinc Telluride Detector”, IEEE Transactions on Medical Imaging, 24(7): 868-877, Jul. 2005.
Jan et al. “Preliminary Results From the AROPET”, IEEE Nuclear Science Symposium Conference Record, Nov. 4-10, 2001, 3: 1607-1610, 2001.
Ohno et al. “Selection of Optimum Projection Angles in Three Dimensional Myocardial SPECT”, IEEE Nuclear Science Symposium Conference Record 2001, 4: 2166-2169, 2001.
Seret et al. “Intrinsic Uniformity Requirements for Pinhole SPECT”, Journal of Nuclear Medicine Technology, 34(1): 43-47, Mar. 2006.
Smither “High Resolution Medical Imaging System for 3-D Imaging of Radioactive Sources With 1 mm FWHM Spatial Resolution”, Proceedings of the SPIE, Medical Imaging 2003: Physics of Medical Imaging, 5030: 1052-1060, Jun. 9, 2003.
Tornai et al. “A 3D Gantry Single Photon Emission Tomograph With Hemispherical Coverage for Dedicated Breast Imaging”, Nuclear Instruments & Methods in Physics Research, Section A, 497: 157-167, 2003.
Official Action Dated Oct. 7, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/750,057.
Official Action Dated Oct. 7, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/932,872.
Notice of Allowance Dated Feb. 2, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/628,074.
Advisory Action Before the Filing of an Appeal Brief Dated Feb. 26, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/514,785.
Notice of Allowance Dated Feb. 21, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/798,017.
Notice of Panel Decision From Pre-Appeal Brief Review Dated Feb. 29, 2012 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Official Action Dated Mar. 1, 2012 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Official Action Dated Dec. 22, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/514,785.
Official Action Dated Feb. 28, 2012 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Notice of Allowance Dated Mar. 7, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 13/726,316.
Official Action Dated Mar. 6, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/792,856.
Restriction Official Action Dated Mar. 9, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/976,852.
Official Action Dated Mar. 11, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 13/345,719.
Communication Pursuant to Article 94(3) EPC Dated May 12, 2010 From the European Patent Office Re. Application No. 06809851.6.
International Preliminary Report on Patentability Dated Apr. 7, 2009 From the International Bureau of WIPO Re. Application No. PCT/IL2007/000918.
International Preliminary Report on Patentability Dated Jan. 13, 2009 From the International Bureau of WIPO Re. Application No. PCT/IL2006/000834.
International Preliminary Report on Patentability Dated May 14, 2008 From the International Bureau of WIPO Re. Application No. PCT/IL2006/001291.
International Preliminary Report on Patentability Dated May 15, 2007 From the International Bureau of WIPO Re. Application No. PCT/IL2005/001173.
International Search Report Dated Jul. 1, 2008 From the International Searching Authority Re. Application No. PCT/IL2006/000834.
International Search Report Dated Jul. 2, 2007 From the International Searching Authority Re. Application No. PCT/IL2006/001291.
International Search Report Dated Aug. 3, 2006 From the International Searching Authority Re. Application No. PCT/IL2005/001173.
International Search Report Dated Oct. 15, 2008 From the International Searching Authority Re. Application No. PCT/2007/000918.
Official Action Dated Nov. 1, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/728,383.
Official Action Dated Mar. 2, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,617.
Official Action Dated Dec. 8, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,690.
Official Action Dated Apr. 9, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/798,017.
Official Action Dated Aug. 13, 2008 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/769,826.
Official Action Dated Feb. 16, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/747,378.
Official Action Dated May 18, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/989,223.
Official Action Dated Apr. 20, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/798,017.
Official Action Dated Dec. 20, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/309,479.
Official Action Dated Sep. 21, 2009 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/798,017.
Official Action Dated May 23, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,690.
Official Action Dated Sep. 30, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/798,017.
Bracco Diagnostics “Cardiotec®: Kit for the Preparation of Technetium Tc 99m Teboroxime. For Diagnostic Use”, Bracco Diagnostics Inc., Product Sheet, 2 P., Jul. 2003.
Bracco Diagnostics “Techneplex®: Kit for the Preparation of Technetium Tc 99m Pentetate Injection. Diagnostic—for Intravenous Use”, Bracco Diagnostics™, Product Sheet, 5 P., Jun. 1995.
Dewaraja et al. “Accurate Dosimetry in 131I Radionuclide Therapy Using Patient-Specific, 3-Dimensional Methods for SPECT Reconstruction and Basorbed Dose Calculation”, The Journal of Nuclear Medicine, 46(5): 840-849, May 2005.
GE Healthcare “Myoview™: Kit for the Preparation of Technetium Tc99m Tetrofosmin for Injection. Diagnostic Radiopharmaceutical. For Intravenous Use Only. Rx Only”, GE Healthcare, Product Sheet, 4 P., Aug. 2006.
Mallinckrodt “Kit for the Preparation of Technetium Tc 99m Sestamibi Injection”, Mallinckrodt Inc., Product Sheet, 2 P., Sep. 8, 2008.
Mallinckrodt “OctreoScan®: Kit for the Preparation of Indium In-111 Pentetreotide. Diagnostic—for Intravenous Use. Rx Only”, Mallinckrodt Inc., Product Sheet, 2 P., Oct. 25, 2006.
Pharmalucence “Kit for the Preparation of Technetium Tc99m Sulfur Colloid Injection for Subcutaneous, Intraperitoneal, Intravenous, and Oral Use”, Pharmalucence Inc., Reference ID: 2977567, Prescribing Information, 10 P., Jul. 2011.
Saltz et al. “Interim Report of Randomized Phase II Trial of Cetuximab/Bevacizumab/Irinotecan (CBI) Versus Cetuximab/Bevacizumab (CB) in Irinotecan-Refractory Colorectal Cancer”, Gastrointestinal Cancer Symposium, Hollywood, FL, USA, Jan. 27-29, 2005, American Society of Clinical Oncology, Abstract 169b, 4P., 2005.
Trikha et al. “Monoclonal Antibodies as Therapeutics in Oncology”, Current Opinion in Biotechnology, 13: 609-614, 2002.
Volkow et al. “Imaging the Living Human Brain: Magnetic Resonance Imaging and Positron Emission Tomography”, Proc. Natl. Acad. Sci. USA, 94: 2787-2788, Apr. 1997.
Notice of Allowance Dated Jun. 14, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 10/616,307.
Notice of Allowance Dated Mar. 14, 2013 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Studen “Compton Camera With Position-Sensitive Silicon Detectors”, Doctoral Thesis, University of Ljubljana, Faculty of Mathematics and Physics, 36 P.
Notice of Allowance Dated Apr. 10, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/976,852.
Notice of Allowance Dated Jun. 21, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,653.
Official Action Dated Sep. 5, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 13/947,198.
Official Action Dated Aug. 31, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,653.
Cancer Medicine “Radiolabeled Monoclonal Antibodies. Historical Perspective”, Cancer Medicine, 5th Ed., Sec.16: Principles of Biotherapeutics, Chap.65: Monoclonal Serotherapy, 2000.
Lange et al. “EM Reconstruction Algorithms for Emission and Transmission Tomography”, Journal of Computer Assisted Tomography, 8(2): 306-316, Apr. 1984.
Ohrvall et al. “Intraoperative Gamma Detection Reveals Abdominal EndocrineTumors More Efficiently Than Somatostatin Receptor Scintigraphy”, 6th Conference on Radioimmunodetection and Radioimmunotherapy of Cancer, Cancer, 80: 2490-2494, 1997.
Rockmore et al. “A Maximum Likelihood Approach to Emission Image Reconstruction From Projections”, IEEE Transactions on Nuclear Science, 23(4): 1428-1432, Aug. 1976.
Shepp et al. “Maximum Likelihood Reconstruction for Emission Tomography”, IEEE Transactions on Medical Imaging, MI-1: 113-122, Oct. 1982.
Sitek et al. “Reconstruction of Dynamic Renal Tomographic Data Acquired by Slow Rotation”, The Journal of Nuclear Medicine, 42(11): 1704-1712, Nov. 2001.
Solanki “The Use of Automation in Radiopharmacy”, Hospital Pharmacist, 7(4): 94-98, Apr. 2000.
Weldon et al. “Quantification of Inflammatory Bowel Disease Activity Using Technetium-99m HMPAO Labelled Leucocyte Single Photon Emission Computerised Tomography (SPECT)”, Gut, 36: 243-250, 1995.
Official Action Dated Oct. 11, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/976,852.
Official Action Dated Oct. 10, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/798,017.
Communication Pursuant to Article 94(3) EPC Dated Nov. 12, 2012 From the European Patent Office Re. Application No. 06756258.7.
Notice of Allowance Dated Nov. 15, 2012 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,683.
Amendment After Allowance Under 37 CFR 1.312 Dated Sep. 13, 2010 to Notice of Allowance of Jul. 22, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/794,799.
Appeal Brief Dated Jan. 19, 2010 to Notice of Appeal of Nov. 16, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Communication Pursuant to Article 93(3) EPC Dated Mar. 8, 2010 From the European Patent Office Re.: Application No. 06832278.3.
Communication Pursuant to Article 94(3) EPC Dated Mar. 2, 2011 From the European Patent Office Re.: Application No. 06756259.5.
Communication Pursuant to Article 94(3) EPC Dated Mar. 8, 2010 From the European Patent Office Re.: Application No. 06832278.3.
Communication Pursuant to Article 94(3) EPC Dated May 12, 2010 From the European Patent Office Re.: Application No. 06809851.6.
Communication Pursuant to Article 94(3) EPC Dated Apr. 16, 2010 From the European Patent Office Re. Application No. 01951883.6.
Communication Pursuant to Article 94(3) EPC Dated Nov. 18, 2011 From the European Patent Office Re. Application No. 05803689.8.
Communication Pursuant to Article 94(3) EPC Dated Oct. 21, 2009 From the European Patent Office Re.: Application No. 02716285.8.
Communication Pursuant to Article 94(3) EPC Dated Jul. 22, 2009 From the European Patent Office Re.: Application No. 06809851.6.
Communication Pursuant to Article 94(3) EPC Dated Sep. 22, 2011 From the European Patent Office Re. Application No. 06756258.7.
Communication Pursuant to Article 96(2) EPC Dated Jun. 19, 2006 From the European Patent Office Re.: Application No. 03810570.6.
Communication Pursuant to Article 96(2) EPC Dated Aug. 30, 2007 From the European Patent Office Re.: Application No. 03810570.6.
Communication Pursuant to Rules 70(2) and 70a(2) EPC Dated Apr. 4, 2011 From the European Patent Office Re. Application No. 05803689.8.
Communication Relating to the Results of the Partial International Search Dated Apr. 18, 2007 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL2006/001291.
Communication Relating to the Results of the Partial international Search Dated May 21, 2008 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL2007/001588.
Examination Report Dated Jun. 22, 2011 From the Government of India, Patent Office, Intellectual Property Building Re. Application No. 2963/CIIENP/2006.
International Preliminary Report on Patentability Dated Apr. 16, 2009 From the International Bureau of WIPO Re.: Applicaiton No. PCT/IL2007/000918.
International Preliminary Report on Patentability Dated Jun. 21, 2007 From the International Bureau of WIPO Re.: Application No. PCT/IL2005/000575.
International Preliminary Report on Patentability Dated Jan. 22, 2009 From the International Bureau of WIPO Re.: Application No. PCT/IL2006/000834.
International Preliminary Report on Patentability Dated Jan. 22, 2009 From the International Bureau of WIPO Re.: Application No. PCT/IL2006/001511.
International Preliminary Report on Patentability Dated May 22, 2007 From the International Preliminary Examining Authority Re.: Application No. PCT/IL06/00059.
International Preliminary Report on Patentability Dated May 22, 2008 From the International Bureau of WIPO Re.: Application No. PCT/IL2006/001291.
International Preliminary Report on Patentability Dated May 24, 2007 From the International Bureau of WIPO Re.: Application No. PCT/IL2005/001173.
International Preliminary Report on Patentability Dated Apr. 26, 2007 From the International Bureau of WIPO Re.: Application No. PCT/IL2005/000394.
International Preliminary Report on Patentability Dated Jan. 31, 2008 From the International Bureau of WIPO Re.: Application No. PCT/IL2006/000840.
International Search Report Dated Oct. 10, 2006 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL06/00059.
International Search Report Dated Jul. 11, 2008 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL06/01511.
International Search Report Dated Jul. 25, 2008 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL2007/001588.
International Search Report Dated Feb. 1, 2006 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/00048.
International Search Report Dated Jul. 1, 2008 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL06/00834.
International Search Report Dated Nov. 1, 2007 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL06/00840.
International Search Report Dated Jul. 2, 2007 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL2006/001291.
International Search Report Dated Aug. 3, 2006 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/001173.
International Search Report Dated May 11, 2006 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/001215.
International Search Report Dated Sep. 11, 2002 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL01/00638.
International Search Report Dated Sep. 12, 2002 From the International Searching Authority of the Patent Cooperation Treaty Re: Application No. PCT/IL02/00057.
International Search Report Dated Oct. 15, 2008 From the International Searching Authority Re.: Application No. PCT/IL07/00918.
International Search Report Dated Mar. 18, 2004 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL03/00917.
International Search Report Dated Mar. 23, 2006 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/00572.
International Search Report Dated May 24, 2007 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/00575.
International Search Report Dated Mar. 26, 2007 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/00394.
Interview Summary Dated Mar. 25, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Interview Summary Dated May 31, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 10/616,301.
Invitation to Pay Additional Fees Dated Jul. 10, 2008 From the International Searching Authority Re.: Application No. PCT/IL06/01511.
Invitation to Pay Additional Fees Dated Feb. 15, 2007 From the International Searching Authority Re.: Application No. PCT/IL05/00575.
Notice of Allowance Dated May 5, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/240,239.
Notice of Allowance Dated May 6, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,617.
Notice of Allowance Dated Oct. 11, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/988,926.
Notice of Allowance Dated Nov. 15, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 10/616,301.
Notice of Allowance Dated Jul. 16, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 12/084,559.
Notice of Allowance Dated Sep. 16, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/607,075.
Notice of Allowance Dated Dec. 17, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,617.
Notice of Allowance Dated Sep. 17, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/533,568. Suppl. IDS VIII in 25855.
Notice of Allowance Dated Jul. 22, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/794,799.
Notice of Allowance Dated Feb. 23, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 12/728,383.
Notice of Allowance Dated Jun. 23, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Notice of Allowance Dated Nov. 23, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 12/084,559.
Notice of Allowance Dated Aug. 25, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,617.
Notice of Allowance Dated Jun. 30, 2010 From the United States Patent and Trademark Office Re.: U.S. Appl. No. 09/727,464.
Notice of Appeal and Pre-Appeal Brief Dated Jan. 4, 2010 to Official Action of Sep. 2, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/343,792.
Notice of Appeal Dated Nov. 16, 2009 to Official Action of Jul. 15, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Notice of Non-Compliant Amendment Dated Feb. 14, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Office Action Dated Dec. 2, 2007 From the Israeli Patent Office Re.: Application No. 158442.
Office Action Dated Jan. 2, 2006 From the Israeli Patent Office Re.: Application No. 154323.
Office Action Dated Sep. 4, 2007 From the Israeli Patent Office Re.: Application No. 157007.
Office Action Dated Jul. 17, 2007 From the Israeli Patent Office Re.: Application No. 154323.
Official Action Dated Jun. 1, 2006 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/686,536.
Official Action Dated Mar. 1, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/794,799.
Official Action Dated Nov. 1, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 12/728,383.
Official Action Dated Sep. 1, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/794,799.
Official Action Dated Jul. 2, 2004 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Official Action Dated Mar. 2, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Official Action Dated Mar. 2, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,617.
Official Action Dated Sep. 2, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/343,792.
Official Action Dated Aug. 3, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Official Action Dated May 3, 2007 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/240,239.
Official Action Dated Sep. 4, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/533,568.
Official Action Dated Sep. 5, 2002 From the US Patent and Trademark Office Re.: U.S. Appl. No. 12/084,559.
Official Action Dated Jan. 7, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Official Action Dated Jul. 7, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/533,568.
Official Action Dated Oct. 7, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated Apr. 8, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Official Action Dated Dec. 8, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/132,320.
Official Action Dated Dec. 8, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Official Action Dated Dec. 8, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Official Action Dated Jan. 8, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Official Action Dated Apr. 9, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Official Action Dated Mar. 9, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated Aug. 10, 2007 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Official Action Dated Nov. 10, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Official Action Dated Nov. 10, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Official Action Dated Aug. 11, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Official Action Dated Jul. 11, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,683.
Official Action Dated Mar. 11, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/607,075.
Official Action Dated Jul. 12, 2007 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated Jul. 12, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Official Action Dated Sep. 12, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,653.
Official Action Dated Sep. 12, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/087,150.
Official Action Dated Dec. 13, 2007 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated May 13, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Official Action Dated Sep. 13, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/976,852.
Official Action Dated May 14, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Official Action Dated Apr. 15, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Official Action Dated Dec. 15, 2006 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated Dec. 15, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/343,792.
Official Action Dated Feb. 15, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/343,792.
Official Action Dated Jul. 15, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Official Action Dated Jul. 15, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated Mar. 15, 2004 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/765,316.
Official Action Dated Mar. 15, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/343,792.
Official Action Dated Sep. 15, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Official Action Dated Sep. 15, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Official Action Dated Dec. 16, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/343,792.
Official Action Dated Sep. 16, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/727,464.
Official Action Dated Jan. 17, 2006 From the United States Patent and Trademark Office Re.: U.S. Appl. No. 11/034,007.
Official Action Dated Jul. 19, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/607,075.
Official Action Dated Jul. 19, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Official Action Dated Mar. 19, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/240,239.
Official Action Dated Apr. 20, 2006 From the United States Patent and Trademark Office Re.: U.S. Appl. No. 10/240,239.
Official Action Dated Apr. 20, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Official Action Dated Jul. 20, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,617.
Official Action Dated Mar. 21, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/533,568.
Official Action Dated Sep. 21, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Official Action Dated Dec. 23, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/727,464.
Official Action Dated Feb. 23, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Official Action Dated Jun. 23, 2006 From the United States Patent and Trademark Office Re.: U.S. Appl. No. 09/727,464.
Official Action Dated May 23, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/667,793.
Official Action Dated May 23, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Official Action Dated May 23, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Official Action Dated Nov. 23, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Official Action Dated Jun. 25, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated Sep. 25, 2006 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated Nov. 26, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/240,239.
Official Action Dated Oct. 26, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Official Action Dated Apr. 27, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Official Action Dated Jul. 27, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Official Action Dated Oct. 27, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Official Action Dated Apr. 28, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated Aug. 28, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/240,239.
Official Action Dated Dec. 28, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/607,075.
Official Action Dated Jan. 28, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Official Action Dated Jun. 28, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/628,074.
Official Action Dated Apr. 29, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Official Action Dated Oct. 30, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Official Action Dated Sep. 30, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Official Action Dated Sep. 30, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Official Action Dated Jan. 31, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/667,793.
Response Dated Jul. 1, 2010 to Official Action of Mar. 2, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Response Dated Jul. 1, 2010 to Official Action of Mar. 2, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,617.
Response Dated Jun. 1, 2010 to Official Action of Mar. 1, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/794,799.
Response Dated Sep. 1, 2010 to Official Action of Aug. 3, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Response Dated Sep. 1, 2011 to Communication Pursuant to Article 94(3) EPC of Mar. 2, 2011 From the European Patent Office Re.: Application No. 06756259.5.
Response Dated Jun. 3, 2010 to Notice of Appeal and Pre-Appeal Brief of Jan. 4, 2010 to Official Action of Sep. 2, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/343,792.
Response Dated Mar. 3, 2011 to Notice of Non-Compliant Amendment of Feb. 14, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Response Dated Apr. 5, 2011 to Official Action of Nov. 10, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Response Dated Oct. 5, 2010 to Official Action of Jul. 19, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Response Dated Apr. 7, 2009 to Official Action of Oct. 7, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Response Dated Jun. 7, 2011 to Official Action of Mar. 9, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Response Dated Dec. 8, 2011 to Restriction Official Action of Nov. 8, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/309,479.
Response Dated Jul. 8, 2010 to Communication Pursuant to Article 94(3) EPC of Mar. 8, 2010 From the European Patent Office Re.: Application No. 06832278.3.
Response Dated Jul. 8, 2010 to Official Action of Apr. 9, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Response Dated Mar. 8, 2011 to Official Action of Dec. 8, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Response Dated Sep. 8, 2010 to Communication Pursuant to Article 94(3) EPC Dated May 12, 2010 From the European Patent Office Re.: Application No. 06809851.6.
Response Dated Dec. 10, 2009 to Official Action of Aug. 11, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Response Dated Feb. 10, 2011 to Notice of Allowance of Nov. 15, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 10/616,301.
Response Dated Feb. 10, 2011 to Official Action of Nov. 10, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Response Dated May 10, 2010 to Official Action of Apr. 8, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Response Dated May 10, 2010 to Official Action of Jan. 8, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Response Dated May 11, 2010 to Official Action of Mar. 11, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/607,075.
Response Dated Oct. 12, 2009 to Notice of Allowance of Jul. 16, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 12/084,559.
Response Dated Sep. 12, 2011 to Official Action of Jul. 11, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,683.
Response Dated Mar. 13, 2008 to Official Action of Dec. 13, 2007 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Response Dated Nov. 13, 2011 to Official Action of Sep. 12, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,653.
Response Dated Aug. 14, 2008 to Official Action of Apr. 15, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/727,464.
Response Dated Jan. 14, 2010 to Official Action of Sep. 15, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Response Dated Jan. 14, 2010 to Official Action of Sep. 15, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Response Dated Jul. 14, 2011 to Official Action of Mar. 15, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/343,792.
Response Dated Nov. 14, 2011 to Official Action of Jul. 12, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Response Dated Nov. 14, 2011 to Official Action of Sep. 12, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/087,150.
Response Dated Oct. 14, 2009 to Official Action of May 14, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Response Dated Oct. 14, 2011 to Communication Pursuant to Rules 70(2) and 70a(2) EPC of Apr. 4, 2011 From the European Patent Office Re. Application No. 05803689.8.
Response Dated Oct. 14, 2011 to Supplementary European Search Report and the European Search Opinion of Mar. 16, 2011 From the European Patent Office Re. Application No. 05803689.8.
Response Dated Dec. 15, 2010 to Official Action of Jul. 19, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Response Dated Mar. 15, 2007 to Official Action of Dec. 15, 2006 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Response Dated Aug. 16, 2010 to Communication Pursuant to Article 94(3) EPC of Apr. 16, 2010 From the European Patent Office Re. Application No. 01951883.6.
Response Dated Nov. 18, 2010 to Official Action of Jul. 19, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/607,075.
Response Dated Sep. 20, 2011 to Official Action of Apr. 20, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Response Dated Jan. 21, 2010 to Official Action of Sep. 21, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Response Dated Feb. 22, 2010 to Communication Pursuant to Article 94(3) EPC of Oct. 21, 2009 From the European Patent Office Re.: Application No. 02716285.8.
Response Dated Sep. 22, 2008 to Official Action of Jun. 25, 2008 From US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Response Dated Jun. 23, 2010 to Official Action of Feb. 23, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Response Dated Nov. 23, 2011 to Official Action of May 23, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Response Dated Nov. 23, 2011 to Official Action of May 23, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Response Dated Mar. 24, 2011 to Official Action of Dec. 8, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Response Dated Oct. 24, 2011 to Official Action of May 23, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/667,793.
Response Dated Aug. 25, 2010 to Official Action of Jul. 27, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Response Dated Nov. 25, 2005 to Office Action of May 13, 2005 From the Patent Office of the People's Republic of China Re.: Application No. 1817689.5.
Response Dated Jul. 26, 2010 to Official Action of Apr. 28, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Response Dated May 26, 2010 to Official Action of Mar. 19, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/240,239.
Response Dated Jan. 27, 2011 to Official Action of Nov. 1, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 12/728,383.
Response Dated Dec. 28, 2009 to Official Action of Aug. 28, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/240,239.
Response Dated Jun. 28, 2011 to Official Action of Dec. 28, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/607,075.
Response Dated Nov. 28, 2011 to Official Action of Jun. 28, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/628,074.
Response Dated Aug. 29, 2011 to Official Action of Apr. 27, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/836,223.
Response Dated Dec. 29, 2011 to Office Action of Jul. 17, 2007 From the Israeli Patent Office Re.: Application No. 154323.
Response Dated Dec. 30, 2009 to Official Action of Sep. 1, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/794,799.
Response Dated Dec. 30, 2009 to Official Action of Oct. 30, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Response Dated Jan. 31, 2011 to Official Action of Sep. 30, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/798,017.
Response Dated Mar. 31, 2011 to Official Action of Jan. 31, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/667,793.
Response Dated Oct. 31, 2007 to Official Action of Jul. 12, 2007 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,301.
Response to the International Search Report and the Written Opinion of Oct. 10, 2006 From the International Searching Authority Re.: Appliction No. PCT/IL06/00059.
Restriction Official Action Dated Nov. 8, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/309,479.
Restriction Official Action Dated Nov. 15, 2011 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,683.
Second International Search Report Dated Jun. 1, 2009 From the International Searching Authority Re.: Application No. PCT/IL07/00918.
Summons to Attend Oral Proceedings Pursuant to Rule 115(1) EPC Dated Jan. 16, 2009 From the European Patent Office Re.: Application No. 03810570.6.
Supplemental Notice of Allowability Dated Oct. 24, 2011 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/607,075.
Supplemental Response After Interview Dated Aug. 4, 2010 to Official Action of Mar. 2, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,617.
Supplemental Response Under 37 C.F.R. § 1.125 Dated Aug. 12, 2010 to Telephonic Interview of Aug. 6, 2010 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,617.
Supplementary European Search Report and the European Search Opinion Dated Mar. 16, 2011 From the European Patent Office Re. Application No. 05803689.8.
Supplementary European Search Report Dated Dec. 12, 2005 From the European Patent Office Re.: Application No. 03810570.6.
Supplementary Partial European Search Report and the European Search Opinion Dated Dec. 15, 2009 From the European Patent Office Re.: Application No. 06832278.3.
Supplementary Partial European Search Report and the European Search Opinion Dated Oct. 16, 2009 From the European Patent Office Re.: Application No. 06756259.5.
Supplementary Partial European Search Report Dated Sep. 4, 2007 From the European Patent Office Re.: Application No. 0 2716285.8.
Supplementary Partial European Search Report Dated Nov. 11, 2008 From the European Patent Office Re.: Application No. 01951883.6.
Supplementary Partial European Search Report Dated Nov. 20, 2007 From the European Patent Office Re.: Application No. 02716285.8.
Translation of Office Action Dated May 13, 2005 From the Patent Office of the People's Republic of China Re.: Application No. 01817689.5.
Written Opinion Dated Feb. 1, 2006 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/00048.
Written Opinion Dated Jul. 1, 2008 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL06/00834.
Written Opinion Dated Jul. 2, 2007 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL2006/001291.
Written Opinion Dated Aug. 3, 2006 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/001173.
Written Opinion Dated Oct. 10, 2006 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL06/00059.
Written Opinion Dated Oct. 15, 2008 From the International Searching Authority Re.: Application No. PCT/IL07/00918.
Written Opinion Dated Mar. 23, 2006 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/00572.
Written Opinion Dated May 24, 2007 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/00575.
Written Opinion Dated Jul. 25, 2008 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/001173.
Written Opinion Dated Mar. 26, 2007 From the International Searching Authority of the Patent Cooperation Treaty Re.: Application No. PCT/IL05/00394.
Aoi et al. “Absolute Quantitation of Regional Myocardial Blood Flow of Rats Using Dynamic Pinhole SPECT”, IEEE Nuclear Science Symposium and Medical Imaging Conference Record, 3: 1780-1783, 2002. Abstract, Figs.
Beekman et al. “Efficient Fully 3-D Iterative SPECT Reconstruction With Monte Carlo-Based Scatter Compensation”, IEEE Transactions on Medical Imaging, 21(8): 867-877, Aug. 2002.
Berman et al. “Dual-Isotope Myocardial Perfusion SPECT With Rest Thallium-201 and Stress Tc-99m Sestamibi”, Cardiology Clinics, 12(2): 261-270, May 1994.
Bloch et al. “Application of Computerized Tomography to Radiation Therapy and Surgical Planning”, Proceedings of the IEEE, 71(3): 351-355, Mar. 1983.
Bromiley et al. “Attenuation Correction in PET Using Consistency Conditions and a Three-Dimensional Template”, IEEE Transactions on Nuclear Science, XP002352920, 48(4): 1371-1377, 2001. p. 1376, col. 2, § 2.
Brown et al. “Method for Segmenting Chest CT Image Data Using an Anatomical Model: Preliminary Results”, IEEE Transactions on Medical Imaging, 16(6): 828-839, Dec. 1997.
Chengazi et al. “Imaging Prostate Cancer With Technetium-99m-7E11-C5.3 (CYT-351)”, Journal of Nuclear Medicine, 38: 675-682, 1997.
Corstens et al. “Nuclear Medicine's Role in Infection and Inflammation”, The Lancet, 354: 765-770, 1999.
Day et al. “Localization of Radioiodinated Rat Fibrogen in Transplanted Rat Tumors”, Journal of the National Cancer Institute, 23(4): 799-812, 1959.
DeGrado et al. “Topics in Integrated Systems Physiology. Tracer Kinetic Modeling in Nuclear Cardiology”, Journal of Nuclear Cardiology, 7: 686-700, 2000.
Del Guerra et al. “An Integrated PET-SPECT Small Animal Imager: Preliminary Results”, Nuclear Science Symposium, IEEE Records, 1: 541-544, 1999.
Ellestad “Stress Testing: Principles and Practice”, XP008143015, 5th Edition, p. 432, Jan. 1, 2003.
Erbil et al. “Use and Limitations of Serum Total and Lipid-Bound Sialic Acid Concentrations as Markers for Colorectal Cancer”, Cancer, 55: 404-409, 1985.
Garcia et al. “Accuracy of Dynamic SPECT Acquisition for Tc-99m Teboroxime Myocardial Perfusion Imaging: Preliminary Results”, American College of Cardiology, 51st Annual Scientific Session, Atlanta, Georgia, USA, 8 P., 2002.
Gilland et al. “A 3D Model of Non-Uniform Attenuation and Detector Response for Efficient Iterative Reconstruction in SPECT”, Physics in Medicine and Biology, XP002558623, 39(3): 547-561, Mar. 1994. p. 549-550, Section 2.3 ‘Active Voxel Reconstruction’, p. 551, Lines 4-8.
Gilland et al. “Long Focal Length, Asymmetric Fan Beam Collimation for Transmission Acquisition With a Triple Camera SPECT System”, IEEE Transactions on Nuclear Science, XP011087666, 44(3): 1191-1196, Jun. 1, 1997.
Gilland et al. “Simultaneous Reconstruction and Motion Estimation for Gated Cardiac ECT”, IEEE Transactions on Nuclear Science, XP011077797, 49(5): 2344-2349, Oct. 1, 2002. p. 2344, Section ‘Introduction’, First §.
Gugnin et al “Radiocapsule for Recording the Ionizing Radiation in the Gastrointestinal Tract”, UDC 615. 417:616.34-005.1-073.916-71 (All-Union Scientific-Research Institute of medical Instrument Design, Moscow. Translated from Meditsinskaya Tekhnika, 1:21-25, Jan.-Feb. 1972).
Hassan et al. “A Radiotelemetry Pill for the Measurement of Ionising Radiation Using a Mercuric Iodide Detector”, Physics in Medicine and Biology, 23(2): 302-308, 1978.
Hayakawa et al. “A PET-MRI Registration Technique for PET Studies of the Rat Brain”, Nuclear Medicine & Biology, 27: 121-125, 2000. p. 121, col. 1.
Herrmann et al. “Mitochondrial Proteome: Altered Cytochtrome C Oxidase Subunit Levels in Prostate Cancer”, Proteomics, XP002625778, 3(9): 1801-1810, Sep. 2003.
Hoffman et al. “Intraoperative Probes and Imaging Probes”, European Journal of Nuclear Medicine, 26(8): 913-935, 1999.
Huesman et al. “Kinetic Parameter Estimation From SPECT Cone-Beam Projection Measurements”, Physics in Medicine and Biology, 43(4): 973-982, 1998.
Jeanguillaume et al. “From the Whole-Body Counting to Imaging: The Computer Aided Collimation Gamma Camera Project (CACAO)”, Radiation Projection Dosimetry, 89(3-4): 349-352, 2000. & RSNA 2000 Infosystem, 87th Scientific Assembly and Annual Meeting, Chicago, Illinois, 2000.
Jessup “Tumor Markers—Prognostic and Therapeutic Implications for Colorectal Carcinoma”, Surgical Oncology, 7: 139-151, 1998.
Kadrmas et al. “Static Versus Dynamic Teboroxime Myocardial Perfusion SPECT in Canines”, IEEE Transactions on Nuclear Science, 47(3): 1112-1117, Jun. 2000.
Kinahan et al. “Attenuation Correction for a Combined 3D PET/CT Scanner”, Medical Physics, 25(10): 2046-2053, Oct. 1998.
Kojima et al. “Quantitative Planar Imaging Method for Measurement of Renal Activity by Using a Conjugate-Emission Image and Transmission Data”, Medical Physics, 27(3): 608-615, 2000. p. 608.
Krieg et al. “Mitochondrial Proteome: Cancer-Altered Metabolism Associated With Cytochrome C Oxidase Subunit Level Variation”, Proteomics, XP002625779, 4(9): 2789-2795, Sep. 2004.
Lavallee et al. “Building a Hybrid Patient's Model for Augmented Reality in Surgery: A Registration Problem”, Computing in Biological Medicine, 25(2): 149-164, 1995.
Li et al. “A HOTLink/Networked PC Data Acquisition and Image Reconstruction System for a High Resolution Whole-Body PET With Respiratory or ECG-Gated Performance”, IEEE Nuclear Sience Symposium and Medical Imaging Conference, Norfolk, VA, USA, Nov. 10-16, 2002, XP010663724, 2: 1135-1139, Nov. 10, 2002. p. 1137, First Col., 2nd §.
Lin et al. “Improved Sensor Pills for Physiological Monitoring”, NASA Technical Brief, JPL New Technology Report, NPO-20652, 25(2), 2000.
Links “Advances in SPECT and PET Imaging”, Annals in Nuclear Medical Science, 13(2): 107-120, Jun. 2000.
Mao et al. “Human Prostatic Carcinoma: An Electron Microscope Study”, Cancer Research, XP002625777, 26(5): 955-973, May 1966.
McJilton et al. “Protein Kinase C? Interacts With Bax and Promotes Survival of Human Prostate Cancer Cells”, Oncogene, 22; 7958-7968, 2003.
Mettler et al. “Legal Requirements and Radiation Safety”, Essentials of Nuclear Medicine Imaging, 2nd Ed., Chap.13: 323-331, 1985.
Meyers et al. “Age, Perfusion Test Results and Dipyridamole Reaction”, Radiologic Technology, XP008142909, 73(5): 409-414, May 1, 2002.
Molinolo et al. “Enhanced Tumor Binding Using Immunohistochemical Analyses by Second Generation Anti-Tumor-Associated Glycoprotein 72 Monoclonal Antibodies versus Monoclonal Antibody B72.3 in Human Tissue”, Cancer Research, 50: 1291-1298, 1990.
Moore et al. “Quantitative Multi-Detector Emission Computerized Tomography Using Iterative Attenuation Compensation”, Journal of Nuclear Medicine, XP002549083, 23(8): 706-714, Aug. 1982. Abstract, p. 707, Section ‘The Multi-Detector Scanner’, First §.
Mori et al. “Overexpression of Matrix Metalloproteinase-7mRNA in Human Colon Carcinomas”, Cancer, 75: 1516-1519, 1995.
Ogawa et al. “Ultra High Resoultion Pinhole SPECT”, IEEE Nuclear Science Symposium, 2: 1600-1604, 1998.
Pardridge et al. “Tracer Kinetic Model of Blood-Brain Barrier Transport of Plasma Protein-Bound Ligands”, Journal of Clinical Investigation, 74: 745-752, 1984. Suppl. IDS in 27480.
Pellegrini et al. “Design of Compact Pinhole SPECT System Based on Flat Panel PMT”, IEEE Nuclear Science Symposium Conference Record, 3: 1828-1832, 2003.
Piperno et al. “Breast Cancer Screening by Impedance Measurements”, Frontiers Med. Biol. Engng., 2(2): 11-17, 1990.
Pluim et al. “Image Registration by Maximization of Combined Mutual Information and Gradient Information”, IEEE Transactions on Medical Imaging, 19(8): 1-6, 2000.
Qi et al. “Resolution and Noise Properties of MAP Reconstruction for Fully 3-D PET”, IEEE Transactions on Medical Imaging, XP002549082, 19(5): 493-506, May 2000. p. 493, col. 2, Lines 10-21, p. 495, col. 1, Last §.
Quartuccia et al. “Computer Assisted Collimation Gama Camera: A New Approach to Imaging Contaminated Tissues”, Radiation Projection Dosimetry, 89(3-4): 343-348, 2000.
Rajshekhar “Continuous Impedence Monitoring During CT-Guided Stereotactic Surgery: Relative Value in Cystic and Solid Lesions”, British Journal of Neurosurgery, 6: 439-444, 1992.
Reutter et al. “Direct Least Squares Estimation of Spatiotemporal Distributions From Dynamic SPECT Projections Using a Spatial Segmentation and Temporal B-Splines”, IEEE Transactions on Medical Imaging, 19(5): 434-450, 2000.
Reutter et al. “Kinetic Parameter Estimation From Attenuated SPECT Projection Measurements”, IEEE Transactions on Nuclear Science, 45(6): 3007-3013, 1998.
Stoddart et al. “New Multi-Dimensional Reconstructions for the 12-Detector, Scanned Focal Point, Single-Photon Tomograph”, Physics in Medicine and Biology, XP020021960, 37(3): 579-586, Mar. 1, 1992. p. 582, § 2-p. 585, § 1.
Storey et al. “Tc-99m Sestamibi Uptake in Metastatic Prostate Carcinoma”, Clinical Nuclear Medicine, XP009145398, 25(2): 133-134, Feb. 2000.
Takahashi et al. “Attenuation Correction of Myocardial SPECT Images With X-Ray CT: Effects of Registration Errors Between X-Ray CT and SPECT”, Annals of Nuclear Medicine, 16(6): 431-435, Sep. 2002.
Wilson et al. “Non-Stationary Noise Characteristics for SPECT Images”, Proceedings of the Nuclear Science Symposium and Medical Imaging Conference, Santa Fe, CA, USA, Nov. 2-9, 1991, XP010058168, p. 1736-1740, Nov. 2, 1991. p. 1736, col. 2, Lines 4-6.
Wu et al. “ECG-Gated Pinhole SPECT in Mice With Millimeter Spatial Resolution”, IEEE Transactions on Nuclear Science, 47(3): 1218-1221, Jun. 2000.
Xu et al. “Quantitative Expression Profile of Androgen-Regulated Genes in Prostate Cancer Cells and Identification of Prostate-Specific Genes”, International Journal of Cancer, 92: 322-328, 2001.
Yu et al. “Using Correlated CT Images in Compensation for Attenuation in PET Image Reconstruction”, Proceedings of the SPIE, Applications of Optical Engineering: Proceedings of OE/Midwest '90, 1396: 56-58, 1991.
Zaidi et al. “Magenetic Resonance Imaging-Guided Attenuation and Scatter Corrections in Three-Dimensional Brain Positron Emission Tomography”, Medical Physics, 30(5): 937-948, May 2003.
Zaidi et al. “MRI-Guided Attenuation Correction in 3D Brain PET”, Neuroimage Human Brain Mapping 2002 Meeting, 16(2): Abstract 504, Jun. 2002.
Zhang et al. “An Innovative High Efficiency and High Resolution Probe for Prostate Imaging”, The Journal of Nuclear Medicine, 68: 18, 2000. Abstract.
Zhang et al. “Potential of a Compton Camera for High Performance Scintimammography”, Physics in Medicine and Biology, XP020024019, 49(4): 617-638, Feb. 21, 2004.
Notice of Allowance Dated Feb. 25, 2013 From the US Patent and Trademark Office Re.: U.S. Appl. No. 09/641,973.
Communication Under Rule 71(3) EPC Dated Feb. 26, 2013 From the European Patent Office Re. Application No. 06756259.5.
Official Action Dated Feb. 22, 2013 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/616,307.
Restriction Official Action Dated Apr. 13, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/989,223.
Applicant-Initiated Interview Summary Dated May 1, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 10/343,792.
Applicant-Initiated Interview Summary Dated May 1, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,653.
Advisory Action before the Filing of an Appeal Brief Dated May 21, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/980,653.
Applicant-Initiated Interview Summary Dated May 9, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/448,473.
Communication Pursuant to Article 94(3) EPC Dated May 29, 2012 From the European Patent Office Re. Application No. 05803689.8.
Communication Under Rule 71(3) EPC Dated May 30, 2012 From the European Patent Office Re.: Application No. 02716285.8.
Dillman “Radiolabeled Anti-CD20 Monoclonal Antibodies for the Treatment of B-Cell Lymphoma”, Journal of Clinical Oncology, 20(16): 3545-3557, Aug. 15, 2002.
Sands et al. “Methods for the Study of the Metabolism of Radiolabeled Monoclonal Antibodies by Liver and Tumor”, The Journal of Nuclear Medicine, 28: 390-398, 1987.
Advisory Action Before the Filing of an Appeal Brief Dated Jul. 12, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/667,793.
Official Action Dated May 18, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/514,785.
Official Action Dated Apr. 19, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/750,057.
Official Action Dated Apr. 23, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/932,987.
Notice of Allowance Dated Jul. 15, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/932,987.
Official Action Dated Jul. 5, 2013 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/656,548.
Official Action Dated Jul. 30, 2012 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,683.
Bowsher et al. “Treatment of Compton Scattering in Maximum-Likelihood, Expectation-Maximization Reconstructions of SPECT Images”, Journal of Nuclear Medicine, 32(6): 1285-1291, 1991.
Official Action Dated Aug. 2, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/087,150.
Handrick et al. “Evaluation of Binning Strategies for Tissue Classification in Computed Tomography Images”, Medical Imaging 2006: Image Processing, Proceedings of the SPIE, 6144: 1476-1486, 2006.
Thorndyke et al. “Reducing Respiratory Motion Artifacts in Positron Emission Tomography Through Retrospective Stacking”, Medical Physics, 33(7): 2632-2641, Jul. 2006.
Notice of Allowance Dated Dec. 26, 2012 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/980,690.
Official Action Dated Dec. 28, 2012 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/343,792.
Charland et al. “The Use of Deconvolution and Total Least Squares in Recovering a Radiation Detector Line Spread Function”, Medical Physics, 25(2): 152-160, Feb. 1998. Abstract Only!
Official Action Dated Jan. 19, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/667,793.
Official Action Dated Jan. 23, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/087,150.
Jin et al. “Reconstruction of Cardiac-Gated Dynamic SPECT Images”, IEEE International Conference on Image Processing 2005, ICIP 2005, Sep. 11-14, 2005, 3: 1-4, 2005.
Toennies et al. “Scatter Segmentation in Dynamic SPECT Images Using Principal Component Analysis”, Progress in Biomedical Optics and Imaging, 4(23): 507-516, 2003.
Office Action Dated Jul. 17, 2007 From the Israeli Patent Office Re.: Application No. 154323 and Its Translation Into English.
Communication Pursuant to Article 94(3) EPC Dated Jun. 11, 2012 From the European Patent Office Re.: Application No. 06756259.5.
Official Action Dated Jul. 30, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 10/343,792.
Official Action Dated Jul. 31, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/667,793.
Supplementary European Search Report and the European Search Opinion Dated Nov. 13, 2012 From the European Patent Office Re. Application No. 06728347.3.
Notice of Allowance Dated Oct. 26, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/514,785.
Action Dated Nov. 30, 2012 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/514,785.
Official Action Dated Feb. 10, 2014 From the US Patent and Trademark Office Re. U.S. Appl. No. 11/667,793.
Official Action Dated Apr. 11, 2014 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/309,479.
Sharir et al. “D-SPECT: High Speed Myocardial Perfusion Imaging: A Comparison With Dual Detector Anger Camera (A-SPECT)”, The Journal of Nuclear Medicine, 48(Suppl.2): 51P, # 169, 2007.
Communication Pursuant to Article 94(3) EPC Dated Nov. 25, 2013 From the European Patent Office Re. Application No. 06756258.7.
Notice of Allowance Dated Dec. 17, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 13/913,804.
Official Action Dated Dec. 16, 2013 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/087,150.
Notice of Allowance Dated Feb. 4, 2014 From the US Patent and Trademark Office Re. U.S. Appl. No. 13/947,198.
Applicant-Initiated Interview Summary Dated Mar. 20, 2014 From the US Patent and Trademark Office Re. U.S. Appl. No. 12/087,150.
Notice of Allowance Dated Apr. 2, 2014 From the US Patent and Trademark Office Re. U.S. Appl. No. 10/343,792.
Notice of Allowance Dated Aug. 20, 2013 From Re. U.S. Appl. No. 11/932,872.
Applicant-Interview Summary Dated Mar. 21, 2012 From Re. U.S. Appl. No. 12/514,785.
Communication Pursuant to Article 94(3) EPC Dated May 8, 2014 From the European Patent Office Re. Application No. 05803689.8.
Official Action Dated May 13, 2014 From Re. U.S. Appl. No. 12/448,473.
Official Action Dated Jun. 17, 2014 From Re. U.S. Appl. No. 12/087,150.
Related Publications (1)
Number Date Country
20120106820 A1 May 2012 US
Provisional Applications (18)
Number Date Country
60816970 Jun 2006 US
60800846 May 2006 US
60800845 May 2006 US
60799688 May 2006 US
60763458 Jan 2006 US
60750597 Dec 2005 US
60750334 Dec 2005 US
60750287 Dec 2005 US
60741440 Dec 2005 US
60720652 Sep 2005 US
60720541 Sep 2005 US
60720034 Sep 2005 US
60702979 Jul 2005 US
60700753 Jul 2005 US
60700752 Jul 2005 US
60700318 Jul 2005 US
60700299 Jul 2005 US
60700317 Jul 2005 US
Continuations (1)
Number Date Country
Parent 11988926 US
Child 13345773 US
Continuation in Parts (4)
Number Date Country
Parent PCT/IL2006/000562 May 2006 US
Child 11988926 US
Parent PCT/IL2006/000059 Jan 2006 US
Child PCT/IL2006/000562 US
Parent PCT/IL2005/001215 Nov 2005 US
Child PCT/IL2006/000059 US
Parent PCT/IL2005/001173 Nov 2005 US
Child PCT/IL2005/001215 US