
Compressed Sensing for Dynamic PETDynamic positron emission tomography (PET) is widely used to measure changes in the biodistribution of radiopharmaceuticals within particular organs of interest over time. However, to retain sufficient temporal resolution, photon counts in each time frame should be limited, so conventional reconstruction algorithms such as the ordered subset expectation maximization (OSEM) produce noisy reconstructions. To address this problem, many advanced reconstruction algorithms have been developed using various spatiotemporal regularizations. In this work, we propose a novel spatiotemporal regularization approach that exploits inherent similarities in inter and intra frames. In such correlations, we are interested in a low rank constraint of spatiotemporal patches that is less sensitive to global intensity variations, which is originated from the matrix completion in compressed sensing (CS) theory. To solve nonconvex problem, we also provide an optimization method using the concaveconvex procedure (CCCP) by exploiting the LegendreFenchel transform. We confirm that the proposed algorithm can improve image quality as well as extract accurate kinetic parameters.
Interior Tomography for Xray CT In xray computed tomography (CT), reconstruction of region of interest (ROI) from local projection data has been called for to reduce radiation exposure in imaging specific organs such as heart, or to reduce size of xray detector for cost saving. Using backprojection filtration (BPF) algorithm, after backprojecting the differentiated projection data on chord lines, the original 2D or 3D reconstruction problem is converted into a set of 1D Hilbert filtration problems. If exact differentiated backprojection (DBP) signal on the chord line can cover the whole object, the object can be reconstructed by 1D finite inverse Hilbert transform of DBP along the chord line. But, in interior problems, enough exact DBP signal is not available, so this Hilbert transform relationship needs to be solved by more complicated method. Now, we are trying to find accurate and fast algorithm to solve this truncated Hilbert transform problem using total variation (TV) penalty and multiscale decomposition, and the simulation result shows the proposed algorithm can produce high quality reconstruction with significantly reduced computational time to practical level even for helical conebeam CT.
