使用2个坐标列和权重列的密度图

Question

我有一个3列的矩阵。前两列是坐标，第三列是重量或强度。

newmat = [ 27.37 -45.69 14.47 
      27.37 -45.68 18.58 
      27.37 -45.67 29.05 
      27.37 -45.66 51.7 
      ...  ...  ... ] 

我已经创建了一个散点图：

不过，我想有类似的密度图（如第二个图here）。我曾尝试使用hist3函数，如here，但我没有弄清楚如何考虑第三列 - 权重。

Answer 1

你可以（使用功能sortrows，unique，并且accumarray）创建从数据newmat矩阵，并绘制它作为image：

newmat = sortrows(newmat, [1 2]); % Sort the first two columns in ascending order 
[x, ~, newmat(:, 1)] = unique(newmat(:, 1));  % Make numeric indices for column 1 
[y, ~, newmat(:, 2)] = unique(newmat(:, 2));  % Make numeric indices for column 2 
M = accumarray(newmat(:, 1:2), newmat(:, 3)).'; % Build the matrix 
imagesc(x, y, M); 

这里是类似于您的格式一些示例数据：

[X, Y] = meshgrid(0:0.1:2, 3:0.1:5); 
Z = peaks(21); 
newmat = [X(:) Y(:) Z(:)]; 

下面是上面的代码从数据产生的情节：