Sheared Interpolation and Gradient Estimation. for Real-Time Volume Rendering. Hanspeter Pster, Frank Wessels, and Arie Kaufman

Size: px
Start display at page:

Download "Sheared Interpolation and Gradient Estimation. for Real-Time Volume Rendering. Hanspeter Pster, Frank Wessels, and Arie Kaufman"

Transcription

1 Sheared Iterpolatio ad Gradiet Estimatio for Real-Time Volume Rederig Haspeter Pster, Frak Wessels, ad Arie Kaufma Departmet of Computer Sciece State Uiversity of New York at Stoy Brook Stoy Brook, NY 11794{4400, U.S.A. Abstract I this paper we preset a techique for the iteractive cotrol ad display of static ad dyamic 3D datasets. We describe ovel ways of tri-liear iterpolatio ad gradiet estimatio for a real-time volume rederig system, usig coherecy betwee rays. We show simulatio results that compare the proposed methods to traditioal algorithms ad preset them i the cotext of Cube-3, a special-purpose architecture capable of rederig bit per voxel datasets at over 20 frames per secod. 1. Itroductio Numerous scietic applicatios, icludig biomedical ad geophysical aalysis, computatioal uid dyamics ad ite elemet models, require the rapid display of dyamically acquired or computer geerated 3D datasets. Real-time visualizatio of dyamic volume data, called 4D (spatial-temporal) visualizatio, permits observatio of 3D data chages, such as the study of uid ow i rocks or the study of a beatig heart. I order to reveal the iteral structure of the data, direct volume rederig methods have to be employed that geerate a image without pre-processig ad allow for the iteractive cotrol of viewig parameters [6]. The massive computatioal resources ecessary to achieve 4D visualizatio at high frame rates place hard to meet requiremets o sequetial implemetatios ad geeral-purpose computers. Oly parallelism amog a dedicated set of processors ca achieve the ecessary high memory badwidth ad arithmetic performace [4, 7, 12, 15] [6, Chapter 6]. While relatively fast algorithms exist for the display of static datasets o massively parallel architectures [17, 18], very little attetio has bee paid to the real-time visualizatio of dyamically chagig high-resolutio 3D data. This is the mai objective of Cube-3, a special-purpose architecture capable of rederig bit per voxel datasets at over 20 frames per secod [16]. Authors' pster@cs.suysb.edu, frakw@cs.suysb.edu, ari@cs.suysb.edu Cube-3 implemets ray-castig, a powerful volume rederig techique that oers high image quality while allowig for algorithmic optimizatios which sigicatly reduce image geeratio times [6, 13, 14]. Rays are cast from the viewig positio ito the volume data. At evely spaced locatios alog each ray, the data is tri-liearly iterpolated usig values of surroudig voxels. Cetral diereces of voxels aroud the sample poit yield a gradiet which is used as a surface ormal approximatio. Usig the gradiet ad the iterpolated sample value, a local shadig model is applied ad a sample opacity is assiged. Fially, ray samples alog the ray are composited ito pixel values to produce a image [11]. A importat problem of ray-castig is the ouiform mappig of samples oto voxels, sice voxels may cotai more tha oe ray sample or may be ivolved i multiple gradiet calculatios. This leads to redudat data accesses ad irregular iterprocessor commuicatio that aect the performace. I Cube-3 we use a ray-castig approach that trasforms the volume ito a itermediate coordiate system for which there is a mappig of ray samples oto the volume that is oe-to-oe. This allows for eciet projectios oto a face of the volume, ad the distorted image is the warped (2D trasformed ad projected) oto the view plae. Usig a similar approach, Yagel ad Kaufma [20] describe a template based ray-castig scheme to simplify path geeratio for rays through the volume, ad Schroder ad Stoll [17] have implemeted this method o a Priceto Egie of 1024 processors ad have achieved sub-secod rederig times for a dataset. Camero ad Uderill [3] ecietly use a itermediate volume trasformatio to reduce data commuicatio i a SIMD parallel processor. Lacroute ad Levoy [10] recetly reported o a fast implemetatio usig a shear-warp trasformatio ad were able to achieve iteractive rederig times for datasets o a graphics workstatio. All these implemetatios require a pre-processig step to calculate the gradiet eld or to geerate color ad opacity volumes ad are therefore ot suitable for 4D visualizatio. Appeared i Eurographics Hardware Workshop, Oslo, September 94 1

2 This paper presets two ew methods that allow for real-time tri-liear iterpolatio ad gradiet estimatio without pre-computatio. They are suitable for 4D visualizatio ad lead to a eciet implemetatio i hardware. Sectio 2 describes the uderlyig realtime ray-castig approach that trasforms the volume ito a itermediate sheared coordiate space. Sectio 3 discusses the problems associated with performig iterpolatio i this sheared space ad itroduces sheared tri-liear iterpolatio as a eective solutio. We the preset a ew way of gradiet approximatio usig coherecy betwee rays i Sectio 4. Sectio 5 describes the mai architectural features of Cube-3 ad Sectio 6 gives results o the proposed iterpolatio ad shadig methods. 2. Real-Time Ray-Castig Our real-time ray-castig algorithm assumes that the volume is sampled o a rectiliear grid. A distorted itermediate image is projected oto the volume face that is most perpedicular to the viewig directio. Usig a term by Yagel ad Kaufma [20] we call this face the base-plae. A 2D warp of the base-plae projectio produces the al image. The rst step is to trasform the volume ito a itermediate coordiate system for which there is a simple mappig of voxels oto base-plae pixels. I a recet approach, Lacroute ad Levoy [10] use a shear-warp factorizatio of the viewig trasform ad project the volume i a slice-parallel fashio oto the base-plae. The volume is treated as a set of 2D slices which are subject to a 2D shear-scale ad resamplig operatio accordig to the viewig trasform. Each slice is treated idepedetly without computig idividual rays, ad the resultig base-plae image is warped oto the viewig plae. Other approaches [20] operate i a ray-parallel fashio, where resamplig ad compositig operatios take place o rays cast from each pixel of the base-plae. I both approaches the 3D volume is traversed oly oce per projectio. The algorithms ivolve oe resamplig of the volume ad a iexpesive 2D image warp. I Cube-3 we adopted the ray-parallel approach because it allows for eciet parallel implemetatios of compositig alog rays. Usig a techique by Yagel ad Kaufma [20], we geerate lookup tables or templates to cast discrete rays from the base-plae ito the volume. Figure 1 shows a example of a parallel ad perspective projectio. 26-coected discrete lies are pre-geerated usig a 3D variatio [8] [6, pp. 280{301] of Breseham's algorithm modied for o-iteger edpoits. This algorithm guaratees costat steppig by a distace of oe alog the major axis (the Z-axis i Figure 1). The steppig alog the two other axes (the X- ad Y-axes i Figure 1) is stored i two templates. For parallel projectios, where eighborig rays follow the exact same path through the volume, the templates store positios for a 3 volume. For perspective projectios they are of size 2 each (see Figure 1). Y-Template Y Z X X-Template Y-Template a) Parallel Projectio b) Perspective Projectio Y Z X X-Template Figure 1: X/Y-Templates for Discrete Rays. Figure 2 schematically shows how the algorithm proceeds. All the discrete rays belogig to the same scalie of the base-plae image reside o the same plae iside the volume, called the Projectio Ray Plae (PRP). By fetchig all voxels o a PRP ad trasformig them accordigly ito a 2D buer, all discrete rays ca be aliged alog a directio parallel to a axis, e.g. horizotal. If we dee beams to be rays parallel to a mai axis of the Cubic Frame Buer (CFB), the for parallel projectios this trasformatio is simply a shear of beams to the left or right (see Figure 2). For perspective projectios each voxel belogig to a discrete ray has to be shifted by a dieret amout. We refer to this process as de-faig, sice divergig rays are stored adjacet to each other i the 2D buer. Orthogoal Beams Image Plae Fetch ad Shear Above PRP Curret Ray Cubic Frame Buffer (CFB) Compositig ad 2D Warp Tri-Liear Iterpolatio Above Curret Below Gradiet Estimatio ad Shadig Figure 2: Real-Time Ray-Castig. As soo as two PRPs are stored i two 2D buers (referred to as the above ad curret buers i Figure 2), a tri-liear iterpolatio is performed to geerate Appeared i Eurographics Hardware Workshop, Oslo, September 94 2

3 sample poits o cotiuous rays usig the voxels of four discrete rays as iput data (see Sectio 3). The two 2D buers geerate oe iterpolated plae of cotiuous rays. Three such plaes, above, below ad curret, are eeded for local gradiet approximatios usig eighborig rays (see Sectio 4). The samples of the rays are shaded ad opacities are assiged usig a user cotrollable trasfer fuctio. The shaded rays are composited ito a al pixel color usig a parallel implemetatio of the frot-to-back (or back-to-frot) compositig: C 0 = C L + (1? L )C R 0 = L + (1? L ) R (1) Here the subscripts L ad R idicate sample color C or opacity from left or right childre of the biary tree, respectively. Other parallel projectio schemes such as rst or last opaque projectio, maximum or miimum voxel value ad weighted summatio ca also be employed. The ext sectio discusses the issues of tri-liear iterpolatio betwee discrete rays to geerate cotiuous rays, ad Sectio 4 shows how to compute the local gradiet at each cotiuous sample poit. 3. Sheared Tri-Liear Iterpolatio Tri-liear iterpolatio geerates a value at oiteger locatios by fetchig the eight surroudig voxels ad iterpolatig as follows: P abc = P 000 (1? a)(1? b)(1? c) + P 100 a(1? b)(1? c) + P 010 (1? a)b(1? c) + P 001 (1? a)(1? b)c + P 101 a(1? b)c + P 011 (1? a)bc + P 110 ab(1? c) + P 111 abc: (2) Here the relative 3D coordiate of a sample poit withi a cube with respect to the corer voxel closest to the origi is ha; b; ci ad the data values associated with the corer voxels of the cube are P ijk, where i, j, k = 0 or 1, ad the iterpolated data value associated with the sample poit is P abc. Dieret optimizatios aim at reducig the arithmetic complexity of this operatio [9, 16], but the arbitrary memory access to fetch eight eighborig voxels for each sample poit makes this oe of the most time cosumig operatios durig volume rederig. By trasformig discrete rays from the PRP so that they are aliged ad storig them i two 2D buers (see Figure 2), we ca greatly reduce this data access ad commuicatio cost. Istead of fetchig the eighteighborhood of each resamplig locatio, four discrete rays are fetched from the buer, two from each of the above ad below plaes. I parallel implemetatios, eighborig rays reside i adjacet iterpolatio modules, requirig oly a local shift operatio of oe voxel uit betwee eighbors. S5 S4 S3 S2 S1 a) Parallel Projectio A B A B b) Perspective Projectio Discrete Ray A Discrete Ray B Missig Voxels Figure 3: Problems with Discrete Ray Iterpolatio. However, there is a problem itrisic to iterpolatio betwee discrete rays. Figure 3 illustrates this i 2D. The samples o the cotiuous ray have to be iterpolated usig bi-liear iterpolatio betwee samples of the discrete rays A (white) ad B (black). Sample S1 ca be correctly iterpolated usig four voxels from A ad B, sice they form a rectagle, i.e., the rays do ot make a discrete step to the left or right. As soo as the discrete rays step to the left or right as is the case for samples S2 ad S4, the eighborig voxels form a parallelogram, ad a straightforward biliear iterpolatio would produce the wrog sample values. The grey shaded square voxels i Figure 3a would be eeded to yield the correct result, but they reside o rays two uits apart from ray B. This problem is exacerbated for perspective projectios (Figure 3b). The discrete rays diverge, ad the correct eighborig voxels are ot eve stored i the 2D plae buers. For example, oly two voxels of ray A cotribute to the correct iterpolatio at sample poit S3. I the 3D case as may as six voxels may be missig i the immediate eighborhood of a sample poit for perspective projectios. The solutio is to perform a sheared tri-liear iterpolatio by factorig it ito four liear ad oe biliear iterpolatio. Istead of specifyig the sample locatio with respect to a corer voxel closest to the origi, each 3D coordiate alog the ray cosists of relative weights for liear iterpolatios alog each axis i possibly sheared voxel eighborhoods. These weights ca be pre-computed ad stored i the X/Y-templates discussed i Sectio 2. Figure 4 shows the ecessary iterpolatio steps i 3D. First we perform four liear iterpolatios i directio of the major axis (the Z-axis i Figure 4) usig Appeared i Eurographics Hardware Workshop, Oslo, September 94 3

