1

我正在使用 ct 扫描原始格式的医学图像。它基本上是一个 3d 体素矩阵（512*512*nb 个切片）。我想将文件的每个切片提取到单独的文件中。

import numpy as np
import matplotlib.pyplot as plt

# reading the raw image into a string. The image files can be found at:
# https://grand-challenge.org/site/anode09/details/
f = open('test01.raw', 'rb')
img_str = f.read()

# converting to a uint16 numpy array
img_arr = np.fromstring(img_str, np.uint16)

# get the first image and plot it
im1 = img_arr[0:512*512]
im1 = np.reshape(im1, (512, 512))
plt.imshow(im1, cmap=plt.cm.gray_r)
plt.show()

结果绝对看起来像胸部 ct 扫描，但图像的纹理很奇怪，好像像素放错了位置。

在此处输入图像描述

一些相关信息可能位于关联的 .mhd 信息文件中，但我不确定在哪里查看：

ObjectType = Image
NDims = 3
BinaryData = True
BinaryDataByteOrderMSB = False
CompressedData = False
TransformMatrix = 1 0 0 0 1 0 0 0 1
Offset = 0 0 0
CenterOfRotation = 0 0 0
AnatomicalOrientation = RPI
ElementSpacing = 0.697266 0.697266 0.7
DimSize = 512 512 459
ElementType = MET_SHORT
ElementDataFile = test01.raw
4

1 回答 1

2

试试这种方式：

Dim_size=np.array((512,512,459),dtype=np.int) #Or read that from your mhd info File

f = open(FileName,'rb') #only opens the file for reading
img_arr=np.fromfile(f,dtype=np.uint16)
img_arr=img_arr.reshape(Dim_size[0],Dim_size[1],Dim_size[2])

如果您的内存有限，请以块的形式读取文件

f = open(FileName,'rb') #only opens the file for reading
for i in range(0,Dim_size[2]):
    img_arr=np.fromfile(f,dtype=np.uint16,count=Dim_size[0]*Dim_size[1])
    img=img.reshape(Dim_size[0],Dim_size[1])
    #Do something with the Slice

显示原始文件中实际内容的一个好方法也是在 ImageJ 中读取它。对于读取此类 ITK 兼容文件，甚至还有可用的插件，但直接原始导入也应该可以工作。 https://imagej.net/ 欢迎 http://ij-plugins.sourceforge.net/plugins/3d-io/

于 2017-03-28T10:41:56.313 回答