给我一段python计算三维图像dice的代码

以下是一段计算三维图像Dice系数的Python代码：

import numpy as np

def dice_coefficient(img1, img2):
    intersection = np.sum(img1 & img2)
    dice = (2. * intersection) / (np.sum(img1) + np.sum(img2))
    return dice

其中，`img1`和`img2`是两个三维图像矩阵，可以使用`numpy`库来表示和处理。这个函数首先计算两个图像的交集，也就是两个图像都为1的像素点的数量，然后用Dice系数公式计算Dice系数。最后返回计算结果。注意，这个函数只适用于二进制图像，即图像中的像素值只能是0或1。