2017-03-22 1 views
1

그래서, CT에서 궤도 (목표 x, y, z 및 접선으로 표시됨)를 따라 밀도 프로파일 (HU)을 계산하려고합니다. 현재 나는 표적을 통과하는 선을 따라 그리고 표적으로부터 일정한 거리 (입구)를 따라 프로파일을 얻을 수있다. 내가하고 싶은 것은 너비 1mm 정도의 볼륨 (이 경우 실린더)에 대한 밀도 프로파일을 얻는 것입니다. 연속 좌표 사이의 간격에 따라 여러 좌표가 동일한 인덱스를 가리킬 수 있기 때문에 보셀을 따라 보간을해야합니다. 예를 들어, 이것이 제가 말하고자하는 것입니다.궤도 볼륨을 따라 HU 값 가져 오기

Trajectories

또한, 나는 예를 들어, 궤도의 끝의 다른 모양의 밀도 프로파일을 좀하고 싶습니다 :

Different shapes

내 생각은 내가 3으로 만들 것입니다 3 행렬을 사용하여 팁의 모양을 나타내고,이를 보셀 값과 결합하여 팁에 해당하는 HU 값을 얻습니다. ITK/VTK를 사용하여 이것을 어떻게 할 수 있습니까?

좀 더 자세한 정보가 필요하면 알려주세요. (이미지가 충분히 선명하기를 바랍니다.)

답변

0

밀도 드릴 팁을 계산하려면, 이미지보다 높은 해상도에서 팁 절단 표면의 마스크를 만드는 것이 가장 쉽습니다. 드릴 이미지를 CT 이미지의 원하는 위치에 놓는 변환 행렬 M을 정의하십시오.

그리고, 보간을 사용하여, 그들 M 변환 마스크에있는 모든 비 - 제로 복셀 물리적 점 transform 지수, apply 통해 그 때의 CT 이미지 값을 샘플 (evaluate)를 반복에 의해 곱 마스크의 불투명도 (비 바이너리 마스크의 경우)를 계산하여 누적 합계에 값을 추가합니다.

결국 합계가 발생 밀도를 나타냅니다. 이 밀도의 합은 팁의 절삭면 마스크의 해상도에 따라 달라집니다. 나는 그것을 당신이 어떤 물리적 양 (뉴턴의 힘에 저항하는 것)과 어떻게 관련시킬 지 모르겠습니다.

+0

정확히 내가 언급 한대로 생각하고 있었던 것이지만 내 이미지보다 더 높은 해상도의 가면을 만들어야하는 이유는 무엇입니까? 또한 이전 답변에서 보간자를 사용하지 않고 이미지를 리샘플링 할 필요가 없다고 언급했습니다.이미지보다 높은 해상도의 가면을 만들면 보간 기가 필요하지 않습니다. –

+0

드릴 팁의 마스크를 이미지보다 거친 해상도로 만들면 마스크가 일부 보셀을 건너 뜁니다. 드릴 팁 마스크 해상도를 이미지 해상도보다 훨씬 높게 설정하면 계산 능력을 낭비하게됩니다. 나는 최적이 2 배 (Nyquist-Shannon theorem) 또는 약간 (2.1x) 더 높다고 가정합니다. 리샘플링은 인터폴 레이터를 사용합니다. 처음 샘플을 리샘플링 한 다음 리샘플링 된 이미지를 계산에 사용하면 실제로 리샘플링이 두 번 반복됩니다. –

+0

예제 코드를 제공 할 수 있습니까? –

0

일부 경로를 따라 프로파일을 얻으려면 resample filter을 사용하십시오. 시작점을 0,0,0으로 변환하고 끝점을 x, 0,0으로 변환하는 변환 행렬을 설정하십시오. 대상 이미지의 크기를 x, 1,1로 설정하고 원본 이미지와 동일한 간격으로 지정하십시오.

두 번째 질문을 이해할 수 없습니다. 팁에서 HU 값을 얻으려면 high quality interpolator (example using linear interpolator)을 사용하여 해당 포인트를 샘플링합니다. 팁의 모양이 왜 중요한지 알 수 없습니다.

+1

답변 해 주셔서 감사합니다. 팁의 모양은 드릴이 만날 뼈/조직의 밀도를 기준으로합니다. 당신이 상상할 수 있듯이, 다른 각도의 원추형 팁 (이미지에서 보여지는 것처럼)은 다른 것과 비교하여 약간 다른 밀도를 만나게 될 것입니다. –

+0

또한 내 용도로 필터를 사용하는 방법에 대한 예제 코드/링크를 제공 할 수 있습니까? 당신이 언급 한 아이디어는 내가 생각했던 것입니다. 그러나 문서가 너무 일반적입니다. –

+0

원래 게시물에 드릴을 언급하지 않았습니다. Drill Tip에서 발생할 밀도를 계산하려면 먼저 CT를 작은 이미지로 다시 샘플링 할 필요가 없습니다. 이는 이중 리샘플링을 희생하면서 모델링을 단순화합니다. –