The Generalized NEQ and Detectability Index for Tomosynthesis and Cone-Beam CT: From Cascaded Systems Analysis to Human Observers M5G2M9;

Size: px
Start display at page:

Download "The Generalized NEQ and Detectability Index for Tomosynthesis and Cone-Beam CT: From Cascaded Systems Analysis to Human Observers M5G2M9;"

Transcription

1 The Generalize NEQ an Detectability Inex for Tomosynthesis an Cone-Beam CT: From Cascae Systems Analysis to Human Observers G. J. Gang, a J. Lee, b J. W. Stayman, c D. J. Twar, c W. Zbijewski, c J. L. Prince, b an J. H. Siewersen a,c a Institute of Biomaterials an Biomeical Engineering, University of Toronto, Toronto ON, Canaa M5GM9; b Electrical an Computer Engineering, Johns Hopkins University, Baltimore MD, USA 118 c Dept. of Biomeical Engineering, Johns Hopkins University, Baltimore MD, USA 105 ABSTRACT Purpose: In the early evelopment of new imaging moalities such as tomosynthesis an cone-beam CT (CBCT) an accurate preictive moel for imaging performance is particularly valuable in ientifying the physical factors that govern image quality an guiing system optimization. In this work, a task-base cascae systems moel for etectability inex is propose that escribes not only the signal an noise propagation in the D (projection) an 3D (reconstruction) imaging chain but also the influence of backgroun anatomical noise. The extent to which generalize etectability inex provies a vali metric for imaging performance was assesse through irect comparison to human observer experiments. Methos: Detectability inex ( ) was generalize to inclue anatomical backgroun noise in the same manner as the generalize noise-equivalent quanta (NEQ) propose by Barrett et al. (Proc. SPIE Me. Imaging, Vol. 1090, 1989). Anatomical backgroun noise was measure from a custom phantom esigne to present power-law spectral ensity comparable to various anatomical sites (e.g., breast an lung). Theoretical calculations of as a function of the sourceetector orbital extent (θ tot ) was obtaine from a 3D cascae systems analysis moel for tomosynthesis an cone-beam CT (CBCT). Four moel observers were consiere in the calculation of : prewhitening (PW), non-prewhitening (NPW), prewhitening with eye filter an internal noise (PWE), an non-prewhitening with eye filter an internal noise (NPWE). Human observer performance was measure from 9AFC tests for a variety of iealize imaging tasks presente within a clutter phantom. Theoretical results ( ) were converte to area uner the ROC curve (A z ) an compare irectly to human observer performance as a function of imaging task an orbital extent. Results: Theoretical results emonstrate reasonable corresponence with human observer response for all tasks across the continuum in θ tot ranging from low-angle tomosynthesis (θ tot ~10 o ) to CBCT (θ tot ~180 o ). Both theoretical an experimental A z were foun to increase with acquisition angle, consistent with increase rejection of out-of-plane clutter for larger tomosynthesis angle. Of the four theoretical moel observers consiere, the prewhitening moels tene to overestimate real observer performance, while the non-prewhitening moels emonstrate reasonable agreement. Conclusions: Generalize etectability inex was shown to provie a meaningful metric for imaging performance, helping to brige the gap between real observer performance an prevalent Fourier-base metrics base in first principles of spatial-frequency-epenent NEQ an imaging task. Keywors: cone-beam CT, anatomical noise, anatomical clutter, cascae systems analysis, noise-power spectrum, noise-equivalent quanta, etectability inex 1. INTRODUCTION Recent years have witnesse a ramatic avancement in the scope of avance 3D imaging applications base on flatpanel etectors, such as tomosynthesis an cone-beam CT (CBCT). Each of these approaches aims to improve the conspicuity of subtle lesions through epth localization an the rejection of out-of-plane anatomical clutter, with a broa spectrum of potential applications ranging from breast an chest imaging to image-guie interventions. The knowlegeable evelopment of such moalities stans to benefit from a general quantitative framework for performance Meical Imaging 010: Physics of Meical Imaging, eite by Ehsan Samei, Norbert J. Pelc, Proc. of SPIE Vol. 76, 760Y 010 SPIE CCC coe: /10/$18 oi: / Proc. of SPIE Vol Y-1

2 moeling, assessment, an system optimization. Two common approaches to image quality assessment inclue physical characterization of system performance [e.g., prevalent Fourier-base metrics of spatial-frequency-epenent signal an noise] an human observer-base characterization [e.g., sensitivity an specificity as assesse in a receiver operating characteristic (ROC) test]. A variety of avantages an limitations can be appreciate for each approach for example, the former is simple an practical, yet subject to numerous assumptions of linearity an shift-invariance, while the latter irectly assesses real observer response, but is time consuming, expensive, an epenent on the generation of real (or highly realistic) images. A general quantitative approach that begins with simple Fourier metrics an briges the gap to real observer performance woul be valuable in the esign an optimization of such imaging systems. Previous work has extene cascae systems analysis of the noise-equivalent quanta (NEQ) for projection raiography 1-4 to avance D an 3D imaging moalities 5-8. Recognizing the importance of imaging task in efining meaningful image quality metrics, a task-riven approach is aopte in which NEQ is combine with the imaging task an moel observers to form the etectability inex ( ). The goal of the current stuy is to investigate the extent to which etectability inex is a vali metric for system performance an optimization through comparison to real human observers. The work below represents preliminary stuies spanning initial analysis an experimentation involving a small number of simple imaging tasks, with a more complete evaluation, incluing a greater variety of imaging conitions an imaging tasks, to be pursue in future work.. GENERALIZED DETECTABILITY INDEX.1 Cascae Systems Analysis (CSA) Moel for Tomosynthesis an CBCT Previous work has emonstrate the utility an accuracy of cascae systems analysis of the performance of D x-ray projection imaging systems base on flat-panel etectors. 3,4 Such has more recently been extene to 3D FPD-base CBCT an tomosynthesis systems, 5-8 proviing a general framework for moeling the 3D moulation transfer function (MTF), noise-power spectrum (NPS), NEQ, an etective quantum efficiency (DQE). The basic approach provies quantitative moeling of the signal an noise propagation in the D an 3D imaging chain as illustrate in Fig. 1. The moel consists of the physical processes of D projection image formation followe by the mathematical processes of 3D filtere backprojection reconstruction. The former has been extensively investigate an evelope over the last 15 years. 1-4,9 The reconstruction cascae inclues 7,8 : log-normalization of projection ata (Stage 8); application of the ramp filter (Stage 9); optional application of apoization filters(stage 10); interpolation of the D projection ata (Stage 11); backprojection of D projections, moele as superposition of the D NPS along vanes in the 3D Fourier omain (Stage 1); an sampling of the 3D reconstruction (Stag 13). Such a moel has emonstrate agreement with measurements of the 3D NPS an NEQ for FPD-CBCT over a broa range of experimental conitions an reconstruction techniques 7,8, proviing a moel for imaging performance of tomosynthesis an CBCT as a contimuum represente by varying extent of acquisition angle (θ tot ). Figure. 1. Cascae systems analysis for 3D tomosynthesis an CBCT. Stages 0-7 moel signal an noise transfer through various physical processes in the etector, while Stages 8 to 13 represent mathematical processes of 3D reconstruction by filtere backprojection as etaile in previous work: 0.) incient x-ray quanta; 1.) interaction of x-rays Proc. of SPIE Vol Y-

