[논문 리뷰] Tomographic X-ray data of a lotus root filled with attenuating objects

Tatiana A. Bubba, Andreas Hauptmann외 3명|arXiv (Cornell University)|2016. 09. 23.

Chromatography in Natural Products인용 수 6

한 줄 요약

이 논문은 희박한 데이터 기반 CT 재구성 알고리즘의 테스트를 가능하게 하기 위해 다양한 감쇠 물체(연필, 흙연기, 세라믹, matches의 머리)로 채워진 연꽃 뿌리의 공개된 2차원 단층촬영 X선 데이터셋을 제시한다. 데이터셋에는 고해상도 시노그램, 측정 행렬, 그리고 기준 FBP 재구성 결과가 포함되어 있어, 제한된 투과각도나 노이즈가 있는 조건에서 정규화 기반 이미지 재구성 기술에 대한 연구를 지원한다.

ABSTRACT

This is the documentation of the tomographic X-ray data of a lotus root, filled with four different attenuating objects, of different sizes. Data are available at www.fips.fi/dataset.php, and can be freely used for scientific purposes with appropriate references to them, and to this document in http://arxiv.org/arXiv. The data set consists of (1) the X-ray sinogram of a single 2D slice of the lotus root with two different resolutions and (2) the corresponding measurement matrices modeling the linear operation of the X-ray transform. Each of these sinograms was obtained from a measured 360-projection fan-beam sinogram by down-sampling and taking logarithms. The original (measured) sinogram is also provided in its original form and resolution.

연구 동기 및 목표

• 희박한 데이터 기반 단층촬영 알고리즘의 벤치마킹을 위한 현실적이고 물리적으로 확보된 CT 데이터셋을 제공함으로써, 복잡한 감쇠 구조를 가진 데이터셋을 확보한다.
• 특히 총변위 또는 티코노프 정규화에 의존하는 정규화 기반 재구성 방법의 개발 및 검증을 지원한다.
• 실제 대상에서 다양한 대비와 물체 기하학적 형태를 가지며 제한된 투과각 또는 낮은 투과수치 조건에서 재구성 알고리즘의 테스트를 가능하게 한다.
• 재현 가능한 알고리즘 평가를 위해 알려진 기준값과 校정된 측정 행렬을 갖춘 표준화된 데이터셋을 제공한다.
• 특정 해상도 및 행렬 노름 특성을 갖는 데이터를 제공함으로써, 이론적 정규화 파rameter 선택 규칙(예: [2]에서 제시된 것들)의 적용을 촉진한다.

제안 방법

• 몰리브덴 X선 소스, 회전 스테이지, 평판 검출기를 갖춘 맞춤형 마이크로-CT 시스템을 사용하여 360개의 프로젝션 팬빔 CT 데이터를 확보하였다.
• 2D 프로젝션의 중심 수평 단면을 추출하여 2221×120의 시노그램을 구성하였으며, 이후 저해상도로 변환하고 로그 변환하여 120개의 프로젝션 시노그램을 생성하였다.
• 128×128 및 256×256 해상도의 두 가지 변형을 생성하였으며, 저해상도 시노그램과 해당하는 희박한 측정 행렬 A(크기 51480×16384 및 51480×65536)를 사용하였다.
• 원본 360개 프로젝션 시노그램(2221×360)과 이 데이터로부터 생성된 고해상도 FBP 재구성 결과(1500×1500)를 기준값으로 제공하였다.
• 시스템 행렬 A의 노름과 상한값을 저장하여 정규화(A/normA)를 가능하게 하여, 단위 연산자 노름을 요구하는 알고리즘에 대한 호환성을 확보하였다.
• 기준값 재구성 결과를 재현하고 데이터 사용법을 시연하기 위한 MATLAB 스크립트(예: FBPgroundtruthRec.m)를 제공하였다.

실험 결과

연구 질문

• RQ1 희박한 프로젝션 수를 가진 실제 CT 데이터셋에서 정규화 기반 재구성 알고리즘이 미세한 구조나 고대비 특징을 얼마나 잘 복원할 수 있는가?
• RQ2 자연적인 생물학적 샘플(다양한 물질을 포함한 연꽃 뿌리)의 복잡하고 비균일한 감쇠 특성이 기존의 재구성 방법에 얼마나 큰 도전을 가하는가?
• RQ3 알려진 시스템 행렬 노름과 해상도 제약 조건을 갖는 실제 CT 데이터에 대해 이론적 정규화 파rameter 선택 규칙을 효과적으로 적용할 수 있는가?
• RQ4 다양한 물체 크기, 형태, 물질 대비를 갖는 데이터셋에 적용했을 때, 티코노프 및 총변위 정규화의 성능은 어떻게 달라지는가?
• RQ5 예를 들어 120개의 프로젝션에서 20개의 프로젝션으로의 데이터 저해상도화 및 프로젝션 수 감소가 물리적으로 확보된 데이터셋에서 재구성 정밀도에 어떤 영향을 미치는가?

주요 결과

• 128×128 및 256×256 해상도의 두 가지 변형이 있으며, 이에 해당하는 측정 행렬과 시노그램을 포함하여, 다양한 공간 해상도에서 재구성 알고리즘의 평가가 가능하다.
• 원본 360개 프로젝션 시노그램(2221×360)은 50 kV 및 1 mA 조건에서 각 프로젝션 당 1초 노출로 확보되었으며, 기준값 생성을 위한 고신호대비노이즈비율을 확보하였다.
• 기준값 FBP 재구성 결과(1500×1500)는 360개 프로젝션 데이터로부터 계산되었으며, 알고리즘 검증을 위한 고정밀 기준으로 기능한다.
• 측정 행렬 A는 알려진 노름과 상한값을 갖는 희박 행렬로 저장되어 있어, 알고리즘 호환성을 위한 단위 연산자 노름 정규화가 가능하다.
• 시노그램은 2D 프로젝션의 중심 행에서 유도되었으며, 픽셀 크기는 128×128 기준 0.627 mm, 256×256 기준 0.315 mm로, 공간 해상도의 물리적 정확성을 확보하였다.
• 시연 스크립트는 120개 프로젝션을 사용한 티코노프 정규화 재구성에서 명확한 구조가 나타나며, 20개 프로젝션의 경우 적절한 정규화를 통해 여전히 복원 가능한 결과를 보여준다.