Singular Value Decomposition of a Large XYZ file

Question

Oguzhan M 2018-8-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/414832-singular-value-decomposition-of-a-large-xyz-file

评论： Oguzhan M 2018-8-21

Hi everyone, I need to perform Singular Value Decomposition on a XYZ file (surface scan file with 3 columns and millions of row) to simplify the surface. In the papers I have checked in the literature, they present a 2D picture (x-z base) colouring with amplitude (y-height), and state that the fluctuating surface elevation of the surface y(x,z) should be decomposed. Any idea how to turn very narrow randomly lined xyz values into a matrix that we can implement SVD? Cheers.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Christine Tobler 2018-8-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/414832-singular-value-decomposition-of-a-large-xyz-file#answer_332743

First thing that occurs to me would be to interpolate the xyz data to lie on a rectangular grid in the x-z plane, and to then apply SVD to the resulting matrix (millions of rows should translate to a number of rows and columns in the thousands?).

Here is a doc example for how to interpolate xyz-data on a grid.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Oguzhan M 2018-8-16

Thanks for the prompt response. When I interpolate xyz data, I get 3 matrices as expected. Should I apply SVD to x and z matrices separately to reduce number of points, if yes, how can I re-combine new x and z matrices with y matrix (unchanged) to achieve the complete simplified surface? Thanks, Oguzhan.

请先登录，再进行评论。

Answer 2

Christine Tobler 2018-8-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/414832-singular-value-decomposition-of-a-large-xyz-file#answer_333022

在 MATLAB Online 中打开

The way I thought of it was to use x and z as the 2d grid on which to interpolate, and y as the function value. Something like this:

 x = randn(10, 1); y = randn(10, 1); z = randn(10, 1);
 [xi,zi] = meshgrid(linspace(-1, 1), linspace(-1, 1));
 yi = griddata(x,z,y,xi,zi); % yi is a matrix with the values of y interpolated on the grid (xi, zi)

In this case, you would do SVD of the yi data, I think?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Oguzhan M 2018-8-21

Sorry for the late response. You were right! SVD should be applied on yi then the final matrix can be re-composed using orthogonal matrices and inverse of the diagonal matrix. I'm now facing "out of memory error", which does not make any sense to me as there is a lot more RAM on my PC than matlab requires. .Anyways, I'll raise it on another topic. Thanks, Oguzhan.

请先登录，再进行评论。