3 in etector;.) generation of seconary quanta; 3.) sprea of seconary quanta; 4.) coupling of quanta to etector apertures; 5.) integration by etector aperture; 6.) sampling of etector pixels; 7.) reaout with aitive noise; 8.) log normalization; 9.) ramp filter; 10.) apoization filter; 11.) interpolation; 1.) 3D backprojection; an 13.) 3D sampling.. Anatomical Backgroun: The Generalize NEQ an Detectability Inex Such moeling has been shown useful in characterizing the physical performance of systems uner the influence of quantum an electronics noise. However, there is consierable evience that image noise (generally interprete as any ranom or semi-ranom component in an image that impees conspicuity) is often ominate by image fluctuations associate not with the Poisson-istribute incient quanta, but rather those associate with the superposition of anatomical clutter in the D image often referre to as anatomical noise or backgroun noise Similarly in tomosynthesis, epening on the acquisition angle, θ tot, out-of-plane structures are superimpose to varying egrees within a slice of interest, with greater orbital extent reucing the contribution of in-plane anatomical noise. In this way, anatomical noise is expecte to be a ominant source of image egraation at low tomosynthesis angles an to be progressively reuce for a larger arc. Potential traeoffs between quantum noise an anatomical noise are immeiately appreciate: For fixe ose per projection: - a smaller number of projections over a shorter arc will increase quantum noise an anatomical noise - a larger number of projections over a shorter arc will reuce quantum noise, reuce view aliasing noise, an increase anatomical noise - a smaller number of projections over a wier arc will increase quantum noise, increase view aliasing noise, an reuce anatomical noise - a larger number of projections over a wier arc will reuce quantum noise an increase anatomical noise For fixe total ose: - a smaller number of projections over a shorter arc will reuce quantum noise an increase anatomical noise - a larger number of projections over a shorter arc will reuce quantum noise, reuce view aliasing noise, an increase anatomical noise - a smaller number of projections over a wier arc will reuce quantum noise, increase view aliasing noise, an reuce anatomical noise - a larger number of projections over a wier arc will increase quantum noise an reuce anatomical noise Clearly, a moel that quantitatively escribes such traeoffs woul be beneficial in unerstaning the factors that limits imaging performance uner various system configurations. Such traeoffs can be quantitatively appreciate through incorporation of anatomical noise (the so-calle anatomical power spectrum, S B ) in the enominator of the NEQ. Barrett et al. 14 terme this a form of generalize NEQ, written for the 3D case as: MTF GNEQ = S + S + S B Q E (1) where S B can take whatever form is appropriate to the anatomical clutter an is often empirically characterize accoring to a power-law form: SB ( f ) κ = () f β where κ escribe the magnitue an β the frequency content of the anatomical clutter. Moeling such image fluctuations in this manner implies treatment of anatomical clutter as a stochastic noise source, which clearly implies numerous assumptions an limitations. Power-law moeling of anatomical noise has shown reasonable agreement with human observer performance measure with real patient backgroun in breast imaging, 15 proviing a convenient basis for Proc. of SPIE Vol Y-3

4 incorporation of a lumpy backgroun in the generalize NEQ. Such approach has been applie extensively in mammography 1,16-0, chest imaging 1-3, an more recently in tomosynthesis an CBCT 4,5. Just as NEQ may be combine with an iealize task function as escribe in ICRU 54 to yiel the etectability inex: 6 MTF (3) ' = ΔH f f f x y z S + S Q E the generalize NEQ can be similarly extene to yiel a generalize etectability inex: MTF (4) ' = ΔH f f f x y z MTF S + S + S B Q E This iealize form for etectability in the presence of anatomical backgroun is equivalent to the prewhitening (PW) observer moel. It has been shown in previous work 1,7-9 to give a useful metric for imaging performance an system optimization. Eq. (4) escribes a fully 3D etectability inex in which the observer can make full use of the volumetric ata. Such may or may not be a reasonable moel for a real observer scrolling slices; rather, it assumes the observer can instantly appreciate the fully 3D statistics of the volume image. The moel can be reuce to escribe the etectability within a ' single slice of 3D image, enote. Since extraction of a single slice from the 3D spatial omain correspons to a slice convolution over the orthogonal irection in the Fourier omain (i.e., convolution with a sinc function corresponing to the slice thickness), a simple form of slice etectability for a prewhitening observer can be obtaine by integrating the numerator an enominator of the 3D etectability inex as: [ T W f ] ' Task y = f f slice x z T S + S + S f B Q E y (5).3 Moel Observers The simple PW moel of Eq. (4) can be extene to slightly more sophisticate moels that have been shown to better approximate real observers uner various conitions. For example, the non-prewhitening (NPW) observer moel is: ' = [ T W ] f f f Task x y z + + B Q E Task x y z ( T S S S ) [ T W ] f f f Similarly, each moel can be extene to inclue a simple escription of a human eye filter an internal noise, referre to respectively as the prewhitening-eye (PWE) an non-prewhitening-eye (NPWE) moels: (6) ' = E [ T W ] Task B Q E i ( ) E T S S S N f f f x y z (7) ' = E [ T W ] f f f Task x y z E ( T S + S + S ) [ T W ] + N f f f 4 B Q E Task i x y z The eye filter an internal noise aopte in this work are consistent previous work by Burgess 30 : (8) n E( f) = f exp( cf ) (9) Proc. of SPIE Vol Y-4

