Medical Image Analysis

 Karen Thompson
 10 months ago
 Views:
Transcription
1 Medicl Imge Anlysis 13 (29) Contents lists ville t ScienceDirect Medicl Imge Anlysis journl homepge: Modelling of imgectheter motion for 3D IVUS Misel Rosles, Peti Rdev, Oriol RodriguezLeor c, Deor Gil d, * Lortorio de Físic Aplicd, Fcultd de Ciencis, Deprtmento de Físic, Universidd de Los Andes Mérid, Venezuel Centre de Visió per Computdor, Universitt de Brcelon, Spin c Universitry Hospitl Germns Tries i Pujol, Bdlon, Spin d Centre de Visió per Computdor, Universitt Autonom de Brcelon, Spin rticle info strct Article history: Received 12 Jnury 27 Received in revised form 13 June 28 Accepted 19 June 28 Aville online 3 July 28 Keywords: Intrvsculr ultrsound (IVUS) Motion estimtion Motion decomposition Fourier Threedimensionl intrvsculr ultrsound (IVUS) llows to visulize nd otin volumetric mesurements of coronry lesions through n explortion of the cross sections nd longitudinl views of rteries. However, the visuliztion nd susequent morphogeometric mesurements in IVUS longitudinl cuts re suject to distortion cused y periodic imge/vessel motion round the IVUS ctheter. Usully, to overcome the imge motion rtifct ECGgting nd imgegted pproches re proposed, leding to slowing the pullck cquisition or disregrding prt of IVUS dt. In this pper, we rgue tht the imge motion is due to 3D vessel geometry s well s crdic dynmics, nd propose dynmic model sed on the trcking of n ellipticl vessel pproximtion to recover the rigid trnsformtion nd lign IVUS imges without loosing ny IVUS dt. We report n extensive vlidtion with synthetic simulted dt nd in vivo IVUS sequences of 3 ptients chieving n verge reduction of the imge rtifct of 97% in synthetic dt nd 79% in reldt. Our study shows tht IVUS lignment improves longitudinl nlysis of the IVUS dt nd is necessry step towrds ccurte reconstruction nd volumetric mesurements of 3D IVUS. Ó 28 Elsevier B.V. All rights reserved. 1. Introduction The introduction of intrvsculr ultrsound (IVUS) in the field of medicl imging (Yock et l., 1988; Grhm et l., 1989; Mintz et l., 21; Whle et l., 1999; Evns et l., 1996) s n explortory technique hs significntly chnged the understnding of rteril diseses of coronry rteries. Ech IVUS plne visulizes the crosssection (Fig. 1) of the rtery llowing, mong others the ssessment of different vessel plques (thromus of the rteriosclerotic plques nd clcium deposits) nd the determintion of morphogeometric prmeters (Metz et l., 1992, 1998; Jumo nd Rimund, 1998). 3D IVUS is otined y ssemling in stck 2D IVUS imges cquired during n utomtic pullck of the ctheter. 3D IVUS is of gret clinicl interest s it llows to evlute nd quntify the effect of therosclerotic plque on rteril distensiility (Ginnttsio et l., 21), to followup the plque nd lumen fter intervention (Nissen, 22), nd to guide the selection of n optiml interventionl procedure (Nissen nd Yock, 21). Crdic dynmics ( Mrio et l., 1995; Delchrtre et l., 1999; Roelndt et l., 1994; Thrush et l., 1997) nd 3D shpe of the vessel (Whle et l., 1999) introduce n imge mislignment tht cn e pprecited in the IVUS sequence s rottion nd displcement of the vessel crosssection. This is min rtifct tht hinders * Corresponding uthor. Emil ddress: (D. Gil). the longitudinl visuliztion of vessel morphology nd ny ccurte ssessment of morphologicl structures long the dt sequence. Their effect in 3D IVUS imges hs een descried (Whle et l., 1999; Ndkrni et l., 25; Di Mrio et l., 1993, 1995; Berry et l., 2; de Korte, 1999; Bruining et l., 1998; Dijkstr et l., 1998) s swtoothshped longitudinl ppernce of the vessel wll tht decreses precision in volumetric mesurements. Intrinsic rottion nd trnsltion of the imging ctheter introduce rtifcts in the longitudinl views mking structures to pper nd dispper s well s provoking toothshped dventiti nd clcified plque, etc. This motion minly ffects the computtion of the vessel dimeter in the Lviews. Fig. 1 shows the sinusoidl wvy pttern of the vessel wll due to the hert dynmics nd lood pulstion tht troules the vessel wll detection nd nlysis of vessel structures long the sequence. The most disturing phenomenon induced y vessel periodic motion is tht the displcement invlidtes ny point correspondence long the sequence s, given fixed ngle in the shortxis IVUS imge, points from different temporl frmes do not correspond to the sme vessel wll segment. On the other hnd, if we project the greylevel pixels in the temporl direction (longitudinl vessel direction) of the IVUS sequences, we cn oserve tht slient fetures usully representing the clcium nd firous plque nd dventitimedi orders descrie trjectory of semirottion of ring round stick (in this cse, the ctheter) (Fig. 1c). Tht is, in the shortxis /$  see front mtter Ó 28 Elsevier B.V. All rights reserved. doi:1.116/j.medi
2 92 M. Rosles et l. / Medicl Imge Anlysis 13 (29) c Shde zone Plque Ctheter Medi Distole Systole Longitudinl cut line reference Distole Intim Ctheter Ctheter Lumen Adventiti Guidewire rtifct Distole Systole Systole Vessel wll evolution Fig. 1. Morphologicl rteril structures nd rtifcts in n IVUS imge (). Sinusoidl shpe in longitudinl cut () nd greylevel shift in cross sectionl view (c). Red nd lue ellipses show the loction of vessel wll t systole nd distole, respectively. view, the rottion nd trnsltion of the ctheter led to high mount of vessel motion. A usul strtegy used to minimize the impct of the ove descried rtifcts is to synchronize sequence cquisition with hert dynmics. This cn e chieved y either specil ECGgted devices (Bruining et l., 1998) or y n imgesed ECGgted 3D IVUS (De Winter et l., 24; O Mlley et l., 27, 28). In the first cse, specil ultrsonic device performing the ctheter pullck t the sme systolic pek is required nd the cquisition time significntly increses. In the second cse, one discrds ll frmes etween two consecutive systolic peks, thus, running the risk of loosing vlule informtion. Keeping ll frmes nd ligning them y suppressing their rottion nd displcement would llow to get relile rdil mesurements during the whole cycle relted to the estimtion of the elstic properties of the vessel. Alterntively, reconstructing the pth using iplne ngiogrphy (Evns et l., 1996; Whle et l., 1999) could llevite the geometric distortions due to the position nd ngulr rottion of the ctheter y trcking the Frenet tringle ccording to the vessel tortuosity, lthough in this cse, imge rottion due to hert dynmics is not considered. The gol of this work is to estimte nd remove the rigid motion effect y using only IVUS imges keeping t the sme time ll the dt. The complex movement of the imging ctheter inside the coronry vessel cuses motion rtifcts tht cn e explined y three phenomen: () rottion of the IVUS imging ctheter with respect to the vessel (i.e. round the xis tngentil to the ctheter), () trnsltion of the imging ctheter in the plne perpendiculr to the ctheter xes, nd (c) trnsltion of the ctheter long its xis forwrd nd ckwrd cusing swinging effect. In this pper, we ddress the first two motion rtifcts in order to help the shortxis s longxis vessel disply, to visulize nd quntify different morphologicl vsculr structures like plque, stents, lumen dimeter, lesions, etc. We present hypothesis tht the vessel wll motion in IVUS imges is due to two min fctors: systemtic contriution cused y vessel wll geometry nd dynmicl periodic contriution due to the hert pulstile contrction influence. This decomposition of the vessel profiles serves to model the vessel dynmics nd determines roust procedure for vessel motion suppression. Aiming to correct rigid trnsformtion the vessel wll is pproximted y n ellipticl model tht is trcked in IVUS frmes to estimte its rottion nd trnsltion. The ellipticl pproximtion (Andrew et l., 1999) is computed on rough segmenttion of slient vessel structures using Neurl Network lgorithm (Ptrick, 1994). By pplying our stiliztion pproch, vessel ppernce in shortxes views llows etter nd esier perception of vessel structures distriution nd quntity. In longitudinl views we remove the toothshpe nd thus ssure tht chnges in vessel dimeter re due to chnges in morphology. In this wy, we ese the nlysis of plque progression long the vessel compred to the originl Lviews. Note tht this step is crucil for further nlysis of vessel iomechnics (plpogrphy, mesuring elstic properties s strin or stress, etc.) since their computtion requires trcking mteril points of the vessel long the crdic cycle. Recently there hs een n incresing interest in compensting motion (Dnilouchkine et l., 26; Leung et l., 26, 25) of the whole stck of IVUS imges in order to improve tissue elstogrphic studies. The existing pproches for motion compenstion in IVUS rely on either registrtion (Leung et l., 26) or trcking (Dnilouchkine et l., 26) of imge intensities. In Leung et l. (25) the uthors report severl techniques s follows: (1) Glol rottion lock mtching (GRBM), (2) contour mpping (CMAP), (3) locl lock mtching (LBM) nd (4) ctheter rottion nd trnsltion (CRT). These methods re sed on similrity estimtion of locl imge fetures defining ROI on frmes of interest. However, when the IVUS segment sequence does not possess ny chrcterizing structure, the similrity etween locks cn ecome very sensitive. In the GRBM cse, the minimizing lgorithm is prone to get stuck t locl minim unless exhustive (computtionlly inefficient) serch is performed. In the cse of trcking lgorithms, lrge displcements cnnot e properly modelled. In ny cse, chnges in imge intensity from one frme to the next one sustntilly ffect the performnce of the lgorithms. Our method is sed on registering the glol representtion of vsculr structures defined y the geometry of the coronry vessel. The proposed pproch does not require ny minimiztion process, since prmeters re given y n explicit formul. It follows tht our computtion of motion prmeters is not ffected y ny intensity chnges etween consecutive frmes. Furthermore, since points on vessel structures re computed seprtely on ech imge, our methodology nturlly hndles lrge displcements. An dvntge of our strtegy is tht it does not need very precise segmenttion of vessel structures in order to correct vessel motion. The set of reported experiments includes vessel motion simultions nd vlidtion of vessel dynmics suppression in reldt of 3 ptients. Our results revel suppression of rigid trnsformtion up to 97% in synthetic dt nd up to 85% in reldt. The pper is orgnized s follows: The theoreticl grounds of the model re given in Section 2: ssumptions on vessel dynmics re presented in Section 2.1 nd the mthemticl formultion in Section 2.2. Our pproch to computtion of prmeters is given in Section 3. Experiments re reported in Section 4 vlidting the model on phntom dt in Section 4.1 nd results on experimentl dt in Section 4.3. The rticle finishes with Discussions nd Conclusions.
3 M. Rosles et l. / Medicl Imge Anlysis 13 (29) A dynmicl model of vessel wll ppernce in IVUS 2.1. Model ssumptions The dynmics of coronry rtery is minly governed y the left ventricle dynmic evolution, lood pressure nd intrinsic geometric vessel properties (Mzumdr, 1992; Young nd Heth, 2; Ndkrni et l., 23; Holzpfel et l., 22). The first order pproximtion to vessel dynmics is given y liner trnsformtion comining trnsltion, rottion nd scling (Wks et l., 1996). The min ssumptions in the model we propose re the following: 1. Vessel trnsltion nd rottion hve two min contriutions: dynmic periodic motion nd systemtic geometric one. Vessel displcement long the sequence is reflected y the vessel wll profile in longitudinl cuts (Gil et l., 1996). Visul inspection of lrge longitudinl cuts (see Fig. 2) shows tht there is periodic wvy profile (lue nd red sinusoidl lines) nd systemtic contriution (green seline). We tke the temporl evolution of the lumen center long the sequence s the min descriptor of vessel displcement. The lumen centers were estimted s the centers of n ellipticl pproximtion to the vessel wll points. The sptil vrition of the lumen center, (Dx, Dy), is min geometric mesure tht provides with relevnt informtion on hert dynmics nd geometric contriutions to the IVUS displcement. Fig. 3 shows the sptil evolution of the lumen center horizontl coordinte, Dx. We note tht such profile decomposes into two min curves: crdic periodic oscilltion (lue line) nd se curve (red thick line) representing vessel wll position. The power spectrl density (Fig. 3) reflects these two min phenomen, s there re two predominnt frequency rnges. Low frequencies (in the rnge (1, 2) = (.2,.33) Hz) correspond to the geometric lumen evolution, while, those frequencies locted t (3,4) = (.78,.91) Hz re hert dynmics contriutions. By the ove considertions, we decompose lumen displcement into geometric, (Dx g,dy g ), nd dynmicl, (Dx d,dy d ), contriution terms: Dx g + Dx d nd Dy = Dy g + Dy d. The decoupled terms re given y the Fourier series: Dx g ðtþ ¼ Xn¼ 2 n¼ 1 ða n cosðnxtþþb n sinðnxtþþ; Using similr rgumenttions, we lso decouple the ngle of rottion, D, into geometric nd dynmicl terms: D g ðtþ Xn¼ 2 n¼ 1 ða n cosðnxtþþb n sinðnxtþþ; D d ðtþ Xn¼ 4 n¼ 3 ðc n cosðnxtþþd n sinðnxtþþ These Fourier coefficients give the hert dynmics nd geometric contriutions mplitudes to vessel motion nd ply centrl role in the simultion of vessel dynmic profiles (see experimentl Section 4). 2. Wll shpe evolution cn e descried y mens of the rottion nd trnsltion of n ellipticl model. In order to correct rigid motion rtifcts due to the reltive displcement of the vessel versus the imging ctheter, we restrict the vessel wll evolution to rigid trnsformtion, tht is, trnsltion followed y rottion. There re two possile wys of computing such trnsformtions: registrtion of IVUS dt sed on imge gry level (Lester nd Arridge, 1999) or correspondence of vessel structures/vessel order (Rosles et l., 25). Using n intensity sed imge registrtion runs the risks of ligning imges tking into ccount the most externl prt of the imge, which contin no structurl informtion on vessel orders ut re minly shdow nd right rtifct res. In order to void such mislignment, we trck the most slient vessel structures locted round the interfce intim, medi nd dventiti. To such purpose, points conforming to the gry vlue chrcteristics defining the trnsition lumentissue re extrcted y mens of neurl network vessel wll detection (see Section 3). Insted of serching for points correspondence etween two consecutive frmes, for the ske of roust nd efficient pproch, we extrct nd trck glol model of the set of points extrcted. Since we focus on glol liner rigid trnsformtion n elliptic model suits our purposes. The motiontrcking lgorithm will succeed s fr s the computed elliptic model keeps stle long the sequence since, in this cse, ny rtifct due to n unsuccessful fitting does not ffect the ellipticl vessel wll pproximtion. Dy g ðtþ Xn¼ 2 ða n cosðnxtþþb n sinðnxtþþ n¼ 1 Dx d ðtþ ¼ Xn¼ 4 ðc n cosðnxtþþd n sinðnxtþþ; n¼ 3 Dy d ðtþ Xn¼ 4 ðc n cosðnxtþþd n sinðnxtþþ n¼ 3 ð1þ ð2þ 2.2. Model formultion The ove ssumptions on vessel dynmics nd the ellipticl shpe pproximtion of vessel cross sections simplify the formultion of vessel motion to descriing rigid trnsformtion of n ellipticl vessel wll shpe. The vessel wll c k (u)=(x k (u),y k (u)) t time t k cn e written s follows: Cut C1 Vessel wll A Ctheter S1 S2 S3 S4 S5 S6 S7 Vessel wll B Fig. 2. Longitudinl cut of IVUS dt cquired during the ctheter pullck. Blue nd red points show the upper nd down vessel wll oundry. The green line shows the systemtic contriution of ctheter motion.
4 94 M. Rosles et l. / Medicl Imge Anlysis 13 (29) Δx in [mm] Amplitude in [u] Time in [sec] 2 1 Vessel Geometry contriutions Vessel Geometry contriutions Hert dynmic contriutions Hert dynmic contriutions Frequency in [Hz] Fig. 3. Temporl evolution of the lumen center Dx coordinte () nd its Fourier spectrum (). Red lines show vessel geometry contriution nd lue line illustrtes the hert dynmic contriution. x k ðuþ ¼ cosðhðuþþd k Þþcx k y k ðuþ ¼ sinðhðuþþd k Þþcy k where < h 6 2p determines the ngulr position of the corresponding point on the ellipse, u is the internl prmeter of the model curve, (,) re the minor nd mjor rdii of the ellipse, d k is its orienttion nd C k =(cx k,cy k ) is its center. We ssume tht the lumen center C k of c k cn e chnged from position (cx k,cy k ) t time t k to C k+1 with position (cx k+1,cy k+1 ) t time t k+1 due to the periodic hert movement. Therefore, the position of the center C k+1 of the new vessel shpe c k+1 t time t k+1 cn e written s cxðt kþ1 Þ¼cxðt k ÞþDx k ; cyðt kþ1 Þ¼cyðt k ÞþDy k where Dx nd Dy give the geometric nd dynmicl chnge induced y the vessel geometry nd hert pulstile contriution on the lumen center sptil nd temporl evolution, given y Eqs. (1) nd (2). When chnging to the new lumen center, new rotted vessel wll c k+1 =(x k+1,y k+1 ) is given, which coordintes hve chnged ccording to rottion of ngle k. Using the imge center s rottionl center, the new vessel wll coordintes cn e written s x kþ1 ðuþ cosð kþ1 Þ sinð kþ1 Þ xk ðuþ ¼ ð5þ y kþ1 ðuþ sinð kþ1 Þ cosð kþ1 Þ y k ðuþ Thus, if we know the rottion ngle k t ech time k, we re le to find the corresponding mteril point position in the next frme. In order to otin the rottion ngle we define s reference point, the point of the vessel model closest to the rottion center pproximted y the center of the imge. The ngle determined y the vectors of the reference point in ech pir of consecutive imges determine the rottion ngle for these frmes. Fig. 4 shows the reltion of geometric prmeters of our model. Without loss of generlity, to illustrte our pproch we use circulr vessel model. Fig. 4 is showing two vessel models c 1 nd c 2 with their corresponding centers C 1 nd C 2 nd rottion ngle D. Note tht k+1 cn e written s k+1 = k + D k, where k+1 = rctn(cy k /cx k ) nd D k ¼ rctnððcy k þ Dy k Þ=ðcx k þ Dx k ÞÞ: 3. A procedure for rigid motion estimtion nd suppression The procedure for rottion suppression splits into three min steps: ð3þ ð4þ ð6þ Ellipticl vessel pproximtion this step consists in vessel detection y Neurl Network trined to clssify slient structures of vessel order nd ellipticl pproximtion of the vessel wll. Motion prmeters computtion this step consists in estimtion of the rottion long the imge sequence using trcking procedure of reference point. Imge motion correction this step consists in correction of imge rottion nd displcement using the extrcted prmeters of the ellipticl pproximtion Ellipticl vessel pproximtion For the computtion of the vessel ellipticl representtion, structures on vessel wll long the sequence must e locted. The detection of the vessel wll is sed on physicl prmeters tht chrcterize the intim, medi nd dventiti lyers. Such descriptors re defined y glol nd locl imge greylevel properties otined from IVUS imges in polr form (Rosles et l., 24) with the origin t the ctheter center. In this coordinte system, the imge intensity, nmely I, depends on the pixel distnce (rdius r in polr coordintes) to the ultrsound trnsducer (Jensen, 21; Vogt et l., 1998). The decrese in the initil em intensity, I, is exponentilly proportionl to the sorption coefficient f, the frequency of the ultrsound, f, nd the prticle size or sctterer numer, N h, locted long the ultrsound em pth: IðrÞ ¼I expð fn h frþ: The sorption coefficient gives the rte of diminution with respect to the distnce long trnsmission pth (Jensen, 21) nd it is loclly otined from the regression line slope (Vogt et l., 1998) of the imge profile. Fig. 5 shows rdil greylevel profile of n IVUS imge, the lumenvessel order is pproximtely represented y the steepest greylevel trnsition tht is used to trin the Neurl Network. The set of chrcteristics chosen to find slient points on the vessel wll re the sorption coefficient f, imge pixel greylevel I(r) nd rdil stndrd men l nd devition r of the dt. The sorption coefficient gives locl informtion out the lumenvessel trnsition cting s n edge detector. The sttistics considered (l nd r) contin regionl informtion on vessel structures: l gives the seline of glol greylevel intensity nd r gives simple descriptor of the regionl texture. In order to find potentil cndidte structures on the vessel wll, Perceptron Multilyer Neurl Network (6:5:6:3) ws trined using stndrd Bck
5 M. Rosles et l. / Medicl Imge Anlysis 13 (29) Vessel wll γ 2 t time t 2 y P 2 (x,y ) 1 y P 2 (x,y ) Δα P 1 (x,y) Δα Δα P 1 (x,y) Δα C 2 C 2 Δy C 1 α 1 C 1 Ctheter (x o,y o ) Δx x (x,y ) Δα (x o,y o ) x C 1 (xc 1,yc 1 ) nd C 2 (xc 2,yc 2 ) Lumen centers position t time t1 nd t2 Vessel wll γ 1 t time t 1 (x,y ) Fig. 4. Geometric prmeters of the model (). Estimtion of the rottion ngle ()..8 Grey level Liner fitting to estimte the sorption coefficient (ζ) Rdil coordinte (R) Fig. 5. Vessel structures loction. Rdil greylevel profile () nd extrction of positive (+) nd negtive ( ) ptterns () to trin the neuronl network. Propgtion Algorithm (Ptrick, 1994). Fig. 5 shows trining exmples of positive (+) corresponding to intim nd negtive ( ) pttern, corresponding to lood, dventiti, shdows, nd rtifct zones. Once pplied the neuronl network to clssify slient points on the vessel order, n ellipticl pproximtion to the vessel wll points is djusted using the men squre ellipse fitting method descried in Andrew et l. (1999) Motion prmeters computtion Once defined the ellipse center s new origin, it is necessry to estimte the ngulr rottion profile of the vessel wll structure. In rottion, it is sufficient to provide the temporl evolution of single point on the vessel wll structure (Kittel et l., 1991; Broucke, 1978), to mesure the ngulr difference etween two consecutive frmes. Since the ultrsound intensity grdully decreses in the rdil direction, it follows tht one of the est visile structure points on the vessel wll is the point nerest to the imge center (Courtney et l., 22; Rosles et l., 24). The sptil loction of this reference point is determined s the positions ðx ; k yþ in k frme k nd ðx kþ1 ; y kþ1þ in frme k + 1 on the vessel structure (see Fig. 4) tht hve the miniml Eucliden distnce to the imge centers (cx k,cy k ) nd (cx k+1,cy k+1 ). The distnce criterion cn e written s ðx k ; y k Þ¼rgmin uððx o x k ðuþþ 2 þðy o y k ðuþþ 2 Þ 1=2 where x k (u) nd y k (u) re the coordintes of the points of the ellipticl vessel model. written s ðt k Þ¼ þ D k where ¼ rctnðy =x Þ is the reference ngle corresponding to the initil frme, nd D k rctnðy k =x kþ is the rottion ngle of frme k Imge motion correction The suppression of the rottion is given y the following liner trnsformtion. The current imge frme I k (x,y) is trnslted y (cx k,cy k ) to center the ellipse on the imge center, following y rottion through n ngle k : x k y k ¼ cosð k Þ sinð k Þ sinð k Þ cosð k Þ x k y k cx k cy k where ðx ; k y kþ nd (x,y) re the new nd old crtesin imge coordintes, k is the rottion ngle nd (cx k,cy k ) is the ellipticl center. By itertively pplying this eqution, we chieve to lign ll frmes within the IVUS sequence. 4. Vlidtion nd results As mentioned ove, our work is to use the geometric model to explin nd remove the discontinuous longitudinl vessel ppernce ð7þ ð8þ
6 96 M. Rosles et l. / Medicl Imge Anlysis 13 (29) in IVUS dt; hence, we pln the following set of experiments to vlidte the motion compenstion procedure: 1. Roustness of the estimtion of motion prmeters we check the impct of inccurcies in the estimting of the ellipsoidl shpe s well s the roustness of the estimtion of motion prmeters y mens of dynmic vessel phntom. 2. A roust mesure of rottion suppression in experimentl dt since in rel cses it is difficult to know in dvnce the expected motion prmeters, it is required to define first quntity properly ssessing the ccurcy in motion correction. The synthetic phntom helps us to define roust mesure of rottion suppression in terms of the vessel wll profile in longitudinl cuts. 3. Assessment of motion suppression in rel sequences we check the efficiency of the proposed technique y study of short sequences of 3 imges ech one extrcted from 3 ptients. 4. Stilized Vessel Model in order to illustrte the improvements provided y IVUS imge lignment we show the results of vessel reconstruction efore nd fter correction Roustness of the estimtion of motion prmeters The vessel phntom for motion correction is generted y pplying motion profile to synthetic IVUS imges. The Fourier decomposition given in Section 2.1 serves to compute the rottion (Eqs. (6) nd (7)) nd the vessel displcement (Eqs. (1) nd (2)) profiles. The result of dding hert dynmics contriutions to vessel geometry is the sinusoidl lue line nd the resulting geometric profile due to the ctheter pullck in tortuous vessel displyed y the red solid line in Fig. 6. We generted the vessel position corresponding to the different temporl frmes nd otined the vessel profiles in the simulted longitudinl IVUS cuts. Two different phntom models hve een used: n ellipticl shpe nd vessel profile extrcted from reldt. The ellipticl model llows to detect rtifcts in the geometric computtion of the rottion ngle. The vessel shpes extrcted from mnul segmenttion of IVUS sequences, ssess the roustness of the method ginst the elliptic modelling of vessel orders. For the ellipticl model we used 73 ellipses with eccentricities uniformly smpled in the rnge [, 1]. For the genertion of the reldt sed phntom, we used, for ech frme k only 25% of the vessel wll dt points rndomly selected. In this wy we cn check the model roustness to lck of informtion in disesed sections of the vessel. The error mesure used to ssess the ccurcy of the estimted rottion profiles is the solute difference etween ngles computed from the true synthetic profile. The correltion coefficients (the slope, m, nd the independent term, ) etween the estimted ngle nd the theoretic one hve een computed. In the idel cse of perfect correltion m = 1 nd =. Fig. 7 shows the regression line otined for n ellipse eccentricity with n = 37.5%, with slope m = 1.3 nd =.23. Fig. 7 shows the ngulr error s function of the eccentricity n. When the ctheter is ner to the lumen center ( 1 6 n 6 1) the liner coefficient m hs singulrity due to the equidistnce etween the ctheter center nd the vessel wll points. In these cses, the reference point s the closest point to the ctheter center, used to estimte the rottion ngle mkes no sense. In the cse of reldt sed phntoms (see Fig. 8), the liner coefficients present stle ehvior long the sequence, with sttisticl vlues for m equl to l 1.2, r.1. This represents pproximtely 1% of error in the estimtion of the rottion ngle which, tking into ccount tht only 25% of the vessel wll dt point hs een used, is n cceptle experimentl result. The liner coefficient gives the se line ngle etween theoreticl nd experimentl ngle. The sttisticl vlues of were l 5, r 4), hence, there is n cceptle se line shift etween theoretic nd rel ngle. The dvntges of incorporting vessel geometric informtion to the computtion of its rigid dynmics re reflected in the ccurcy chieved in synthetic phntoms. According to (Dnilouchkine et l., 26) only motion elow 3.5 cn e successfully modelled. For those rottions the men solute ngulr errors did not exceed 7.8. Our methodology supports lrge rottions in the rnge [ 25,25] degrees with n solute error of 5 ± 4. Errors sttisticl rnges etween theoretic nd estimted prmeters where extrcted from the regression coefficients (the slope m nd the independent term ) of the errors otined for ech vessel dynmic phntom. Since the verge slope is 1.2, the rnge for the independent term reflects the rnge for the solute error. Figs. 7 nd 8 show the regression lines for n elliptic (Fig. 7) nd rel vessel profile (Fig. 8) phntom models A roust mesure of rottion suppression Numericlly checking to wht extend the vessel wll rigid motion hs een correctly computed cn only e done in synthetic models, where the true displcement prmeters re known. In the cse of rel imges extrction of quntittive mesurement it is not possile nd one hs to se on motion correction in vessel wll visul ppernce. Vessel profiles in longitudinl cuts nd lignment of vessel slient structures re good cndidtes. After proper trnsltion nd rottion suppression (Pluim et l., 23), the vessel shpe ppernce of longitudinl cut is stright line, s oth swlike periodic ptterns nd vessel curvture due to n origin different from lumen center hve een removed (Seemntini et l., 25; De Winter et l., 24). We propose using the Eucliden distnce, Dd, etween the vessel profile nd its liner fitting s ojective mesure of the mount of motion removed. Fig. 9 illustrtes the mechnism followed to compute Dd nd the difference etween the dynmic synthetic profile (Fig. 9) nd its corrected counterprt (Fig. 9). The sttisticl rnge (given y the men ± stndrd devition) of Dd long the segment mesures the reduction in vessel movement. The men vlue of Dd reflects the overll performnce of the proposed suppression nd the stndrd devition mesures its roustness. Fig. 1 displys Dd histogrm efore (Fig. 1) nd fter (Fig. 1) rottion suppression for the reldt sed α in [degrees] 1 5 Vessel geometry contriutions Hert dynmics contriutions Time in [sec] Fig. 6. Rottion profile for synthetic simultion of vessel dynmics. Blue line shows the hert dynmics contriution nd the red line the vessel geometry contriution.
7 M. Rosles et l. / Medicl Imge Anlysis 13 (29) Estimted ngle (α) in [degrees] Liner coefficient m = 1.3 =.23 Error in [degress] Rel ngle (α) in [degrees] Ctheter eccentricity (ξ) in [%] Fig. 7. Rottion estimtion for elliptic phntom. Theoreticl vs. estimted rottion profile for n eccentricity n = 37.5%, (), nd error in ngle estimtion s function of the ellipse eccentricity, (). Fig. 8. Rottion estimtion for reldt phntoms. Rel (in lue) nd estimted (in red) rottion profile () induced to rel vessel wll points nd its liner correltion (). Vessel wll points Liner fitting Distnce (Δd) etween liner fitting nd vessel wll points Vessel wll points Liner fitting Distnce (Δd) etween liner fitting nd vessel wll points Fig. 9. Longitudinl cut efore () nd fter () rottion suppression. phntom. Before movement suppression, the distnce distriution, Dd A, is in the rnge (1.1 ±.5) mm, while fter correction the distnce, Dd B, is within (.4 ±.1) mm. This represents n verge movement reduction of 97% nd strighten of the sinusoidl vessel wll shpe induced y rottion nd trnsltion. The rtio, R, etween the distnce fter, Dd A, nd efore, Dd B, indictes the mount of movement suppressed. A vlue of zero is chieved when the vessel profile is stright line nd vlue ove 1 suggests tht the lgorithm hs fils to reduce motion. In experimentl dt, rdil diltion nd vessel curvture prevent the rtio R chieving its optiml lower ound even if rigid motion hs een properly removed. In fct, in short segments, rdil diltion is the min source of devition from stright pttern. According to clinicl studies (Willims et l., 1999; Dodge et l., 1992) lood pressure introduces n verge diltion of 14% of the vessel dimeter. In order to tke into ccount this seline elstic deformtion, we define rigid rtio s Rr ¼ jdd A :14 VesselDimj Dd B for VesselDim the men dimeter of the vessel segment under study. ð9þ
8 98 M. Rosles et l. / Medicl Imge Anlysis 13 (29) Count Distnce Δd in [mm] Count Distnce Δd in [mm] Fig. 1. Histogrm of Dd efore () nd fter () rottion suppression Assessment of motion suppression in rel sequences In order to study the ehvior of the reported motion suppression strtegy we nlyzed set of sequences of 3 imges ech one for 3 ptients independent from the set used to trin the Neurl Network. The sequences were cquired using Boston Sci. equipment Boston Scientific Corportion, 1998 t 4 MHz t constnt pullck speed of.5 mm/sec. The processing time cn e divided in four steps (see Section 3): 1. Neuronl network trining is reltively fst, since the min chrcteristics difference etween vessel wll nd the rest is notly lrge, nerly 5 min (using PC Pentium t 2.4 GHz). 2. Vessel wll detection y Neuronl Network (.8 s/imge) 3. Vessel wll ellipticl fitting (.1 s/imge). 4. Correction of imge rottion nd displcement using the extrcted prmeters of the ellipticl pproximtion (.1 s/ imge). Longitudinl cuts efore nd fter rottion elimintion re shown in Fig. 11. One cn notice the smoothing effect of the rottion correcting on the ppernce of the lumen order (red line in Fig. 11) in the longitudinl view. Still complete smooth ppernce of the vessel in the longitudinl imges is not expected due to the different degree of rdil expnsion/contrction of the vessel (Fig. 11d nd f) ccording to its elstic properties (since we remove/correct only rigid vessel motion effects). We should lso tke into ccount tht y ligning vessel lumens, we hve n opticl rtifct due to the externl dventiti right tissue nd the ctheter which dopt wvy swshpe ppernce (see Fig. 11). The structures in longitudinl views tht reflect correct lignment of vessel IVUS plnes, re the ones ordered y the red lines tht show the trnsition etween lood nd tissue. Correction of center position strightens the vessel profile (see Fig. 11 c), removing most of the swshpe effect if present (see Fig. 11d f). Rottion correction produces n lignment of structures which is etter pprecited in the presence of slient plque such s clcium (ottom of Fig. 11). The periodic vessel rottion produces shdowright nded profile t clcium sectors which disppers converting it into uniform right or drk nd (s it is the cse of Fig. 11) menwhile the rottion is suppressed. As for the geometric systemtic rottion induced y the vessel geometry, it results in continuous vrition of slient structures which do not suddenly dispper ut smoothly chnge long the longitudinl cut (see Fig. 11e). Our sttisticl nlysis focused on checking the roustness in prmeter estimtion s well s reduction in movement. Roustness in the estimtion procedure is ssessed y computing prmeters of originl nd corrected imges. Intuitively, we cn ssume tht the less vrition in prmeters extrcted from utomticlly computed still imges, the more relile the reported geometric motion estimtion is. The stndrd devition, r, long ech sequence indictes such vriility for ech prmeter. The stndrd devition in the rdil position of the lumen centers efore trnsformtion correction presents significtive decrement for ll cses from Dr (2.3 ± 3) mm to Dr (.5 ±.3) mm fter imge lignment. In the cse of rottion estimtion,wehvereductionfromd (69 ± 2) efore to d (11 ± 8) degrees fter rottion suppression tht supposes 85% of rottion suppression. Fig. 12 shows the stndrd devition of motion prmeters for the 3 ptients nlyzed, lumen center position in Fig. 12 nd rottion ngle in Fig. 12. The quntittive ssessment of motion reduction is given y the rigid rtio Rr defined y (9) nd using tht the verge vessel dimeter in our dt is pproximtely 2.5 mm. Fig. 13 shows the rigid rtio Rr nd its histogrm in percentges (Fig. 13) for the 3 cses. The rnge for Rr in percentges is 79.3 ± nd in 76% of the cses is over 72% Stilized vessel model In order to otin relistic model of the vessel shpe, we must seprte the geometric vessel properties from hert dynmics contriutions. Let c(t) =(x(t), y(t), z(t)) e the vessel order s function of time t written s n ellipticl pproximtion: xðtþ ¼ðtÞ cosðh þ dðtþþ þ cxðtþ; yðtþ ¼ðtÞ sinðh þ dðtþþ þ cyðtþ; zðtþ ¼czðtÞ ð1þ where ( < h 6 2p), ((t),(t)) re the minor nd mjor rdii of the ellipse, d(t) is its orienttion nd C =(cx(t), cy(t), cz(t)) its center. In 2.5D IVUS reconstruction cz(t) is unknown, ut in our model we cn tke it s the ctheter position in longitudinl direction, cz = vt, where v is the ctheter velocity. Tking into ccount the decoupling into geometricl nd hert dynmics contriutions given in Section 2.1, the vessel wll prmeters cn e written s ðtþ ¼ g ðtþþ d ðtþ; ðtþ ¼ g ðtþþ d ðtþdðtþ ¼d g ðtþþd d ðtþ cxðtþ ¼cx g ðtþþcx d ðtþ; cyðtþ ¼cy g ðtþþcy d ðtþ where the suscript g stnds for the geometric component nd the d for the periodic hert movement. In order to reconstruct the vessel wll from ctheter point of view only the Fourier coefficient tht corresponds to geometricl contriution must e tken into ccount, so tht Eq. (1) reduces to
9 M. Rosles et l. / Medicl Imge Anlysis 13 (29) Fig. 11. Longitudinl cut efore (( f) up) nd fter (( f) down) rottion correction Ellipses center efore Before rottion suppression std in [mm] Ellipses center fter std in [degrees] After rottion suppression Ptient numer Ptient numer Fig. 12. Stndrd devition of motion prmeters efore (lue) nd fter (red) rottion suppression: ellipses centers () nd rottion ngle (). xðtþ ¼ g ðtþ cosðh þ d g ðtþþ þ c xg ðtþ; yðtþ ¼ g ðtþ sinðh þ d g ðtþþ þ c yg ðtþ; zðtþ ¼vt Fig. 14 shows sttic vessel reconstruction of n IVUS pullck of 19 imges corresponding to 25 mm long vessel segment. Fig. 14 shows the 2.5D grphicl representtion of the segment efore
10 1 M. Rosles et l. / Medicl Imge Anlysis 13 (29) Rr rtio Ptient numer Fig. 13. Movement reduction in experimentl dt: reduction rtio for the 3 ptients () nd histogrm of reduction percentge (). Fig. 14. Sttic vessel model: two views of 2.5D vessel wll reconstruction efore () nd fter () dynmic suppression nd the corresponding longitudinl IVUS cuts efore (c) nd fter (d) motion suppression. () nd fter () dynmics suppression. Fig. 14c nd d illustrte the corresponding longitudinl IVUS cuts efore nd fter motion compenstion. 5. Discussions 5.1. Ellipticl fitting of vessel orders When considering the prolem of imge registrtion nd lignment, nturl question is which generl strtegy to use: imgesed (where ll pixels of the imges contriute to the estimte of the imge trnsformtion) or lndmrksed (only some, usully hving specil mening, points of the imges re tken into ccount in order to lign the imges). Our pproch follows the second strtegy ligning the elliptic pproximtions to set of points corresponding to the lumentissue trnsition. Although n ellipsoidl shpe cn not exctly pproximte the vessel wll shpe, such inccurcies do not ffect the performnce of the proposed pproch since its purpose is to estimte the rigid trnsformtion of the vessel. By using n ellipticl pproximtion of the vessel represents glol pproch to detect motion nd rottion of the vessel; t the sme time the ellipticl pproximtion is le to etter cope with smll errors in the clssifiction of the contour points y the neurl network (compred to flexile snke model or other locl deformle models). Empiriclly, we found tht in 88% of segmented imges the ellipticl shpe model could pproximte the vessel eing the rest of imges minly cses of ifurctions to min vessels nd lrge NURD rtifcts. In order to evlute the qulity of the vessel pproximtion, we estimted distnce mp of the found contour points y the neurl network nd pproximted the vessel points y the ellipticl model. We consider tht the ellipticl model pproximtion is correct if the vlue of its points on the distnce mp is less thn given epsilon. In our cse, epsilon ws set empiriclly to.2 mm. Fig. 16 shows the men error nd stndrd devition of estimted distnces s well s their histogrms. In cse the distnce surpsses this threshold, hypothesis for segmenttion filure is generted nd the frme corresponding to this prt of the vessel is disregrded from nlysis nd motion estimtion is
11 M. Rosles et l. / Medicl Imge Anlysis 13 (29) Fig. 15. Segmented vessel wll y neurl network (in yellow) nd fitted ellipse (in lue). Fig. 16. Men error ( k ) () its stndrd devition () etween segmented points y the neuronl network nd the djusted ellipse. The corresponding histogrm re given in (c) nd (d). computed ligning the next vlid frme. Finlly, dynmic prmeters of the filed frme re computed y interpoltion of the neighor vlid frmes. Ellipticl fitting qulity cn e visully checked in the imges of Fig. 15, where we show the points extrcted from the neurl network (in yellow) nd the fitted ellipse (in lue) for six consecutive frmes. From our tests on the 3 ptients IVUS videos we conclude the following: () Our methodology is not ffected y clcified plque presence s it presents n rupt trnsition from lumen to tissue which is esily identified y the neurl network (see Fig. 17).
12 12 M. Rosles et l. / Medicl Imge Anlysis 13 (29) Fig. 17. Performnce in the presence of rtifcts: clcifiction (), guide rtifct () nd ifurctions (c). () Imging rtifcts, such s sensor guide, constitute sprse phenomen leding to locl responses in the shortxes imge view tht do not ffect the ellipse pproximtion (see Fig. 17). (c) In the cse of NURDS, such sequences re of limited clinicl interest y their low qulity (impossiility to quntify vessel structures) nd re out of the scope of our nlysis, neither recommended to e quntified in clinicl studies. Still, the procedure for ellipticl fitting is vlid provided tht the NURD rtifct does not ffect more thn 3% of the vessel Motion rtifct removl The complex motion of the imging ctheter inside the coronry vessel cuses motion rtifcts tht cn e pproximted y three phenomen: () rottion of the IVUS imging ctheter with respect to the vessel (i.e., round the xis tngentil to the ctheter), () trnsltion of the imging ctheter in the plne perpendiculr to the ctheter xes, nd (c) trnsltion of the ctheter long its xis forwrd nd ckwrd cusing swinging effect. In this pper, we ddress the first two motion rtifcts in order to help the shortxis s longxis vessel disply, to visulize nd quntify different morphologicl vsculr structures like plque, stents, lumen dimeter, lesions, etc. The intrinsic rottion nd trnsltion of the imging ctheter (relted to type () nd ()) introduce rtifcts in the longitudinl views mking structures to pper nd dispper s well s provoking toothshped dventiti nd clcified plque, etc. This motion minly ffects the computtion of the vessel dimeter in the Lviews. In the shortxis view, the rottion nd trnsltion of the ctheter led to high mount of vessel motion. By pplying our stiliztion pproch, vessel ppernce in shortxes views llows etter nd esier perception of vessel structures distriution nd quntity. In longitudinl views we remove the toothshpe nd thus ssure tht chnges in vessel dimeter re due to chnges in morphology. In this wy, we improve the nlysis of plque progression long the vessel compred to the originl Lviews. Note tht this step is crucil for further nlysis of vessel iomechnics (plpogrphy, mesuring elstic properties s strin or stress, etc.) since their computtion requires trcking mteril points of the vessel long the crdic cycle. Relile length mesurements long the vessel must either ccount (i.e. correct) for ctheter swinging or simulte imgesed ECGgting. We note tht the crdic components of the rottionl ngle cn serve to simulte n ECGgting y smpling the sequence t the right (crdic) rte. Therefore, lthough we do not explicitly ddress it, our estimtion of dynmics enles relile volumetric mesurements s well s exploring vessel (continuous) elstic properties y keeping the originl informtion (complete set of dt O Mlley et l., 26) Our pproch nd IVUS gting In this pper we del with the prolem of rottion suppression of IVUS imges tht cn e oserved in non ECGgted s well s ECGgted sequences. Thus, our pproch is not depending nd completely pplicle to ECGgted procedures (sed on hrdwre or imgegting). ECGgting hs the dvntge to llow more precise longitudinl mesurements; the price to e pid consists of incresing the cquisition time in up to 7 times with respect nongted motorized pullcks t.5 mm/seg. (von Birgelen et l., 1997), which results in n increse of the intervention time of the ptient. Besides the specil cquisition devices re not commonly ville t ll medicl centers. The ove resons hve recently motivted the development of imgesed ECGgted IVUS (Seemntini et l., 25; De Winter et l., 24). Such techniques discrd prt of the informtion (frmes) ville in n IVUS sequence. Since in our pproch to lign imges we do not put ny constrint on the pir of imges, our pproch llows using the whole stck of IVUS imges provided y the cquisition device s well s the ECGgted suset of imges in order to overcome the mislignment rtifcts introduced y vessel rottion nd trnsltion. Current ECGgting s well imgegting pproches still cn suffer from the imge rottion rtifcts mking difficult to ddress prolems s plpogrphy estimtion, plque following long the vessel, volumetric mesurements, etc. Moreover, lthough we do not explicitly ddress it, we note tht the crdic components of the rottionl ngle cn serve to simulte n ECGgting y smpling the sequence t the right (crdic) rte. Therefore, our estimtion of dynmics enles relile volumetric mesurements s well s exploring vessel (continuous) elstic properties y keeping the originl informtion (complete set of dto Mlley et l., 26) Ctheter oscillting oliquity Crdic dynmics introduce two min rtifcts ffecting the ctheter trjectory. Firstly, due to the chnge in vessel curvture during the crdic cycle, the ctheter moves long the xis of the vessel (longitudinl displcement). This forwrd nd ckwrd displcement implies tht, in mechnicl pullck, IVUS imges re not uniformly distriuted long the sequence zxis. Secondly, tilting of the ctheter introduces n nisotropic scling of IVUS cross sections if the ctheter is not coxil within the vessel. Such phenomenon produces more ellipticl shped vessel, which might result in n overestimtion of the re nd misinterprettion of the vessel shpe (Mintz et l., 21).
13 M. Rosles et l. / Medicl Imge Anlysis 13 (29) Correction of ctheter longitudinl movement would require dynmic model of the ctheter 3D trjectory nd is out of the scope of the presented work. Tilting of the ctheter crosssectionl plne cn not e properly corrected unless complete 3D reconstruction nd modelling of the vessel nd the ctheter inside the vessel re ville. Since we im t improving 3D IVUS without ffecting vessel mesurements such s length or re, ctheter oliquity will e disregrded Rdil deformtion of the coronry vessel Rdil deformtion of crdic rteries is minly due to lood pressure nd vessel wll elstic properties. By Hooke s lw (Mzumdr, 1992; Ndkrni et l., 23; Holzpfel nd Weizsäcker, 1998; Humphrey, 1995), the rdil increment, rr, is proportionl to the grdient of lood pressure, rp vi the reltion (Mzumdr, 1992; Ndkrni et l., 23; Holzpfel nd Weizsäcker, 1998; Humphrey, 1995): rr ¼ðjDP=pÞ 1=2 ð11þ where j is the elsticity coefficient. At vessel sections presenting non homogenous plque, the elsticity coefficient vries for ech pixel. In order to estimte the rdil deformtion, it suffices to use constnt elstic coefficient (Hook s lw for uniform medi) descriing the verge deformtion rnges of the vessel wll. According to clinicl studies (Willims et l., 1999), j is in the rnge j = (.1 ±.2) mm 2 /mmhg nd DP 4 mmhg. By Hook s lw (11), it follows tht rr.35 mm. Tking into ccount tht the rdii of coronry segments (Dodge et l., 1992) is in the rnge r = 2.64 ±.3 mm, we hve tht the reltive rdil deformtion induced y lood pressure scling is dr =(rr/r) 1 13%. Applying scling of the vessel in the IVUS imges might remove the deformtion due to the lood pumping tht is directly relted to the elsticity of the vessel. Besides the existing rdil scling does not ffect the computtion of the trnsltion nd rottion prmeters. Therefore, rdil scling is not included into the model formultion Clinicl pplicility The methodology proposed nd its implementtion improve the visuliztion in shortview IVUS videos s well s in the longitudinl views with success of correctly ligned sequences. Currently, the methodology is to e implemented in pckge for IVUS nlysis nd plque chrcteriztion to e instlled in clinicl conditions. On the other hnd, the proposed methodology for rottion suppression hs different possile clinicl pplictions: 1. Generl mesures extrction from stilized IVUS sequence. We consider tht this ppliction is first nd necessry step towrds: () Estimtion of geometric prmeters such s vessel length using longitudinl cuts views sed on methodology pplying the ECGgting to ddress the swinging rtifcts of the IVUS ctheter. () Vessel wll strin estimtion y plpogrphy. In order to otin vlid strin estimtion, we consider tht plpogrphy extrction needs: first, lignment of imges nd second, determining the frmes corresponding to endsystole nd enddistole peks (eing estimted y ECGsignl or imge nlysis). 2. Hert dynmic estimtion. Being possile to seprte the geometric contriution from dynmics contriution, we developed the sis to 2.5D vessel reconstruction only using IVUS dt. This method cn give promising evidences for vessel wll dynmics estimtion such s the introduction of n lterntive technique to estimte locl hert dynmics. In this wy, we provide new possiility of studying the vessel dynmics nd geometry s well s estlishing new dignostic tools. For exmple, n extensive study of IVUS rottion profile versus pumping hert efficiency otined experimentlly using distolic nd systolic ngiogrphy hert views, cn hint positive correltion etween hert rottion nd pumping hert efficiency. 6. Conclusions In this rticle we developed geometric nd kinemtic model in order to study the evolution of coronry rtery wll. The model is sed on the ssumption tht the evolution of the rteril wll cn e modeled ssuming two principl contriutions tht come from different physicl resons. The first one, systemtic contriution cused y geometric intrinsic rteril properties, nd the second one, n oscillting contriution tht comes from ventricle dynmics. These contriutions govern in mjor degree the profiles ppernce of rteril wll in longitudinl views. Using this ssumption, we propose methodologicl strtegy in order to estimte nd suppress rigid IVUS dynmicl distortions. IVUS imge lignment is very importnt in prolems tht study vessel deformtion long the crdic cycle, such s plpogrphy which compres deformtions t distole nd systole. Exploring vessel wll kinemtics needs vessel order points trcking in order to extrct vessel deformtion nd to judge out coronry elsticity. Definitely, imge rottion oserved in current IVUS sequences hinders the trcking of the corresponding points/zones of the vessel order. On the other hnd, the ccurte rdil deformtion of the vessel on the imge not only depends on the mutul imge/ vessel rottion, ut lso on the rtifct produced y the oscillting longitudinl oliquity tht is induced y the ventricle pulstile dynmics. Usully, oscillting longitudinl ctheter motion is overcome y ECGgted or imgegted sequence construction tht improves significntly the longitudinl view of the vessel, ut cquires picture of the vessel shpe only in distole (De Winter et l., 24). Looking for lterntive technicl pproches to solve the prolem of imge sptil loction covering the whole crdic cycle is nother reserch opportunity to dvnce in the more ccurte longitudinl visuliztion, vessel rdil deformtion nlysis, vessel/plque elsticity estimte, prediction of therosclerotic ccumultion nd possile plque ruptures from IVUS dt. Seprting the geometric contriution from dynmicl contriution cn e the fundmentl sis to 2.5D vessel reconstruction only using IVUS dt nd could e n importnt dvnce in vessel wll dynmics estimtion tht introduces n lterntive technique to estimte locl hert dynmics. In this wy, we provide new possiility of studying roustly the vessel dynmics, creting sttisticl models of vessel dynmics nd estlishing new dignostic tools. Acknowledgement We would like to thnk Crlos Rodríguez from the Mthemticl Dep. of the Universitt Autonom de Brcelon for his support nd dvise. This work ws supported in prt y reserch grnt from projects TIN C2, FISPI6129, PI7454 nd PI71188 nd CONSOLIDERINGENIO 21 (CSD2718). The lst uthor hs een supported y The Rmón y Cjl Progrm. References Andrew, F., Murizio, P., Roert, B., Direct lest squre fitting of ellipses. IEEE Trns. Pttern Anl. Mch. Intell. 21 (5), Berry, E., Kelly, S., Hutton, J., Lindsy, H., et l., 2. Intrvsculr ultrsoundguided interventions in coronry rtery disese, NHS RD HTA Progrmme.
Chapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationA027 Uncertainties in Local Anisotropy Estimation from Multioffset VSP Data
A07 Uncertinties in Locl Anisotropy Estimtion from Multioffset VSP Dt M. Asghrzdeh* (Curtin University), A. Bon (Curtin University), R. Pevzner (Curtin University), M. Urosevic (Curtin University) & B.
More information4 VECTORS. 4.0 Introduction. Objectives. Activity 1
4 VECTRS Chpter 4 Vectors jectives fter studying this chpter you should understnd the difference etween vectors nd sclrs; e le to find the mgnitude nd direction of vector; e le to dd vectors, nd multiply
More informationTHERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION
XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es
More informationAnalytically, vectors will be represented by lowercase boldface Latin letters, e.g. a, r, q.
1.1 Vector Alger 1.1.1 Sclrs A physicl quntity which is completely descried y single rel numer is clled sclr. Physiclly, it is something which hs mgnitude, nd is completely descried y this mgnitude. Exmples
More information2. VECTORS AND MATRICES IN 3 DIMENSIONS
2 VECTORS AND MATRICES IN 3 DIMENSIONS 21 Extending the Theory of 2dimensionl Vectors x A point in 3dimensionl spce cn e represented y column vector of the form y z zxis yxis z x y xxis Most of the
More informationa * a (2,1) 1,1 0,1 1,1 2,1 hkl 1,0 1,0 2,0 O 2,1 0,1 1,1 0,2 1,2 2,2
18 34.3 The Reciprocl Lttice The inverse of the intersections of plne with the unit cell xes is used to find the Miller indices of the plne. The inverse of the dspcing etween plnes ppers in expressions
More informationSection 4: Integration ECO4112F 2011
Reding: Ching Chpter Section : Integrtion ECOF Note: These notes do not fully cover the mteril in Ching, ut re ment to supplement your reding in Ching. Thus fr the optimistion you hve covered hs een sttic
More informationu( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 218, pp 4448): Determine the equation of the following graph.
nlyzing Dmped Oscilltions Prolem (Medor, exmple 218, pp 4448): Determine the eqution of the following grph. The eqution is ssumed to e of the following form f ( t) = K 1 u( t) + K 2 e!"t sin (#t + $
More informationP 1 (x 1, y 1 ) is given by,.
MA00 Clculus nd Bsic Liner Alger I Chpter Coordinte Geometr nd Conic Sections Review In the rectngulr/crtesin coordintes sstem, we descrie the loction of points using coordintes. P (, ) P(, ) O The distnce
More informationLinear Systems with Constant Coefficients
Liner Systems with Constnt Coefficients 4305 Here is system of n differentil equtions in n unknowns: x x + + n x n, x x + + n x n, x n n x + + nn x n This is constnt coefficient liner homogeneous system
More informationSection 14.3 Arc Length and Curvature
Section 4.3 Arc Length nd Curvture Clculus on Curves in Spce In this section, we ly the foundtions for describing the movement of n object in spce.. Vector Function Bsics In Clc, formul for rc length in
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More information#6A&B Magnetic Field Mapping
#6A& Mgnetic Field Mpping Gol y performing this lb experiment, you will: 1. use mgnetic field mesurement technique bsed on Frdy s Lw (see the previous experiment),. study the mgnetic fields generted by
More informationTests for the Ratio of Two Poisson Rates
Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06073 HIKARI Ltd, www.mhikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationMA123, Chapter 10: Formulas for integrals: integrals, antiderivatives, and the Fundamental Theorem of Calculus (pp.
MA123, Chpter 1: Formuls for integrls: integrls, ntiderivtives, nd the Fundmentl Theorem of Clculus (pp. 27233, Gootmn) Chpter Gols: Assignments: Understnd the sttement of the Fundmentl Theorem of Clculus.
More informationIntensity transformations
Intensity trnsformtions Stefno Ferrri Università degli Studi di Milno stefno.ferrri@unimi.it Methods for Imge Processing cdemic yer 2017 2018 Sptil domin The sptil domin of n imge is the plne tht contins
More informationexpression simply by forming an OR of the ANDs of all input variables for which the output is
2.4 Logic Minimiztion nd Krnugh Mps As we found ove, given truth tle, it is lwys possile to write down correct logic expression simply y forming n OR of the ANDs of ll input vriles for which the output
More informationMACsolutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences  Finite Elements  Finite Volumes  Boundry Elements MACsolutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More information5: The Definite Integral
5: The Definite Integrl 5.: Estimting with Finite Sums Consider moving oject its velocity (meters per second) t ny time (seconds) is given y v t = t+. Cn we use this informtion to determine the distnce
More informationdx dt dy = G(t, x, y), dt where the functions are defined on I Ω, and are locally Lipschitz w.r.t. variable (x, y) Ω.
Chpter 8 Stility theory We discuss properties of solutions of first order two dimensionl system, nd stility theory for specil clss of liner systems. We denote the independent vrile y t in plce of x, nd
More information1 Nondeterministic Finite Automata
1 Nondeterministic Finite Automt Suppose in life, whenever you hd choice, you could try oth possiilities nd live your life. At the end, you would go ck nd choose the one tht worked out the est. Then you
More informationConducting Ellipsoid and Circular Disk
1 Problem Conducting Ellipsoid nd Circulr Disk Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 (September 1, 00) Show tht the surfce chrge density σ on conducting ellipsoid,
More informationThis chapter will show you. What you should already know. 1 Write down the value of each of the following. a 5 2
1 Direct vrition 2 Inverse vrition This chpter will show you how to solve prolems where two vriles re connected y reltionship tht vries in direct or inverse proportion Direct proportion Inverse proportion
More information5.2 Volumes: Disks and Washers
4 pplictions of definite integrls 5. Volumes: Disks nd Wshers In the previous section, we computed volumes of solids for which we could determine the re of crosssection or slice. In this section, we restrict
More informationChapter 9 Definite Integrals
Chpter 9 Definite Integrls In the previous chpter we found how to tke n ntiderivtive nd investigted the indefinite integrl. In this chpter the connection etween ntiderivtives nd definite integrls is estlished
More informationFORM FIVE ADDITIONAL MATHEMATIC NOTE. ar 3 = (1) ar 5 = = (2) (2) (1) a = T 8 = 81
FORM FIVE ADDITIONAL MATHEMATIC NOTE CHAPTER : PROGRESSION Arithmetic Progression T n = + (n ) d S n = n [ + (n )d] = n [ + Tn ] S = T = T = S S Emple : The th term of n A.P. is 86 nd the sum of the first
More informationPrecalculus Spring 2017
Preclculus Spring 2017 Exm 3 Summry (Section 4.1 through 5.2, nd 9.4) Section P.5 Find domins of lgebric expressions Simplify rtionl expressions Add, subtrct, multiply, & divide rtionl expressions Simplify
More informationSynoptic Meteorology I: Finite Differences September Partial Derivatives (or, Why Do We Care About Finite Differences?
Synoptic Meteorology I: Finite Differences 1618 September 2014 Prtil Derivtives (or, Why Do We Cre About Finite Differences?) With the exception of the idel gs lw, the equtions tht govern the evolution
More information4.1. Probability Density Functions
STT 1 4.14. 4.1. Proility Density Functions Ojectives. Continuous rndom vrile  vers  discrete rndom vrile. Proility density function. Uniform distriution nd its properties. Expected vlue nd vrince of
More informationTriangles The following examples explore aspects of triangles:
Tringles The following exmples explore spects of tringles: xmple 1: ltitude of right ngled tringle + xmple : tringle ltitude of the symmetricl ltitude of n isosceles x x  4 +x xmple 3: ltitude of the
More informationAcceptance Sampling by Attributes
Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire
More informationMeasuring Electron Work Function in Metal
n experiment of the Electron topic Mesuring Electron Work Function in Metl Instructor: 梁生 Office: 7318 Emil: shling@bjtu.edu.cn Purposes 1. To understnd the concept of electron work function in metl nd
More informationOrdinary differential equations
Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview IntroductionModelling Bsic concepts to understnd n ODE. Description nd properties
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationAlg. Sheet (1) Department : Math Form : 3 rd prep. Sheet
Ciro Governorte Nozh Directorte of Eduction Nozh Lnguge Schools Ismili Rod Deprtment : Mth Form : rd prep. Sheet Alg. Sheet () [] Find the vlues of nd in ech of the following if : ) (, ) ( 5, 9 ) ) (,
More informationPurpose of the experiment
Newton s Lws II PES 6 Advnced Physics Lb I Purpose of the experiment Exmine two cses using Newton s Lws. Sttic ( = 0) Dynmic ( 0) fyi fyi Did you know tht the longest recorded flight of chicken is thirteen
More informationUsing QM for Windows. Using QM for Windows. Using QM for Windows LEARNING OBJECTIVES. Solving Flair Furniture s LP Problem
LEARNING OBJECTIVES Vlu%on nd pricing (November 5, 2013) Lecture 11 Liner Progrmming (prt 2) 10/8/16, 2:46 AM Olivier J. de Jong, LL.M., MM., MBA, CFD, CFFA, AA www.olivierdejong.com Solving Flir Furniture
More informationThermal Diffusivity. Paul Hughes. Department of Physics and Astronomy The University of Manchester Manchester M13 9PL. Second Year Laboratory Report
Therml iffusivity Pul Hughes eprtment of Physics nd Astronomy The University of nchester nchester 3 9PL Second Yer Lbortory Report Nov 4 Abstrct We investigted the therml diffusivity of cylindricl block
More informationPhys 7221, Fall 2006: Homework # 6
Phys 7221, Fll 2006: Homework # 6 Gbriel González October 29, 2006 Problem 37 In the lbortory system, the scttering ngle of the incident prticle is ϑ, nd tht of the initilly sttionry trget prticle, which
More informationMASKING OF FERROMAGNETIC ELLIPTICAL SHELL IN TRANSVERSE MAGNETIC FIELD
POZNAN UNVE RSTY OF TE HNOLOGY AADE M JOURNALS No 7 Electricl Engineering Kzimierz JAKUUK* Mirosł WOŁOSZYN* Peł ZMNY* MASKNG OF FERROMAGNET ELLPTAL SHELL N TRANSVERSE MAGNET FELD A ferromgnetic oject,
More informationPartial Derivatives. Limits. For a single variable function f (x), the limit lim
Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the righthnd side limit equls to the lefthnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles
More informationVisual motion. Many slides adapted from S. Seitz, R. Szeliski, M. Pollefeys
Visul motion Mn slides dpted from S. Seitz, R. Szeliski, M. Pollefes Outline Applictions of segmenttion to video Motion nd perceptul orgniztion Motion field Opticl flow Motion segmenttion with lers Video
More informationA SignalLevel Fusion Model for ImageBased Change Detection in DARPA's Dynamic Database System
SPIE Aerosense 001 Conference on Signl Processing, Sensor Fusion, nd Trget Recognition X, April 160, Orlndo FL. (Minor errors in published version corrected.) A SignlLevel Fusion Model for ImgeBsed
More informationI. Equations of a Circle a. At the origin center= r= b. Standard from: center= r=
11.: Circle & Ellipse I cn Write the eqution of circle given specific informtion Grph circle in coordinte plne. Grph n ellipse nd determine ll criticl informtion. Write the eqution of n ellipse from rel
More informationOn the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou
On the Uncertinty of Sensors Bsed on Mgnetic Effects E. ristoforou, E. Kyfs, A. Kten, DM Kepptsoglou Ntionl Technicl University of Athens, Zogrfou Cmpus, Athens 1578, Greece Tel: +3177178, Fx: +3177119,
More informationMTH 505: Number Theory Spring 2017
MTH 505: Numer Theory Spring 207 Homework 2 Drew Armstrong The Froenius Coin Prolem. Consider the eqution x ` y c where,, c, x, y re nturl numers. We cn think of $ nd $ s two denomintions of coins nd $c
More informationSection 6.1 Definite Integral
Section 6.1 Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot e determined
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
More informationOn the application of explicit spatial filtering to the variables or fluxes of linear equations
Journl of Computtionl Physics 225 (27) 2 27 www.elsevier.com/locte/jcp Short Note On the ppliction of explicit sptil filtering to the vriles or fluxes of liner equtions Christophe Bogey *, Christophe Billy
More informationHow to simulate Turing machines by invertible onedimensional cellular automata
How to simulte Turing mchines by invertible onedimensionl cellulr utomt JenChristophe Dubcq Déprtement de Mthémtiques et d Informtique, École Normle Supérieure de Lyon, 46, llée d Itlie, 69364 Lyon Cedex
More information1 Online Learning and Regret Minimization
2.997 DecisionMking in LrgeScle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in
More informationAN020. a a a. cos. cos. cos. Orientations and Rotations. Introduction. Orientations
AN020 Orienttions nd Rottions Introduction The fct tht ccelerometers re sensitive to the grvittionl force on the device llows them to be used to determine the ttitude of the sensor with respect to the
More informationFactors affecting the phonation threshold pressure and frequency
3SC Fctors ffecting the phontion threshold pressure nd frequency Zhoyn Zhng School of Medicine, University of Cliforni Los Angeles, CA, USA My, 9 57 th ASA Meeting, Portlnd, Oregon Acknowledgment: Reserch
More informationMinimum Energy State of Plasmas with an Internal Transport Barrier
Minimum Energy Stte of Plsms with n Internl Trnsport Brrier T. Tmno ), I. Ktnum ), Y. Skmoto ) ) Formerly, Plsm Reserch Center, University of Tsukub, Tsukub, Ibrki, Jpn ) Plsm Reserch Center, University
More informationCharacterization of Impact Test Response of PCCP with System Identification Approaches
7th World Conference on Nondestructive Testing, 2528 Oct 2008, Shnghi, Chin Chrcteriztion of Impct Test Response of PCCP with System Identifiction Approches Astrct Zheng LIU, Alex WANG, nd Dennis KRYS
More informationCOUPLING OF DAMAGE MECHANICS AND PROBABILISTIC APPROACH FOR LIFETIME PREDICTION OF COMPOSITE STRUCTURES
ORAL/POSTER REFERENCE : COUPLING OF DAMAGE MECHANICS AND PROBABILISTIC APPROACH FOR LIFETIME PREDICTION OF COMPOSITE STRUCTURES Y. Bruner, J. Renrd, D. Jeulin nd A. Thionnet Centre des Mtériux PierreMrie
More informationProblem Set 9. Figure 1: Diagram. This picture is a rough sketch of the 4 parabolas that give us the area that we need to find. The equations are:
(x + y ) = y + (x + y ) = x + Problem Set 9 Discussion: Nov., Nov. 8, Nov. (on probbility nd binomil coefficients) The nme fter the problem is the designted writer of the solution of tht problem. (No one
More informationSome basic concepts of fluid dynamics derived from ECE theory
Some sic concepts of fluid dynmics 363 Journl of Foundtions of Physics nd Chemistry, 2, vol. (4) 363 374 Some sic concepts of fluid dynmics derived from ECE theory M.W. Evns Alph Institute for Advnced
More informationShape and measurement
C H A P T E R 5 Shpe nd mesurement Wht is Pythgors theorem? How do we use Pythgors theorem? How do we find the perimeter of shpe? How do we find the re of shpe? How do we find the volume of shpe? How do
More information10. AREAS BETWEEN CURVES
. AREAS BETWEEN CURVES.. Ares etween curves So res ove the xxis re positive nd res elow re negtive, right? Wrong! We lied! Well, when you first lern out integrtion it s convenient fiction tht s true in
More information[ ( ) ( )] Section 6.1 Area of Regions between two Curves. Goals: 1. To find the area between two curves
Gols: 1. To find the re etween two curves Section 6.1 Are of Regions etween two Curves I. Are of Region Between Two Curves A. Grphicl Represention = _ B. Integrl Represention [ ( ) ( )] f x g x dx = C.
More informationLINEAR ALGEBRA APPLIED
5.5 Applictions of Inner Product Spces 5.5 Applictions of Inner Product Spces 7 Find the cross product of two vectors in R. Find the liner or qudrtic lest squres pproimtion of function. Find the nthorder
More information1.2 What is a vector? (Section 2.2) Two properties (attributes) of a vector are and.
Homework 1. Chpters 2. Bsis independent vectors nd their properties Show work except for fillinlnksprolems (print.pdf from www.motiongenesis.com Textooks Resources). 1.1 Solving prolems wht engineers
More informationMath 113 Exam 1Review
Mth 113 Exm 1Review September 26, 2016 Exm 1 covers 6.17.3 in the textbook. It is dvisble to lso review the mteril from 5.3 nd 5.5 s this will be helpful in solving some of the problems. 6.1 Are Between
More informationCalculus Module C21. Areas by Integration. Copyright This publication The Northern Alberta Institute of Technology All Rights Reserved.
Clculus Module C Ares Integrtion Copright This puliction The Northern Alert Institute of Technolog 7. All Rights Reserved. LAST REVISED Mrch, 9 Introduction to Ares Integrtion Sttement of Prerequisite
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk out solving systems of liner equtions. These re prolems tht give couple of equtions with couple of unknowns, like: 6= x + x 7=
More informationHomework Assignment 3 Solution Set
Homework Assignment 3 Solution Set PHYCS 44 6 Ferury, 4 Prolem 1 (Griffiths.5(c The potentil due to ny continuous chrge distriution is the sum of the contriutions from ech infinitesiml chrge in the distriution.
More informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
More informationTechnical Appendix: Childhood Family Structure and Schooling Outcomes: Evidence for Germany
Technicl Appendix: Childhood Fmily Structure nd Schooling Outcomes: Evidence for Germny Mrco Frncesconi* Stephen P. Jenkins Thoms Siedler Universy of Essex nd Universy of Essex Universy of Essex Instute
More informationSimulation of Eclipsing Binary Star Systems. Abstract
Simultion of Eclipsing Binry Str Systems Boris Yim 1, Kenny Chn 1, Rphel Hui 1 Wh Yn College Kowloon Diocesn Boys School Abstrct This report briefly introduces the informtion on eclipsing binry str systems.
More informationWhat Is Calculus? 42 CHAPTER 1 Limits and Their Properties
60_00.qd //0 : PM Pge CHAPTER Limits nd Their Properties The Mistress Fellows, Girton College, Cmridge Section. STUDY TIP As ou progress through this course, rememer tht lerning clculus is just one of
More informationMORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)
MORE FUNCTION GRAPHING; OPTIMIZATION FRI, OCT 25, 203 (Lst edited October 28, 203 t :09pm.) Exercise. Let n be n rbitrry positive integer. Give n exmple of function with exctly n verticl symptotes. Give
More informationThe Basic Functional 2 1
2 The Bsic Functionl 2 1 Chpter 2: THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 2.1 Introduction..................... 2 3 2.2 The First Vrition.................. 2 3 2.3 The Euler Eqution..................
More informationHomework Solution  Set 5 Due: Friday 10/03/08
CE 96 Introduction to the Theory of Computtion ll 2008 Homework olution  et 5 Due: ridy 10/0/08 1. Textook, Pge 86, Exercise 1.21. () 1 2 Add new strt stte nd finl stte. Mke originl finl stte nonfinl.
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationData Assimilation. Alan O Neill Data Assimilation Research Centre University of Reading
Dt Assimiltion Aln O Neill Dt Assimiltion Reserch Centre University of Reding Contents Motivtion Univrite sclr dt ssimiltion Multivrite vector dt ssimiltion Optiml Interpoltion BLUE 3dVritionl Method
More informationIntro to Nuclear and Particle Physics (5110)
Intro to Nucler nd Prticle Physics (5110) Feb, 009 The Nucler Mss Spectrum The Liquid Drop Model //009 1 E(MeV) n n(n1)/ E/[ n(n1)/] (MeV/pir) 1 C 16 O 0 Ne 4 Mg 7.7 14.44 19.17 8.48 4 5 6 6 10 15.4.41
More informationFarey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University
U.U.D.M. Project Report 07:4 Frey Frctions Rickrd Fernström Exmensrete i mtemtik, 5 hp Hledre: Andres Strömergsson Exmintor: Jörgen Östensson Juni 07 Deprtment of Mthemtics Uppsl University Frey Frctions
More information1 Error Analysis of Simple Rules for Numerical Integration
cs41: introduction to numericl nlysis 11/16/10 Lecture 19: Numericl Integrtion II Instructor: Professor Amos Ron Scries: Mrk Cowlishw, Nthnel Fillmore 1 Error Anlysis of Simple Rules for Numericl Integrtion
More information12.1 Nondeterminism Nondeterministic Finite Automata. a a b ε. CS125 Lecture 12 Fall 2016
CS125 Lecture 12 Fll 2016 12.1 Nondeterminism The ide of nondeterministic computtions is to llow our lgorithms to mke guesses, nd only require tht they ccept when the guesses re correct. For exmple, simple
More informationMath 100 Review Sheet
Mth 100 Review Sheet Joseph H. Silvermn December 2010 This outline of Mth 100 is summry of the mteril covered in the course. It is designed to be study id, but it is only n outline nd should be used s
More information9.4. The Vector Product. Introduction. Prerequisites. Learning Outcomes
The Vector Product 9.4 Introduction In this section we descrie how to find the vector product of two vectors. Like the sclr product its definition my seem strnge when first met ut the definition is chosen
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be relvlues nd smooth The pproximtion of n integrl by numericl
More informationUnit #10 De+inite Integration & The Fundamental Theorem Of Calculus
Unit # De+inite Integrtion & The Fundmentl Theorem Of Clculus. Find the re of the shded region ove nd explin the mening of your nswer. (squres re y units) ) The grph to the right is f(x) = x + 8x )Use
More informationPART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point.
PART MULTIPLE CHOICE Circle the pproprite response to ech of the questions below. Ech question hs vlue of point.. If in sequence the second level difference is constnt, thn the sequence is:. rithmetic
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
More informationu(t)dt + i a f(t)dt f(t) dt b f(t) dt (2) With this preliminary step in place, we are ready to define integration on a general curve in C.
Lecture 4 Complex Integrtion MATHGA 2451.001 Complex Vriles 1 Construction 1.1 Integrting complex function over curve in C A nturl wy to construct the integrl of complex function over curve in the complex
More informationMore precisely, given the collection fx g, with Eucliden distnces between pirs (; b) of ptterns: = p (x? x b ) ; one hs to nd mp, ' : R n distnceerro
Improved Multidimensionl Scling Anlysis Using Neurl Networks with DistnceError Bckpropgtion Llus Grrido (), Sergio Gomez () nd Jume Roc () () Deprtment d'estructur i Constituents de l Mteri/IFAE Universitt
More informationA. Limits  L Hopital s Rule ( ) How to find it: Try and find limits by traditional methods (plugging in). If you get 0 0 or!!, apply C.! 1 6 C.
A. Limits  L Hopitl s Rule Wht you re finding: L Hopitl s Rule is used to find limits of the form f ( x) lim where lim f x x! c g x ( ) = or lim f ( x) = limg( x) = ". ( ) x! c limg( x) = 0 x! c x! c
More informationLine and Surface Integrals: An Intuitive Understanding
Line nd Surfce Integrls: An Intuitive Understnding Joseph Breen Introduction Multivrible clculus is ll bout bstrcting the ides of differentition nd integrtion from the fmilir single vrible cse to tht of
More informationSection 4.8. D v(t j 1 ) t. (4.8.1) j=1
Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions
More informationA new algorithm for generating Pythagorean triples 1
A new lgorithm for generting Pythgoren triples 1 RH Dye 2 nd RWD Nicklls 3 The Mthemticl Gzette (1998; 82 (Mrch, No. 493, pp. 86 91 http://www.nicklls.org/dick/ppers/mths/pythgtriples1998.pdf 1 Introduction
More informationSTRAND J: TRANSFORMATIONS, VECTORS and MATRICES
Mthemtics SKE: STRN J STRN J: TRNSFORMTIONS, VETORS nd MTRIES J3 Vectors Text ontents Section J3.1 Vectors nd Sclrs * J3. Vectors nd Geometry Mthemtics SKE: STRN J J3 Vectors J3.1 Vectors nd Sclrs Vectors
More informationA Convergence Theorem for the Improper Riemann Integral of Banach Spacevalued Functions
Interntionl Journl of Mthemticl Anlysis Vol. 8, 2014, no. 50, 24512460 HIKARI Ltd, www.mhikri.com http://dx.doi.org/10.12988/ijm.2014.49294 A Convergence Theorem for the Improper Riemnn Integrl of Bnch
More informationpivot F 2 F 3 F 1 AP Physics 1 Practice Exam #3 (2/11/16)
AP Physics 1 Prctice Exm #3 (/11/16) Directions: Ech questions or incomplete sttements below is followed by four suggested nswers or completions. Select one tht is best in ech cse nd n enter pproprite
More informationset is not closed under matrix [ multiplication, ] and does not form a group.
Prolem 2.3: Which of the following collections of 2 2 mtrices with rel entries form groups under [ mtrix ] multipliction? i) Those of the form for which c d 2 Answer: The set of such mtrices is not closed
More information