给定GMM时，立方体体积中的预期点数

Question

我有一个三维空间，我将其离散为体素（体积的立方体）。我在这样的空间里也有一组三维点。我想知道给定体素中的预期点数。为此，我选择了GMM作为模型，但我不知道如何从μ、西格玛和每个高斯的权重开始计算我想要的。

到目前为止，我设法适应了GMM（简单）：

obj = gmdistribution.fit(points', 20);

我通过绘制

figure(1);
hold on;
for i = 1:k
    plot_gaussian_ellipsoid(obj.mu(i,:), obj.Sigma(:,:,i));
end
axis equal;

这就是我所期望的，那是一张地图，颜色告诉我点的浓度。

问题是，给定体素的中心(x,y,z)和侧面s，我如何提取体素中的预期点数？

Answer 1

您可以使用（参见此处的示例http://www.mathworks.nl/help/stats/gmdistribution.cluster.html)

idx = cluster(gm,points);