2017-02-22 82 views
1

我相對較新與dicom files.Thanks提前工作。DICOM:對齊

我有不同的時間間隔拍攝的同一個病人的2個dicom文件。 它們的尺寸不完全相同。

第一個是:cube1104X163X140dimesions,第二個是CUBE2107X164X140dimesions。我想在原點對齊這兩個立方體並進行比較。

所述第一文件的ImagePositionPatient爲:[-207.4748,-151.3715 -198.7500]

所述第二文件的ImagePositionPatient爲:[-207.4500,-156.3500 -198.7500]

兩個文件有相同的ImageOrientationPatient - [1 0 0 0 1 0]

有人可以請任何機會給我看一個例子嗎?我不知道如何將物理平面映射回圖像平面?

非常感謝提前,

=============================== ================================

補充:23/2/17

我已經使用下面的矩陣公式基於我的情況下的鏈接:

立方體1的IPP(Sxyz) = [-2 07.4748,-151.3715-198.7500]

Xxyz(IOP) = [1,0,0]

Yxyz(IOP) = [1,0,0]

delta_i = 2.5

delta_j = 2.5

enter image description here

因此,對於cube1的i = 0:103和j = 0:162的值,我應該計算出Pxyz的值。

下一步是什麼?對不起,我不明白這將如何幫助我將具有不同IPP的兩個立方體對準圖像平面?

很抱歉的新手問題...

+0

在這裏看到:http://stackoverflow.com/questions/36996353/finding-the-coordinates-mm-of-identical- slice-locations-for-mr-datasets-ac/36997168#36997168 and here: http://stackoverflow.com/questions/40008507/how-to-spatial-order-files-in-a- dicom-data-sequence/40014740#40014740 –

+1

請注意,上述建議的鏈接中沒有一個(截至2017年2月22日)描述了「引用UID框架」,因此基於IOP和IPP的對齊只有在兩者的屬性相同時纔有意義系列。 – malat

+0

@kritzel_sw任何機會我可以得到一個例子?我是否需要計算矩陣的逆矩陣 - 我不知道如何將物理平面映射回像平面? – Ash

回答

1

我沒有驗證您構建的矩陣。但是,如果計算正確,則可以在體積座標系(VCS)(x1,y1,z1)(其中x1 = column,y1 = row且z1 =切片號)之間轉換爲患者座標系(PCS)(x2, y2,z2) - 這些座標以毫米爲單位定義患者體內的點。

通過反轉矩陣,您可以從PCS轉換回VCS。假設,體積1的變換矩陣:= M1,體積2的變換矩陣PCS - > VCS:= M2。然後,可以通過使用M1將其轉換爲PCS,並使用M2'(反轉的M2)將PCS轉換爲第2卷,將第1卷中的點p1轉換爲第2卷中的相應點p2。

通過乘以M1和M2',可以計算直接從volume1轉換爲volume2的矩陣。

所以:

P2 =(M1 * M2' )* P1