in divergent beam ramp fbp reconstruction algorithms , it is known that the properties such as noise , resolution and divergent beam artifacts of reconstructed image are influenced by the position - dependent weight in the backprojection . the backprojection weight ( which is related inversely to the distance of the source to the point being reconstructed ) causes spatially nonuniform distribution of noise and resolution , and the divergent beam artifacts are more pronounced away from the rotation axis [ 13 ] . pan [ 2 , 4 ] used shift variant ramp filtering approach while wang et al .   proposed spatially variant weighting in the backprojection . recently , divergent beam fbp reconstruction algorithms have been proposed with inverse distance weight   and no backprojection weight [ 7 , 8 ] . the reconstructed images with these algorithms have better properties in terms of reduced noise level and uniformity in noise distribution and uniformity in resolution as compared to the classical ramp fbp reconstruction algorithm . the fan - beam reconstruction algorithms with inverse distance backprojection weight have been extended to obtain 3d cone - beam ( cb ) fdk - type approximate reconstruction algorithms with improved noise and resolution properties . on the other hand , algorithms have been developed to improve the noise performance with the modified redundancy weight in the projection domain for circular and helical cb ct [ 1012 ] . the fan - beam reconstruction algorithms with no backprojection weight can be extended to obtain 3d cb fdk - type approximate reconstruction algorithms to improve noise and resolution properties . the most widely used practical method for cb reconstruction has been the approximate algorithm of feldkamp , davis and kress ( fdk )   originally proposed for circular scan cb reconstruction , and fdk has been extended to helical scan cb geometry by wang et al . . the circular scan trajectory geometry does not satisfy tuy 's data sufficiency condition on the projections for exact and stable reconstruction [ 15 , 16 ] permitting only approximate reconstruction and also leads to intensity drop in the reconstructed images . many methods have been suggested to give improvements in the image quality when reconstructing from circular cb data . some of them use the iteration schemes to reconstruct the image with circular scan fdk as an intermediate result [ 17 , 18 ] . another approach to reduce the axial intensity drop is given by hu   which adds , to fdk reconstruction , a correction term obtained from unused circular scan cb data . propose to further add a  missing term  to obtain a more accurate fdk reconstruction . the missing term is obtained by estimating the missing data , required by the data sufficiency condition , from the cb scan projection data .   in this paper , two exact hilbert fbp reconstruction algorithms with no backprojection weight for 2d fan - beam with equispace linear array detector geometry ( complement of the equi - angular curved array detector ) have been derived . using these fan - beam hilbert fbp inversion algorithms , feldkamp - type reconstruction algorithms with no backprojection weight for 3d reconstruction for circular and helical further , we present efficient algorithms to reduce the axial intensity drop for circular cb ct using the proposed fdk reconstruction algorithms with no backprojection weight . simulation studies demonstrate the properties of the proposed reconstruction algorithms with no backprojection weight and the efficient algorithms for reducing axial intensity drop in 3d cb ct . in this section , we revisit the relation between the parallel beam projections and the fan - beam projections through the hilbert kernel . using this relation , we derive the fan - beam reconstruction formulae with no backprojection weight for equispace linear array detector , which are counterpart of the algorithms for the curved array detector [ 7 , 8 ] . let f(x ) denote the 2d object density function to be reconstructed where x=(x , y)t . the parallel - beam projections of f are given by the expression 
 ( 1)p(,s)=2f(x)(xs)dx , 
 where s   , =(-sin,cos ) is any unit vector and  is dirac delta function . the delta function in ( 1 ) selects the value of x that satisfies x=s , so that p(,s ) is the integral of f on the line perpendicular to  at a distance s from the origin x=(0,0 ) as shown in figure 1 . the inversion of f from its projection p using the fbp algorithm involves two steps . the first step is to filter the projection p(,s ) by applying 1d ramp filter in s , 
 ( 2)pr(,s)=(hrp)(,s)=p(,s)hr(ss)ds , 
 where h 
 r =  
  |  | exp(j2s)d is the impulse response of || filter . the ramp filter is seen as equivalent to the successive application of a derivative in s and hilbert transform ,  
 ( 3)pr(,s)=(hhsp)(,s)=12sph(,s ) , 
 where h 
 h(s ) = 1/2 
 j  sgn()exp(j2s)d. the second step in the inversion is the backprojection of filtered projections resulting from the application of hilbert filter . mathematically , 
 ( 4)f(x)=  1402sph(,s)|s = xd. 
  the fan - beam projections are obtained by using the fan - beam sampling by an equi - angular curved or equispaced collinear detector arrays . the fan - beam projection are measured by moving the x - ray source along a circular trajectory as shown in figure 2 . the source trajectory can be parameterized by an angular parameter  and is given by a()=(r0cos ,  r0sin ) where r 
 0 is the radius of the circle . the fan - beam projections can be represented by 
 ( 5)g(,)=0f(a()+t)dt ,  s1 , 
 where s 
  is the set of all possible unit vectors in 2d space . recalling that p 
 h denotes the hilbert transform of the 2d radon transform p of f. the line m corresponds to parallel beam ray and fan - beam ray as shown in figure 2 . mathematical representation of this line satisfies s  =  a() , where a( ) is a point on this line . the value of p 
 h corresponding to this line can be obtained from gh(, ) , and the relation between p 
 h and gh(, ) is given by 
 ( 6)ph(,s)|s = a()=gh(, ) , 
 where 
 ( 7)gh(,)=s1hh()g(,)d. 
 this relation is given by hamaker et al . . the use of this relation has been explained through the parameterization of =cosew+sineu and =-sinew+coseu , where eu and ew are unit vectors and these vectors are illustrated in figure 3 . the proof of ( 6 ) can be found in [ 6 , 21 ] . using ( 6 ) , we derive two reconstruction algorithms with no backprojection weight for fan - beam equispace linear array detector . let g( , u ) be the fan - beam projections for equispace detector and p( , s ) be the parallel - beam projections . using change of variable  = tan(u / r 
 0 ) in ( 7 ) , we define hilbert transform of g( , u)(8)gh(,u)=1r02+u2umumg(,u)uu1r02+u2du . 
 from the relation of hamaker ( 6 ) , =   tan(u/r 
 0 ) , and s = r0u/r02+u2 , we have g 
 h( , u ) = p 
 h( tan(u/r 
 0 ) , r 
 0 
 u/(r 
 0 
  + u 
 ) ) . the derivative of g 
 h( , u ) with respect to u and  using chain rule is given by 
 ( 9)ugh(,u)=[sus+u]ph(,s),gh(,u)=[]ph(,s ) . since the parallel beam parameters are functions of the fan - beam parameters ,  s/u , /u , and / are given by 
 ( 10)su=r03r02+u21r02+u2 ,  u=r0r02+u2 ,  =1 . substituting ( 10 ) in ( 9 ) , we get 
 ( 11)ugh(,u)=[r03r02+u21r02+u2s+r0r02+u2 ]  ph(tan1ur0,r0ur02+u2 ) , 
 
 ( 12)gh(,u)=ph(tan1ur0,r0ur02+u2 ) . 
  substituting ( 12 ) in ( 11 ) , ( 11 ) can be written as 
 ( 13)ugh(,u)=r0r02+u2[sph(tan1ur0,r0ur02+u2 )  r02r02+u2gh(,u) ] . 
 in ( 4 ) , ( /s)p 
 h is needed to reconstruct the function f , therefore ( 13 ) can be rewritten as 
 ( 14)sph(tan1ur0,r0ur02+u2 )  = ( ugh(,u)r02+u2r0+gh(,u))r02+u2r02 . this result is used in ( 4 ) to get the first fan - beam inversion formula with no backprojection weight for equispace linear array detector . the inversion formula is given by 
 ( 15)f(x)=1402[(ugh(,u)r02+u2r0+gh(,u ) )  r02+u2r02]d , 
 where  =   tan(u/r 
 0 ) and u=xr0/r02-(x)2 . the integration in this formula being over  , it is a parallel - beam backprojection . equation ( 15 ) is applicable for full scan , short scan , and very short scan . for the short scan and over scan , ( 15 ) degenerates to the 180 parallel - beam reconstruction . in the case of super short scan , you and zeng   suggested the redundancy weight similar to noo 's redundancy weighting . if the scanning range is [ 0   ] , the redundancy weight can be calculated as follows : 
 ( 16)w(,u)=wi(10.5wa ) , 
 where wi is equal to 1 if  = (  + tan(u/r 
 0 ) )  floor(( + tan ( u/r 
 0))/2)2   otherwise wi = 0 and wa is equal to 1 if  = (   tan(u/r 
 0 ) +  )  floor(( tan(u/r 
 0 ) + )/2)2   otherwise wa is equal to 0 . the reconstruction formula for this scanning range is given by 
 ( 17)f(x)=120w(,u)[(ugh(,u)r02+u2r0  + gh(,u))r02+u2r02]d. 
  noo et al . have proposed hilbert fbp fan - beam inversion formula which can also reconstruct with less than short scan . the inversion formula is given by 
 ( 18)f(x)=12021l[w(,x)gh(,u)]u=u(,x)d , 
 where 
 ( 19)gh(,u)=umumhh(uu)r0r02+u2  (+r02+u2r0u)g(,u)du , 
 ( 20)u(x , y,)=r0xeur0+xew , l = r0+xew . x. l is the distance from source position obtained by projecting the point x on the ray passing through the origin . illustration of u and l are shown figure 3 . since the redundancy weighting w( , u ) in ( 18 ) applies after the hilbert filtering , it can be applied in the image domain as part of the backprojection step . this is not feasible with the standard fbp formula because redundancy weight is applied before filtering . have suggested the valid redundancy weight for full scan , 
 ( 21)w(,x)=||xa()||2r0cos   , 
 where  = arctan(u/r 
 0 ) . from the geometry shown in figures 3 and 4 , it can be proved that w( ,  ) + w( 
 c ,  
 c ) = 1 , where w( ,  ) and w( 
 c ,  
 c ) are the redundancy weights corresponding to the redundant line k which passes through the point x. equation ( 21 ) , in terms of equispace array geometry parameters , can be written as 
 ( 22)w(,x)=lr02+u22r03 , 
 where l=||x-a()||cos=r0+xew . the second fan - beam inversion formula with no backprojection weight for equispace linear array detector is obtained by substituting ( 22 ) in ( 18 ) and is given as 
 ( 23)f(x)=1202r02+u22r03gh(,u)u=u(,x)d , 
 where g 
 h is given by ( 19 ) and u=r0xeu/(r0+xew ) . thus , the backprojection weight ( l ) in ( 18 ) can be eliminated for each ray that corresponds to a line integral that is measured twice by weighting each ray with function w(,x ) of ( 21 ) during the backprojection for short scan or super short scan . however , it does not allow elimination of the weight for singly measured rays ( i.e , in case of short scan or a super short scan ) . have developed 3d cb ct reconstruction algorithms for circular and helical trajectories , respectively , from the conventional fan - beam reconstruction formula . similarly , in this section , we propose two new fdk - type inversion formulae for 3d cb reconstruction with the flat panel array detector cb geometry which are obtained by heuristic extension of the fan - beam inversion formulae with no backprojection weight . these formulae are the counterpart of the formulae for cylindrical array detector given by narasimhadhan et al . let f(x , y , z ) be the object density to be reconstructed and g( , u , v ) be the flat panel detector cb projection data collected with a helical scan trajectory given by a()=[r0cos ,  r0sin,h]t where r 
 0 is the radius of the source trajectory , h = p/2 , and p represents the helical pitch . to describe the cb data for a flat detector , we introduce the local detector coordinates ( u , w , v ) with unit vectors : eu()=[-sin ,  cos,0]t , ew()=[-cos ,  -sin,0]t , and ev()=t . here ew( ) points from the source point a( ) to the center of the detector and the unit vectors eu( ) and ev( ) span the detector . the cb projection for flat panel detector geometry is given by 
 ( 24)g(,u , v )  = 0f(a()+t(ueu()+dew()+vev()u2+d2+v2))dt , 
 where d is the source to detector distance . to simplify the formulae given below , we assume that the flat panel virtual detector is parallel to the actual detector panel at a distance r 
 0 from a( ) and hence contains the z - axis . scaling all equations by substituting u r 
 0 
 u / d and v  r 
 0 
 v / d is sufficient to describe the case of a flat panel array detector located at a distance d from the vertex . the hfdkw-1 formula is obtained from ( 15 ) by the standard fdk extension technique for a particular slice z , 
 ( 25)f^(x , y , z )  = 14(z / h)(z / h)+r02+u2r03(g^(,u , v)hh(uu)du)d , 
 where 
 ( 26)g^(,u , v )  = r0r02+u2+v2(+r02+u2r0u+uvr0v)g(,u , v ) , 
 (x , y ,  ) =  + (x , y ,  ) , and (x , y ,  ) = arcsin  ( ( xsin + ycos)/r 
 0 ) . the variables u and v are computed by the following equations : 
 ( 27)u(x , y,)=r0(xsin(x , y,)+ycos(x , y,))r0xcosysin(x , y,),v(x , y,)=r0(zh(x , y,))r0xcos(x , y,)ysin(x , y, ) . the hfdkw-2 formula is obtained from the ( 23 ) by standard fdk extension technique 
 ( 28)f^(x , y , z )  = 14(z / h)(z / h)+r02+u2r03g^(,u , v)hh(uu)du  d , 
 where 
 ( 29)u(x , y,)=r0(xsin+ycos)r0xcosysin,v(x , y , z,)=r0(zh)r0xcosysin , 
 and g^(,u , v ) is given by ( 26 ) . note that the integration in second formula in ( 28 ) is over  and therefore it is a fan - beam backprojection while the integration in first hfdkw formula ( 25 ) being over  results in a parallel - beam backprojection . the first hfdkw formula ( 25 ) can be modified to include redundancy weight for partial scan data as given by you and zeng   ( 25 ) . the weight l 
  in the backprojection , where l is the distance from source position obtained by projecting the point x on the ray passing through the origin , is therefore position - dependent . the weight ( r 
 0 
  + u 
 ) /r 
 0 
  required in the hfdkw-2 can be evaluated and multiplied with filtered projection data before the backprojection . since there is no position dependent weight in the backprojection step , hfdkw-2 is computationally more efficient than the conventional fdk algorithm . hfdkw-1 algorithm performs an additional interpolation in the backprojection , hence the computational complexity of the hfdkw-1 is more than the conventional hfdk and hfdkw-2 algorithms . the hfdkw-1 and hfdkw-2 algorithms can be reduced to circular fdk algorithms with no backprojection weight by setting the helical pitch p = 0 . it is seen from the figure 7 that the cfdkw algorithms , hu 's algorithm , and t - fdk algorithm give similar performance in terms of axial intensity drop . this motivates us to develop an efficient approach using the estimate of missing term as given by zhu et al . . it has been shown by reformulating grangeat 's algorithm for circular cb trajectory that the reconstruction of the true object f can be written as the sum of three terms : 
 ( 30)f^=f^fdk+f^h+f^n , 
 where f^fdk is the reconstructed object using fdk algorithm , f^h is hu 's correction term , and f^n is the missing term . hu 's term represents the information contained in a circular cb scan but not utilized in the fdk reconstruction . one approach is to estimate the missing term from the available circular cb scan data as proposed by zhu et al . , which includes all the three terms in ( 30 ) . in this section , we propose to use the circular fdk inversion with no backprojection weight ( cfdkw algorithms ) and the estimate of missing term f^n using zhu et al . thus , the reconstructed object f^ of the true object f , with the two algorithms , cfdkw-1 and cfdkw-2 , respectively , can be written as : 
 ( 31)f^=f^cfdkw-1+f^n , f^=f^cfdkw-2+f^n . these two efficient algorithms are called em-1 and em-2 , respectively , for further reference . the estimation of missing term is based on hu 's theory   and grangeat 's theory . the formula for estimating the missing term   is given by 
 ( 32)f^n(x)=142r02+z2r02(1r02z2r0)02g^(,z)d , 
 where g^(,z)=(/z2)-(r0/r02+u2+z2)g(,u , z)du . cone - beam circular and helical x - ray source trajectory scanner is used in the computer simulations . the radius of the circular and helical trajectory of the source r 
 0 is 2.4  m and the object radius is 1  m for the simulations . circular cb projection data is simulated with 450 uniformly spaced source positions over 2 interval , using 283 detector elements , and 283 detector rows . helical cb projection data is simulated with 450 uniformly spaced source positions over 2 interval , using 283 detector elements , and 29 ( 111 ) detector rows with pitch value p = 0.125  m ( 0.5  m ) . image matrix for a slice of the 3d phantom is 256    256 . simulations are carried out on full , short , and very short scan data obtained from 3d shepp - logan low contrast phantom   for hfdkw algorithms . algorithm is based on fan - beam backprojection over 2 scan and therefore it can only handle full - scan reconstruction . hfdkw-1 is a more general algorithm which can accommodate both short scan as well as very short scan data . figure 5(a ) shows the reconstructed slice at z = 0.25 of shepp - logan low - contrast phantom using hfdkw-1 for short scan ( scan angle of  + fan angle ) and figure 5(b ) shows the reconstruction with very short scan case ( scan angle =  ) . the 3d shepp - logan low - contrast phantom , forbild head phantom and defrise disk phantom , are used to illustrate the performance of the proposed algorithms for 3d reconstruction with circular scan trajectories under different conditions . each disc has a uniform attenuation coefficient of 0.7 , and the ellipsoid has a diameter of 1.2  m and a thickness of 0.12  m , with a distance of 0.2  m between discs . we evaluate the performance of the cfdkw algorithms in terms of longitudinal or axial intensity drop . we compare the cfdkw algorithms with three different algorithms : the fdk algorithm which is only the first term of ( 30 ) , hu 's algorithm ( consists of first two terms of ( 30 ) ) , and t - fdk algorithm which has been developed heuristically with a structure of shift - invariant filtering . we also illustrate the performance of the em-1 and em-2 in terms of reduction in the axial intensity drop , noise level , and computational complexity . 's   algorithm ( i.e , all the three terms in ( 30 ) ) . the axial intensity drop and variation become significant in conventional helical fdk algorithms due to the approximate nature of the algorithm when large helical pitch is used . we reconstructed the slices of 3d shepp - logan low - contrast phantom with helical pitch 0.5  m using fdk , hfdkw-1 , and hfdkw-2 algorithms . it is observed from figures 6(a ) , 6(b ) , and 6(c ) that the reconstructed image using the helical fdk shows the axial intensity drop and variation compared to the reconstructed images using hfdkw algorithms . we study the axial intensity drop associated with the circular fdk algorithms due to data incompleteness of the circular trajectory in 3d reconstruction . the shepp - logan low - contrast phantom , forbild phantom , and the defrise disk phantom are used to illustrate the performance of the cfdkw algorithms , in comparison with fdk , hu 's , and t - fdk algorithms . the first column of the figure 7 shows the reconstruction of the shepp - logan low - contrast phantom with full scan circular trajectory using fdk , hu 's , t - fdk and cfdkw algorithms . it is observed that the drop in density in the case of hu 's , t - fdk , cfdkw-1 , and cfdkw-2 in figure 7 is same and significantly less than the circular fdk . the second column of the figure 7 shows the reconstruction of forbild phantom which has high - density object in the central region leading to artifacts . it is seen in figure 7 that the reconstructed images of forbild head phantom have white and black streaks using the fdk and hu 's algorithms . these streaks are significantly reduced in the reconstructed images using the cfdkw and t - fdk algorithms . the defrise disk phantom is used to show the performance of the algorithms to reconstruct sharp change in the axial intensity . the strong high - frequency components in the axial direction due to sharp change in density represents the most challenging aspect of the reconstruction using circular cb data . figure 9 shows the 1d axial intensity profile of defrise disk phantom illustrating the observations made earlier in figure 7 . it is seen from the axial intensity profile in figure 9 that only the central disk is reconstructed properly by all the algorithms . we reconstruct the three phantoms using the em algorithms and compare them with the reconstructions obtained by zhu et al . 's algorithm , em-1 , and em-2 algorithms of shepp - logan phantom in the first column , forbild phantom in the second column and defrise disk in the third column , respectively . 's , em-1 , and em-2 reconstruction algorithms give the same performance in terms of the axial intensity drop and significant improvement over reconstructions shown in figure 7 . this is also evident quantitatively from the axial intensity profiles shown in figures 11 and 12 of the shepp - logan and defrise disk phantoms in figure 10 , respectively . we simulate noisy projection data using the shepp - logan low - contrast phantom with poisson noise based on an emission of n 
 0 = 3  10 photons . the number of photons detected n was calculated by ln(n / n 
 0 ) . with this as the mean to a poisson random process , n for the case of noisy projection is computed , corresponding noisy projection is then given by ln(n / n 
 0 ) . the image is reconstructed using the eight different algorithms to make a qualitative and quantitative evaluation of the cfdkw and the em algorithms in comparison to other fdk algorithms as given in table 1 . figure 13 shows the reconstructed noisy images for the eight algorithms in the same order as given in table 1 . to compute the error variances of the images , the error images are obtained by taking the difference of reconstructed noisy images ( figure 13 ) and corresponding noise - free images as shown in figures 7 and 10 . the variances of error images are measured and the variances are given in the table 1 . it is seen from table 1 that the noise level is reduced in the cfdkw and em algorithms . this is due to the additional numerical differentiation step , which is usually implemented as two- or three - point difference , and these algorithms are based on hilbert filtering which avoids the approximation introduced by the nonuniform cutoff frequency required in the ramp filter - based fbp algorithm . since cfdkw-1 and em-1 algorithms perform additional interpolation in the source direction in the backprojection , these algorithms give better noise performance among all algorithms . we compute the execution time for fdk algorithm , hu 's term , missing term , hu 's algorithm , zhu et al . 's algorithm , cfdkw-2 , and em-2 algorithms in matlab 7.0 on a intel p4 desktop computer with a dual core 2.4  ghz processor platform . the algorithm modules considered for estimating the execution time are ( 1 ) 1d filtering of the cb projections and ( 2 ) backprojection of a single slice . the execution time of the reconstruction for fdk , hu 's term , missing term , zhu et al . 's algorithm , cfdk-2 , and em-2 algorithms are computed using these modules . the execution times are given in the table 2 . cfdkw-2 does not have position - dependent weight in the backprojection , hence cfdkw-2 algorithm has less computational complexity than fdk and hu 's algorithms . the em-2 algorithm consists of the reconstruction of the image using cfdkw-2 algorithm and estimating the missing term , whereas zhu et al . 's algorithm consists of the reconstruction of the image using conventional fdk , estimation of hu 's term , and estimation of missing term . we have given two fan - beam algorithms with no backprojection weight for equispace linear array detector , which are counterparts of the algorithms with no backprojection weight for equi - angular detector geometry . two new variants of fdk algorithms for circular and helical cb scan geometry with planar detector geometry are proposed . the cfdkw-1 and hfdkw-1 are applicable to full , short and very short scan data . the cfdkw and hfdkw algorithms are based on hilbert filtering and have no position - dependent weights in the backprojection step and hence reduction in the axial intensity drop and variation when compared with the classical fdk algorithm . the cfdkw , t - fdk , and hu 's algorithms give identical performance in the reduction of axial intensity drop in the circular cb ct . cfdkw-2 algorithm is more efficient than the fdk and hu 's algorithm in terms of computational complexity . we have given efficient algorithms to reduce the axial intensity drop in the circular scan cb ct 's algorithm give a similar performance in terms of the axial intensity drop . the em-2 algorithm noise performance is better in case of cfdkw and em algorithms than conventional fdk , hu 's algorithm , t - fdk and zhu et al .