4 Discrete Ray A Discrete Ray B Liearly Iterpolated Samples Ray Sample Y Z X a) Parallel Projectio Back Plaes Frot Plaes b) Perspective Projectio Figure 4: Sheared Tri-Liear Iterpolatio. eight voxels of four eighborig discrete rays iside the 2D buers. These eight voxels are the vertices of a oblique parallelepiped for parallel projectios (see Figure 4a) or of a frustum of a pyramid for perspective projectios (see Figure 4b). Four voxels each reside o two separate plaes oe uit apart, which we call the frot or the back plae depedig o whe it is ecoutered durig ray traversal i the directio of the major axis. Therefore, oly oe weight factor has to be stored, correspodig to the distace betwee the frot plae ad the positio of the ray sample poit. The resultig four iterpolated values form a rectagle ad ca be bi-liearly iterpolated to yield the al sample value. We split this bi-liear iterpolatio ito two liear iterpolatios betwee the corer values ad a al liear iterpolatio betwee the edge values. At the bottom of Figure 4 this is show as two iterpolatios i X- directio followed by oe iterpolatio i Y-directio. Discrete Ray A Discrete Ray B Out-of-Rage Samples a) No Offset b) Offset i Rage c) Offset out of Rage Figure 5: Variable Ray Osets i Major Directio. The sample poits correspodig to the cotiuous rays have to be iside the polyhedro deed by the voxels o the four surroudig discrete rays. Whe costructig the discrete rays, all cotiuous rays start at iteger positios of the base plae, i.e., they coicide with voxels of the rst slice of the volume dataset. However, as Figure 5a shows, usig these rays durig raycastig eectively reduces the tri-liear iterpolatio to a bi-liear iterpolatio, because all sample poits alog the ray fall oto the frot plaes of the parallelepipeds or pyramid frustum. Usig X ad Y iteger positios o the base-plae we ca allow a oset from the base-plae i major directio as a degree of freedom ad are able to perform sheared tri-liear iterpolatios (Figure 5b). But for osets i major directio that are too big, as show i Figure 5c), some of the samples alog the rays may fall outside the boudig box deed by the discrete rays. I order to get a upper boud for admissible osets we have to uderstad how steps i o-major directio alog discrete rays occur. Figure 6 shows the situatio i 2D. The view vector is split ito a dx compoet alog the X-axis (dx ad dy i 3D) ad a uit vector i directio of the major axis (the Y-axis i Figure 6). Steppig i directio of the major axis, we add the viewig vector to the curret sample positio at S i order to get the ew sample positio at S. dy=1 dx View-Vector o dx 1 - dx dx I-Rage Samples Figure 6: Maximum Oset Estimatio. Suppose that the additio of dx at poit S leads to a step of the discrete rays i x directio. This step ca oly occur if S has a relative x oset with respect to the lower left corer voxel of more tha 1? dx for positive dx (or less tha 1 + dx for egative dx). I other words, sample S was iside the rectagle of size dx by 1 show i Figure 6. However, oly the shaded regio of this rectagle cotais sample positios iside the parallelepiped deed by the corer voxels. Takig the smallest side i major axis as the worst-case, this meas that i-rage samples have a maximal relative y oset of o more tha 1? dx for positive dx (o less tha 1 + dx for egative dx). Sice we step with a uit vector i the directio of the major axis, all relative osets alog the ray are determied by the osets of the rst ray samples from the base-plae. The above argumet easily exteds to 3D, makig the maximum allowed oset i directio of the major axis: mi(1? dx; 1? dy); dx; dy 0 Appeared i Eurographics Hardware Workshop, Oslo, September 94 4

5 mi(1 + dx; 1? dy); dx < 0; dy 0 mi(1? dx; 1 + dy); dx 0; dy < 0 mi(1 + dx; 1 + dy); dx; dy < 0; (3) where dx ad dy are the compoets of the viewig vector i x ad y directio, respectively. Notice that for 45 o viewig agle dx ad dy are 1, yieldig a oset of 0 ad bi-liear iterpolatio as i Figure 5a. This fact will be of importace whe discussig the results i Sectio 6. I our implemetatio we cast a sigle ray from the origi of the image plae oto the base-plae usig uiform distace betwee samples ad choose the oset i major directio of the rst sample after peetratio of the base-plae. If ecessary the oset is iteratively reduced util it satises the above coditio. This leads to view depedet osets i major directio ad to varyig resamplig of the dataset. The variatio of resamplig poits accordig to viewig directio is a advatage for iteractive visualizatio, because more of the iteral data structure ca be revealed. Each discrete ray cosists of voxels, idepedet of the viewig directio. Sice the maximum viewig agle dierece with the major axis is ot more tha 45 degrees, the volume sample rate is deed by the diagoal through the cube ad is by a factor of p 3 higher for orthographic viewig. We foud that for raycompositig this is ot a importat cosideratio due to the averagig ature of the compositig operator. A more severe problem is the varyig size of the sample eighborhood (see Figure 4). For parallel projectios, the eight voxels surroudig the sample poit either form a cube with sides of legth oe or a oblique parallelepiped as i Figure 4a. For perspective projectios, however, the surroudig voxels may form the frustum of a pyramid with parallel frot ad back plaes as i Figure 4b. Due to the divergece of rays towards the back of the dataset, the volume spaed by this frustum icreases, thereby reducig the precisio of the tri-liear iterpolatio. However, we foud that the distace betwee eighborig discrete rays at the ed of the volume ever exceeded two voxels for a dataset while still achievig a high amout of perspectivity. Furthermore, i typical datasets the samples at the back of the volume have little iuece o the al pixel color due to compositig alog the ray. The ceter of projectio C ad the eld-of-view (FOV) i perspective projectios also iuece the samplig rate (see Figure 7). The discrete lie algorithm casts exactly oe ray per pixel of the base-plae, or a maximum of 2 rays per scalie. I cases where the FOV exteds across the the dataset (Figure 7a) this guaratees better samplig tha regular image order ray-castig, which would cast rays spaig the FOV ad sed wasteful rays that miss the dataset. However, for a small FOV the discrete lie steppig yields udersamplig i the active regios of the base-plae (Figure C a) Correct Samplig FOV FOV C b) Udersamplig C c) Two Base-Plae Projectios Figure 7: Samplig for Perspective Projectios. 7b). Figure 7c shows a case where two base-plae images cotribute to the al view image. The worst case i 3D is the geeratio of three base-plae projectios for a sigle perspective image. Sectio 6 presets comparisos betwee image order ray-castig usig a view idepedet samplig rate alog the rays, tri-liear iterpolatio employig equatio 2 usig the correct voxels, ad the proposed sheared tri-liear iterpolatio amog discrete rays. The ext sectio describes methods for gradiet estimatio usig samples o eighborig rays. 4. ABC Gradiet Estimatio To approximate the surface ormals ecessary for shadig ad classicatio we use the gray-level gradiet which is computed by the diereces betwee the values of the curret sample ad its immediate eighbors [5]. I order to evaluate the gradiet at a particular poit, we form cetral diereces betwee the tri-liearly iterpolated values of rays o the immediate left, right, above ad below, as well as the values of the curret ray. Sice this amouts to storig three cosecutive plaes of ray samples, we call this method ABC gradiet estimatio for the above, below, ad curret ray sample buers. The simplest approach, show i Figure 8 for 2D, is to use the 6-eighborhood gradiet, which uses the diereces of eighborig sample values alog the ray, P (;m+1)? P (;m?1) i base-plae directio ad P (;m?1)? P (?1;m+1) i the ray directio. Although the left, right, above ad below ray samples are i the same plae ad orthogoal to each other, the samples i the ray directio may be slated. A more critical problem occurs durig a switch of base-plae. Figure 8a shows the situatio for almost 45 o viewig directio, where a image is projected oto the horizotal baseplae. For ay agle greater tha 45 o a switch of baseplaes occurs, ad the values of P (;m)?p (?1;m) are used istead to calculate the gradiet i the base-plae directio. This leads to itolerable temporal aliasig. We also simulated the use of a 26-eighborhood gradiet (Figure 9). Istead of fetchig sample values from Appeared i Eurographics Hardware Workshop, Oslo, September 94 5