5 where n = 1.3 an c = 3 gives an eye filter, E(f), that peaks at 4 cycles/eg for a viewing istance of 50 cm. Internal noise is assume to be uncorrelate an is set to a fraction (0.0) of the DC component of the noise-power spectrum at a viewing istance of 100 cm: D N = 0.0 NPS(0,0) i 100 where the viewing istance, D, is set to 50 cm. (10).4 Relating to Human Observer Performance Equations (4) (8) above yiel the etectability inex for a given task an observer moel, with the funamental signal an noise properties escribe by a cascae systems moel for the imaging chain. These forms are almost purely theoretical, with the main semi-empirical components of the moel being: the incient x-ray spectrum; the Poisson excess (Swank noise) of the x-ray converter; an the eye filter parameters (taken from Ref. 30). Relation of erive theoretically from cascae systems analysis calculations can be converte to area-uner the ROC curve (A z ) by the relation: 1 ' A = (1 + erf ( )) (11) Z uner assumption that istribution of cases uner two hypotheses are Gaussian an homosceastic IMAGING TASK 3.1 Task Phantoms To investigate task functions presenting stimuli in a power-law backgroun (anatomical noise), a custom phantom was evelope specifically to present power-law noise. The phantom esign principles an associate backgroun power spectra are etaile in Ref. 3. To summarize, the phantom consists of a ranom assortment of acrylic spheres of iameter ranging mm, with proportions of each sphere selecte to give a power-law backgroun NPS with β ~ 3. For the small-sphere etection task, an assortment of six 3.18 mm iameter Teflon spheres was place within the central coronal (x-z) slice of the clutter phantom as shown in Fig. 3(b). For the large sphere etection task, six 10.3 mm polypropylene spheres were use as stimuli. For the encapsulate-sphere iscrimination task, each stimulus consists of a 10 mm acrylic sphere coate with an 3.18mm outer shell of paraffin wax. Example images of each task phantom are illustrate in Fig. 3(b). Images of each task phantom were acquire over a range in θ tot as escribe above, an 3D images were reconstructe using the Felkamp algorithm for 3D filtere backprojection. A central coronal slice was extracte from the volume image through the central plane of each sphere, giving the slice images use in human observer tests escribe below. The six stimuli give six inepenent trials for each task. 3. Imaging Tasks A brige between simple Fourier metrics of imaging performance (e.g., MTF, NPS, an NEQ) an task-epenent performance metrics can be obtaine through a quantitative escription of the spatial frequencies of interest in performing a given task i.e., through a spatial-frequency-epenent task function or template. A simple form for such task functions is given by the Fourier transform of the ifference in two hypotheses H 1 (x, y, z) an H (x, y, z) where H 1 an H represent the spatial omain representations of, for example, the signal-present an noise-only hypotheses. For example, etection of a given stimulus [H 1 (x,y,z)] against a uniform backgroun [H (x,y,z) = constant] is represente by the task function: WTask =Δμ FT[ H ] (1) 1 where FT[H 1 ] is the Fourier transform of the stimulus, an Δμ are the ifference in attenuation coefficients of the stimulus an backgroun, respectively. Proc. of SPIE Vol Y-5

6 For etection of a given stimulus against the power-law backgroun escribe previously, images presenting H1 contains a stimulus (e.g., a sphere of a ifferent contrast) at the center of the ROI, whereas images presenting H have an acrylic sphere of the same size place at the center. Thus, the two hypotheses can be written as: FT [ H1 ] = μ1 FT [ H1 ] + clutter (13a) FT [ H ] = μ FT [ H1 ] + clutter (13b) Therefore, the task function is given by: WTask = Δμ FT [ H1 ] (13c) Finally, iscrimination of two structures the first ientifie as the stimulus an the secon the normal case is given by: WTask = μ1 FT [ H1 ] μ FT [ H ] (14) where H1 an H represent the two structures to be iscriminate. These three simple task functions represent ifferent regions of the spatial frequency omain an therefore imply ifferent susceptibility to various noise sources. High frequency tasks, such as iscrimination tasks or etection of a small object, are prone to egraation by high frequency quantum an electronics noise, an are therefore more epenent on factors of etector esign an reconstruction technique. Low frequency tasks, such as etection of a large, low-contrast lesion, are more epenent on low-frequency noise characteristics, such as anatomical backgroun. Aitional task functions representing higher-orer tasks an representing a varie richness of spatial frequency content are uner investigation in ongoing work. Example task functions are shown in Fig., along with real images acquire as a function of θtot. Figure. Hypotheses an task functions. Task 1: Detection of a small sphere in clutter. Task : Detection of a large sphere in clutter. Task 3. Discrimination of an encapsulate sphere vs. a uniform sphere. 3.3 Image Acquisition An experimental imaging bench for tomosynthesis an cone-beam CT is illustrate in Fig. 3(a) an escribe in previous work7,3. The main components of the bench inclue: an x-ray tube (Ra 94 in a Sapphire housing; Varian Meical Proc. of SPIE Vol Y-6

7 Systems, Salt Lake City, UT); a flat-panel etector (RID-1640A, Perkin Elmer Optoelectronics, Santa Clara CA), an a motion control system in which a phantom place at isocenter is rotate through an arc corresponing to θtot (6K series translation stages, Parker Daeal, Harrison PA, an Dynaserv rotation motor, Parker Hannifin, Rohnert Park, CA). Images were acquire over a range of θtot spanning a continuum from low-angle tomosynthesis (θtot = 10o) to fully 3D CBCT (180o+fan an 360o). In the experiments reporte below, the angular increment (Δθ) was hel fixe at 0.45o, an the number of projections (Nproj) was varie such that θtot range 10o to 360o. In all cases, the imaging technique was 10 kvp (1.53 mm Al mm Cu ae filtration), an 0.63 mas per projection. The range in Nproj was 3 (at 10o) to 800 (at 360o), corresponing to a range in ose of mgy. Figure 3. (a) Experimental bench showing components an coorinate systems for CBCT an tomosynthesis. (b) Experimental phantom presenting power-law backgroun esigne using principles of fractal self-similarity. Six stimuli were inserte to the central coronal slice. 4. EXPERIMENTAL VALIDATION: HUMAN OBSERVER STUDY Human observer tests were conucte to measure task performance as a function of θtot for each of the imaging tasks escribe above. An efficient test well suite to simple phantom stuies is offere by the multiple-alternative forcechoice (MAFC) methoology. In the stuies escribe below, 9AFC tests were use to measure the proportion correct (Pcorr) which may in turn be relate to as: Pcorr ( ', M ) = 1 π ( x ' ) M 1 [φ ( x) ] x exp (15) where M is the number of alternatives, an Φ is the cumulative Gaussian istribution. Note that for a AFC test, Pcorr = Az uner the assumption of an ieal test in which observer response oes not vary over the course of the test. For the 9AFC stuy, Az corresponing to experimental value of Pcorr is convert to Az via a Pcorr - - Az lookup table given by Eq. (11) an (15). Proc. of SPIE Vol Y-7

