Ожидаемое количество точек в кубическом объеме с учетом GMM

Question

У меня есть трехмерное пространство, которое я дискретизировал в вокселях (кубах объема). У меня также есть набор 3D-точек в таком пространстве. Я хочу знать ожидаемое количество очков в данном вокселе. Я выбрал 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);