How can I interpolate or spline a dataset into a curve?

Question

FsC 2021-11-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1583529-how-can-i-interpolate-or-spline-a-dataset-into-a-curve

评论： Mathieu NOE 2021-11-18

采纳的回答： Mathieu NOE

Data1.mat

在 MATLAB Online 中打开

Hello,

I have a dataset that I need to interpolate or create a curve from; however, the datasets include different curves and I'm unsure how to interpolate or create an average line from them. I've attached an examples from one of the datasets to this post. When plotting the points it looks like this:

% for Data1 dataset
plot(Data1(:,1),Data1(:,2),'b.')

Ideally, when connecting the points or plotting a line, the curve would look something like this:

However, When I connect the points in a line this is the result:

I thought perhaps an interpolation approach might work. If I use interp1, in the following code the result is:

xq = 0:1:length(Data1(:,1))-1; 
x = Data1(:,1);
y = Data1(:,2);
vq = interp1(x,y,xq,'nearest');
plot(xq,vq)

Am I applying the interp1() function incorrectly or should I be using another method? or is there a way to rearrange the matrix to get the line shown above?

Thank you for your help!

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

dpb 2021-11-10

By combining the multi-valued curves, the interpolation routines have fit the input data in strict ascending order of the points globally, not sequentially in the piecewise segments.

MATLAB doesn't have builtin tools designed for such problems, unfortunately; to use interp1 or spline you'll have to treat the segments as non-overlapping segments individually and then paste the results together. That likely won't be all that satisfactory at the disjoint locations as you have overlap there.

One trick that might help would be to reverse the sense of the x-y coordinates for the LH center section -- it looks to be univalued in the plotted y-direction so that one could smooth it in x- instead.

请先登录，再进行评论。

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

Answer 1

Mathieu NOE 2021-11-10

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1583529-how-can-i-interpolate-or-spline-a-dataset-into-a-curve#answer_828564

在 MATLAB Online 中打开

hello

this would be my suggestion

clc
clearvars
%% load('');
load('Data1.mat')
x = Data1(:,1);
y = Data1(:,2);
centroid_x = mean(x);
centroid_y = mean(y);
[theta,r] = cart2pol(x-centroid_x,y-centroid_y);
% sort theta in ascending order
[theta_sorted,ind] = sort(theta);
r_sorted = r(ind);
% remove duplicates before interpolation
[theta_unique,IA,IC] = unique(theta_sorted);
r_unique = r_sorted(IA);
% put angle in range 0 - 2pi
ind_neg = theta_unique<=0;
ind_pos = theta_unique>0;
theta_new = [theta_unique(ind_pos); theta_unique(ind_neg)+2*pi];
r_new = [r_unique(ind_pos); r_unique(ind_neg)];
% robust average ()
N = 15;  % adjust smoothing factor 
rr = smoothdata(r_new,'sgolay', N,'Degree',1);
[u,v] = pol2cart(theta_new,rr);
u = u + centroid_x;
v = v + centroid_y;
plot(x,y,'.b',u,v,'r');grid