8 Figure 4. (a) Observer performing a 9AFC test in a arkene reaing room. (b) Nine images were shown on a 3x3 gri as isplaye by Matlab-base OPTEx software in a ranomize orer. The example shows the encapsulate sphere task, with the stimulus highlighte in the lower-left corner. All tests were performe in a controlle image reaing environment with subue lighting an on a iagnostic quality monochrome isplay (Fig. 4). Reaing orer was ranomize (using the Matlab-base OPTEx utility for human observer stuies), an each test was precee by a training set compose of images rawn from the same task phantom image as a function of θtot with images istinct from those in the actual test. For these simple phantom stuies, physicists / engineers were sufficiently expert to perform each imaging task following the training set. A total of 6 observers were use, giving 5 inepenent trials for each task at each setting (angle, case), an a total time of ~60 min require for each test. Responses from each observer for each setting of a task are assume to be inepenent, i.e. no inter- an intraobserver variability were assume. Mean proportion correct (Pcorr) was therefore calculate from the number of correct responses out of a total of 30 (=5 trials x 6 observers), an converte to Az from the lookup table as mentione previously. A 95% confience interval was calculate for Az an shown in the results below. 5. BRIDGING THE GAP: THEORETICAL AND EXPERIMENTAL OBSERVER PERFORMANCE Theoretical calculations of Az for all four observer moels are plotte in Figure 5 in comparison to Az measure from real observer stuies Pcorr. In each case, an upwar tren with θtot is observe for both theoretical an experimental results, consistent with improve rejection of clutter observe in Figure. Overall, the prewhitening moels (PW an PWE) are seen to overestimate real observer performance in each case, whereas the non-prewhitening (NPW an NPWE) observers show reasonable agreement with experimental results. The large sphere etection (iscrimination of polyethylene an acrylic spheres) task appears to be most ifficult, with observer performance challenge up to θtot ~10o, an with the NPW an NPWE moels bouning the human observer measurements. The small sphere etection an encapsulate sphere iscrimination tasks are a bit more conspicuous, each challenging observers up to θtot ~70-100o an escribe reasonably well by the NPW or NPWE moels. Overall, agreement of theory an measurement is fairly goo, proviing evience that simple Fourier metrics such as NEQ combine with imaging task can offer meaningful escriptions of real observer performance. Proc. of SPIE Vol Y-8

9 Figure 5. Theoretical an experimental imaging performance (A z ) plotte as a function of θ tot for 3 imaging tasks. Reasonable agreement between theory an experiment is observe, with prewhitening observer moels (PW an PWE) appearing to overestimate human response, while non-prewhitening observer moels (NPW an NPWE) emonstrating closer agreement. 6. DISCUSSION AND CONCLUSION Fourier-base physical characterization of imaging systems has been commonly applie in etector esign. In esigning an evaluating meical imaging systems with increasing complexity as in tomosynthesis an CBCT this approach can be extene to provie simple, elegant performance metrics, but correlation of theoretical preictions with human observer response has not been previously emonstrate. On the other han, human observer-base characterization realistically escribes imaging performance, but can be expensive an time-consuming. This work presente generalize etectability inex ( ) as a quantitative metric combining NEQ, anatomical backgroun, imaging task, an moel observers. Theoretical calculations of were valiate with real human observer performance, showing to be a meaningful metric in moeling system performance. It is important to acknowlege a variety of assumptions an limitations associate with in this stuy. Cascae systems analysis of imaging systems assumes linearity, stationarity, an shift invariance. While tomosynthesis an CBCT are recognize to be neither strictly stationary nor shift invariant, we assume the first an secon orer statistics are at least locally stationary an that system response is shift-invariant to the extent that the stuies involve tasks greater than the pixel / voxel size analogous to the D shift-invariance escribe in Ref. 33. In calculating observer performance metrics A z, the istribution of the two hypotheses are assume to be Gaussian an homosceestic. Finally, the evaluation involve a small number of simple, iealize imaging tasks, an the clutter phantom, though measure to present powerlaw backgroun, oes not truly simulate the complexity of structures in human organs. Future work involves the investigation of a broaer range of acquisition techniques, especially in the low-ose limit where etection transitions between anatomical noise-limite an quantum-limite regimes. Other reconstruction techniques will also be stuie, incluing ifferent reconstruction filters an sampling, which have been shown to impart significant effects on system NEQ. More complex an higher-orer tasks will also be investigate as a fuller examination of the extent to which generalize etectability inex emonstrates agreement with real observer performance. 7. ACKNOWLEDGEMENTS The authors acknowlege numerous stimulating conversations an collaboration with Dr. Rebecca Fahrig (Stanfor University), Dr. Sungwon Yoon (Stanfor University), Dr. Angel Pinea (California State University Fullerton), an Dr. Art Burgess (Harvar University, retire). This work was fune by National Institute of Health Grant No. R01-CA Proc. of SPIE Vol Y-9

