Question

我有一个带有（x，y）坐标的点集及其在矩阵a中的相应权重，其中第1，第2和第3列分别是x，y和重量。我想将此点集划分为网格单元格，并计算每个网格中的点数和每个网格的总重量。

我尝试了下面的小例子，但它没有用。在这里，我尝试将此数据集划分为2x2小网格，并尝试计算点数及其权重总和。此外，我有大数据集，所以当我需要不同的网格步长时，我不能进一步扩展这种方法。

有人可以帮助我开发一种更简单的方法吗？

function dataTree
count=zeros(9,1);
avg=zeros(9,1);
data=[1 3 100; 2 1 120; 3 5 110; 4 2 100; 5 3 150; 6 2 100];

for i=1:6
    if data(i,1)<=2
        for j=1:6
            if data(j,2)<=2
                count(1) = count(1) + 1;
                avg(1) = avg(1) + data(j,3);
            elseif data(j,2)<=4
                    count(2) = count(2) + 1;
                    avg(2) = avg(2) + data(j,3);
             elseif data(j,2)<=6
                    count(3) = count(3) + 1;
                    avg(3) = avg(3) + data(j,3);
            end
        end
    elseif data(i,1)<=4
        for j=1:6
            if data(j,2)<=2
                count(4) = count(4) + 1;
                avg(4) = avg(4) + data(j,3);
            elseif data(j,2)<=4
                    count(5) = count(5) + 1;
                    avg(5) = avg(5) + data(j,3);
             elseif data(j,2)<=6
                    count(6) = count(6) + 1;
                    avg(6) = avg(6) + data(j,3);
            end
        end
    elseif data(i,1)<=6
        for j=1:6
            if data(j,2)<=2
                count(7) = count(7) + 1;
                avg(7) = avg(7) + data(j,3);
            elseif data(j,2)<=4
                    count(8) = count(8) + 1;
                    avg(8) = avg(8) + data(j,3);
             elseif data(j,2)<=6
                    count(9) = count(9) + 1;
                    avg(9) = avg(9) + data(j,3);
            end
        end

    end
end
count'
avg'

Answer 1

如果您的x和y尚未四舍五入到任意单位，请先执行此操作：

x = round((x - min(x))/edgelength+1);

这确保您获得具有edgelength大小的正方形的网格，其由非零整数表示。对y执行相同的操作。

然后，您可以使用sparse或accumarray来获取总重量。 sparse更快，但适用范围较小：

gridWeight = sparse(x,y,weight);

如果您想获得平均权重，请为每个条目添加1并除以该矩阵：

NumEntries = sparse(x,y,1);
MeanWeights = gridWeight./NumEntries;

accumarray可以一次完成这两项操作：

gridWeight = accumarray([x y],weight);
MeanWeights = accumarray([x y], weight,[],@mean); %//add ,[], 'issparse' for sparse matrix

请注意，通过设置sparse，accumarary是accumarray=([i,j],val,[],@sum,[],'issparse')的子功能。唯一可以处理sparse的函数是@sum，它的唯一填充值是0，而对于accumarray，可以使用其他函数和值。

用于分析网格上数据的Matlab代码

1 个答案: