Functional tomography using a time-gated ICCD camera
|
|
- Elmer Wilson
- 6 years ago
- Views:
Transcription
1 Functionl tomogrphy using time-gted ICCD cmer Qing Zho, 1 Lorenzo Spinelli, 2 Andre Bssi, 3 Ginluc Vlentini, 2,3 Dvide Contini, 3 Alessndro Torricelli, 3 Rinldo Cubeddu, 2,3 Giovnni Zccnti, 4 Fbrizio Mrtelli, 4 nd Antonio Pifferi 2,3,* 1 Dept Robotics Brin nd Cognitive Sciences, Istituto Itlino di Tecnologi, vi Morego 30, I Genov, Itly 2 Istituto di Fotonic e Nnotecnologie, Sezione di Milno, Consiglio Nzionle delle Ricerche, pizz Leonrdo d Vinci 32, I Miln, Itly 3 Diprtimento di Fisic, Politecnico di Milno, pizz Leonrdo d Vinci 32, Miln, Itly 4 Diprtimento di Fisic e Astronomi, Università degli Studi di Firenze, vi G. Snsone 1, Sesto Fiorentino, Firenze, Itly *ntonio.pifferi@polimi.it Abstrct: We present system for ner infrred functionl tomogrphy bsed on single pulsed source nd time-gted cmer, for non-contct collection over lrge re. The men penetrtion depth of diffusely reflected photons is dependent on the rrivl time of photons, but not on the source detector distnce. Thus, time-encoded dt cn be used to recover depth informtion while photon exiting point is exploited for lterl locliztion. This pproch ws tested ginst simultions, demonstrting both detection nd locliztion cpbilities. Preliminry mesurements on inhomogeneous phntoms showed good detection sensibility, even for low opticl perturbtion, nd locliztion cpbilities, yet with decresing sptil resolution for incresing depths. Potentil ppliction of this method to in vivo functionl studies on the brin is discussed Opticl Society of Americ OCIS codes: ( ) Photon migrtion; ( ) Medicl optics instrumenttion; ( ) Time-resolved imging. References nd links 1. M. S. Ptterson, B. Chnce, nd B. C. Wilson, Time resolved reflectnce nd trnsmittnce for the non-invsive mesurement of tissue opticl properties, Appl. Opt. 28(12), (1989). 2. J. Steinbrink, H. Wbnitz, H. Obrig, A. Villringer, nd H. Rinneberg, Determining chnges in NIR bsorption using lyered model of the humn hed, Phys. Med. Biol. 46(3), (2001). 3. A. Liebert, H. Wbnitz, J. Steinbrink, H. Obrig, M. Möller, R. Mcdonld, A. Villringer, nd H. Rinneberg, Time-resolved multidistnce ner-infrred spectroscopy of the dult hed: intrcerebrl nd extrcerebrl bsorption chnges from moments of distribution of times of flight of photons, Appl. Opt. 43(15), (2004). 4. D. Contini, A. Torricelli, A. Pifferi, L. Spinelli, F. Pgli, nd R. Cubeddu, Multi-chnnel time-resolved system for functionl ner infrred spectroscopy, Opt. Express 14(12), (2006). 5. J. Selb, J. J. Stott, M. A. Frnceschini, A. G. Sorensen, nd D. A. Bos, Improved sensitivity to cerebrl hemodynmics during brin ctivtion with time-gted opticl system: nlyticl model nd experimentl vlidtion, J. Biomed. Opt. 10(1), (2005). 6. J. C. Hebden, A. Gibson, R. M. Yusof, N. Everdell, E. M. Hillmn, D. T. Delpy, S. R. Arridge, T. Austin, J. H. Meek, nd J. S. Wytt, Three-dimensionl opticl tomogrphy of the premture infnt brin, Phys. Med. Biol. 47(23), (2002). 7. Y. Hoshi, I. Od, Y. Wd, Y. Ito, M. Yutk Ymshit, K. Od, Y. Oht, Ymd, nd Mmoru Tmur, Visuosptil imgery is fruitful strtegy for the digit spn bckwrd tsk: study with ner-infrred opticl tomogrphy, Brin Res. Cogn. Brin Res. 9(3), (2000). 8. J. Selb, D. K. Joseph, nd D. A. Bos, Time-gted opticl system for depth-resolved functionl brin imging, J. Biomed. Opt. 11(4), (2006). 9. J. Selb, E. M. C. Hillmn, D. Joseph, nd D. A. Bos, Discrimintion between superficil nd cerebrl signls during functionl brin imging with time-gted system, presented t the Europen Conferences on Biomedicl Optics (ECBO), Munich, Germny, June 13 16, J. Selb, A. M. Dle, nd D. A. Bos, Liner 3D reconstruction of time-domin diffuse opticl imging differentil dt: improved depth locliztion nd lterl resolution, Opt. Express 15(25), (2007). # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 705
2 11. W. Becker, Advnced TCSPC with Advnced Time-Correlted Single Photon Counting Techniques (Springer, Berlin, Germny 2006). 12. A. Torricelli, A. Pifferi, L. Spinelli, R. Cubeddu, F. Mrtelli, S. Del Binco, nd G. Zccnti, Time-resolved reflectnce t null source-detector seprtion: improving contrst nd resolution in diffuse opticl imging, Phys. Rev. Lett. 95(7), (2005). 13. A. Pifferi, A. Torricelli, L. Spinelli, D. Contini, R. Cubeddu, F. Mrtelli, G. Zccnti, A. Tosi, A. D. Mor, F. Zpp, nd S. Cov, Time-resolved diffuse reflectnce t null source-detector seprtion using fst gted single-photon vlnche diode, Phys. Rev. Lett. 100(138), 101 (2008). 14. P. Swosz, M. Kcprzk, A. Liebert, nd R. Mniewski, Appliction of time-gted, intensified CCD cmer for imging of bsorption chnges in non-homogenous medium, in 11th Mediterrnen Conference on Medicl nd Biomedicl Engineering nd Computing 2007, T. Jrm, P. Krmr, nd A. Zupnic, eds., vol. 16 of IFMBE Proceedings (Interntionl Federtion for Medicl nd Biologicl Engineering, 2007), pp R. B. Schulz, J. Peter, W. Semmler, C. D Andre, G. Vlentini, nd R. Cubeddu, Comprison of noncontct nd fiber-bsed fluorescence-medited tomogrphy, Opt. Lett. 31(6), (2006). 16. S. Del Binco, F. Mrtelli, nd G. Zccnti, Penetrtion depth of light re-emitted by diffusive medium: theoreticl nd experimentl investigtion, Phys. Med. Biol. 47(23), (2002). 17. S. Crrresi, T. S. M. Shtir, F. Mrtelli, nd G. Zccnti, Accurcy of perturbtion model to predict the effect of scttering nd bsorbing inhomogeneities on photon migrtion, Appl. Opt. 40(25), (2001). 18. S. R. Arridge, Opticl tomogrphy in medicl imging, Inverse Probl. 15(2), R41 R93 (1999). 19. H. W. Engl, M. Hnke, nd A. Neubuer, Regulriztion of Inverse Problems (Kluwer, Dordrecht, 1996). 20. D. Contini, F. Mrtelli, nd G. Zccnti, Photon migrtion through turbid slb described by model bsed on diffusion pproximtion. I. Theory, Appl. Opt. 36(19), (1997). 21. F. Mrtelli, S. Del Binco, A. Ismelli, nd G. Zccnti, Light Propgtion Through Biologicl Tissue nd Other Diffusive Medi (SPIE, Bellinghm, Wshington, 2010), Chps. 4 nd L. Azizi, K. Zrycht, D. Ettori, E. Tinet, nd J.-M. Tulle, Ultimte sptil resolution with diffuse opticl tomogrphy, Opt. Express 17(14), (2009). 23. A. Sssroli, F. Mrtelli, nd S. Fntini, Perturbtion theory for the diffusion eqution by use of the moments of the generlized temporl point-spred function. III. Frequency-domin nd time-domin results, J. Opt. Soc. Am. A 27(7), (2010). 24. L. Spinelli, F. Mrtelli, A. Frin, A. Pifferi, A. Torricelli, R. Cubeddu, nd G. Zccnti, Clibrtion of scttering nd bsorption properties of liquid diffusive medium t NIR wvelengths. Time-resolved method, Opt. Express 15(11), (2007). 1. Introduction Light is powerful tool for in-vivo investigtion of biologicl tissues. In the nm rnge it is not hrmful for biologicl medi t low power densities (few mw/mm 2 ), permitting design of in vivo non-invsive dignostics. The low tissue bsorption in this wvelength rnge mkes it possible to look into the body t few cm of depth (e.g. the brin cortex) or through more thn 6 cm of tissue in trnsmission (e.g. compressed brest). Furthermore, light crries multiple useful informtion from the visited tissues linked either to bsorption (e.g. tissue composition, oxygention), or to scttering (e.g. tissue microstructure) or even by exploiting selective fluorescent mrkers to specific biochemicl trgets. In the sme spectrl rnge most tissues re highly scttering. This is most severe chllenge cusing scttering to bsorption coupling, strong ttenution of the remitted light with exponentil dependence on the visited depth, blurring effects impiring sptil resolution. An ttrctive pproch to study photon propgtion in highly scttering medi (photon migrtion) is the doption of time-domin scheme [1]. Seprting photons propgted through the medium s function of their trveling time permits to uncouple bsorption from scttering contributions, to probe the medium t incresing men penetrtion depths by collecting longer lived photons, nd to improve sptil resolution using erly less dispersed photons. In prticulr, in the cse of brin imging, ddressed to detection of hemoglobin content nd oxygention, time-domin mesurements hve been proposed s mens to rech deeper structures (brin cortex) uncoupling for the msking effects cused by more superficil tissues [2 7]. Different pproches hve been proposed to exploit time-encoded depth informtion, such s selection of lte time-gtes [4,5,8], subtrction of erly gtes from lte gtes [4,9], derivtion of the vrince of the Distribution of Time of Flight (DTOF) tht is lmost exclusively sensitive to bsorption chnges in deep tissue lyers [3], use of the whole DTOF interpreted with Monte Crlo bsed sensitivity fctors [2], s well s the doption of full 3D reconstruction scheme [6,10]. # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 706
3 At present, in most cses, only tiny frction of the potentilly vilble informtion is retrieved in typicl time-resolved mesurement. The ultimte physicl limittion would correspond to uniform illumintion over the hed with power density corresponding to the sfety limits (few mw/mm 2 ), combined with corresponding collection of the emitted signl from the sme re with n ccepting solid ngle of 2π. Current time-resolved systems re quite fr from this limit. Just smll frction of the vilble signl is collected due to the very poor (<1%) sptil coverge of the hed, both for sources nd detectors. Furthermore, in the cse of Time-Correlted Single-Photon counting (TCSPC) systems, the mximum detectble signl is limited by the mximum count rte of the detection electronics tht t 6 present is t best bout 5 10 counts/s for bord running t 50% of the sturted count rte [11]. This limit is esily reched with source power of bout 1 mw, nd set of collecting bundles with n interfiber distnce of few cm. Thus, it is not possible with clssicl schemes to fully exploit the future vilbility of compct pulsed fiber sources with t lest 100 times more power thn diode lsers, s well s to foresee dense coverge of the hed, or even to tke dvntge of the higher contrst nd sptil resolution offered by very short interfiber distnces [12]. A possibility to overcome these limittions is to use Single-Photon Avlnche Diode (SPAD) operted in time-gted mode [13], even though the smll detector re is t present limiting fctor if single SPAD is employed. Alterntively, it is possible to dopt time-gted Intensified CCD cmer (ICCD) tht hs lower limittions in term of mximum signl level, even though the signl qulity is possibly worse thn for TCSPC systems. This pproch hs been widely studied by Selb et l. [5,8,10] empowering set of bundles of different lengths to smple different propgtion times with single time-gte. Alterntively, Swosz et l. [14] hve dopted completely non-contct scheme with scnning source nd wide field collection. In this pper we hve implemented system cpble to cquire lrge dt set from wide region surrounding single injection point, exploiting the time informtion to recover 3D tomogrphy. The system is bsed on single injection source nd on gted cmer for cquisition. Compred to previous work of Selb et l. [8] here we explore lrge re with non-contct pproch tht could permit better light hrvesting together with higher number of source-detector pirs. In the cse of moleculr imging, wide-field non-contct pproch hs been proven superior to fiber-bsed collection schemes for tomogrphic reconstruction both in terms of sptil resolution nd imge qulity [15]. Also, sequentil cquisition of different time windows permits to chnge light ttenution for ech imge so to properly fill the whole dynmic rnge of the CCD, nd rech better equliztion of the signl t different time gtes. Conversely, with respect to the work of Swosz et l. [14], here we present complete 3D tomogrphic pproch. In ddition, the use of single fixed lunching fiber mkes it possible to recover full 3D tomogrphy in short time, comptible with functionl imging studies. Compred to our current implementtion, the work of Swosz et l. is more ggressive, iming t performing mesurements even t null interfiber distnce, tht permits, lso, to chieve fully non-contct scheme. Here we follow more conservtive pproch, using shield to msk smll circle round the lunching fiber tht prevents to expose the cmer to the high photon flux rising from smll distnces, nd most of ll to reduce rdildependent vritions in signl level. Nevertheless, it will be possible to implement in future sptil filter on the imging optics so to fcilitte fully non-contct mesurements. In the following, we will first revise two concepts tht constitute the bsis of the present work, nd then we will test the vlidity of the method for 3D functionl tomogrphy both on simultions nd on phntom mesurements. Finlly, the possible ppliction to in vivo functionl imging of brin ctivity will be discussed. 2. Concepts The key ides underpinning this work re depicted in Fig. 1, nd Fig. 2. A first concept (Fig. 1) is the use of time to reconstruct 3D tomogrphy. Considering n homogeneous diffusive medium, upon incresing the rrivl time of photons, the men depth of visited tissue increses s well. Most noticebly, the men depth is not ffected by the source-detector # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 707
4 distnce, being encoded solely in the time. This opens the possibility for 3D tomogrphic scheme with single injection point, where depth informtion is provided by the photon trnsit time, while lterl informtion is gined from the different collection points. The grph in Fig. 1 shows the men penetrtion depth s function of the photon rrivl time, clculted s described in Ref [16]. for homogeneous medium with reduced scttering coefficient ( µ s ' ) of 10 cm 1. The function does not depend neither on the bsorption coefficient ( µ ), nor on the source-detector seprtion ( ρ ). () men depth (cm) (b) time (ps) Fig. 1. Use of time to explore different depths in the medium. Left: upon incresing the photon rrivl time, the probbility function of photon pths (bnn shpe) gets deeper, in sme wy for ll source-detector couples. Right: men depth of photon pths s function of the photon trveling time, clculted for µ ' = 10 cm 1 s described in Ref [16]. The plot is the sme for s ny source-detector seprtion. The second concept is to fully exploit light exiting from the tissue (Fig. 2). In typicl setup for opticl functionl imging, certin number of bundles re used to couple light to nd from the tissue. Sptil coverge is rther poor (bout 0.25%) since reltively lrge sourcedetector distnce is dopted nd only limited number of bundles cn be positioned on the hed. Here we wnt to exploit much lrger filling fctor (bout 90%) by collecting most of the exiting light strting from lower distnce of 1 cm from the injection point up to mximum distnce of 3 cm. The inner region is excluded (shielded) to void sturtion of the ICCD due to the overwhelming number of erly photons t short distnces. Conversely, t lrger distnces, the rdil dependence of signl becomes smoother, prticulrly t lter times. The grph in Fig. 2 shows the rdil dependence of the emitted signl for 5 different delys (1, 2, 3, 4, nd 5 ns). All curves re normlized to the vlue derived for ρ = 1 cm. By pplying these two concepts we wnt to ttin enough informtion to get 3D reconstruction from single injection point nd in short time. These two concepts were shown for homogeneous medium. Similr trends re expected lso for n inhomogeneous cse, with smll perturbtion, since the effect both on depth sensitivity (Fig. 1) nd signl level (Fig. 2) is smll. Out of the smll perturbtion regime, those concepts re still vlid qulittively, lthough the trends could be quite distorted. 3. Reconstruction strtegy 3.1. Forwrd model The geometry of the problem is depicted in Fig. 3. One single source ws used both in simultions nd in phntom experiments. The position of the source ws defined s the origin of the Crtesin coordinte system. Collection points on the medium surfce re identified by r, while r * loctes ny point (voxel) within the diffusive medium. The time t represents the rrivl time of photons t the detector point. # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 708
5 () detection re useful re (φ 6 cm) injection fiber collection bundle (φ 0.3 cm) filling fctor = 0.25% signl counts (.u.) (b) dynmic rnge of signl ns 2 ns ns ns 5 ns rdil distnce (cm) Fig. 2. Use of lrge collection re. Left: ssuming source-detector pir with detector rdius of 0.15 cm, compred to the useful re with rdius of 3 cm, the sptil coverge (detector re over totl re) is only 0.25% of the physicl limit. Right: signl level clculted s function of distnce from the source for different vlues of the photon rrivl time. All curves re normlized to the vlue t 1 cm, which is the minimum distnce used in this pper. Upon incresing time, the signl gets more uniform over the whole collection re. Φ0 r,t reference photon fluence for the unperturbed medium (for n in vivo functionl imging experiment, this cn correspond to n initil rest stte), nd with Φ r,t the perturbed photon fluence obtined when the bsorption coefficient in the medium We indicte with ( ) ( ) µ ( r *) t loction r* is chnged compred to the reference initil stte 0 ( *) mount µ ( r* ) = µ ( r* ) µ ( r *). 0 µ r by the For smll bsorption chnges the problem cn be simplified ssuming linerized perturbtion pproch: smll loclized perturbtion µ ( r *) yields perturbed photon fluence given by [2,3,17]: where ( rr, *,t ) int ( µ ( r ) ( rr t) ) Φ = Φ exp *, *, (1) 0 int is the time-resolved internl men pthlength, tht photons, detected t time t, hve trveled inside the inhomogeneity [17]. Thus, the observed chnge in signl ttenution cn be expressed s: ( µ, rr, *, ) ln ( ) µ ( r* ) ( rr, *, ) A t = ΦΦ 0 = int t (2) Moving from the continuous to the discretized cse, the reconstructed region is volume ( cm 3 ) centered beneth the injection source, nd divided into n voxels cubic voxels with side 0.3 cm. The xy plne, prllel to the surfce (detector plne), ws discretized t steps of 0.3 cm both in the x nd y directions. Then, ll pixels included in the field of view of the imging system (5.5 7 cm 2 ), nd out of the circulr blck shield ( 2 cm) were selected s detectors, yielding totl of n SD source-detector pirs. For every collection point, set of n gtes mesurements re derived, corresponding to integrtion from n initil time t G up to t G + t, where t is the width of the time-gted window. The width t is kept fixed, while t G is moved t incresed delys, yielding mesurements with incresing verge propgtion times. Mesurements re stored in the rry ΔA, with dimensions n SD n gtes. A globl chnge in the bsorption properties of the medium cn be expressed by the rry Δμ, n elements, defined for every voxel of the medium. Assuming tht the time-resolved with voxels # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 709
6 fluence resulting from distributed chnge in the bsorption coefficient Δμ cn be obtined s the superposition of loclized perturbtions, from Eq. (2) we cn derive where ( i, j, t G) ΔA = L Δμ (3) Lrr is mtrix with indices ddressing the detector position r i, the voxel position r j, nd the time window G nsd ngtes nvoxels. Mtrix L represents the men internl pthlengths spent within ech voxel by photons mesured t ech detector position nd ech time gte. Since the longer the men pthlength in given voxel, the lrger the effect produced on the mesurements, the men pthlength mtrix L is lso clled sensitivity mtrix [10]. We clculted it using the nlyticl expression for int ( r, r*,t) provided in Ref. [17], Eq. (6), for homogeneous medium with known opticl properties, evluted t time t G, tht is t the beginning of ech time window. This pproximtion cn be justified for lrger t G tking into ccount tht the mjority of the photons collected in ech time window rrives to the detector t the opening of the gte, being the photon fluence rpidly decying in time. Yet, some lrger discrepncy is expected t erly times. t, nd dimensions ( ) 5.5 cm detector plne 7x5.5 cm 2 shield source O x r detector y 7 cm r* 3.6 cm reconstructed volume 6x6x3.6 cm 3 z 6 cm Fig. 3. Geometry of the problem. A single injection source is set on the origin of the Crtesin system. The volume is divided into n voxels cubic voxels, ddressed by r*. The surfce not covered by the blck circulr shield round the injection source is divided into n SD squre detectors, ddressed by r Inverse problem The estimted imge of [18,19]: Δμ in Eq. (3) ws reconstructed using the Tikhonov regulriztion T ( ) 1 T Δμ = L LL +αi ΔOD (4) 6 cm The prmeter α in Eq. (4) ws used to restrin the sensitivity of Δμ to the noise level of the dt. Let s go bck now to the two leding concepts formulted in the previous prgrph. Regrding Fig. 1, we observe tht the men photon internl pthlength int in deeper voxels increses upon incresing the photon rrivl times t, resulting in higher sensitivity of the mesurement to the locl perturbtion (see Fig. 10 in Ref. [17]). Thus the L mtrix provides higher weight for deeper voxels upon incresing t G. In the CW cse, for ech detector position # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 710
7 single mesurement is chieved, nd the L mtrix contins the internl men pthlengths corresponding to n verge propgtion time t. Thus, deeper regions re reched exploiting lrge source-detector distnces tht yield lrger verge propgtion time. Conversely, the use of multiple time gtes permits, for every source-detector couple, to hve multiple mesurements with incresing depth sensitivity. In wy, lterl nd depth sensitivity re seprted, being the former encoded in the detector position, nd the ltter in the time gte t G. Furthermore, depth informtion is reched just by incresing t G, without chnging source-detector position, thus keeping exctly the sme collection efficiency. Regrding Fig. 2, we observe tht incresing the detection coverge we increse both the totl collected signl nd the sptil coverge of the mesurement ΔA. Since depth informtion is lredy provided by time, this should help in mending the ill-posed nture of the problem. The determintion of the optiml sectioning of the detection surfce nd the optiml choice of the time gtes is not obvious nd depends on deeper nlysis of the mtrix L nd of the vilble SNR, nd goes beyond the scope of the present pper. Conversely, we hve dopted simple choice by sectioning the detection surfce with the sme size of the voxels, nd tking number of time-gtes comptible with rel-time in vivo cquisition. Anyhow, it is resonble tht the mximum informtion is gined when the entire signl exiting the surfce is collected, nd tht spre coverge, s the one depicted in the left of Fig. 2 in the cse of the bundle is quite fr from the ultimte limit. 4. Mterils nd methods 4.1. Simultions The photon fluence Φ, simulting locl perturbtion, ws clculted using n higher order perturbtive pproch yielding better ccurcy thn the purely liner eqution [17]: ( ) ( ) int int 0 ( ) ( ) µ ( ) ( ( ) µ ( )) Φ r, t = 1 rr,, t r exp rr,, t r Φ r, t (5) where the photon fluence for the unperturbed stte Φ 0 ws clculted using the timedependent diffusion eqution solved for homogeneous slb pplying the extrpolted boundry conditions [20,21], while int ws clculted gin using Eq. (6) of Ref 17. Then, Φ ( r,t) ws integrted between t G nd t G + t for ech dely. Since, both for simultions nd phntom mesurements the slb thickness ws rther lrge, identicl results would be obtined using the semi-infinite model. The perturbtions were simulted by sphere (rdius 0.6 cm) embedded in homogeneous -1-1 medium with µ = 0.1 cm, µ s ' = 10 cm n medium = 1.33, n externl = 1 nd slb thickness of 15 cm (lmost semi-infinite). The reduced scttering coefficient of the inclusion ws kept the INC sme s the bckground, while 5 different vlues for the bsorption coefficient ( µ ) were considered, nmely 0.1 cm 1 (1 µ ), sme s bckground, 0.2 cm 1 ( 2 µ ), 0.4 cm 1 ( 4 µ ), 0.8 cm 1 ( 8 µ ), nd 2.0 cm 1 ( 20 µ ), to mimic completely bsorbing inclusion. Position in the xy plne ws (0,1) cm. Along the direction z, four depths were tested: 0.5, 1.0, 1.5, nd 2.0 cm. Eight time-delys t G were used (2.25, 2.5, 2.75, 3, 3.25, 3.5, 3.75, 4 ns), while keeping fixed the width t = 500 ps. Poisson noise ws dded to the simultion in ccordnce to wht expected for the mesurements with the ICCD. In detil, simulted imges were multiplied by proper fctor so to yield given signl intensity. At ech detector position, the signl-to-noise rtio (SNR) ws clculted s the squre root of the number of counts, thus ssuming purely Poisson noise. Then, the globl SNR ws derived s the medin of the SNR t ech detector position nd ech dely. The ctul dynmic rnge of the ICCD ws effectively tken into ccount by # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 711
8 limiting the pek vlue in ech simulted imge to the mximum counts, nd normlizing the entire imge ccordingly Experimentl set-up A scheme of the experimentl set-up is depicted in Fig. 4. A pulsed diode lser (PDL800B, PicoQunt, Germny), operting t 690 nm, with pulse width 200 ps, nd repetition rte of 80 MHz ws used s the illuminting source. Light ws delivered to the tissue by mens of 50 μm core grded index fiber, nd expnded to few mm 2 using spcer. The mximum power impinging onto the phntom or onto the tissue ws 1 mw. The lunching fiber ws inserted into blck cylinder ( = 2 cm) so to msk the region too close to the injection point tht could sturte the cmer t n erly gte. The cylinder ws kept over the surfce, with the PVC in contct with the wter solution. Detection ws chieved using time-gted intensifier tube (HRI, Kentech, Didcot, UK), coupled to low-noise, Peltier cooled CCD cmer (PCO GmbH, Gottingen, Germny), set t distnce of 20 cm over the surfce nd collecting diffusely remitted light from n re centered in the injection point. The cmer ws operted with 500 ps wide gte nd rising edge of 120 ps set t given progrmmble dely with respect to the lser pulse. A hrdwre binning of 2x2 pixels ws pplied to the ICCD so to reduce the incidence of the red-out noise. All opertions involved in the mesurement were utomted, tht re: i) moving the opticl inhomogeneity using trnslting stge; ii) djusting the gte dely; iii) rotting the circulr vrible ttenutor to reduce injected power t erly delys; iv) setting the intensifier gin; v) cquiring nd sving the imges. The imging region of the ICCD ws cm 2. Before processing, softwre binning procedure ws pplied by dding pixels in cm 2 blocks (28 28 pixels) to reduce noise level. To eliminte the defocused shding effect of the lunching fiber crossing the imged field, prt of the rw imge ws removed. The remining 205 pixels were selected s detectors. The imging region for reconstruction ws 6 6 cm 2, s indicted by the red squre block in Fig. 4(b). The blck circle indictes the position of the blck cylinder. The region enclosed by the blck lines nd the bottom line of the rectngle ws removed to void influence of the opticl fiber. 5. Results nd Discussion 5.1. Vlidtion on simulted dt sets The proposed method ws tested first on simultions. Figure 5 shows verticl yz sections of the 3D reconstructions of µ t x = 0 for incresing vlues of µ INC (rows) nd of the perturbtion depth z (columns) recovered from simulted dt. For the sme depth, the reconstructed µ ugments for incresing vlues of µ INC (ech column hs the sme colorbr). The perturbtion cn be detected t ll depths, lthough with decresing contrst nd worse sptil resolution. Also, the lterl nd xil positions of the perturbtion re recovered with resonble greement. The SNR of simultions in Fig. 5 is 1000, leding to choice of α = To investigte the relevnce of the signl level on the outcome of the reconstruction, we hve repeted the exercise using SNR = 40. The reconstructed µ is reported in Fig. 6. Although detection is still possible, compred to higher SNR simultions (Fig. 5), sptil locliztion is worsened, since the perturbed re is quite lrge, nd depth dependence is lost. This confirms tht sptil resolution in diffuse opticl tomogrphy is ultimtely limited by the vilble SNR [22]. # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 712
9 Fig. 4. () Experimentl set-up. Light source: pulsed diode lser, operted t 690 nm with 80 MHz repetition rte. Detection: time-gted intensified CCD cmer, gte width 1 ns, risetime 120 ps. The lunching fiber is embedded into blck cylinder (1 cm rdius) to shield the cmer from erly photons t short source-detector distnces. (b) One rw imge cptured by the ICCD cmer. The red rectngle shows the imging region in X-Y plne for reconstruction. The blck circle indictes the position of the blck cylinder. The region enclosed by the blck lines nd the bottom line of the rectngle ws removed to void influence of the opticl fiber. It must be reclled tht the reconstruction model is bsed on the ssumption of smll perturbtion. It hs been shown [17], tht for sphere with 0.5 cm rdius set t depth of 2 cm, this pproximtion strts deviting from the rel cse for µ = 0.2 cm 1. Lrger perturbtions could be better described using higher order perturbtion pproch [23]. Yet, this would led to non-liner problem, to be solved with n itertive pproch Phntom mesurements For the homogeneous medium we used wter solution of Intrlipid nd blck ink with -1-1 µ = 0.1 cm nd µ s ' = 10 cm [24], poured in lrge tnk (surfce 15x20 cm 2, height 19 cm). A totl of 5 blck PVC cylindricl (height equl to dimeter) inclusions, with incresing volume (0.021, 0.050, 0.098, 0.270, cm 3 ) were hold beneth the point (0,1.5) cm of the xy plne t 4 incresing depths z : 0.5, 1.0, 1.5 nd 2.0 cm (mesured from the upper surfce of the inclusion). The gting window pplied to the ICCD ws shifted in time t 8 incresing time-delys t G (2.25, 2.5, 2.75, 3, 3.25, 3.5, 3.75, 4 ns), while the gte width t ws fixed t 500 ps. The cquisition time of the CCD cmer ws 2.5 s for ech dely. Figure 7 shows verticl yz sections of the 3D reconstructions of µ t x = 0 for incresing volumes of the PVC cylinder (rows) nd of depth z (columns). The perturbtion cn be detected with enough contrst up to the lrger depth (2.0 cm). Furthermore, there is resonble dependence of the reconstructed µ on the mount of perturbtion (PVC volume). Locliztion properties re lso dequte both for the y coordinte, nd for the depth z dependence. Sptil resolution suffers for the lower contrst t lrger depths, with swelling nd decresing of mximum µ upon incresing z. # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 713
10 Fig. 5. Verticl section long the yz plne of 3D reconstructions of µ from simulted dt, using α = A loclized perturbtion (sphere with rdius 0.6 cm) with incresing vlues of INC µ compred to µ (rows) is set t incresing depths z (columns). The first nd lst row INC represent completely homogeneous cse ( µ = 1 µ ), nd very high (completely INC bsorbing) inhomogeneity ( µ = 20 µ ), respectively. Axis dimensions re in cm, while colorbr is in cm 1. The use of blck PVC cylinders ws motivted by the need to produce controllble perturbtions. Different solutions hve been proposed to simulte opticl inhomogeneities, such s glss/plstic tubes or gloves filled with liquid solution, yet some rrngements re prone to light-guiding effects, others re difficult to hndle or to control in shpe nd position. Conversely, we hve checked vi Monte Crlo simultions tht smll totlly bsorbing inhomogeneities produce perturbtion on the DTOF of the sme shpe nd mplitude of relistic bsorption chnges provided tht the depth of the inhomogeneity is equl to or lrger thn 1 cm. Although phntom mesurements hve been crried out lso for the depth z = 0.5 cm, in this sitution, for which the inhomogeneity my be very close to source or detector, the perturbtions of blck PVC cylinders is not equivlent to tht of relistic bsorption chnges. Conversely, for z 1cm, the equivlence holds true. For instnce, totl bsorber with volume of cm 3 yields temporl perturbtion quite similr to n opticl inhomogeneity with volume of 1 cm 3-1 nd µ = 0.5 cm (dt not shown). Aprt from the centrl blck disk, the medium boundry ws ir (92% of the imged re), thus n ext = 1 ws used for the refrctive index of the externl medium. Yet, some perturbtions cused by the different refrctive index of the blck holder cnnot be ruled out. The vilble source power in this experiment ws rther low (1 mw). Results presented on simultions show tht gin in SNR could indeed led to better sptil resolution nd locliztion. Thus, there is lrge spce for improvement, considering tht round 100 mw of source power could be used in vivo, provided tht the source dimeter is enlrged (e.g cm) so to sty within the sfety limits. Yet, the power stbility of the source needs to be well controlled so to permit detection of time-gted intensity chnges s low s 1%. # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 714
11 Fig. 6. Verticl section long the yz plne of 3D reconstructions of µ from simulted dt, generted s in Fig. 5, prt from lower SNR = 40 nd vlue of α = 0.1. A loclized perturbtion (sphere with rdius 0.6 cm) with incresing vlues of µ INC compred to µ (rows) is set t incresing depths z (columns). The first nd lst row represent completely INC homogeneous cse ( µ = 1 µ ), nd very high (completely bsorbing) inhomogeneity INC ( µ = 20 µ ), respectively. Axis dimensions re in cm, colorbr is in cm Potentil ppliction to functionl imging of the brin One of the key dvntges of the proposed pproch is the exploittion of both significnt light hrvesting s well s optiml collection of vilble informtion. As consequence full 3D cquisition is fesible in reltively short time. Phntom mesurements hve been demonstrted with totl cquisition time per reconstructed imge of 15 s. The vilble power ws just 1 mw, with room for relistic gin by 2 orders of mgnitude, possibly resulting in better imge qulity nd/or reduced cquisition time. This mkes the technique prticulrly ppeling for in vivo functionl imging, involving the reconstruction of bsorption chnges compred to n initil reference stte. The study of the pplicbility of this pproch to in vivo functionl imging will be ddressed in the ner future with the id of more powerful source. By now we just performed some preliminry tests on shved subject leding to 3 observtions: i) when the subject is lying on bed, it is rther esy to keep the sme working distnce of the cmer from the hed (20 cm) nd field of view (5.5x7.0 cm 2 ) s for phntom mesurements; ii) the signl level t 690 nm is similr to wht observed on phntoms; iii) movement rtifcts seem not to be so criticl, with shift in the position of the lunching fiber, identified by the blck cylinder fiber holder, less thn the size of one voxel (<3 mm) during 15 min. Clerly, the presence of hir could hmper prcticl pplicbility of the proposed pproch. In ny cse, ll mesurements performed on the forehed (e.g. those relted to cognitive tsks) could be esily mnged. Also, on regions covered by hir, proper combing could crete rows or zones of exposed skin. We will study whether this kind of reduced dt set is still vluble to provide 3D reconstruction. # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 715
12 Fig. 7. Verticl section long the yz plne with x = 0 cm of tomogrphies of µ from phntom mesurements, using α = 0.1. From left to right the depth (Z) of the inclusion is incresed. From top to bottom, the volume of the PVC cylinder is incresed. Axis dimensions re in cm, while colorbr is in cm Conclusion We hve proposed novel scheme to functionl tomogrphy bsed on single source injection point nd on time-gted ICCD cmer for detection. A non-contct modlity is dopted for collection, with direct imging of the medium surfce vi the cmer lens. Shielding of smll region round the lunching fiber helps preventing reduction of the useful dynmic rnge. Simultions of loclized perturbtion of the bsorption coefficient demonstrte the fesibility of the proposed pproch both for detection nd locliztion of smll perturbtion, with n imge qulity clerly dependent on the SNR. Phntom mesurements using blck PVC cylinders demonstrte detection of cm 3 perturbtion down to depth of 2 cm within diffusive tissue-like medium. Also, firly good locliztion is chieved both in depth nd lterl direction, while sptil resolution nd contrst is degrded for incresing depth. Totl cquisition time to cquire full 3D dt set ws 15 s for n injected power of 1 mw. Work is in progress to study the pplicbility of this method to in vivo functionl imging of the brin nd to exploit more powerful (100 mw) sources. Acknowledgments The reserch leding to these results hs been prtilly funded by the Europen Community s Seventh Frmework Progrmme under the neuropt project (FP7-HEALTH ). # $15.00 USD Received 23 Dec 2010; revised 18 Feb 2011; ccepted 19 Feb 2011; published 25 Feb 2011 (C) 2011 OSA 1 Mrch 2011 / Vol. 2, No. 3 / BIOMEDICAL OPTICS EXPRESS 716
MATH 2530: WORKSHEET 7. x 2 y dz dy dx =
MATH 253: WORKSHT 7 () Wrm-up: () Review: polr coordintes, integrls involving polr coordintes, triple Riemnn sums, triple integrls, the pplictions of triple integrls (especilly to volume), nd cylindricl
More informationCHAPTER III IMAGE DEWARPING (CALIBRATION) PROCEDURE
CHAPTER III IMAGE DEWARPING (CALIBRATION) PROCEDURE 3.1 Scheimpflug Configurtion nd Perspective Distortion Scheimpflug criterion were found out to be the best lyout configurtion for Stereoscopic PIV, becuse
More informationFig.1. Let a source of monochromatic light be incident on a slit of finite width a, as shown in Fig. 1.
Answer on Question #5692, Physics, Optics Stte slient fetures of single slit Frunhofer diffrction pttern. The slit is verticl nd illuminted by point source. Also, obtin n expression for intensity distribution
More informationa < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1
Mth 33 Volume Stewrt 5.2 Geometry of integrls. In this section, we will lern how to compute volumes using integrls defined by slice nlysis. First, we recll from Clculus I how to compute res. Given the
More informationDigital approximation to extended depth of field in no telecentric imaging systems
Journl of Physics: Conference Series Digitl pproximtion to extended depth of field in no telecentric imging systems To cite this rticle: J E Meneses nd C R Contrers 0 J Phys: Conf Ser 74 0040 View the
More informationSection 10.4 Hyperbolas
66 Section 10.4 Hyperbols Objective : Definition of hyperbol & hyperbols centered t (0, 0). The third type of conic we will study is the hyperbol. It is defined in the sme mnner tht we defined the prbol
More informationInternational Conference on Mechanics, Materials and Structural Engineering (ICMMSE 2016)
\ Interntionl Conference on Mechnics, Mterils nd tructurl Engineering (ICMME 2016) Reserch on the Method to Clibrte tructure Prmeters of Line tructured Light Vision ensor Mingng Niu1,, Kngnin Zho1, b,
More informationStained Glass Design. Teaching Goals:
Stined Glss Design Time required 45-90 minutes Teching Gols: 1. Students pply grphic methods to design vrious shpes on the plne.. Students pply geometric trnsformtions of grphs of functions in order to
More informationSUPPLEMENTARY INFORMATION
Supplementry Figure y (m) x (m) prllel perpendiculr Distnce (m) Bird Stndrd devition for distnce (m) c 6 prllel perpendiculr 4 doi:.8/nture99 SUPPLEMENTARY FIGURE Confirmtion tht movement within the flock
More informationChapter 2 Sensitivity Analysis: Differential Calculus of Models
Chpter 2 Sensitivity Anlysis: Differentil Clculus of Models Abstrct Models in remote sensing nd in science nd engineering, in generl re, essentilly, functions of discrete model input prmeters, nd/or functionls
More informationComplete Coverage Path Planning of Mobile Robot Based on Dynamic Programming Algorithm Peng Zhou, Zhong-min Wang, Zhen-nan Li, Yang Li
2nd Interntionl Conference on Electronic & Mechnicl Engineering nd Informtion Technology (EMEIT-212) Complete Coverge Pth Plnning of Mobile Robot Bsed on Dynmic Progrmming Algorithm Peng Zhou, Zhong-min
More informationParallel Square and Cube Computations
Prllel Squre nd Cube Computtions Albert A. Liddicot nd Michel J. Flynn Computer Systems Lbortory, Deprtment of Electricl Engineering Stnford University Gtes Building 5 Serr Mll, Stnford, CA 945, USA liddicot@stnford.edu
More information3.5.1 Single slit diffraction
3..1 Single slit diffrction ves pssing through single slit will lso diffrct nd produce n interference pttern. The reson for this is to do with the finite width of the slit. e will consider this lter. Tke
More informationOn the Detection of Step Edges in Algorithms Based on Gradient Vector Analysis
On the Detection of Step Edges in Algorithms Bsed on Grdient Vector Anlysis A. Lrr6, E. Montseny Computer Engineering Dept. Universitt Rovir i Virgili Crreter de Slou sin 43006 Trrgon, Spin Emil: lrre@etse.urv.es
More information6.3 Volumes. Just as area is always positive, so is volume and our attitudes towards finding it.
6.3 Volumes Just s re is lwys positive, so is volume nd our ttitudes towrds finding it. Let s review how to find the volume of regulr geometric prism, tht is, 3-dimensionl oject with two regulr fces seprted
More informationOptics and Optical design Problems
Optics nd Opticl design 0 Problems Sven-Görn Pettersson / Cord Arnold 0-09-06 4:3 This mteril is tken from severl sources. Some problems re from the book Våglär och Optik by Görn Jönsson nd Elisbeth Nilsson.
More information4452 Mathematical Modeling Lecture 4: Lagrange Multipliers
Mth Modeling Lecture 4: Lgrnge Multipliers Pge 4452 Mthemticl Modeling Lecture 4: Lgrnge Multipliers Lgrnge multipliers re high powered mthemticl technique to find the mximum nd minimum of multidimensionl
More information6.2 Volumes of Revolution: The Disk Method
mth ppliction: volumes by disks: volume prt ii 6 6 Volumes of Revolution: The Disk Method One of the simplest pplictions of integrtion (Theorem 6) nd the ccumultion process is to determine so-clled volumes
More informationIf f(x, y) is a surface that lies above r(t), we can think about the area between the surface and the curve.
Line Integrls The ide of line integrl is very similr to tht of single integrls. If the function f(x) is bove the x-xis on the intervl [, b], then the integrl of f(x) over [, b] is the re under f over the
More informationOPTICS. (b) 3 3. (d) (c) , A small piece
AQB-07-P-106 641. If the refrctive indices of crown glss for red, yellow nd violet colours re 1.5140, 1.5170 nd 1.518 respectively nd for flint glss re 1.644, 1.6499 nd 1.685 respectively, then the dispersive
More informationThe Nature of Light. Light is a propagating electromagnetic waves
The Nture of Light Light is propgting electromgnetic wves Index of Refrction n: In mterils, light intercts with toms/molecules nd trvels slower thn it cn in vcuum, e.g., vwter The opticl property of trnsprent
More informationEngineer To Engineer Note
Engineer To Engineer Note EE-169 Technicl Notes on using Anlog Devices' DSP components nd development tools Contct our technicl support by phone: (800) ANALOG-D or e-mil: dsp.support@nlog.com Or visit
More information3.5.1 Single slit diffraction
3.5.1 Single slit diffrction Wves pssing through single slit will lso diffrct nd produce n interference pttern. The reson for this is to do with the finite width of the slit. We will consider this lter.
More information1. SEQUENCES INVOLVING EXPONENTIAL GROWTH (GEOMETRIC SEQUENCES)
Numbers nd Opertions, Algebr, nd Functions 45. SEQUENCES INVOLVING EXPONENTIAL GROWTH (GEOMETRIC SEQUENCES) In sequence of terms involving eponentil growth, which the testing service lso clls geometric
More informationPower Transmittance of a Laterally Shifted Gaussian Beam through a Circular Aperture
Poer Trnsmittnce of Lterlly Shifted Gussin Bem through Circulr Aperture Triq Shmim Khj 1 nd Syed Azer Rez 1 1. Deprtment of Electricl Engineering, Lhore University of Mngement Sciences, DHA, Lhore 5479,
More informationBefore We Begin. Introduction to Spatial Domain Filtering. Introduction to Digital Image Processing. Overview (1): Administrative Details (1):
Overview (): Before We Begin Administrtive detils Review some questions to consider Winter 2006 Imge Enhncement in the Sptil Domin: Bsics of Sptil Filtering, Smoothing Sptil Filters, Order Sttistics Filters
More informationMath 35 Review Sheet, Spring 2014
Mth 35 Review heet, pring 2014 For the finl exm, do ny 12 of the 15 questions in 3 hours. They re worth 8 points ech, mking 96, with 4 more points for netness! Put ll your work nd nswers in the provided
More informationL. Yaroslavsky. Fundamentals of Digital Image Processing. Course
L. Yroslvsky. Fundmentls of Digitl Imge Processing. Course 0555.330 Lecture. Imge enhncement.. Imge enhncement s n imge processing tsk. Clssifiction of imge enhncement methods Imge enhncement is processing
More informationUSING HOUGH TRANSFORM IN LINE EXTRACTION
Stylinidis, Efstrtios USING HOUGH TRANSFORM IN LINE EXTRACTION Efstrtios STYLIANIDIS, Petros PATIAS The Aristotle University of Thessloniki, Deprtment of Cdstre Photogrmmetry nd Crtogrphy Univ. Box 473,
More information2 Computing all Intersections of a Set of Segments Line Segment Intersection
15-451/651: Design & Anlysis of Algorithms Novemer 14, 2016 Lecture #21 Sweep-Line nd Segment Intersection lst chnged: Novemer 8, 2017 1 Preliminries The sweep-line prdigm is very powerful lgorithmic design
More informationRepresentation of Numbers. Number Representation. Representation of Numbers. 32-bit Unsigned Integers 3/24/2014. Fixed point Integer Representation
Representtion of Numbers Number Representtion Computer represent ll numbers, other thn integers nd some frctions with imprecision. Numbers re stored in some pproximtion which cn be represented by fixed
More informationAn Integrated Simulation System for Human Factors Study
An Integrted Simultion System for Humn Fctors Study Ying Wng, Wei Zhng Deprtment of Industril Engineering, Tsinghu University, Beijing 100084, Chin Foud Bennis, Dmien Chblt IRCCyN, Ecole Centrle de Nntes,
More information10.5 Graphing Quadratic Functions
0.5 Grphing Qudrtic Functions Now tht we cn solve qudrtic equtions, we wnt to lern how to grph the function ssocited with the qudrtic eqution. We cll this the qudrtic function. Grphs of Qudrtic Functions
More informationprisms Prisms Specifications Catalogue number BK7 Wedge, Beam Deviation, deg
Cotings Wedge Steer bems in opticl systems Cn be used in pirs for continuous ngulr djustment T Hving selected n pproprite wedge, it is esy to crete precise bem devition without ffecting other bem prmeters.
More informationUnit 5 Vocabulary. A function is a special relationship where each input has a single output.
MODULE 3 Terms Definition Picture/Exmple/Nottion 1 Function Nottion Function nottion is n efficient nd effective wy to write functions of ll types. This nottion llows you to identify the input vlue with
More information3 Talk to Us First. Reasons You Should. Non-Contact Temperature Measurement Solutions
3 Tlk to Us First Resons You Should Non-Contct Temperture Mesurement Solutions Reson # 1 Fixed System Spot Thermometers FACTS Verstile technologies re must for the ever expnding list of pplictions Lnd
More informationChapter 2. 3/28/2004 H133 Spring
Chpter 2 Newton believe tht light ws me up of smll prticles. This point ws ebte by scientists for mny yers n it ws not until the 1800 s when series of experiments emonstrte wve nture of light. (But be
More informationSolutions to Math 41 Final Exam December 12, 2011
Solutions to Mth Finl Em December,. ( points) Find ech of the following its, with justifiction. If there is n infinite it, then eplin whether it is or. ( ) / ln() () (5 points) First we compute the it:
More informationMA1008. Calculus and Linear Algebra for Engineers. Course Notes for Section B. Stephen Wills. Department of Mathematics. University College Cork
MA1008 Clculus nd Liner Algebr for Engineers Course Notes for Section B Stephen Wills Deprtment of Mthemtics University College Cork s.wills@ucc.ie http://euclid.ucc.ie/pges/stff/wills/teching/m1008/ma1008.html
More informationAgilent Mass Hunter Software
Agilent Mss Hunter Softwre Quick Strt Guide Use this guide to get strted with the Mss Hunter softwre. Wht is Mss Hunter Softwre? Mss Hunter is n integrl prt of Agilent TOF softwre (version A.02.00). Mss
More informationA New Learning Algorithm for the MAXQ Hierarchical Reinforcement Learning Method
A New Lerning Algorithm for the MAXQ Hierrchicl Reinforcement Lerning Method Frzneh Mirzzdeh 1, Bbk Behsz 2, nd Hmid Beigy 1 1 Deprtment of Computer Engineering, Shrif University of Technology, Tehrn,
More information2. What are the types of diffraction and give the differences between them? (June 2005, June 2011)
UNIT-1 b DIFFRACTION Diffrction:A) Distinction between Fresnel nd Frunhofer diffrction, B) diffrction due to single slit, N-slits,C) Diffrction grting experiment. 1 A) Distinction between Fresnel nd Frunhofer
More informationEngineer-to-Engineer Note
Engineer-to-Engineer Note EE-245 Technicl notes on using Anlog Devices DSPs, processors nd development tools Contct our technicl support t dsp.support@nlog.com nd t dsptools.support@nlog.com Or visit our
More informationLecture 10 Evolutionary Computation: Evolution strategies and genetic programming
Lecture 10 Evolutionry Computtion: Evolution strtegies nd genetic progrmming Evolution strtegies Genetic progrmming Summry Negnevitsky, Person Eduction, 2011 1 Evolution Strtegies Another pproch to simulting
More informationsuch that the S i cover S, or equivalently S
MATH 55 Triple Integrls Fll 16 1. Definition Given solid in spce, prtition of consists of finite set of solis = { 1,, n } such tht the i cover, or equivlently n i. Furthermore, for ech i, intersects i
More informationWhat do all those bits mean now? Number Systems and Arithmetic. Introduction to Binary Numbers. Questions About Numbers
Wht do ll those bits men now? bits (...) Number Systems nd Arithmetic or Computers go to elementry school instruction R-formt I-formt... integer dt number text chrs... floting point signed unsigned single
More informationA Fast Imaging Algorithm for Near Field SAR
Journl of Computing nd Electronic Informtion Mngement ISSN: 2413-1660 A Fst Imging Algorithm for Ner Field SAR Guoping Chen, Lin Zhng, * College of Optoelectronic Engineering, Chongqing University of Posts
More informationIn the last lecture, we discussed how valid tokens may be specified by regular expressions.
LECTURE 5 Scnning SYNTAX ANALYSIS We know from our previous lectures tht the process of verifying the syntx of the progrm is performed in two stges: Scnning: Identifying nd verifying tokens in progrm.
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd business. Introducing technology
More informationGeometric transformations
Geometric trnsformtions Computer Grphics Some slides re bsed on Shy Shlom slides from TAU mn n n m m T A,,,,,, 2 1 2 22 12 1 21 11 Rows become columns nd columns become rows nm n n m m A,,,,,, 1 1 2 22
More informationHW Stereotactic Targeting
HW Stereotctic Trgeting We re bout to perform stereotctic rdiosurgery with the Gmm Knife under CT guidnce. We instrument the ptient with bse ring nd for CT scnning we ttch fiducil cge (FC). Above: bse
More information12-B FRACTIONS AND DECIMALS
-B Frctions nd Decimls. () If ll four integers were negtive, their product would be positive, nd so could not equl one of them. If ll four integers were positive, their product would be much greter thn
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd business. Introducing technology
More informationWhat do all those bits mean now? Number Systems and Arithmetic. Introduction to Binary Numbers. Questions About Numbers
Wht do ll those bits men now? bits (...) Number Systems nd Arithmetic or Computers go to elementry school instruction R-formt I-formt... integer dt number text chrs... floting point signed unsigned single
More informationUnit #9 : Definite Integral Properties, Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties, Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More informationTilt-Sensing with Kionix MEMS Accelerometers
Tilt-Sensing with Kionix MEMS Accelerometers Introduction Tilt/Inclintion sensing is common ppliction for low-g ccelerometers. This ppliction note describes how to use Kionix MEMS low-g ccelerometers to
More informationEngineer-to-Engineer Note
Engineer-to-Engineer Note EE-295 Technicl notes on using Anlog Devices DSPs, processors nd development tools Visit our Web resources http://www.nlog.com/ee-notes nd http://www.nlog.com/processors or e-mil
More informationClass-XI Mathematics Conic Sections Chapter-11 Chapter Notes Key Concepts
Clss-XI Mthemtics Conic Sections Chpter-11 Chpter Notes Key Concepts 1. Let be fixed verticl line nd m be nother line intersecting it t fixed point V nd inclined to it t nd ngle On rotting the line m round
More informationMENSURATION-IV
MENSURATION-IV Theory: A solid is figure bounded by one or more surfce. Hence solid hs length, bredth nd height. The plne surfces tht bind solid re clled its fces. The fundmentl difference between plne
More informationQuestions About Numbers. Number Systems and Arithmetic. Introduction to Binary Numbers. Negative Numbers?
Questions About Numbers Number Systems nd Arithmetic or Computers go to elementry school How do you represent negtive numbers? frctions? relly lrge numbers? relly smll numbers? How do you do rithmetic?
More informationLETKF compared to 4DVAR for assimilation of surface pressure observations in IFS
LETKF compred to 4DVAR for ssimiltion of surfce pressure oservtions in IFS Pu Escrià, Mssimo Bonvit, Mts Hmrud, Lrs Isksen nd Pul Poli Interntionl Conference on Ensemle Methods in Geophysicl Sciences Toulouse,
More informationAn Efficient Divide and Conquer Algorithm for Exact Hazard Free Logic Minimization
An Efficient Divide nd Conquer Algorithm for Exct Hzrd Free Logic Minimiztion J.W.J.M. Rutten, M.R.C.M. Berkelr, C.A.J. vn Eijk, M.A.J. Kolsteren Eindhoven University of Technology Informtion nd Communiction
More informationMATH 25 CLASS 5 NOTES, SEP
MATH 25 CLASS 5 NOTES, SEP 30 2011 Contents 1. A brief diversion: reltively prime numbers 1 2. Lest common multiples 3 3. Finding ll solutions to x + by = c 4 Quick links to definitions/theorems Euclid
More informationP(r)dr = probability of generating a random number in the interval dr near r. For this probability idea to make sense we must have
Rndom Numers nd Monte Crlo Methods Rndom Numer Methods The integrtion methods discussed so fr ll re sed upon mking polynomil pproximtions to the integrnd. Another clss of numericl methods relies upon using
More informationA REINFORCEMENT LEARNING APPROACH TO SCHEDULING DUAL-ARMED CLUSTER TOOLS WITH TIME VARIATIONS
A REINFORCEMENT LEARNING APPROACH TO SCHEDULING DUAL-ARMED CLUSTER TOOLS WITH TIME VARIATIONS Ji-Eun Roh (), Te-Eog Lee (b) (),(b) Deprtment of Industril nd Systems Engineering, Kore Advnced Institute
More informationAccelerating 3D convolution using streaming architectures on FPGAs
Accelerting 3D convolution using streming rchitectures on FPGAs Hohun Fu, Robert G. Clpp, Oskr Mencer, nd Oliver Pell ABSTRACT We investigte FPGA rchitectures for ccelerting pplictions whose dominnt cost
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd business. Introducing technology
More informationModeling and Simulation of Short Range 3D Triangulation-Based Laser Scanning System
Modeling nd Simultion of Short Rnge 3D Tringultion-Bsed Lser Scnning System Theodor Borngiu Anmri Dogr Alexndru Dumitrche April 14, 2008 Abstrct In this pper, simultion environment for short rnge 3D lser
More informationNATIONAL UNIVERSITY OF SINGAPORE. Imaging FCS. A plugin for ImageJ. Thorsten Wohland August 17, 2015
NATIONAL UNIVERSITY OF SINGAPORE Imging FCS A plugin for ImgeJ Thorsten Wohlnd August 17, 015 Imging_FCS 1.45 is bsic ImgeJ plugin to clculte nd view sptio-temporl correltion functions from 16 bit GRAY
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd processes. Introducing technology
More informationIllumination and Shading
Illumintion nd hding In order to produce relistic imges, we must simulte the ppernce of surfces under vrious lighting conditions. Illumintion models: given the illumintion incident t point on surfce, wht
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd processes. Introducing technology
More informationSlides for Data Mining by I. H. Witten and E. Frank
Slides for Dt Mining y I. H. Witten nd E. Frnk Simplicity first Simple lgorithms often work very well! There re mny kinds of simple structure, eg: One ttriute does ll the work All ttriutes contriute eqully
More information9 Graph Cutting Procedures
9 Grph Cutting Procedures Lst clss we begn looking t how to embed rbitrry metrics into distributions of trees, nd proved the following theorem due to Brtl (1996): Theorem 9.1 (Brtl (1996)) Given metric
More information1 Drawing 3D Objects in Adobe Illustrator
Drwing 3D Objects in Adobe Illustrtor 1 1 Drwing 3D Objects in Adobe Illustrtor This Tutoril will show you how to drw simple objects with three-dimensionl ppernce. At first we will drw rrows indicting
More informationZZ - Advanced Math Review 2017
ZZ - Advnced Mth Review Mtrix Multipliction Given! nd! find the sum of the elements of the product BA First, rewrite the mtrices in the correct order to multiply The product is BA hs order x since B is
More informationFault injection attacks on cryptographic devices and countermeasures Part 2
Fult injection ttcks on cryptogrphic devices nd countermesures Prt Isrel Koren Deprtment of Electricl nd Computer Engineering University of Msschusetts Amherst, MA Countermesures - Exmples Must first detect
More informationEngineer To Engineer Note
Engineer To Engineer Note EE-186 Technicl Notes on using Anlog Devices' DSP components nd development tools Contct our technicl support by phone: (800) ANALOG-D or e-mil: dsp.support@nlog.com Or visit
More informationA Heuristic Approach for Discovering Reference Models by Mining Process Model Variants
A Heuristic Approch for Discovering Reference Models by Mining Process Model Vrints Chen Li 1, Mnfred Reichert 2, nd Andres Wombcher 3 1 Informtion System Group, University of Twente, The Netherlnds lic@cs.utwente.nl
More informationMidterm 2 Sample solution
Nme: Instructions Midterm 2 Smple solution CMSC 430 Introduction to Compilers Fll 2012 November 28, 2012 This exm contins 9 pges, including this one. Mke sure you hve ll the pges. Write your nme on the
More informationFall 2018 Midterm 1 October 11, ˆ You may not ask questions about the exam except for language clarifications.
15-112 Fll 2018 Midterm 1 October 11, 2018 Nme: Andrew ID: Recittion Section: ˆ You my not use ny books, notes, extr pper, or electronic devices during this exm. There should be nothing on your desk or
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd business. Introducing technology
More informationCSCI 446: Artificial Intelligence
CSCI 446: Artificil Intelligence Serch Instructor: Michele Vn Dyne [These slides were creted by Dn Klein nd Pieter Abbeel for CS188 Intro to AI t UC Berkeley. All CS188 mterils re vilble t http://i.berkeley.edu.]
More informationMath 17 - Review. Review for Chapter 12
Mth 17 - eview Ying Wu eview for hpter 12 1. Given prmetric plnr curve x = f(t), y = g(t), where t b, how to eliminte the prmeter? (Use substitutions, or use trigonometry identities, etc). How to prmeterize
More informationE/ECE/324/Rev.2/Add.112/Rev.3/Amend.1 E/ECE/TRANS/505/Rev.2/Add.112/Rev.3/Amend.1
6 August 2013 Agreement Concerning the Adoption of Uniform Technicl rescriptions for Wheeled Vehicles, Equipment nd rts which cn be itted nd/or be Used on Wheeled Vehicles nd the Conditions for Reciprocl
More informationA COLOUR CORRECTION PREPROCESSING METHOD FOR MULTIVIEW VIDEO CODING
A COLOR CORRECTO REROCESSG METHOD FOR MLTEW DEO CODG Colin Doutre nd nos siopoulos Deprtment of Electricl nd Computer Engineering, niversity of British Columbi 66 Min Mll, 6T Z4, ncouver, BC, Cnd emil:
More informationII. THE ALGORITHM. A. Depth Map Processing
Lerning Plnr Geometric Scene Context Using Stereo Vision Pul G. Bumstrck, Bryn D. Brudevold, nd Pul D. Reynolds {pbumstrck,brynb,pulr2}@stnford.edu CS229 Finl Project Report December 15, 2006 Abstrct A
More informationComputing offsets of freeform curves using quadratic trigonometric splines
Computing offsets of freeform curves using qudrtic trigonometric splines JIULONG GU, JAE-DEUK YUN, YOONG-HO JUNG*, TAE-GYEONG KIM,JEONG-WOON LEE, BONG-JUN KIM School of Mechnicl Engineering Pusn Ntionl
More informationMultiview hyperspectral topography of tissue structural and functional characteristics
Multiview hyperspectrl topogrphy of tissue structurl nd functionl chrcteristics Shiwu Zhng,b, Peng Liu, Jiwei Hung b, Ronld Xu,b Deprtment of Precision Mchinery nd Precision Instrumenttion, University
More informationApproximation by NURBS with free knots
pproximtion by NURBS with free knots M Rndrinrivony G Brunnett echnicl University of Chemnitz Fculty of Computer Science Computer Grphics nd Visuliztion Strße der Ntionen 6 97 Chemnitz Germny Emil: mhrvo@informtiktu-chemnitzde
More informationEngineer To Engineer Note
Engineer To Engineer Note EE-208 Technicl Notes on using Anlog Devices' DSP components nd development tools Contct our technicl support by phone: (800) ANALOG-D or e-mil: dsp.support@nlog.com Or visit
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd processes. Introducing technology
More informationToday. Search Problems. Uninformed Search Methods. Depth-First Search Breadth-First Search Uniform-Cost Search
Uninformed Serch [These slides were creted by Dn Klein nd Pieter Abbeel for CS188 Intro to AI t UC Berkeley. All CS188 mterils re vilble t http://i.berkeley.edu.] Tody Serch Problems Uninformed Serch Methods
More informationDynamic Programming. Andreas Klappenecker. [partially based on slides by Prof. Welch] Monday, September 24, 2012
Dynmic Progrmming Andres Klppenecker [prtilly bsed on slides by Prof. Welch] 1 Dynmic Progrmming Optiml substructure An optiml solution to the problem contins within it optiml solutions to subproblems.
More information1 Quad-Edge Construction Operators
CS48: Computer Grphics Hndout # Geometric Modeling Originl Hndout #5 Stnford University Tuesdy, 8 December 99 Originl Lecture #5: 9 November 99 Topics: Mnipultions with Qud-Edge Dt Structures Scribe: Mike
More informationThe Distributed Data Access Schemes in Lambda Grid Networks
The Distributed Dt Access Schemes in Lmbd Grid Networks Ryot Usui, Hiroyuki Miygi, Yutk Arkw, Storu Okmoto, nd Noki Ymnk Grdute School of Science for Open nd Environmentl Systems, Keio University, Jpn
More informationImproper Integrals. October 4, 2017
Improper Integrls October 4, 7 Introduction We hve seen how to clculte definite integrl when the it is rel number. However, there re times when we re interested to compute the integrl sy for emple 3. Here
More informationModel of the Human Eye Based on ABCD Matrix
Model of the Humn Eye Bsed on ABCD Mtrix G. Díz González nd M. Dvid Iturbe Cstillo Cittion: AIP Conf. Proc. 992, 108 (2008); doi: 10.1063/1.2926797 View online: http://dx.doi.org/10.1063/1.2926797 View
More informationMath 464 Fall 2012 Notes on Marginal and Conditional Densities October 18, 2012
Mth 464 Fll 2012 Notes on Mrginl nd Conditionl Densities klin@mth.rizon.edu October 18, 2012 Mrginl densities. Suppose you hve 3 continuous rndom vribles X, Y, nd Z, with joint density f(x,y,z. The mrginl
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd processes. Introducing technology
More information