Convert x y coordinates and z data to Matrix

36 次查看(过去 30 天)
I have three vectors
>> whos x y VAR1
Name Size Bytes Class Attributes
x 21242x1 679744 double
y 21242x1 679744 double
var1 21242x1 679744 double
where X and Y represent the location of the corresponding value of var1. What Im trying to do is to convert all this data to a matrix.
What I've done until now is:
a1 = 21242;
m = zeros(a1,a1).*NaN;
indexes = sub2ind([a1,a1],x,y);
m(indexes) = var1
Based on this question but I'm getting an error
Error using sub2ind (line 43)
Out of range subscript.
Which, according to another question about this error, is because I'm trying access elements which does not exist on the matrix.
Any help?
Thanks!
  5 个评论
user20912
user20912 2022-9-15
Thank for the input. However, it's still not working.
I have tried to element by element with a for cycle:
for i1=1:a1
indexes = sub2ind([a1 a1], abs(x(i1)), y(i1));
i1
end
and the problem is in the element 3897:
>> x(3896:3898)
-1.0021
-0.9993
-1.0026
what's the problem with that numbers? I can not find the reason why keep failing.

请先登录,再进行评论。

采纳的回答

Star Strider
Star Strider 2022-9-15
It would help to have the actual vectors.
It may be possible to use the reshape function to create matrices from the vectors, if the vectors are from a gridded original matrix. The easiest way to determine the dimensions of the arrays is to use the unique function and look at the second output. The elements of the second output should be regularly-spaced, and that difference is one dimension to specify in the reshape call, with the other one left as the empty vector [].
Otherwise, use linspace, ndgrid and griddata to create appropriate matrices from scattered data.
.
  6 个评论
user20912
user20912 2022-9-19
Thank you!
This was helpfull to do what I needed to do.

请先登录,再进行评论。

更多回答(1 个)

Robert Daly
Robert Daly 2023-8-11
I would use interpolation to convert to a regular grid.
scatteredInterpolant is the funtion to do this with your irregularly spaced data.
[X,Y] = meshgrid(linspace(min(x),max(x),500),linspace(min(y),max(y),500));% make 500x500 grid
f = scatteredInterpolant(x,y,var1); % set up the interpolater with your data
Grid_Var1 = f(X,Y); % Do the interpolation onto the grid set up above

类别

Help CenterFile Exchange 中查找有关 Matrices and Arrays 的更多信息

产品


版本

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by