we develop two feldkamp - type reconstruction algorithms with no backprojection weight for circular and helical trajectory with planar detector geometry . 
 advances in solid - state electronic detector technologies lend importance to ct systems with the equispaced linear array , the planar ( flat panel ) detectors , and the corresponding algorithms . 
 we derive two exact hilbert filtered backprojection ( fbp ) reconstruction algorithms with no backprojection weight for 2d fan - beam equispace linear array detector geometry ( complement of the equi - angular curved array detector ) . 
 based on these algorithms , the feldkamp - type algorithms with no backprojection weight for 3d reconstruction are developed using the standard heuristic extension of the divergent beam fbp algorithm . 
 the simulation results show that the axial intensity drop in the reconstructed image using the fdk algorithms with no backprojection weight with circular trajectory is similar to that obtained by using hu 's and t - fdk , algorithms . 
 further , we present efficient algorithms to reduce the axial intensity drop encountered in the standard fdk reconstructions in circular cone - beam ct . the proposed algorithms consist of mainly two steps : reconstruction of the object using fdk algorithm with no backprojection weight and estimation of the missing term . 
 the efficient algorithms are compared with the fdk algorithm , hu 's algorithm , t - fdk , and zhu et al . 
 's algorithm in terms of axial intensity drop and noise . 
 simulation shows that the efficient algorithms give similar performance in axial intensity drop as that of zhu et al . 
 's algorithm while one of the efficient algorithms outperforms zhu et al . 
 's algorithm in terms of computational complexity .