def point_to_index(point, origin, spacing):
    """Transform image data point coordinates to voxel."""
    i = (point[0] - origin[0]) / spacing[0]
    j = (point[1] - origin[1]) / spacing[1]
    k = (point[2] - origin[2]) / spacing[2]
    return (i, j, k)