0

我有两个向量恰好是从形状文件中提取的经纬度。我想对这些值进行插值,因此我的分辨率与我从中找到距离的点相同。

我发现了这个问题另一个几乎解决了我的问题的问题。不同之处在于我不知道我的矢量点之间的间距(它们会有所不同),并且我想确保我的点之间始终最多有 0.0042 度。下面是我的 lat/long 向量的前 10 个点。

latlim=interp1(latlim,1:0.001:numel(latlim));得到了我想要避免任何混叠的东西,但为了提高效率,我想在点之间精确地得到 0.0042 而不仅仅是 1000 个新点。我有很多要点要迭代。

latlim=[78.1530 77.9963 77.6970 77.4092  77.7322 78.0511 78.1530 NaN 78.8044 78.6020];
latlim=interp1(latlim,1:0.001:numel(latlim));

lonlim=[-111.2644 -109.8545 -110.1869 -112.0512 -113.5343 -112.7246 -111.2644 NaN -110.9637 -109.6631];
lonlim=interp1(lonlim,1:0.001:numel(lonlim));

谢谢!

4

1 回答 1

0

我希望我理解你的问题:

如果您想要 0.042 的恒定步长,请尝试:

lonlim = interp1( lonlim, lonlim(1): 0.042 :lonlim(end) );

这将创建一个向量,从 lonlim 的第一个值 lonlim(1) 开始,步长为 0.042,直到 lonlim 的最后一个值,即 lonlim(end)。

希望它有效!:)

于 2013-03-13T18:13:47.973 回答