6 m-2 m-1 m m+1 m m-2 m-1 m m+1 m+2 (a) 6-eighborhood,, (b) 6-eighborhood, Horizotal Base-Plae Vertical Base-Plae Figure 8: 6-eighborhood Gradiet. four eighborig rays, 26 iterpolated samples from 8 eighborig rays are fetched. Each sample is assiged a weight factor correspodig to the iverse Mahatta distace i the iterpolated buer to the ceter sample. For example, sample P (;m?1) i Figure 9a has a weight of 1, whereas sample P (;m?2) has a weight of 1. I 2 3D we also get weight factors of 1 for the corer samples of the 26-eighborhood. However, to simplify the 3 arithmetic we use powers of 2, so that these samples are multiplied by a weight of 1. The gradiet is estimated 4 by takig weighted sums of ray samples ad diereces betwee opposite sample plaes. For the 2D example i Figure 9a this correspods to: G base = [ 1 2 P (;m) + P (;m+1) P (?1;m+2) ]? G ray = [ 1 [ 1 2 P (;m?2) + P (;m?1) P (?1;m) ] 2 P (;m?2) + P (;m?1) P (;m) ]? [ 1 2 P (?1;m) + P (?1;m+1) P (?1;m+2) ] (4) This method leads to better overall image quality whe compared to the 6-eighborhood gradiet, but the switchig of major axis is still oticeable (compare Figure 9a ad 9b). m-2 m-1 m m+1 (a) 26-eighborhood, Horizotal Base-Plae m+2 m-2 m-1 m m+1 (b) 26-eighborhood, Vertical Base-Plae Figure 9: 26-eighborhood Gradiet m+2 To circumvet this problem we take a similar approach to the 6-eighborhood method but use a additioal liear iterpolatio step to resample the rays o correct orthogoal positios. Figure 10 shows how the roud samples o the left ad right ray are used to liearly iterpolate the correct square samples. We call this approach the 10-eighborhood gradiet estimatio for the 3D case, sice 10 voxels participate i the computatio. It adequately solves the problem of switchig the major axis durig object rotatios ad yields high image quality. The liear iterpolatio weights are costat alog a ray ad correspod to a shift of all samples i the viewig directio. Sectio 6 presets a direct compariso betwee the 6-, 10- ad 26-eighborhood gradiet methods. -2 m-2 m-1 m m+1 m+2 (a) 10-eighborhood, Horizotal Base-Plae m-2 m-1 m m+1 m+2 (b) 10-eighborhood, Vertical Base-Plae Figure 10: 10-eighborhood Gradiet. I the case of perspective projectios, the frot of each PRP is uiformly sampled with rays oe uit apart. As the rays diverge towards the back of the volume, the distace betwee rays icreases, ad the gradiet estimatio becomes less accurate. However, because of the usually small distace betwee rays ad due to the averagig ature of shadig, classicatio ad compositig, these eects do ot iuece image quality for typical datasets. With the gradiet estimatio ad light vector directios, the sample itesity ca be geerated usig a variety of shadig methods (e.g., usig lookup tables [10]). Opacity values for compositig are geerated usig a trasfer fuctio represeted as a 2D lookup table idexed by sample desity ad gradiet magitude [11]. The ext sectio shows how the preseted sheared tri-liear iterpolatio ad ABC gradiet estimatio are supported i the Cube-3 architecture i order to achieve real-time 4D visualizatio. 5. Cube-3 Architecture Cube-3 is a special-purpose real-time volume visualizatio system that allows for the display of highresolutio bit per voxel datasets at frames rates over 20 Hz. It cotais a large CFB memory to hold the volumetric dataset ad performs base-plae projectios accordig to user cotrolled parameters. A host computer, coected to Cube-3 ad cotaiig the frame buer for the al image display, rus the user iter- Appeared i Eurographics Hardware Workshop, Oslo, September 94 6

7 face software ad performs the al 2D image warp oto the viewig plae. Real-time acquisitio devices such as a cofocal microscope, microtomograph, ultrasoud, or a computer ruig a simulatio model are tightly coupled to the Cube-3 memory usig high-badwidth optical liks for the iput of dyamically chagig 3D datasets. The Cube-3 architecture is highly-parallel ad pipelied [16]. Figure 11 shows a block diagram of the overall dataow. The CFB is a 3D memory orgaized i dual-access memory modules, each storig 2 voxels. A special 3D skewed orgaizatio eables the coictfree access to ay beam of voxels [7]. PRPs are fetched as a sequece of voxel beams ad stored i cosecutive 2D Skewed Buers (2DSB). A high-badwidth itercoectio etwork, the Fast Bus, allows the aligmet of the discrete rays o the PRP parallel to a mai axis i the 2DSB modules. PRP Cubic Frame Buffer (CFB) Frame Buffer Parallel Beam Fetch Fast Bus Projectio 2D Warpig 2D Skewed Buffer (2DSB) Ray Projectio Coe (RPC) Discrete Ray Fetch TRILIN Tri-Liear Iterpolatio Shadig Figure 11: Cube-3 System Overview. ABC Shadig Uits Three 2DSBs are used i a pipelied fashio to support sheared tri-liear iterpolatio. Aliged discrete rays from 2DSBs are fetched coict-free ad placed ito special purpose Tri-Liear Iterpolatio (TRILIN) uits. The resultig cotiuous projectio rays are placed oto ABC Shadig Uits, where the gradiets are estimated ad each ray sample is coverted ito both a itesity ad a associated opacity value accordig to lightig ad data segmetatio parameters. These itesity/opacity ray samples are fed ito the leaves of a Ray Projectio Coe (RPC). The RPC is a folded biary tree that geerates i parallel ad i a pipelied fashio the al pixel value usig a variety of projectio schemes o the coe odes. The resultig base-plae pixel is trasmitted to the host where it is post-processed (e.g., post-shaded or splatted) ad 2D trasformed (warped) oto the viewig plae. The result is stored i the 2D frame-buer. The parallel coict-free memory architecture of Cube-3 reduces the memory access bottleeck from O( 3 ) per projectio to O( 2 ) ad allows for very high data throughput. For a dataset size of bit voxels we estimate a performace of up to 30 frames per secod. Such a system would require 8 boards ad a custom fabricated backplae. Cube-3 is a scalable ad exible architecture that allows the user to iteractively cotrol the followig parameters: viewig agle from ay parallel ad perspective directio, cotrol over shadig ad projectio (e.g., rst opaque, maximum value, x-ray, compositig), color segmetatio ad thresholdig, cotrol over traslucecy, sectioig ad slicig. It will provide a rederig performace that is a order of magitude higher tha that of previously reported systems ad thereby revolutioize the way scietists coduct their studies. 6. Results We implemeted the dieret iterpolatio ad gradiet estimatio methods i software ad coducted several experimets. The rst program, VolRe implemets traditioal image order volume rederig. Rays are cast from the image plae ito the volume ad sampled at uiform steps. The tri-liear iterpolatio is performed accordig to Equatio 2 usig the correct 8- eighborhood aroud sample poits. The gradiet is estimated usig cetral diereces of tri-liear iterpolated values i a 6-eighborhood aroud each sample. The secod program, True3D, uses our real-time discrete ray-castig method, but istead of performig sheared tri-liear iterpolatio it fetches the exact 8- eighborhood aroud each sample poit. The last program, Sheared3D, implemets the same algorithm but with the proposed sheared tri-liear iterpolatio. Both True3D ad Sheared3D ca use ay of the 6-, 26- or 10-eighborhood gradiet methods for compariso purposes. For the implemetatio of these algorithms we used the VolVis volume visualizatio system, developed at the State Uiversity of New York at Stoy Brook [2, 1]. (The source code of VolVis is freely available by sedig to volvis@cs.suysb.edu.) 6.1. Tri-Liear Iterpolatio Compariso First we compare images resultig from Sheared3D to results obtaied from VolRe ad True3D. The gradiet approximatio method used for Sheared3D ad True3D was the proposed 10-eighborhood gradiet estimatio. The dataset, a CT study of a cadaver head of size voxels at 8-bit per voxel, was take o a Geeral Electric CT Scaer ad provided courtesy of North Carolia Memorial Hospital. All programs use the same shadig model ad a opacity trasfer fuctio that maps voxel values below 80 to = 0, has a liear ramp for from 0 to 0:75 for values betwee 80 ad 100, ad assigs = 0:75 to values above 100. We chose this particular trasfer fuctio to classify boe Appeared i Eurographics Hardware Workshop, Oslo, September 94 7

