• Journal of Biomedical Engineering Research
• /
• v.20 no.6
• /
• pp.523-529
• /
• 1999

방출컴퓨터단층영상술을 위한 영상재구성법에 있어서 기대값 최대화(EM)를 사용한 maximum likeihood 방법이 기존의 filtered backprojection 방법에 비해 현저한 장점을 지니고 있다는 점에서 지속적으로 그 가치가 인정되어 왔다. 그러나, 이러한 방법은 projection 및 backprojection 의 반복계산을 요하므로 영상재구성을 위한 총 계산시간이 projector 및 backprojector 의 성능에 크게 좌우된다. 본 논문에서는 EM에 근거한 영상재구성 알고리즘의 계산량을 감소시키는 방법에 관하여 논한다. 특히, projection 및 backprojection 계산을 위한 행렬의 원소중 중요한 양들을 구하는 방법과 이들을 미리 계산하여 적절한 양의 메모리에 저장하는 방법에 관하여 고찰한다. 실험에서 제안된 방법을 사용할 경우 EM 알고리즘의 계산시간을 92%까지 현저히 감소시킬 수있음을 보였다.

• Journal of The Korean Radiological Technologist Association
• /
• v.29 no.1
• /
• pp.150-150
• /
• 2003

• KIPS Transactions on Computer and Communication Systems
• /
• v.12 no.1
• /
• pp.41-46
• /
• 2023

Iterative reconstruction of CT takes a long time because projection and back-projection are alternatively repeated until taking a good image. To reduce the reconstruction time, we need a fast algorithm for calculating the projection which is a time-consuming step. In this paper, we proposed a new algorithm to calculate the line integral and the algorithm is approximately 10% faster than the well-known Siddon method (Jacobs version) and has a good image quality. Although the algorithm has been investigated for the case of parallel beams, it can be extended to the case of fan and cone beam geometries in the future.

• Journal of Biomedical Engineering Research
• /
• v.24 no.5
• /
• pp.473-480
• /
• 2003

Iterative reconstruction methods have played a prominent role in emission computed tomography due to their remarkable advantages over the conventional filtered backprojection method. However, since iterative reconstructions typically are comprised of repeatedly projecting and backprojecting the data, the computational load required for reconstructing an image depends highly on the performance of the projector-backprojector pair used in the algorithm. In this work we compare quantitative performance of representative methods for implementing projector-backprojector pairs. To reduce the overall cost for the projection-backprojection operations for each method, we investigate how previously computed results can be reused so that the number of redundant calculations can be minimized. Our experimental results demonstrate that the ray tracing method not only outperforms other methods in computation time, but also provides improved reconstructions with good accuracy.

• Proceedings of the Korea Information Processing Society Conference
• /
• 2022.05a
• /
• pp.312-315
• /
• 2022

CT의 반복재구성방법은 투영과 역투역을 번갈아가며 최적의 단면영상을 얻을 때까지 반복한다. 영상재구성 시간을 단축하기 위하여 시간이 많이 소요되는 투영을 빠르게 수행할 수 있는 알고리즘이 필요하다. 본 논문은 Siddon 알고리즘을 개선한 Jacobs 버전보다 대략 10% 빠른 알고리즘을 제안한다. 평행빔의 경우에 대해 조사되었지만 향후 부채살빔 및 콘빔의 경우로 확장이 가능하다.

• Journal of Biomedical Engineering Research
• /
• v.21 no.5
• /
• pp.441-448
• /
• 2000

방출단층영상 재구성을 위한 최대우도 기대값최대화(maximum likelihood expectation maximization, MLEM) 방법은 영상 획득과정을 통계학적으로 모델링하여 영상을 재구성한다. MLEM은 일반적으로 사용하여 여과후역투사(filtered backprojection)방법에 비해 많은 장점을 가지고 있으나 반복횟수 증가에 따른 발산과 재구성 시간이 오래 걸리는 단점을 가지고 있다. 이 논문에서는 이러한 단점을 보완하기 위해 계산시간을 현저히 단축시킨 배열된부분집합 기대값최대화(ordered subsets expectation maximization. OSEM)에 Gibbs 선행치인 membrance (MM) 또는 thin plate(TP)을 첨가한 OSEM-MAP (maximum a posteriori)을 구현함으로써 알고리즘의 안정성 및 재구성된 영상의 질을 향상시키고자 g나다. 실험에서 알고리즘의 수렴시간을 가속화하기 위해 투사 데이터를 16개의 부분집합으로 분할하여 반복연산을 수행하였으며, 알고리즘의 성능을 비교하기 위해 소프트웨어 모형(원숭이 뇌 자가방사선, 수학적심장흉부)을 사용한 영상재구성 결과를 제곱오차로 비교하였다. 또한 알고리즘의 사용 가능성을 평가하기 위해 물리모형을 사용하여 PET 기기로부터 획득한 실제 투사 데이터를 사용하였다.

