3

提供了来自 DICOM 标头的以下信息，我如何计算体素大小的第三个值？我假设前两个值是 0.515625 和 0.515625。

BitsAllocated: "16"
BitsStored: "12"
Columns: 512
HighBit: "11"
ImageOrientation: "1\0\0\0\-1\0"
ImagePosition: "-144\-34.7242241\925.599976"
ImageType: "ORIGINAL\PRIMARY\AXIAL\HELIX"
InstanceNumber: "456"
Modality: "CT"
PhotometricInterpretation: "MONOCHROME2"
PixelRepresentation: "0"
PixelSpacing: "0.515625\0.515625"
RescaleIntercept: "1"
Rows: 512
SamplesPerPixel: "1"
SeriesDescription: "CERVEAU SANS IV"
SeriesNumber: "3"
SliceThickness: "1.50"
WindowCenter: "00040\00040"
WindowWidth: "00120\00120"
imagesFormat: "jpg"
modality: "CT"
name: "soft tissue"
nodeId: "557621"
pixelHeight: "0.515625"
pixelWidth: "0.515625"

注意：我收到的是 JPEG 图像堆栈，而不是 DICOM，它带有一个文件，其中包含我在上面发布的值。如果需要，我可以返回并要求提供文件中的其他信息。

4

2 回答 2

5

仅给定一个切片的标签，您必须将其SliceThickness用作第三维，尽管我建议不要这样做，因为这不能保证给出切片之间的距离。有SpacingBetweenSlices提供此信息的标签，但在您的情况下似乎不存在。

最好的方法是使用ImagePositionPatient相邻切片之间的差异。为此，您当然还需要下一个切片的标签。附带说明：在您的列表中，ImageOrientation应该ImagePosition更好地阅读ImageOrientationPatient和ImagePositionPatient，ImageOrientation以及ImagePosition其他标签（CT 图像中不存在）。

ImagePositionPatient给出 DICOM 患者坐标中切片左上角的位置，计算距离时必须考虑切片在该坐标系中的方向。这由 给出ImageOrientationPatient，其中包含 DICOM 坐标中切片的归一化行和列方向余弦向量。您可以在DICOM 标准中阅读。

方向矩阵的前两个分量由ImageOrientationPatient（例如第一个和第二个三个数字）提供，第三个分量可以通过取这两个分量的叉积来计算。

所以，在伪代码中，这看起来像这样：

orient1 = vector(ImageOrientationPatient[0], ImageOrientationPatient[1], ImageOrientationPatient[2])
orient2 = vector(ImageOrientationPatient[3], ImageOrientationPatient[4], ImageOrientationPatient[5])
orient3 = orient1 x orient2 // cross product
orient_matrix = matrix(orient1, orient2, orient3)

pos1 = vector(ImagePositionPatient[0], ImagePositionPatient[1], ImagePositionPatient[2]) // from current slice
pos2 = vector(ImagePositionPatient[0], ImagePositionPatient[1], ImagePositionPatient[2]) // from adjacent slice
diff_pos = pos2 - pos1

image_pos = orient_matrix o diff_pos / length(orient3) // normalized dot product

voxel_z = image_pos.z

更新：正如@gofal 所指出的，第一个版本不正确。我还包括了规范化（例如 delete by length(orient3)），尽管严格来说这些值应该已经规范化了。

于 2020-12-21T09:38:08.767 回答
4

PixelSpacing (0028,0030) 为您提供图像中 Columns 和 Row 方向的两种尺寸。您要求的第三个值是切片厚度 (0018,0050)。

在某些重建（MIP 或类似）的情况下，切片厚度的值可能高于两个相邻图像的距离。因此，如果您将此单个图像用作自己的图像，那么您将不得不使用切片厚度。但是，如果您有一堆图像并且想要进行 3d 渲染或重建或类似操作，那么我建议您通过两个相邻图像的距离来计算第三个体素大小。

参考文献： https ://dicom.innolitics.com/ciods/ct-image/image-plane http://dicom.nema.org/medical/dicom/current/output/chtml/part03/sect_C.34.12.html

于 2020-12-21T09:18:33.833 回答