8 Figure 12: Dataset redered with sheared tri-liear iterpolatio (left) ad the dierece image to traditioal volume rederig (right) for 45 o rotatio agle. This is the worst case for sheared tri-liear iterpolatio. i the dataset as opaque i order to try to maximize the display of aliasig eects o the forehead of the CT skull. For the experimets we rotated the dataset by 70 o aroud the horizotal axis with respect to the world coordiate system, ad durig aimatios we rotated it aroud a vertical axis betwee 0 o ad 90 o i steps of 5 o. As error measure betwee the resultig images we use the average Euclidea distace of RGB values betwee correspodig pixels. Figure 12 shows the dataset rotated by 45 o aroud the vertical axis. The left image was geerated usig Sheared3D ad the image o the right is the dierece image, mapped to gray-scale, comparig the correspodig Sheared3D ad VolRe images for this rotatio agle. Figure 13 shows the relative Euclidea error i percetage betwee images from Sheared3D ad VolRe ad betwee Sheared3D ad True3D, respectively. The compariso with VolRe (top curve) shows how the error raises towards 45 o rotatio agle ad reaches a miimum at 0 o ad 90 o. The peak at 45 o is due to the differet samplig distace alog the ray, which is by p 3 bigger for discrete lie steppig (see Sectio 3). Furthermore, due to the oset cosideratios explaied i Sectio 3, our algorithm performs oly bi-liear iterpolatio as opposed to the the tri-liear iterpolatio i VolRe. The compariso to True3D shows zero error for 45 o because both algorithms perform bi-liear iterpolatio ad use the same gradiet estimatio techique. The relative error i percet compared to VolRe stays below 1:3%, ad compared to True3D it stays below 0:3%. Percet Error Sheared Tri-Liear Iterpolatio VolRe True3D Degrees Figure 13: Sheared Iterpolatio Percetage Error ABC Gradiet Estimatio Compariso For the compariso of the dieret ABC gradiet estimatio techiques we use a voxelized model of a sphere as dataset. The sphere is sca-coverted usig the volume samplig method described i [19]. The surface itersectio poits are obtaied by thresholdig, i.e., as soo as a certai voxel value is exceeded we calculate the gradiet at that poit. Each gradiet is compared to the true geometric surface ormal. As error measure we use the magitude of agular dierece betwee two vectors. All diereces are accumulated ad averaged over all surface itersectio poits. Appeared i Eurographics Hardware Workshop, Oslo, September 94 8

9 Figure 14: Error magitude of comparig surface ormals of 10- (Top) ad 26-eighborhood gradiets (Bottom) to the true aalytic ormal of the voxelized sphere. Notice the jump of regios of high error for the 26-eighborhood gradiet betwee 45 o ad 50 o rotatio agle. Dark: 0 o jej < 8:5 o, Medium: 8:5 o jej < 20 o, Light: 20 o jej < 31:5 o, White: jej 31:5 o. Rotatio agles (left to right): 30 o ; 35 o ; 40 o ; 45 o ; 50 o ; 55 o ; 60 o. Figure 15 shows the results of rotatig the sphere aroud a vertical axis betwee 0 o ad 90 o i steps of 5 o. The top two curves compare the aalytic ormal with the 26- ad the 6-eighborhood gradiet, respectively. The error icreases towards 45 o rotatio agle due to the o-orthogoality of the gradiet directios which reaches a maximum at 45 o. Although the 26-gradiet shows a little higher error magitude, the dierece betwee these two methods is ot sigicat. Average Error Magitude (Degrees) Gradiet Average Error 26-Grad 6-Grad 10-Grad Degrees Rotatio Figure 15: Average Error Magitude for ABC Gradiet Estimatios Compared to the Aalytic Normal. The curve o the bottom i Figure 15 shows the compariso of the aalytic ormal with the 10- eighborhood gradiet estimatio. The error magitude is sigicatly smaller tha for the other gradiet methods. The error also icreases towards 45 o rotatio agle. This is due to the dieret distaces betwee samples that are used for the gradiet calculatios i the three orthogoal directios. Figure 14 shows how the error propagates aroud the sphere for rotatio agles from 30 o to 60 o i steps of 5 o. Dark shaded regios idicate regios of low error magitude, light shaded regios idicate higher error magitudes. The top row shows the 10-eighborhood gradiet method with a fairly regular error trasitio from left to right durig a switch of base-plaes at 45 o (ceter sphere). The bottom row, depictig the 26- eighborhood gradiet method, shows a geerally larger error magitude. Additioally, the regio of largest error jumps from the right side of the sphere to the left durig the switch of base-plaes. This jump leads to oticeable chages i image itesity durig object rotatio, a eect that we described as temporal aliasig i Sectio Coclusios I order to achieve the goal of real-time visualizatio of dyamic datasets we developed Cube-3, a scalable architecture that exploits parallelism ad pipeliig. I this paper we preseted the uderlyig real-time ray-castig approach that allows for a mappig of raysamples oto voxels that is oe-to-oe. Usig templates ad shearig/de-faig of beams, we fetch 2D plaes from the volume dataset ad perform sheared tri-liear iterpolatio betwee discrete eighborig rays. Usig the resultig iterpolated ray samples from above, curret ad below plaes, we described ovel ways of gradiet estimatio usig coherecy betwee rays. Usig software simulatios we compared the proposed methods to traditioal image order ray-castig. The error of usig sheared tri-liear iterpolatio istead of performig image order ray-castig is below 1:3% relative dierece i Euclidea distace of the resultig image pixels. We showed that use of the proposed 10-eighborhood istead of a 6- or 26- eighborhood gradiet approach reduces both the average error compared to aalytically computed ormals ad the temporal aliasig that arises from switchig base-plaes durig object rotatios. We preseted both Appeared i Eurographics Hardware Workshop, Oslo, September 94 9

