การวัดในมิติเวลา

การวิเคราะห์ภาพทางชีวการแพทย์ด้วย Python

Stephen Bailey

Instructor

Ejection fraction

$$ Ejection\ Fraction = \frac{LV_{max} - LV_{min}}{LV_{max}} $$

การเต้นของหัวใจ

Ejection fraction

ขั้นตอน

แบ่งส่วนช่องซ้ายของหัวใจ
คำนวณปริมาตรสำหรับแต่ละ 3D volume ในอนุกรมเวลา
เลือกค่าต่ำสุดและสูงสุด
คำนวณ ejection fraction

คำนวณปริมาตรในแต่ละจุดเวลา

# Stored in (t,z,x,y) format
vol_ts.shape
labels.shape

(20, 12, 256, 256)
(20, 12, 256, 256)

กราฟอนุกรมเวลา

# Calculate voxel volume in mm^3
d0,d1,d2,d3=vol_ts.meta['sampling']
dvoxel = d1 * d2 * d3

# Instantiate empty list
ts = np.zeros(20)

# Loop through volume time series
for t in range(20):

    nvoxels=ndi.sum(1,
                    labels[t], 
                    index=1)

    ts[t] = nvoxels * dvoxel

plt.plot(ts)
plt.show()

คำนวณ ejection fraction

min_vol = ts.min()
max_vol = ts.max()

ejec_frac = (max_vol - min_vol) / max_vol

ejec_frac

0.58672

ปริมาตรสูงสุดและต่ำสุดแบบเปรียบเทียบ

มาฝึกกันเถอะ!

การวิเคราะห์ภาพทางชีวการแพทย์ด้วย Python