[논문 리뷰] Tomographic X-ray data of a walnut

K. Hämäläinen, Lauri Harhanen외 4명|arXiv (Cornell University)|2015. 02. 11.

Nuclear Physics and Applications인용 수 19

한 줄 요약

이 논문은 고해상도 시노그램, 측정 행렬, 그리고 기준 FBP 재구성 결과를 포함하여 공개된 호두의 단층촬영 X선 데이터 세트를 제시한다. 이는 희박한 데이터 기반 및 정규화 기반 CT 재구성 알고리즘의 벤치마킹을 가능하게 하며, 파동수 기반 희박성 방법, 총 변동량 정규화, 그리고 저투사 수준에서의 S-곡선 파rameter 선택법 테스트를 지원한다.

ABSTRACT

This is the documentation of the tomographic X-ray data of a walnut made available at http://www.fips.fi/dataset.php . The data can be freely used for scientific purposes with appropriate references to the data and to this document in arXiv. The data set consists of (1) the X-ray sinogram of a single 2D slice of the walnut with three 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 120-projection fan-beam sinogram by down-sampling and taking logarithms. The original (measured) sinogram is also provided in its original form and resolution. In addition, a larger set of 1200 projections of the same walnut was measured and a high-resolution filtered back-projection reconstruction was computed from this data; both the sinogram and the FBP reconstruction are included in the data set, the latter serving as a ground truth reconstruction.

연구 동기 및 목표

• 의료 및 산업 영상에서 희박한 데이터 재구성 알고리즘 테스트를 위한 실제 세계의 CT 데이터를 제공하기 위해.
• 단층촬영에서의 제한된 투사 데이터 문제를 해결하기 위해 고대비 및 비볼록한 구조를 지닌 실제적이고 복잡한 생물 샘플을 제공하기 위해.
• 낮은 신호 영역에서 Tikhonov 및 총 변동량 정규화와 같은 정규화 기법의 개발 및 검증을 지원하기 위해.
• 1200개의 투사로부터 기준 재구성 결과를 제공하여 재구성 정확도의 정량적 평가를 가능하게 하기 위해.
• 측정 기하학 및 해상도 수준이 명확히 정의된 표준화되고 재현 가능한 데이터 세트를 제공함으로써 역문제 연구의 방법론적 발전을 촉진하기 위해.

제안 방법

• 80 kV, 200 µA, 50 µm 픽셀 해상도를 갖춘 맞춤형 µCT 시스템을 사용하여 호두의 120개 투사 패널-빔 CT 데이터를 확보하였다.
• 2차원 투사 이미지의 중심 행을 추출하여 중심 2차원 단면을 나타내는 2296×120 패널-빔 시노그램을 구성하였다.
• 원본 시노그램을 빈닝을 통해 저해상도로 변환하고 로그 변환을 적용하여 세 가지 해상도 수준(82×82, 164×164, 328×328)을 생성하였다.
• 선형 X선 변환을 모델링하기 위해 크기가 N²×M인 희박한 측정 행렬 A를 구성하였다. 여기서 M은 투사 수, N은 재구성 픽셀 수이다.
• 더 높은 정확도를 확보하기 위해 0.3° 각도 간격으로 1200개 투사의 두 번째 고해상도 데이터 세트를 확보하고, 필터링된 역백산법(FBP)을 사용하여 고해상도 기준 재구성 결과를 계산하였다.
• MATLAB 호환 데이터 파일(시노그램, 측정 행렬, FBP 재구성 결과 포함)을 제공하여 알고리즘 테스트 및 검증에 직접 활용할 수 있도록 하였다.

실험 결과

연구 질문

• RQ1 희박성 촉진 재구성 방법은 제한된 투사 수를 가진 실제 복잡한 생물학적 데이터에서 얼마나 잘 성능을 발휘할 수 있는가?
• RQ2 희박한 데이터 영역에서 정규화 파rameter 선택이 재구성 품질에 어떤 영향을 미치는가?
• RQ3 측정 행렬의 해상도 선택이 반복적 재구성 알고리즘의 정확도와 안정성에 어떤 영향을 미치는가?
• RQ4 1200개 투사로부터 유도된 기준 재구성 결과가 희박한 데이터 방법의 평가에 얼마나 신뢰할 수 있는 기준이 되는가?
• RQ5 층상 껍질과 비볼록한 내부 구조를 지닌 호두의 구조적 복잡성은 임상적 또는 산업적 CT 과제를 효과적으로 시뮬레이션할 수 있는가?

주요 결과

• 해상도 수준 82×82, 164×164, 328×328에서 해당 측정 행렬 및 시노그램 데이터를 포함하는 세 개의 저해상도로 변환된 시노그램이 포함되어 있다.
• 원본 120개 투사 시노그램의 치수는 2296×120이며, 세 가지 해상도 수준에서 픽셀 크기는 각각 0.513 mm, 0.257 mm, 0.128 mm이다.
• 0.3° 각도 간격으로 1200개 투사로부터 고해상도 FBP 재구성 결과를 계산하여 2296×2296 기준 이미지를 생성하였다.
• 측정 행렬 A는 희박하며 선형 시스템 A*x = m(:)를 표현하므로 Tikhonov 및 총 변동량 정규화 프레임워크에서 활용할 수 있다.
• 이 데이터 세트는 파동수 기반 희박성 방법 및 정규화 파rameter 선택을 위한 S-곡선 방법 테스트에 성공적으로 사용되었다.
• 정확한 인용을 포함한 과학적 목적을 위한 무료 공개로, 역문제 연구의 재현 가능성 및 벤치마킹을 지원한다.