10 methods i the cotext of Cube-3, a special purpose architecture aimed at real-time 4D visualizatio of highresolutio volumetric datasets. 8. Ackowledgmets This work has bee supported by the Natioal Sciece Foudatio uder grat CCR We would like to thak Lisa Sobierajski ad Rick Avila for their helpful suggestios durig the developmet of these methods. Sidey Wag provided us with the sphere dataset ad helped with the geeratio of Figure 14. A discussio with Claudio Silva gave us the isight ito the various error metrics we used. We also thak Patrick Tora for helpful system admiistratio durig the more hectic momets i the developmet of this project. Refereces 1. Avila, R., He, T., Hog, L., Kaufma, A., Pfister, H., Silva, C., Sobierajski, L., ad Wag, S. VolVis: A diversied system for volume visualizatio research ad developmet. To appear i Proceedigs of Visualizatio '94 (Washigto, DC, Oct. 1994). 2. Avila, R., Sobierajski, L., ad Kaufma, A. Towards a comprehesive volume visualizatio system. I Proceedigs of Visualizatio '92 (Bosto, MA, Oct. 1992), IEEE Computer Society Press, pp. 13{ Camero, G., ad Uderill, P. E. Rederig volumetric medical image data o a SIMD architecture computer. I Proceedigs of Third Eurographics Workshop o Rederig (May 1992). 4. Fuchs, H., Poulto, J., Eyles, J., Greer, T., Goldfeather, J., Ellsworth, D., Molar, S., Turk, G., Tebbs, B., ad Israel, L. Pixel- Plaes 5: A heterogeeous multiprocessor graphics system usig processor-ehaced memories. Computer Graphics 23, No. 3 (July 1989), 79{ Hohe, K. H., ad Berstei, R. Shadig 3Dimages from CT usig gray-level gradiets. IEEE Trasactios o Medical Imagig MI-5, 1 (Mar. 1986), 45{ Kaufma, A. Volume Visualizatio. IEEE CS Press Tutorial, Los Alamitos, CA, Kaufma, A., ad Bakalash, R. Memory ad processig architecture for 3D voxel-based imagery. IEEE Computer Graphics & Applicatios 8, 6 (Nov. 1988), 10{ Kaufma, A., ad Shimoy, E. 3D scacoversio algorithms for voxel-based graphics. I ACM Workshop o Iteractive 3D Graphics (Chapel Hill, NC, Oct. 1986), pp. 45{ Kittel, G. VERVE: Voxel egie for realtime visualizatio ad examiatio. I Computer Graphics Forum (September 1993), vol. 12, No. 3, pp. C{37 { C{ Lacroute, P., ad Levoy, M. Fast volume rederig usig a shear-warp factorizatio of the viewig trasform. Computer Graphics, Proceedigs of SIGGRAPH '94 (July 1994), 451{ Levoy, M. Display of surfaces from volume data. IEEE Computer Graphics & Applicatios 8, 5 (May 1988), 29{ Levoy, M. Desig for real-time high-quality volume rederig workstatio. I 1989 Workshop o Volume Visualizatio (Chapel Hill, NC, May 1989), pp. 85{ Levoy, M. Eciet ray tracig of volume data. ACM Trasactios o Graphics 9, 3 (July 1990), 245{ Levoy, M. Volume rederig by adaptive reemet. The Visual Computer (July 1990), 2{ Molar, S., Eyles, J., ad Poulto, J. Pixelow: High-speed rederig usig image compositio. Computer Graphics 26, 2 (July 1992), 231{ Pfister, H., Kaufma, A., ad Chiueh, T. Cube-3: A Real-Time Architecture for High- Resolutio Volume Visualizatio. To appear i 1994 Workshop o Volume Visualizatio (Washigto, DC, Oct. 1994). 17. Schroder, P., ad Stoll, G. Data parallel volume rederig as lie drawig. I 1992 Workshop o Volume Visualizatio (Bosto, MA, Oct. 1992), pp. 25{ Vezia, G., Fletcher, P., ad Robertso, P. Volume rederig o the MasPar MP-1. I 1992 Workshop o Volume Visualizatio (Bosto, MA, Oct. 1992), pp. 3{ Wag, S., ad Kaufma, A. Volume sampled voxelizatio of geometric primitives. I Proceedigs of Visualizatio '93 (Sa Jose, CA, Oct. 1993), IEEE Computer Society Press, pp. 78{ Yagel, R., ad Kaufma, A. Template-based volume viewig. Computer Graphics Forum 11, 3 (Sept. 1992), 153{167. Appeared i Eurographics Hardware Workshop, Oslo, September 94 10

Sheared Interpolation and Gradient Estimation. State University of New York at Stony Brook. sample value, a local shading model is applied and a

Sheared Interpolation and Gradient Estimation. State University of New York at Stony Brook. sample value, a local shading model is applied and a Sheared Iterpolatio ad Gradiet Estimatio for Real-Time Volume Rederig Haspeter Pæster, Frak Wessels, ad Arie Kaufma Departmet of Computer Sciece State Uiversity of New York at Stoy Brook Stoy Brook, NY

More information

Bezier curves. Figure 2 shows cubic Bezier curves for various control points. In a Bezier curve, only

Bezier curves. Figure 2 shows cubic Bezier curves for various control points. In a Bezier curve, only Edited: Yeh-Liag Hsu (998--; recommeded: Yeh-Liag Hsu (--9; last updated: Yeh-Liag Hsu (9--7. Note: This is the course material for ME55 Geometric modelig ad computer graphics, Yua Ze Uiversity. art of

More information

Accuracy Improvement in Camera Calibration

Accuracy Improvement in Camera Calibration Accuracy Improvemet i Camera Calibratio FaJie L Qi Zag ad Reihard Klette CITR, Computer Sciece Departmet The Uiversity of Aucklad Tamaki Campus, Aucklad, New Zealad fli006, qza001@ec.aucklad.ac.z r.klette@aucklad.ac.z

More information

Improvement of the Orthogonal Code Convolution Capabilities Using FPGA Implementation

Improvement of the Orthogonal Code Convolution Capabilities Using FPGA Implementation Improvemet of the Orthogoal Code Covolutio Capabilities Usig FPGA Implemetatio Naima Kaabouch, Member, IEEE, Apara Dhirde, Member, IEEE, Saleh Faruque, Member, IEEE Departmet of Electrical Egieerig, Uiversity

More information

Elementary Educational Computer

Elementary Educational Computer Chapter 5 Elemetary Educatioal Computer. Geeral structure of the Elemetary Educatioal Computer (EEC) The EEC coforms to the 5 uits structure defied by vo Neuma's model (.) All uits are preseted i a simplified

More information

Ones Assignment Method for Solving Traveling Salesman Problem

Ones Assignment Method for Solving Traveling Salesman Problem Joural of mathematics ad computer sciece 0 (0), 58-65 Oes Assigmet Method for Solvig Travelig Salesma Problem Hadi Basirzadeh Departmet of Mathematics, Shahid Chamra Uiversity, Ahvaz, Ira Article history:

More information

A New Morphological 3D Shape Decomposition: Grayscale Interframe Interpolation Method

A New Morphological 3D Shape Decomposition: Grayscale Interframe Interpolation Method A ew Morphological 3D Shape Decompositio: Grayscale Iterframe Iterpolatio Method D.. Vizireau Politehica Uiversity Bucharest, Romaia ae@comm.pub.ro R. M. Udrea Politehica Uiversity Bucharest, Romaia mihea@comm.pub.ro

More information

Computer Graphics Hardware An Overview

Computer Graphics Hardware An Overview Computer Graphics Hardware A Overview Graphics System Moitor Iput devices CPU/Memory GPU Raster Graphics System Raster: A array of picture elemets Based o raster-sca TV techology The scree (ad a picture)

More information

Pattern Recognition Systems Lab 1 Least Mean Squares

Pattern Recognition Systems Lab 1 Least Mean Squares Patter Recogitio Systems Lab 1 Least Mea Squares 1. Objectives This laboratory work itroduces the OpeCV-based framework used throughout the course. I this assigmet a lie is fitted to a set of poits usig

More information

Normals. In OpenGL the normal vector is part of the state Set by glnormal*()

Normals. In OpenGL the normal vector is part of the state Set by glnormal*() Ray Tracig 1 Normals OpeG the ormal vector is part of the state Set by glnormal*() -glnormal3f(x, y, z); -glnormal3fv(p); Usually we wat to set the ormal to have uit legth so cosie calculatios are correct

More information

Alpha Individual Solutions MAΘ National Convention 2013

Alpha Individual Solutions MAΘ National Convention 2013 Alpha Idividual Solutios MAΘ Natioal Covetio 0 Aswers:. D. A. C 4. D 5. C 6. B 7. A 8. C 9. D 0. B. B. A. D 4. C 5. A 6. C 7. B 8. A 9. A 0. C. E. B. D 4. C 5. A 6. D 7. B 8. C 9. D 0. B TB. 570 TB. 5

More information

A Note on Least-norm Solution of Global WireWarping

A Note on Least-norm Solution of Global WireWarping A Note o Least-orm Solutio of Global WireWarpig Charlie C. L. Wag Departmet of Mechaical ad Automatio Egieerig The Chiese Uiversity of Hog Kog Shati, N.T., Hog Kog E-mail: cwag@mae.cuhk.edu.hk Abstract

More information

IMP: Superposer Integrated Morphometrics Package Superposition Tool

IMP: Superposer Integrated Morphometrics Package Superposition Tool IMP: Superposer Itegrated Morphometrics Package Superpositio Tool Programmig by: David Lieber ( 03) Caisius College 200 Mai St. Buffalo, NY 4208 Cocept by: H. David Sheets, Dept. of Physics, Caisius College

More information

Image Segmentation EEE 508

Image Segmentation EEE 508 Image Segmetatio Objective: to determie (etract) object boudaries. It is a process of partitioig a image ito distict regios by groupig together eighborig piels based o some predefied similarity criterio.

More information

Fundamentals of Media Processing. Shin'ichi Satoh Kazuya Kodama Hiroshi Mo Duy-Dinh Le

Fundamentals of Media Processing. Shin'ichi Satoh Kazuya Kodama Hiroshi Mo Duy-Dinh Le Fudametals of Media Processig Shi'ichi Satoh Kazuya Kodama Hiroshi Mo Duy-Dih Le Today's topics Noparametric Methods Parze Widow k-nearest Neighbor Estimatio Clusterig Techiques k-meas Agglomerative Hierarchical

More information

The Closest Line to a Data Set in the Plane. David Gurney Southeastern Louisiana University Hammond, Louisiana

The Closest Line to a Data Set in the Plane. David Gurney Southeastern Louisiana University Hammond, Louisiana The Closest Lie to a Data Set i the Plae David Gurey Southeaster Louisiaa Uiversity Hammod, Louisiaa ABSTRACT This paper looks at three differet measures of distace betwee a lie ad a data set i the plae:

More information

Creating Exact Bezier Representations of CST Shapes. David D. Marshall. California Polytechnic State University, San Luis Obispo, CA , USA

Creating Exact Bezier Representations of CST Shapes. David D. Marshall. California Polytechnic State University, San Luis Obispo, CA , USA Creatig Exact Bezier Represetatios of CST Shapes David D. Marshall Califoria Polytechic State Uiversity, Sa Luis Obispo, CA 93407-035, USA The paper presets a method of expressig CST shapes pioeered by

More information

Dynamic Programming and Curve Fitting Based Road Boundary Detection

Dynamic Programming and Curve Fitting Based Road Boundary Detection Dyamic Programmig ad Curve Fittig Based Road Boudary Detectio SHYAM PRASAD ADHIKARI, HYONGSUK KIM, Divisio of Electroics ad Iformatio Egieerig Chobuk Natioal Uiversity 664-4 Ga Deokji-Dog Jeoju-City Jeobuk

More information

CS Polygon Scan Conversion. Slide 1

CS Polygon Scan Conversion. Slide 1 CS 112 - Polygo Sca Coversio Slide 1 Polygo Classificatio Covex All iterior agles are less tha 180 degrees Cocave Iterior agles ca be greater tha 180 degrees Degeerate polygos If all vertices are colliear

More information

Cubic Polynomial Curves with a Shape Parameter

Cubic Polynomial Curves with a Shape Parameter roceedigs of the th WSEAS Iteratioal Coferece o Robotics Cotrol ad Maufacturig Techology Hagzhou Chia April -8 00 (pp5-70) Cubic olyomial Curves with a Shape arameter MO GUOLIANG ZHAO YANAN Iformatio ad

More information

Parallel Polygon Approximation Algorithm Targeted at Reconfigurable Multi-Ring Hardware

Parallel Polygon Approximation Algorithm Targeted at Reconfigurable Multi-Ring Hardware Parallel Polygo Approximatio Algorithm Targeted at Recofigurable Multi-Rig Hardware M. Arif Wai* ad Hamid R. Arabia** *Califoria State Uiversity Bakersfield, Califoria, USA **Uiversity of Georgia, Georgia,

More information

Fast Fourier Transform (FFT) Algorithms

Fast Fourier Transform (FFT) Algorithms Fast Fourier Trasform FFT Algorithms Relatio to the z-trasform elsewhere, ozero, z x z X x [ ] 2 ~ elsewhere,, ~ e j x X x x π j e z z X X π 2 ~ The DFS X represets evely spaced samples of the z- trasform

More information

Analysis Metrics. Intro to Algorithm Analysis. Slides. 12. Alg Analysis. 12. Alg Analysis

Analysis Metrics. Intro to Algorithm Analysis. Slides. 12. Alg Analysis. 12. Alg Analysis Itro to Algorithm Aalysis Aalysis Metrics Slides. Table of Cotets. Aalysis Metrics 3. Exact Aalysis Rules 4. Simple Summatio 5. Summatio Formulas 6. Order of Magitude 7. Big-O otatio 8. Big-O Theorems

More information

Polymorph: Morphing Among Multiple Images

Polymorph: Morphing Among Multiple Images Feature Article Polymorph: Morphig Amog Multiple Images Image metamorphosis has prove to be a powerful visual effects tool. May breathtakig examples ow appear i film ad televisio, depictig the fluid trasformatio

More information

CSC 220: Computer Organization Unit 11 Basic Computer Organization and Design

CSC 220: Computer Organization Unit 11 Basic Computer Organization and Design College of Computer ad Iformatio Scieces Departmet of Computer Sciece CSC 220: Computer Orgaizatio Uit 11 Basic Computer Orgaizatio ad Desig 1 For the rest of the semester, we ll focus o computer architecture:

More information

3D Model Retrieval Method Based on Sample Prediction

3D Model Retrieval Method Based on Sample Prediction 20 Iteratioal Coferece o Computer Commuicatio ad Maagemet Proc.of CSIT vol.5 (20) (20) IACSIT Press, Sigapore 3D Model Retrieval Method Based o Sample Predictio Qigche Zhag, Ya Tag* School of Computer

More information

Copyright 2016 Ramez Elmasri and Shamkant B. Navathe

Copyright 2016 Ramez Elmasri and Shamkant B. Navathe Copyright 2016 Ramez Elmasri ad Shamkat B. Navathe CHAPTER 19 Query Optimizatio Copyright 2016 Ramez Elmasri ad Shamkat B. Navathe Itroductio Query optimizatio Coducted by a query optimizer i a DBMS Goal:

More information

Chapter 3 Classification of FFT Processor Algorithms

Chapter 3 Classification of FFT Processor Algorithms Chapter Classificatio of FFT Processor Algorithms The computatioal complexity of the Discrete Fourier trasform (DFT) is very high. It requires () 2 complex multiplicatios ad () complex additios [5]. As

More information

Assignment 5; Due Friday, February 10

Assignment 5; Due Friday, February 10 Assigmet 5; Due Friday, February 10 17.9b The set X is just two circles joied at a poit, ad the set X is a grid i the plae, without the iteriors of the small squares. The picture below shows that the iteriors

More information

. Written in factored form it is easy to see that the roots are 2, 2, i,

. Written in factored form it is easy to see that the roots are 2, 2, i, CMPS A Itroductio to Programmig Programmig Assigmet 4 I this assigmet you will write a java program that determies the real roots of a polyomial that lie withi a specified rage. Recall that the roots (or

More information

Pseudocode ( 1.1) Analysis of Algorithms. Primitive Operations. Pseudocode Details. Running Time ( 1.1) Estimating performance

Pseudocode ( 1.1) Analysis of Algorithms. Primitive Operations. Pseudocode Details. Running Time ( 1.1) Estimating performance Aalysis of Algorithms Iput Algorithm Output A algorithm is a step-by-step procedure for solvig a problem i a fiite amout of time. Pseudocode ( 1.1) High-level descriptio of a algorithm More structured

More information

Copyright 2016 Ramez Elmasri and Shamkant B. Navathe

Copyright 2016 Ramez Elmasri and Shamkant B. Navathe Copyright 2016 Ramez Elmasri ad Shamkat B. Navathe CHAPTER 18 Strategies for Query Processig Copyright 2016 Ramez Elmasri ad Shamkat B. Navathe Itroductio DBMS techiques to process a query Scaer idetifies

More information

Lecture 5. Counting Sort / Radix Sort

Lecture 5. Counting Sort / Radix Sort Lecture 5. Coutig Sort / Radix Sort T. H. Corme, C. E. Leiserso ad R. L. Rivest Itroductio to Algorithms, 3rd Editio, MIT Press, 2009 Sugkyukwa Uiversity Hyuseug Choo choo@skku.edu Copyright 2000-2018

More information

APPLICATION NOTE PACE1750AE BUILT-IN FUNCTIONS

APPLICATION NOTE PACE1750AE BUILT-IN FUNCTIONS APPLICATION NOTE PACE175AE BUILT-IN UNCTIONS About This Note This applicatio brief is iteded to explai ad demostrate the use of the special fuctios that are built ito the PACE175AE processor. These powerful

More information

An array based design for Real-Time Volume Rendering

An array based design for Real-Time Volume Rendering A array based desig for Real-Time Volume Rederig Michael Doggett* School of Computer Sciece ad Egieerig The Uiversity of New South Wales 11 Abstract This paper describes a ew algorithm ad hardware desig

More information

Running Time. Analysis of Algorithms. Experimental Studies. Limitations of Experiments

Running Time. Analysis of Algorithms. Experimental Studies. Limitations of Experiments Ruig Time Aalysis of Algorithms Iput Algorithm Output A algorithm is a step-by-step procedure for solvig a problem i a fiite amout of time. Most algorithms trasform iput objects ito output objects. The

More information

Texture Mapping. Jian Huang. This set of slides references the ones used at Ohio State for instruction.

Texture Mapping. Jian Huang. This set of slides references the ones used at Ohio State for instruction. Texture Mappig Jia Huag This set of slides refereces the oes used at Ohio State for istructio. Ca you do this What Dreams May Come Texture Mappig Of course, oe ca model the exact micro-geometry + material

More information

EE123 Digital Signal Processing

EE123 Digital Signal Processing Last Time EE Digital Sigal Processig Lecture 7 Block Covolutio, Overlap ad Add, FFT Discrete Fourier Trasform Properties of the Liear covolutio through circular Today Liear covolutio with Overlap ad add

More information

Running Time ( 3.1) Analysis of Algorithms. Experimental Studies. Limitations of Experiments

Running Time ( 3.1) Analysis of Algorithms. Experimental Studies. Limitations of Experiments Ruig Time ( 3.1) Aalysis of Algorithms Iput Algorithm Output A algorithm is a step- by- step procedure for solvig a problem i a fiite amout of time. Most algorithms trasform iput objects ito output objects.

More information

Analysis of Algorithms

Analysis of Algorithms Aalysis of Algorithms Iput Algorithm Output A algorithm is a step-by-step procedure for solvig a problem i a fiite amout of time. Ruig Time Most algorithms trasform iput objects ito output objects. The

More information

Math 10C Long Range Plans

Math 10C Long Range Plans Math 10C Log Rage Plas Uits: Evaluatio: Homework, projects ad assigmets 10% Uit Tests. 70% Fial Examiatio.. 20% Ay Uit Test may be rewritte for a higher mark. If the retest mark is higher, that mark will

More information

condition w i B i S maximum u i

condition w i B i S maximum u i ecture 10 Dyamic Programmig 10.1 Kapsack Problem November 1, 2004 ecturer: Kamal Jai Notes: Tobias Holgers We are give a set of items U = {a 1, a 2,..., a }. Each item has a weight w i Z + ad a utility

More information

1. Introduction o Microscopic property responsible for MRI Show and discuss graphics that go from macro to H nucleus with N-S pole

1. Introduction o Microscopic property responsible for MRI Show and discuss graphics that go from macro to H nucleus with N-S pole Page 1 Very Quick Itroductio to MRI The poit of this itroductio is to give the studet a sufficietly accurate metal picture of MRI to help uderstad its impact o image registratio. The two major aspects

More information

Performance Plus Software Parameter Definitions

Performance Plus Software Parameter Definitions Performace Plus+ Software Parameter Defiitios/ Performace Plus Software Parameter Defiitios Chapma Techical Note-TG-5 paramete.doc ev-0-03 Performace Plus+ Software Parameter Defiitios/2 Backgroud ad Defiitios

More information

Data Structures and Algorithms. Analysis of Algorithms

Data Structures and Algorithms. Analysis of Algorithms Data Structures ad Algorithms Aalysis of Algorithms Outlie Ruig time Pseudo-code Big-oh otatio Big-theta otatio Big-omega otatio Asymptotic algorithm aalysis Aalysis of Algorithms Iput Algorithm Output

More information

1. SWITCHING FUNDAMENTALS

1. SWITCHING FUNDAMENTALS . SWITCING FUNDMENTLS Switchig is the provisio of a o-demad coectio betwee two ed poits. Two distict switchig techiques are employed i commuicatio etwors-- circuit switchig ad pacet switchig. Circuit switchig

More information

SD vs. SD + One of the most important uses of sample statistics is to estimate the corresponding population parameters.

SD vs. SD + One of the most important uses of sample statistics is to estimate the corresponding population parameters. SD vs. SD + Oe of the most importat uses of sample statistics is to estimate the correspodig populatio parameters. The mea of a represetative sample is a good estimate of the mea of the populatio that

More information

Harris Corner Detection Algorithm at Sub-pixel Level and Its Application Yuanfeng Han a, Peijiang Chen b * and Tian Meng c

Harris Corner Detection Algorithm at Sub-pixel Level and Its Application Yuanfeng Han a, Peijiang Chen b * and Tian Meng c Iteratioal Coferece o Computatioal Sciece ad Egieerig (ICCSE 015) Harris Corer Detectio Algorithm at Sub-pixel Level ad Its Applicatio Yuafeg Ha a, Peijiag Che b * ad Tia Meg c School of Automobile, Liyi

More information

COMP 558 lecture 6 Sept. 27, 2010

COMP 558 lecture 6 Sept. 27, 2010 Radiometry We have discussed how light travels i straight lies through space. We would like to be able to talk about how bright differet light rays are. Imagie a thi cylidrical tube ad cosider the amout

More information

Consider the following population data for the state of California. Year Population

Consider the following population data for the state of California. Year Population Assigmets for Bradie Fall 2016 for Chapter 5 Assigmet sheet for Sectios 5.1, 5.3, 5.5, 5.6, 5.7, 5.8 Read Pages 341-349 Exercises for Sectio 5.1 Lagrage Iterpolatio #1, #4, #7, #13, #14 For #1 use MATLAB

More information

RESEARCH ON AUTOMATIC INSPECTION TECHNIQUE OF REAL-TIME RADIOGRAPHY FOR TURBINE-BLADE

RESEARCH ON AUTOMATIC INSPECTION TECHNIQUE OF REAL-TIME RADIOGRAPHY FOR TURBINE-BLADE RESEARCH ON AUTOMATIC INSPECTION TECHNIQUE OF REAL-TIME RADIOGRAPHY FOR TURBINE-BLADE Z.G. Zhou, S. Zhao, ad Z.G. A School of Mechaical Egieerig ad Automatio, Beijig Uiversity of Aeroautics ad Astroautics,

More information

CIS 121 Data Structures and Algorithms with Java Spring Stacks, Queues, and Heaps Monday, February 18 / Tuesday, February 19

CIS 121 Data Structures and Algorithms with Java Spring Stacks, Queues, and Heaps Monday, February 18 / Tuesday, February 19 CIS Data Structures ad Algorithms with Java Sprig 09 Stacks, Queues, ad Heaps Moday, February 8 / Tuesday, February 9 Stacks ad Queues Recall the stack ad queue ADTs (abstract data types from lecture.

More information

Improving Template Based Spike Detection

Improving Template Based Spike Detection Improvig Template Based Spike Detectio Kirk Smith, Member - IEEE Portlad State Uiversity petra@ee.pdx.edu Abstract Template matchig algorithms like SSE, Covolutio ad Maximum Likelihood are well kow for

More information

Algorithms for Disk Covering Problems with the Most Points

Algorithms for Disk Covering Problems with the Most Points Algorithms for Disk Coverig Problems with the Most Poits Bi Xiao Departmet of Computig Hog Kog Polytechic Uiversity Hug Hom, Kowloo, Hog Kog csbxiao@comp.polyu.edu.hk Qigfeg Zhuge, Yi He, Zili Shao, Edwi

More information

6.854J / J Advanced Algorithms Fall 2008

6.854J / J Advanced Algorithms Fall 2008 MIT OpeCourseWare http://ocw.mit.edu 6.854J / 18.415J Advaced Algorithms Fall 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 18.415/6.854 Advaced Algorithms

More information

Appendix D. Controller Implementation

Appendix D. Controller Implementation COMPUTER ORGANIZATION AND DESIGN The Hardware/Software Iterface 5 th Editio Appedix D Cotroller Implemetatio Cotroller Implemetatios Combiatioal logic (sigle-cycle); Fiite state machie (multi-cycle, pipelied);

More information

Analysis of Algorithms

Analysis of Algorithms Aalysis of Algorithms Ruig Time of a algorithm Ruig Time Upper Bouds Lower Bouds Examples Mathematical facts Iput Algorithm Output A algorithm is a step-by-step procedure for solvig a problem i a fiite

More information

South Slave Divisional Education Council. Math 10C

South Slave Divisional Education Council. Math 10C South Slave Divisioal Educatio Coucil Math 10C Curriculum Package February 2012 12 Strad: Measuremet Geeral Outcome: Develop spatial sese ad proportioal reasoig It is expected that studets will: 1. Solve

More information

GPUMP: a Multiple-Precision Integer Library for GPUs

GPUMP: a Multiple-Precision Integer Library for GPUs GPUMP: a Multiple-Precisio Iteger Library for GPUs Kaiyog Zhao ad Xiaowe Chu Departmet of Computer Sciece, Hog Kog Baptist Uiversity Hog Kog, P. R. Chia Email: {kyzhao, chxw}@comp.hkbu.edu.hk Abstract

More information

Combination Labelings Of Graphs

Combination Labelings Of Graphs Applied Mathematics E-Notes, (0), - c ISSN 0-0 Available free at mirror sites of http://wwwmaththuedutw/ame/ Combiatio Labeligs Of Graphs Pak Chig Li y Received February 0 Abstract Suppose G = (V; E) is

More information

ECE4050 Data Structures and Algorithms. Lecture 6: Searching

ECE4050 Data Structures and Algorithms. Lecture 6: Searching ECE4050 Data Structures ad Algorithms Lecture 6: Searchig 1 Search Give: Distict keys k 1, k 2,, k ad collectio L of records of the form (k 1, I 1 ), (k 2, I 2 ),, (k, I ) where I j is the iformatio associated

More information

Course Site: Copyright 2012, Elsevier Inc. All rights reserved.

Course Site:   Copyright 2012, Elsevier Inc. All rights reserved. Course Site: http://cc.sjtu.edu.c/g2s/site/aca.html 1 Computer Architecture A Quatitative Approach, Fifth Editio Chapter 2 Memory Hierarchy Desig 2 Outlie Memory Hierarchy Cache Desig Basic Cache Optimizatios

More information

Octahedral Graph Scaling

Octahedral Graph Scaling Octahedral Graph Scalig Peter Russell Jauary 1, 2015 Abstract There is presetly o strog iterpretatio for the otio of -vertex graph scalig. This paper presets a ew defiitio for the term i the cotext of

More information

COMPUTER ORGANIZATION AND DESIGN The Hardware/Software Interface. Chapter 4. The Processor. Part A Datapath Design

COMPUTER ORGANIZATION AND DESIGN The Hardware/Software Interface. Chapter 4. The Processor. Part A Datapath Design COMPUTER ORGANIZATION AND DESIGN The Hardware/Software Iterface 5 th Editio Chapter The Processor Part A path Desig Itroductio CPU performace factors Istructio cout Determied by ISA ad compiler. CPI ad

More information

Chapter 1. Introduction to Computers and C++ Programming. Copyright 2015 Pearson Education, Ltd.. All rights reserved.

Chapter 1. Introduction to Computers and C++ Programming. Copyright 2015 Pearson Education, Ltd.. All rights reserved. Chapter 1 Itroductio to Computers ad C++ Programmig Copyright 2015 Pearso Educatio, Ltd.. All rights reserved. Overview 1.1 Computer Systems 1.2 Programmig ad Problem Solvig 1.3 Itroductio to C++ 1.4 Testig

More information

Parabolic Path to a Best Best-Fit Line:

Parabolic Path to a Best Best-Fit Line: Studet Activity : Fidig the Least Squares Regressio Lie By Explorig the Relatioship betwee Slope ad Residuals Objective: How does oe determie a best best-fit lie for a set of data? Eyeballig it may be

More information

arxiv: v2 [cs.ds] 24 Mar 2018

arxiv: v2 [cs.ds] 24 Mar 2018 Similar Elemets ad Metric Labelig o Complete Graphs arxiv:1803.08037v [cs.ds] 4 Mar 018 Pedro F. Felzeszwalb Brow Uiversity Providece, RI, USA pff@brow.edu March 8, 018 We cosider a problem that ivolves

More information

EE 459/500 HDL Based Digital Design with Programmable Logic. Lecture 13 Control and Sequencing: Hardwired and Microprogrammed Control

EE 459/500 HDL Based Digital Design with Programmable Logic. Lecture 13 Control and Sequencing: Hardwired and Microprogrammed Control EE 459/500 HDL Based Digital Desig with Programmable Logic Lecture 13 Cotrol ad Sequecig: Hardwired ad Microprogrammed Cotrol Refereces: Chapter s 4,5 from textbook Chapter 7 of M.M. Mao ad C.R. Kime,

More information

A General Framework for Accurate Statistical Timing Analysis Considering Correlations

A General Framework for Accurate Statistical Timing Analysis Considering Correlations A Geeral Framework for Accurate Statistical Timig Aalysis Cosiderig Correlatios 7.4 Vishal Khadelwal Departmet of ECE Uiversity of Marylad-College Park vishalk@glue.umd.edu Akur Srivastava Departmet of

More information

Operating System Concepts. Operating System Concepts

Operating System Concepts. Operating System Concepts Chapter 4: Mass-Storage Systems Logical Disk Structure Logical Disk Structure Disk Schedulig Disk Maagemet RAID Structure Disk drives are addressed as large -dimesioal arrays of logical blocks, where the

More information

Lighting and Shading. Outline. Raytracing Example. Global Illumination. Local Illumination. Radiosity Example

Lighting and Shading. Outline. Raytracing Example. Global Illumination. Local Illumination. Radiosity Example CSCI 480 Computer Graphics Lecture 9 Lightig ad Shadig Light Sources Phog Illumiatio Model Normal Vectors [Agel Ch. 6.1-6.4] February 13, 2013 Jerej Barbic Uiversity of Souther Califoria http://www-bcf.usc.edu/~jbarbic/cs480-s13/

More information

Evaluation of Different Fitness Functions for the Evolutionary Testing of an Autonomous Parking System

Evaluation of Different Fitness Functions for the Evolutionary Testing of an Autonomous Parking System Evaluatio of Differet Fitess Fuctios for the Evolutioary Testig of a Autoomous Parkig System Joachim Wegeer 1 ad Oliver Bühler 2 1 DaimlerChrysler AG, Research ad Techology, Alt-Moabit 96 a, D-10559 Berli,

More information

Mobile terminal 3D image reconstruction program development based on Android Lin Qinhua

Mobile terminal 3D image reconstruction program development based on Android Lin Qinhua Iteratioal Coferece o Automatio, Mechaical Cotrol ad Computatioal Egieerig (AMCCE 05) Mobile termial 3D image recostructio program developmet based o Adroid Li Qihua Sichua Iformatio Techology College

More information

Force Network Analysis using Complementary Energy

Force Network Analysis using Complementary Energy orce Network Aalysis usig Complemetary Eergy Adrew BORGART Assistat Professor Delft Uiversity of Techology Delft, The Netherlads A.Borgart@tudelft.l Yaick LIEM Studet Delft Uiversity of Techology Delft,

More information

On (K t e)-saturated Graphs

On (K t e)-saturated Graphs Noame mauscript No. (will be iserted by the editor O (K t e-saturated Graphs Jessica Fuller Roald J. Gould the date of receipt ad acceptace should be iserted later Abstract Give a graph H, we say a graph

More information

Computer Systems - HS

Computer Systems - HS What have we leared so far? Computer Systems High Level ENGG1203 2d Semester, 2017-18 Applicatios Sigals Systems & Cotrol Systems Computer & Embedded Systems Digital Logic Combiatioal Logic Sequetial Logic

More information

A SOFTWARE MODEL FOR THE MULTILAYER PERCEPTRON

A SOFTWARE MODEL FOR THE MULTILAYER PERCEPTRON A SOFTWARE MODEL FOR THE MULTILAYER PERCEPTRON Roberto Lopez ad Eugeio Oñate Iteratioal Ceter for Numerical Methods i Egieerig (CIMNE) Edificio C1, Gra Capitá s/, 08034 Barceloa, Spai ABSTRACT I this work

More information

performance to the performance they can experience when they use the services from a xed location.

performance to the performance they can experience when they use the services from a xed location. I the Proceedigs of The First Aual Iteratioal Coferece o Mobile Computig ad Networkig (MobiCom 9) November -, 99, Berkeley, Califoria USA Performace Compariso of Mobile Support Strategies Rieko Kadobayashi

More information

ISSN (Print) Research Article. *Corresponding author Nengfa Hu

ISSN (Print) Research Article. *Corresponding author Nengfa Hu Scholars Joural of Egieerig ad Techology (SJET) Sch. J. Eg. Tech., 2016; 4(5):249-253 Scholars Academic ad Scietific Publisher (A Iteratioal Publisher for Academic ad Scietific Resources) www.saspublisher.com

More information

Filter design. 1 Design considerations: a framework. 2 Finite impulse response (FIR) filter design

Filter design. 1 Design considerations: a framework. 2 Finite impulse response (FIR) filter design Filter desig Desig cosideratios: a framework C ı p ı p H(f) Aalysis of fiite wordlegth effects: I practice oe should check that the quatisatio used i the implemetatio does ot degrade the performace of

More information

Fast Interpolation of Grid Data at a Non-Grid Point

Fast Interpolation of Grid Data at a Non-Grid Point Fast Iterpolatio of Grid Data at a No-Grid Poit Hiroshi Ioue IBM Research - Tokyo Tokyo, Japa iouehrs@jp.ibm.com Abstract Defiig data at a o-grid poit by iterpolatig grid data is a commo operatio i may

More information

The Magma Database file formats

The Magma Database file formats The Magma Database file formats Adrew Gaylard, Bret Pikey, ad Mart-Mari Breedt Johaesburg, South Africa 15th May 2006 1 Summary Magma is a ope-source object database created by Chris Muller, of Kasas City,

More information

. Perform a geometric (ray-optics) construction (i.e., draw in the rays on the diagram) to show where the final image is formed.

. Perform a geometric (ray-optics) construction (i.e., draw in the rays on the diagram) to show where the final image is formed. MASSACHUSETTS INSTITUTE of TECHNOLOGY Departmet of Electrical Egieerig ad Computer Sciece 6.161 Moder Optics Project Laboratory 6.637 Optical Sigals, Devices & Systems Problem Set No. 1 Geometric optics

More information

Descriptive Statistics Summary Lists

Descriptive Statistics Summary Lists Chapter 209 Descriptive Statistics Summary Lists Itroductio This procedure is used to summarize cotiuous data. Large volumes of such data may be easily summarized i statistical lists of meas, couts, stadard

More information

WestminsterResearch

WestminsterResearch WestmisterResearch http://www.westmister.ac.uk/westmisterresearch Parallelizig the Chambolle Algorithm for Performace-Optimized Mappig o FPGA Devices Beretta, I., Raa, V., Aki, A., Nacci, A. A., Sciuto,

More information

An Efficient Algorithm for Graph Bisection of Triangularizations

An Efficient Algorithm for Graph Bisection of Triangularizations Applied Mathematical Scieces, Vol. 1, 2007, o. 25, 1203-1215 A Efficiet Algorithm for Graph Bisectio of Triagularizatios Gerold Jäger Departmet of Computer Sciece Washigto Uiversity Campus Box 1045, Oe

More information

A Resource for Free-standing Mathematics Qualifications

A Resource for Free-standing Mathematics Qualifications Ope.ls The first sheet is show elow. It is set up to show graphs with equatios of the form = m + c At preset the values of m ad c are oth zero. You ca chage these values usig the scroll ars. Leave the

More information

The Nature of Light. Chapter 22. Geometric Optics Using a Ray Approximation. Ray Approximation

The Nature of Light. Chapter 22. Geometric Optics Using a Ray Approximation. Ray Approximation The Nature of Light Chapter Reflectio ad Refractio of Light Sectios: 5, 8 Problems: 6, 7, 4, 30, 34, 38 Particles of light are called photos Each photo has a particular eergy E = h ƒ h is Plack s costat

More information

Designing a learning system

Designing a learning system CS 75 Machie Learig Lecture Desigig a learig system Milos Hauskrecht milos@cs.pitt.edu 539 Seott Square, x-5 people.cs.pitt.edu/~milos/courses/cs75/ Admiistrivia No homework assigmet this week Please try

More information

A Study on the Performance of Cholesky-Factorization using MPI

A Study on the Performance of Cholesky-Factorization using MPI A Study o the Performace of Cholesky-Factorizatio usig MPI Ha S. Kim Scott B. Bade Departmet of Computer Sciece ad Egieerig Uiversity of Califoria Sa Diego {hskim, bade}@cs.ucsd.edu Abstract Cholesky-factorizatio

More information

Load balanced Parallel Prime Number Generator with Sieve of Eratosthenes on Cluster Computers *

Load balanced Parallel Prime Number Generator with Sieve of Eratosthenes on Cluster Computers * Load balaced Parallel Prime umber Geerator with Sieve of Eratosthees o luster omputers * Soowook Hwag*, Kyusik hug**, ad Dogseug Kim* *Departmet of Electrical Egieerig Korea Uiversity Seoul, -, Rep. of

More information

Sorting in Linear Time. Data Structures and Algorithms Andrei Bulatov

Sorting in Linear Time. Data Structures and Algorithms Andrei Bulatov Sortig i Liear Time Data Structures ad Algorithms Adrei Bulatov Algorithms Sortig i Liear Time 7-2 Compariso Sorts The oly test that all the algorithms we have cosidered so far is compariso The oly iformatio

More information

An Efficient Algorithm for Graph Bisection of Triangularizations

An Efficient Algorithm for Graph Bisection of Triangularizations A Efficiet Algorithm for Graph Bisectio of Triagularizatios Gerold Jäger Departmet of Computer Sciece Washigto Uiversity Campus Box 1045 Oe Brookigs Drive St. Louis, Missouri 63130-4899, USA jaegerg@cse.wustl.edu

More information

Lecture Notes 6 Introduction to algorithm analysis CSS 501 Data Structures and Object-Oriented Programming

Lecture Notes 6 Introduction to algorithm analysis CSS 501 Data Structures and Object-Oriented Programming Lecture Notes 6 Itroductio to algorithm aalysis CSS 501 Data Structures ad Object-Orieted Programmig Readig for this lecture: Carrao, Chapter 10 To be covered i this lecture: Itroductio to algorithm aalysis

More information

n n B. How many subsets of C are there of cardinality n. We are selecting elements for such a

n n B. How many subsets of C are there of cardinality n. We are selecting elements for such a 4. [10] Usig a combiatorial argumet, prove that for 1: = 0 = Let A ad B be disjoit sets of cardiality each ad C = A B. How may subsets of C are there of cardiality. We are selectig elemets for such a subset

More information

Stone Images Retrieval Based on Color Histogram

Stone Images Retrieval Based on Color Histogram Stoe Images Retrieval Based o Color Histogram Qiag Zhao, Jie Yag, Jigyi Yag, Hogxig Liu School of Iformatio Egieerig, Wuha Uiversity of Techology Wuha, Chia Abstract Stoe images color features are chose

More information

How do we evaluate algorithms?

How do we evaluate algorithms? F2 Readig referece: chapter 2 + slides Algorithm complexity Big O ad big Ω To calculate ruig time Aalysis of recursive Algorithms Next time: Litterature: slides mostly The first Algorithm desig methods:

More information

Fast algorithm for skew detection. Adnan Amin, Stephen Fischer, Tony Parkinson, and Ricky Shiu

Fast algorithm for skew detection. Adnan Amin, Stephen Fischer, Tony Parkinson, and Ricky Shiu Fast algorithm for skew detectio Ada Ami, Stephe Fischer, Toy Parkiso, ad Ricky Shiu School of Computer Sciece ad Egieerig Uiversity of New South Wales, Sydey NSW, 2052 Australia ABSTRACT Documet image

More information

End Semester Examination CSE, III Yr. (I Sem), 30002: Computer Organization

End Semester Examination CSE, III Yr. (I Sem), 30002: Computer Organization Ed Semester Examiatio 2013-14 CSE, III Yr. (I Sem), 30002: Computer Orgaizatio Istructios: GROUP -A 1. Write the questio paper group (A, B, C, D), o frot page top of aswer book, as per what is metioed

More information

Graphics (Output) Primitives. Chapters 3 & 4

Graphics (Output) Primitives. Chapters 3 & 4 Graphics (Output) Primitives Chapters 3 & 4 Graphic Output ad Iput Pipelie Sca coversio coverts primitives such as lies, circles, etc. ito pixel values geometric descriptio Þ a fiite scee area Clippig

More information