• Journal of Korea Multimedia Society
• /
• v.12 no.3
• /
• pp.372-382
• /
• 2009

When the range of projection angles is limited, tomographic reconstruction suffers from artifacts caused by incomplete data. One can consider a data completion technique, which estimates projection data at unobserved angles using a prior knowledge or mathematical exploration, but the result is often not improved; the improvement by the data completion often undermined by the artifacts by inaccurate estimation, In this paper, we propose an iterative method, which computes projection data at unobserved angles by using the current estimate on the image, links the computed projection data to the observed ones by using the consistence condition of Radon transform, and reconstruct the next estimate on the image by filtered backprojection. The proposed method does not require a prior knowledge on the image, and has much faster approximation rate than the expectation maximization method. The performance of the proposed method was tested through several simulation studies.

• Journal of the Korean Society of Radiology
• /
• v.17 no.1
• /
• pp.141-148
• /
• 2023

The algebraic reconstruction technique is one of the reconstruction methods in CT and shows good image quality against noise-dominant conditions. The number of iteration is one of the key factors determining the execution time for the algebraic reconstruction technique. However, there are some rules for determining the number of iterations that result in more than a few hundred iterations. Thus, the rules are difficult to apply in practice. In this study, we proposed a method to determine the number of iterations for practical applications. The reconstructed image quality shows slow convergence as the number of iterations increases. Image quality 𝜖 < 0.001 was used to determine the optimal number of iteration. The Shepp-Logan head phantom was used to obtain noise-free projection and projections with noise for 360, 720, and 1440 views were obtained using Geant4 Monte Carlo simulation that has the same geometry dimension as a clinic CT system. Images reconstructed by around 10 iterations within the stop condition showed good quality. The method for determining the iteration number is an efficient way of replacing the best image-quality-based method, which brings over a few hundred iterations.

• Proceedings of the KIEE Conference
• /
• 2006.11a
• /
• pp.272-274
• /
• 2006

본 논문에서는 한국전력공사의 배전계통 운영계획시스템(Distribution Line PLANning system)에서 배전계통 재구성 최적화에 적용하고 있는 선로간의 부하 평준화 알고리즘에 대하여 연구하였다. 기존의 알고리즘은 계통 내 임의의 연계 개폐기 혹은 구분 개폐기를 조작한 후 방사상의 운전조건을 만족하면 선로별 부하를 계산한다. 그리고 최대부하 선로와 최소부하 선로간의 부하 편차를 구하고 직전의 최적해와 비교하여 최적의 해를 찾아내는 방법이다. 그러나 계통이 커질 경우 시험해야할 경우의 수가 많아져서 시간이 많이 소요되고 또한 수렴속도를 빠르게 하기위하여 반복횟수를 제한하기 때문에 절대적인 최적의 해를 보장하지 못하였다. 따라서 본 연구에서는 탐색시간을 만족시키면서 전수조사하는 정도의 부하 평준화 최적해를 찾아서 배전계통의 재구성을 최적화하는 알고리즘을 구현하고자 한다.

• Nuclear Medicine and Molecular Imaging
• /
• v.43 no.5
• /
• pp.459-467
• /
• 2009

Purpose: The maximum likelihood-expectation maximization (ML-EM) is the statistical reconstruction algorithm derived from probabilistic model of the emission and detection processes. Although the ML-EM has many advantages in accuracy and utility, the use of the ML-EM is limited due to the computational burden of iterating processing on a CPU (central processing unit). In this study, we developed a parallel computing technique on GPU (graphic processing unit) for ML-EM algorithm. Materials and Methods: Using Geforce 9800 GTX+ graphic card and CUDA (compute unified device architecture) the projection and backprojection in ML-EM algorithm were parallelized by NVIDIA's technology. The time delay on computations for projection, errors between measured and estimated data and backprojection in an iteration were measured. Total time included the latency in data transmission between RAM and GPU memory. Results: The total computation time of the CPU- and GPU-based ML-EM with 32 iterations were 3.83 and 0.26 see, respectively. In this case, the computing speed was improved about 15 times on GPU. When the number of iterations increased into 1024, the CPU- and GPU-based computing took totally 18 min and 8 see, respectively. The improvement was about 135 times and was caused by delay on CPU-based computing after certain iterations. On the other hand, the GPU-based computation provided very small variation on time delay per iteration due to use of shared memory. Conclusion: The GPU-based parallel computation for ML-EM improved significantly the computing speed and stability. The developed GPU-based ML-EM algorithm could be easily modified for some other imaging geometries.