Error using interp1

1 次查看(过去 30 天)
David E.S.
David E.S. 2021-2-26
评论: David E.S. 2021-2-26
Hi,
I want to run the following code in order to interpolate missing values of latitudes and longitudes between two dates.
% location of the known latitudes and longitudes
indKnownLATITUDE = ~isnan(arch.latitude);
indKnownLONGITUDE = ~isnan(arch.longitude);
% interpolate to obtain values of unknown values
arch.latitude(~indKnownLATITUDE) = interp1(arch.FECHATXT(indKnownLATITUDE),arch.latitude(indKnownLATITUDE),arch.FECHATXT(~indKnownLATITUDE)) ;
arch.longitude(~indKnownLONGITUDE) = interp1(arch.FECHATXT(indKnownLONGITUDE),arch.longitude(indKnownLONGITUDE),arch.FECHATXT(~indKnownLONGITUDE)) ;
The error is:
Error using griddedInterpolant
The grid vectors must contain unique points.
F = griddedInterpolant(X,V,method);
I know I have to use the unique( ) function, so I've tried to introduce it in the code given above but it still gives me errors. Can you help me?
Thanks!
  2 个评论
Walter Roberson
Walter Roberson 2021-2-26
Is your original data truly scattered with unique latitude and longitude, or is it really grid based?
If you have two locations with the same latitude and different latitude, are you expecting they would have the same value?
David E.S.
David E.S. 2021-2-26
编辑:David E.S. 2021-2-26
This is an example of my data:
Although the same date value is always shown, the decimal places are different. Latitudes and longitudes can be repeated and they're not sorted.
If you have two locations with the same latitude and different latitude, are you expecting they would have the same value?
No. For example a point can be (28.47, -17.21) and another (28.47, -16.14).

请先登录,再进行评论。

采纳的回答

Walter Roberson
Walter Roberson 2021-2-26
Convert the Date + Time into datetimes . Use those as the independent variable and lat as the dependent variable. Remove the entries corresponding to unknown lat/long, then interp1() on what remains at the times the missing records appear.
I was going to suggest that you use spline interpolation, but as the times are a few minutes apart, it is no longer clear to me that smooth interpolation is warranted. It depends what is being tracked. Looks like vehicle location, in which case motion could probably not be assumed to be smooth.

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Interpolation 的更多信息

产品


版本

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by