10 REFERENCES 1. Cunningham, I. A., Westmore, M. S., an Fenster, A. "A Spatial-Frequency Depenent Quantum Accounting Diagram an Detective Quantum Efficiency Moel of Signal an Noise Propagation in Cascae Imaging Systems," Me Phys. 1[3]: (1994).. Cunningham, I. A., Yao, J, an Subotic, V. "Cascae Moels an the DQE of Flat-Penel Imagers: Noise Aliasing, Seconary Quantum Noise an Reabsorption," Physics of Meical Imaging. Proc.SPIE Meical Imaging 468: 61-7 (00). 3. Siewersen, J. H., Antonuk, L. E., el-mohri, Y., Yorkston, J., Huang, W., Boury, J. M., an Cunningham, I. A. "Empirical an Theoretical Investigation of the Noise Performance of Inirect Detection, Active Matrix Flat-Panel Imagers (AMFPIs) for Diagnostic Raiology," Me Phys. 4[1]: (1997). 4. Siewersen, J. H., Antonuk, L. E., el-mohri, Y., Yorkston, J., Huang, W., an Cunningham, I. A. "Signal, Noise Power Spectrum, an Detective Quantum Efficiency of Inirect-Detection Flat-Panel Imagers for Diagnostic Raiology," Me Phys. 5[5]: (1998). 5. Siewersen, J. H. an Jaffray, D. A. "Cone-Beam Compute Tomography With a Flat-Panel Imager: Noise Consierations for Fully 3-D Imaging," Proc.SPIE 3977: (000). 6. Siewersen, J. H. an Jaffray, D. A. "Three-Dimensional NEQ Transfer Characteristics of Volume CT Using Direct an Inirect-Detection Flat-Panel Imagers," Proc.SPIE 5030: 9-10 (003). 7. Twar, D. J. an Siewersen, J. H. "Cascae Systems Analysis of the 3D Noise Transfer Characteristics of Flat- Panel Cone-Beam CT," Me Phys. 35[1]: (008). 8. Twar, D. J., Siewersen, J. H., Fahrig, R. A., an Pinea, A. R. "Cascae Systems Analysis of the 3D NEQ for Cone-Beam CT an Tomosynthesis," Proc.SPIE 6913: 69131S-69131S-1 (008). 9. Zhao, W. an Rowlans, J. A. "Digital Raiology Using Active Matrix Reaout of Amorphous Selenium: Theoretical Analysis of Detective Quantum Efficiency," Me Phys. 4[1]: (1997). 10. Bochu, F. O., Verun, F. R., Valley, J. F., Hessler, C., an Moeckli, R. "The Importance of Anatomical Noise in Mammography," Proc.SPIE 3036: (1997). 11. Kunel, H. L., Noine, C. F., Thickman, D., Carmoy, D., an Toto, L. "Noule Detection With an Without a Chest Image," Investigative raiology 0[1]: (1985). 1. Samei, E., Flynn, M. J., an Eyler, W. R. "Detection of Subtle Lung Noules: Relative Influence of Quantum an Anatomic Noise on Chest Raiographs," Raiology 13[3]: (1999). 13. Samei, E., Eyler, W. R., an Baron, L. "Effects of Anatomical Structure on Signal Detection," Beutel, J., Kunel, H. L., an Van Metter, R. L. Hanbook of Meical Imaging: Physics an Psychophysics. 1[1]: (000). Washington, SPIE Publications. 14. Barrett, H. H., Rollan, J. P., Wagner, R. F., an Myers, K. J. "Detection an Discrimination of Known Signals in Inhomogenous, Ranom Backgrouns," Proc.SPIE Meical Imaging III: Image formation 1090: (1989). 15. Burgess, A. E., Jacobson, F. L., an Juy, P. F. "Human Observer Detection Experiments With Mammograms an Power-Law Noise," Me.Phys. 8[4]: (001). Proc. of SPIE Vol Y-10

11 16. Burgess, A. E. "Mammographic Structure: Data Preparation an Spatial Statistics Analysis," Proc.SPIE 3661: 64 (1999). 17. Glick, S. J., Veantham, S., an Karellas, A. "Investigation of Optimal KVp Setting for CT Mammography Using a Flat-Panel Imager," Proc.SPIE 468: (00). 18. Gong, X., Glick, S. J., Liu, B., Veula, A. A., an Thacker, S. "A Computer Simulation Stuy Comparing Lesion Detection Accuracy With Digital Mammography, Breast Tomosynthesis, an Cone-Beam CT Breast Imaging," Me.Phys. 33[4]: (006). 19. Heine, J. J., Deans, S. R., Velthuizen, R. P., an Clarke, L. P. "On the Statistical Nature of Mammograms," Me.Phys. 6[11]: (1999). 0. Heine, J. J. an Velthuizen, R. P. "Spectral Analysis of Full Fiel Digital Mammography Data," Me.Phys. 9[5]: (00). 1. Richar, S., Siewersen, J. H., Jaffray, D. A., Moseley, D. J., an Bakhtiar, B. "Generalize DQE Analysis of Raiographic an Dual-Energy Imaging Using Flat-Panel Detectors," Me.Phys. 3[5]: (005).. Richar, S. an Siewersen, J. H. "Optimization of Dual-Energy Imaging Systems Using Generalize NEQ an Imaging Task," Me.Phys. 34[1]: (007). 3. Yoon, S., Gang, J. G., Twar, D. J., Siewersen, J. H., an Fahrig, R. "Analysis of Lung Noule Detectability an Anatomical Clutter in Tomosynthesis Imaging of the Chest," Proc.SPIE 758: 7581M-7581M-11 (009). 4. Engstrom, E., Reiser, I., an Nishikawa, R. "Comparison of Power Spectra for Tomosynthesis Projections an Reconstructe Images," Me.Phys. 36[5]: (009). 5. Metheany, K. G., Abbey, C. K., Packar, N., an Boone, J. M. "Characterizing Anatomical Variability in Breast CT Images," Me.Phys. 35[10]: (008). 6. ICRU Report 54. "Meical Imaging - The Assessment of Image Quality," (1996). Bethesa, MD, International Commission on Raiation Units an Measurements. 7. Siewersen, J. H. an Jaffray, D. A. "Optimization of X-Ray Imaging Geometry (With Specific Application to Flat- Panel Cone-Beam Compute Tomography)," Me Phys. 7[8]: (000). 8. Chawla, A. S., Samei, E., Sauners, R. S., Lo, J. Y., an Baker, J. A. "A Mathematical Moel Platform for Optimizing a Multiprojection Breast Imaging System," Me Phys. 35[4]: (008). 9. Tapiovaara, M. J. "Efficiency of Low-Contrast Detail Detectability in Fluoroscopic Imaging," Me Phys. 4[5]: (1997). 30. Burgess, A. E., Li, X., an Abbey, C. K. "Visual Signal Detectability With Two Noise Components: Anomalous Masking Effects," J.Opt.Soc.Am.A Opt.Image Sci.Vis. 14[9]: (1997). 31. Beutel, J., Kunel, H. L., an Van Metter, R. L. "Hanbook of Meical Imaging, Volume 1. Physics an Psychophysics," [1] (000). SPIE Press. 3. Gang, G. J., Twar, D. J., Lee, J, an Siewersen, J. H. "Generalize NEQ an Detectability Inex in Cone-Beam CT: Influence of Anatomical Backgroun Clutter," Me Phys. In press. (010). 33. Albert, M. an Maiment, A. D. "Linear Response Theory for Detectors Consisting of Discrete Arrays," Me Phys. 7[10]: (000). Proc. of SPIE Vol Y-11

The Generalized NEQ and Detectability Index for Tomosynthesis and Cone-Beam CT

The Generalized NEQ and Detectability Index for Tomosynthesis and Cone-Beam CT The Generalized NEQ and Detectability Inde for Tomosynthesis and Cone-Beam CT Grace J. Gang 1, Junghoon Lee 2, J. Webster Stayman 3, Daniel Tward 3,W. Zbijewski 3, Jerry L. Prince 2, and Jeffrey H. Siewerdsen

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Ke Li and Guang-Hong Chen

Ke Li and Guang-Hong Chen Ke Li and Guang-Hong Chen Brief introduction of x-ray differential phase contrast (DPC) imaging Intrinsic noise relationship between DPC imaging and absorption imaging Task-based model observer studies

More information

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In the usual geometric derivation of Bragg s Law one assumes that crystalline Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this

More information

On the Surprising Behavior of Distance Metrics in High Dimensional Space

On the Surprising Behavior of Distance Metrics in High Dimensional Space On the Surprising Behavior of Distance Metrics in High Dimensional Space Charu C. Aggarwal, Alexaner Hinneburg 2, an Daniel A. Keim 2 IBM T. J. Watson Research Center Yortown Heights, NY 0598, USA. charu@watson.ibm.com

More information

3.2 Shot peening - modeling 3 PROCEEDINGS

3.2 Shot peening - modeling 3 PROCEEDINGS 3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

Analysis of Image Noise in 3D Cone-Beam CT: Spatial and Fourier Domain Approaches under Conditions of Varying Stationarity

Analysis of Image Noise in 3D Cone-Beam CT: Spatial and Fourier Domain Approaches under Conditions of Varying Stationarity Analysis of Image Noise in 3D Cone-Beam CT: Spatial and Fourier Domain Approaches under Conditions of Varying Stationarity Angel R. Pineda *a, Jeffrey H. Siewerdsen b,c, Daniel J. Tward b a Department

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

The Press-Schechter mass function

The Press-Schechter mass function The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for

More information

Improving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers

Improving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers International Journal of Statistics an Probability; Vol 6, No 5; September 207 ISSN 927-7032 E-ISSN 927-7040 Publishe by Canaian Center of Science an Eucation Improving Estimation Accuracy in Nonranomize

More information

Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N.

Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N. Submitte to the Journal of Rheology Moeling the effects of polyispersity on the viscosity of noncolloial har sphere suspensions Paul M. Mwasame, Norman J. Wagner, Antony N. Beris a) epartment of Chemical

More information

Evaporating droplets tracking by holographic high speed video in turbulent flow

Evaporating droplets tracking by holographic high speed video in turbulent flow Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

A PAC-Bayesian Approach to Spectrally-Normalized Margin Bounds for Neural Networks

A PAC-Bayesian Approach to Spectrally-Normalized Margin Bounds for Neural Networks A PAC-Bayesian Approach to Spectrally-Normalize Margin Bouns for Neural Networks Behnam Neyshabur, Srinah Bhojanapalli, Davi McAllester, Nathan Srebro Toyota Technological Institute at Chicago {bneyshabur,

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2 Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z

More information

Necessary and Sufficient Conditions for Sketched Subspace Clustering

Necessary and Sufficient Conditions for Sketched Subspace Clustering Necessary an Sufficient Conitions for Sketche Subspace Clustering Daniel Pimentel-Alarcón, Laura Balzano 2, Robert Nowak University of Wisconsin-Maison, 2 University of Michigan-Ann Arbor Abstract This

More information

Robustness and Perturbations of Minimal Bases

Robustness and Perturbations of Minimal Bases Robustness an Perturbations of Minimal Bases Paul Van Dooren an Froilán M Dopico December 9, 2016 Abstract Polynomial minimal bases of rational vector subspaces are a classical concept that plays an important

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Approaches for Predicting Collection Efficiency of Fibrous Filters

Approaches for Predicting Collection Efficiency of Fibrous Filters Volume 5, Issue, Summer006 Approaches for Preicting Collection Efficiency of Fibrous Filters Q. Wang, B. Maze, H. Vahei Tafreshi, an B. Poureyhimi Nonwovens Cooperative esearch Center, North Carolina State

More information

Both the ASME B and the draft VDI/VDE 2617 have strengths and

Both the ASME B and the draft VDI/VDE 2617 have strengths and Choosing Test Positions for Laser Tracker Evaluation an Future Stanars Development ala Muralikrishnan 1, Daniel Sawyer 1, Christopher lackburn 1, Steven Phillips 1, Craig Shakarji 1, E Morse 2, an Robert

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz A note on asymptotic formulae for one-imensional network flow problems Carlos F. Daganzo an Karen R. Smilowitz (to appear in Annals of Operations Research) Abstract This note evelops asymptotic formulae

More information

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads IJR International Journal of Railway Vol. 6, No. 4 / December 13, pp. 155-159 The Korean Society for Railway Vibration Analysis of Railway Tracks Force by Distribute Moving Loas Sinyeob Lee*, Dongkyu Kim*,

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

One-dimensional I test and direction vector I test with array references by induction variable

One-dimensional I test and direction vector I test with array references by induction variable Int. J. High Performance Computing an Networking, Vol. 3, No. 4, 2005 219 One-imensional I test an irection vector I test with array references by inuction variable Minyi Guo School of Computer Science

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

A Review of Multiple Try MCMC algorithms for Signal Processing

A Review of Multiple Try MCMC algorithms for Signal Processing A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications

More information

WESD - Weighted Spectral Distance for Measuring Shape Dissimilarity

WESD - Weighted Spectral Distance for Measuring Shape Dissimilarity 1 WESD - Weighte Spectral Distance for Measuring Shape Dissimilarity Ener Konukoglu, Ben Glocker, Antonio Criminisi an Kilian M. Pohl Abstract This article presents a new istance for measuring shape issimilarity

More information

KNN Particle Filters for Dynamic Hybrid Bayesian Networks

KNN Particle Filters for Dynamic Hybrid Bayesian Networks KNN Particle Filters for Dynamic Hybri Bayesian Networs H. D. Chen an K. C. Chang Dept. of Systems Engineering an Operations Research George Mason University MS 4A6, 4400 University Dr. Fairfax, VA 22030

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Hybrid Fusion for Biometrics: Combining Score-level and Decision-level Fusion

Hybrid Fusion for Biometrics: Combining Score-level and Decision-level Fusion Hybri Fusion for Biometrics: Combining Score-level an Decision-level Fusion Qian Tao Raymon Velhuis Signals an Systems Group, University of Twente Postbus 217, 7500AE Enschee, the Netherlans {q.tao,r.n.j.velhuis}@ewi.utwente.nl

More information

Evaluation of Column Breakpoint and Trajectory for a Plain Liquid Jet Injected into a Crossflow

Evaluation of Column Breakpoint and Trajectory for a Plain Liquid Jet Injected into a Crossflow ILASS Americas, 1 st Annual Conference on Liqui Atomization an Spray Systems, Orlano, Floria, May 008 Evaluation of Column Breakpoint an Trajectory for a Plain Liqui Jet Injecte into a Crossflow S.M. Thawley,

More information

Capacity Analysis of MIMO Systems with Unknown Channel State Information

Capacity Analysis of MIMO Systems with Unknown Channel State Information Capacity Analysis of MIMO Systems with Unknown Channel State Information Jun Zheng an Bhaskar D. Rao Dept. of Electrical an Computer Engineering University of California at San Diego e-mail: juzheng@ucs.eu,

More information

Search for Long-Lived Particles and Lepton-Jets with the ATLAS detector

Search for Long-Lived Particles and Lepton-Jets with the ATLAS detector EPJ Web of Conferences 6, 177 (213) DOI:.51/epjconf/2136177 Owne by the authors, publishe by EDP Sciences, 213 Search for Long-Live Particles an Lepton-Jets with the ALAS etector D. Salvatore 1,a On behalf

More information

Optical wire-grid polarizers at oblique angles of incidence

Optical wire-grid polarizers at oblique angles of incidence JOURNAL OF APPLIED PHYSICS VOLUME 93, NUMBER 8 15 APRIL 003 Optical wire-gri polarizers at oblique angles of incience X. J. Yu an H. S. Kwok a) Center for Display Research, Department of Electrical an

More information

Cascaded redundancy reduction

Cascaded redundancy reduction Network: Comput. Neural Syst. 9 (1998) 73 84. Printe in the UK PII: S0954-898X(98)88342-5 Cascae reunancy reuction Virginia R e Sa an Geoffrey E Hinton Department of Computer Science, University of Toronto,

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Shape Effect on Blind Frequency for Depth Inversion in Pulsed Thermography

Shape Effect on Blind Frequency for Depth Inversion in Pulsed Thermography Shape Effect on Blin Frequency for Depth Inversion in Pulse Thermography M. Genest 1, E. Grinzato 2, P. Bison 2, S. Marinetti 2 C. Ibarra-Castaneo 1, X. Malague 1 1 Electrical an Computing Eng. Dept.,

More information

Lower Bounds for Local Monotonicity Reconstruction from Transitive-Closure Spanners

Lower Bounds for Local Monotonicity Reconstruction from Transitive-Closure Spanners Lower Bouns for Local Monotonicity Reconstruction from Transitive-Closure Spanners Arnab Bhattacharyya Elena Grigorescu Mahav Jha Kyomin Jung Sofya Raskhonikova Davi P. Wooruff Abstract Given a irecte

More information

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 29 Sep 1997

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 29 Sep 1997 Dopant Spatial Distributions: Sample Inepenent Response Function An Maximum Entropy Reconstruction D. P. Chu an M. G. Dowsett Department of Physics, University of Warwick, Coventry CV4 7AL, UK arxiv:con-mat/970933v

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

Formulation of statistical mechanics for chaotic systems

Formulation of statistical mechanics for chaotic systems PRAMANA c Inian Acaemy of Sciences Vol. 72, No. 2 journal of February 29 physics pp. 315 323 Formulation of statistical mechanics for chaotic systems VISHNU M BANNUR 1, an RAMESH BABU THAYYULLATHIL 2 1

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an

More information

Delay Limited Capacity of Ad hoc Networks: Asymptotically Optimal Transmission and Relaying Strategy

Delay Limited Capacity of Ad hoc Networks: Asymptotically Optimal Transmission and Relaying Strategy Delay Limite Capacity of A hoc Networks: Asymptotically Optimal Transmission an Relaying Strategy Eugene Perevalov Lehigh University Bethlehem, PA 85 Email: eup2@lehigh.eu Rick Blum Lehigh University Bethlehem,

More information

DEGREE DISTRIBUTION OF SHORTEST PATH TREES AND BIAS OF NETWORK SAMPLING ALGORITHMS

DEGREE DISTRIBUTION OF SHORTEST PATH TREES AND BIAS OF NETWORK SAMPLING ALGORITHMS DEGREE DISTRIBUTION OF SHORTEST PATH TREES AND BIAS OF NETWORK SAMPLING ALGORITHMS SHANKAR BHAMIDI 1, JESSE GOODMAN 2, REMCO VAN DER HOFSTAD 3, AND JÚLIA KOMJÁTHY3 Abstract. In this article, we explicitly

More information

Optimal Cooperative Spectrum Sensing in Cognitive Sensor Networks

Optimal Cooperative Spectrum Sensing in Cognitive Sensor Networks Optimal Cooperative Spectrum Sensing in Cognitive Sensor Networks Hai Ngoc Pham, an Zhang, Paal E. Engelsta,,3, Tor Skeie,, Frank Eliassen, Department of Informatics, University of Oslo, Norway Simula

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is:

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is: Noll 3 CT Notes : Page Compute Tomograph Notes Part Challenges with Projection X-ra Sstems The equation that governs the image intensit in projection imaging is: z I I ep µ z Projection -ra sstems are

More information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How the potentials in different gauges yield the same retarded electric and magnetic fields How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department

More information

Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets

Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets Proceeings of the 4th East-European Conference on Avances in Databases an Information Systems ADBIS) 200 Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets Eleftherios Tiakas, Apostolos.

More information

Estimation of District Level Poor Households in the State of. Uttar Pradesh in India by Combining NSSO Survey and

Estimation of District Level Poor Households in the State of. Uttar Pradesh in India by Combining NSSO Survey and Int. Statistical Inst.: Proc. 58th Worl Statistical Congress, 2011, Dublin (Session CPS039) p.6567 Estimation of District Level Poor Househols in the State of Uttar Praesh in Inia by Combining NSSO Survey

More information

Sparse Reconstruction of Systems of Ordinary Differential Equations

Sparse Reconstruction of Systems of Ordinary Differential Equations Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA

More information

Journal of Engineering Science and Technology Review 7 (1) (2014) Research Article

Journal of Engineering Science and Technology Review 7 (1) (2014) Research Article Jestr Journal of Engineering Science an Technology Review 7 () (04) 8 Research Article JOURNA OF Engineering Science an Technology Review www.jestr.org Research on the Measurement Error of MWIR Average

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Topic Modeling: Beyond Bag-of-Words

Topic Modeling: Beyond Bag-of-Words Hanna M. Wallach Cavenish Laboratory, University of Cambrige, Cambrige CB3 0HE, UK hmw26@cam.ac.u Abstract Some moels of textual corpora employ text generation methos involving n-gram statistics, while

More information

To understand how scrubbers work, we must first define some terms.

To understand how scrubbers work, we must first define some terms. SRUBBERS FOR PARTIE OETION Backgroun To unerstan how scrubbers work, we must first efine some terms. Single roplet efficiency, η, is similar to single fiber efficiency. It is the fraction of particles

More information

3.0 DETECTION THEORY. Figure 3-1 IF Receiver/Signal Processor Representation

3.0 DETECTION THEORY. Figure 3-1 IF Receiver/Signal Processor Representation 3.0 DEECON HEORY 3. NRODUCON n some of our raar range equation problems we looke at fining the etection range base on NRs of 3 an 0 B. We now want to evelop some of the theory that explains the use of

More information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

American Society of Agricultural Engineers PAPER NO PRAIRIE RAINFALL,CHARACTERISTICS

American Society of Agricultural Engineers PAPER NO PRAIRIE RAINFALL,CHARACTERISTICS - PAPER NO. 79-2108 PRAIRIE RAINFALL,CHARACTERISTICS G.E. Dyck an D.M. Gray Research Engineer an Chairman Division of Hyrology University of Saskatchewan Saskatoon, Saskatchewan, Canaa For presentation

More information

Perturbation Analysis and Optimization of Stochastic Flow Networks

Perturbation Analysis and Optimization of Stochastic Flow Networks IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,

More information

A Modification of the Jarque-Bera Test. for Normality

A Modification of the Jarque-Bera Test. for Normality Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

Power law behavior in diffraction from fractal surfaces

Power law behavior in diffraction from fractal surfaces Surface Science 409 (1998 L703 L708 Surface Science Letters Power law behavior in iffraction from fractal surfaces Y.-P. Zhao *, C.-F. Cheng 1, G.-C. Wang, T.-M. Lu Department of Physics, Applie Physics

More information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,

More information

Error Floors in LDPC Codes: Fast Simulation, Bounds and Hardware Emulation

Error Floors in LDPC Codes: Fast Simulation, Bounds and Hardware Emulation Error Floors in LDPC Coes: Fast Simulation, Bouns an Harware Emulation Pamela Lee, Lara Dolecek, Zhengya Zhang, Venkat Anantharam, Borivoje Nikolic, an Martin J. Wainwright EECS Department University of

More information

Collective optical effect in complement left-handed material

Collective optical effect in complement left-handed material Collective optical effect in complement left-hane material S.-C. Wu, C.-F. Chen, W. C. Chao, W.-J. Huang an H. L. Chen, National Nano Device Laboratories, 1001-1 Ta-Hsueh Roa, Hsinchu, Taiwan R.O.C A.-C.

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Wavelet-based Local Tomography

Wavelet-based Local Tomography Inian Society for Non-Destructive Testing Hyeraba Chapter Proc. National Seminar on Non-Destructive Evaluation Dec. 7-9, 006, Hyeraba Wavelet-base Local Tomography M.V. Gopala Rao, E.M.L. Tanua an S. Vathsal

More information

Non-Linear Bayesian CBRN Source Term Estimation

Non-Linear Bayesian CBRN Source Term Estimation Non-Linear Bayesian CBRN Source Term Estimation Peter Robins Hazar Assessment, Simulation an Preiction Group Dstl Porton Down, UK. probins@stl.gov.uk Paul Thomas Hazar Assessment, Simulation an Preiction

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Parameter estimation: A new approach to weighting a priori information

Parameter estimation: A new approach to weighting a priori information Parameter estimation: A new approach to weighting a priori information J.L. Mea Department of Mathematics, Boise State University, Boise, ID 83725-555 E-mail: jmea@boisestate.eu Abstract. We propose a

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Generalized DQE analysis of radiographic and dual-energy imaging using flat-panel detectors

Generalized DQE analysis of radiographic and dual-energy imaging using flat-panel detectors Generalized DQE analysis of radiographic and dual-energy imaging using flat-panel detectors S. Richard Department of Medical Biophysics, University of Toronto, Ontario, Canada M5G M9 J. H. Siewerdsen a

More information

Asymptotics of a Small Liquid Drop on a Cone and Plate Rheometer

Asymptotics of a Small Liquid Drop on a Cone and Plate Rheometer Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity

More information

Thermal-Photon and Residual-Gas Scattering in the NLC Beam Delivery

Thermal-Photon and Residual-Gas Scattering in the NLC Beam Delivery Thermal-Photon an Resiual-Gas Scattering in the NLC Beam Delivery I. Reichel, F. Zimmermann, P. Tenenbaum, T.O. Raubenheimer Stanfor Linear Accelerator Center, Stanfor University, California 94309 Abstract

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

Robust Low Rank Kernel Embeddings of Multivariate Distributions

Robust Low Rank Kernel Embeddings of Multivariate Distributions Robust Low Rank Kernel Embeings of Multivariate Distributions Le Song, Bo Dai College of Computing, Georgia Institute of Technology lsong@cc.gatech.eu, boai@gatech.eu Abstract Kernel embeing of istributions

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information