我有两个经度值矩阵。我想找出两个列表中距离很近的点。我的问题是列表有不同的长度
其中一个名为Lon_pair,另一个名为Lon_epa。
Lon_pair的大小
size(Lon_pair)
12635 1Lon_epa的大小
size(Lon_epa)
20560 1如何在Lon_pair中找到与Lon_epa中的点接近的点?还是找出它们之间的区别?
编辑:添加代码
clear all
close all
filename_pair = 'C:\Users\tadams15\Desktop\ALL_PAIR.csv';
data = xlsread(filename_pair);
Lon_pair = data(:,2);
Lat_pair = data(:,3);
Label_pair = data(:,4);
filename_epa = 'C:\Users\tadams15\Desktop\aqs.csv'
data2 = xlsread(filename_epa);
Lon_epa = data2(:,4);
Lat_epa = data2(:,5);
%% Find Lat/Lon Combos
for j = 1:20560
for i = 2:12636
d_lon(j) = Lon_epa(:) - Lon_pair(i);
d_lon_m(j) = d_lon/(9e-06);
if d_lon(:) < 300
[row,col] = find(d_lon)
end
end
end发布于 2019-03-24 09:04:30
对于非常大的数据集,这不是最有效的,但它可能对您有效(或者至少是一个起点)。下面的双循环遍历两个数据集,找到点之间的距离,保持与最小距离相对应的伴生数据集的最小和行数。
% set up dummy data to mimic UTM x,y (northing,easting) coordinates
x1 = randi([0,50000],20,1);
y1 = randi([-230000,420000],20,1);
x2 = randi([0,50000],12,1);
y2 = randi([-230000,420000],12,1);
% end set up
mydist = NaN(length(x1), 2); % pre-allocate
for ii = 1:length(mydist)
xa = x1(ii);
ya = y1(ii);
temp_dist = NaN(length(x2), 2);
for jj = 1:length(x2)
xb = x2(jj);
yb = y2(jj);
temp_dist(jj, 1) = sqrt((xa-xb)^2 + (ya-yb)^2); % if you have the Statistics Toolbox, you can just use pdist
temp_dist(jj, 2) = jj; % this is the row number in set2
[value, index] = min(temp_dist(:,1));
mydist(ii, 1) = value;
mydist(ii, 2) = temp_dist(index, 2);
end
end
mydist_table = array2table(mydist);
mydist_table.Properties.VariableNames = {'Clostest_Distance', 'Set2_RowNumber'};https://stackoverflow.com/questions/55281731
复制相似问题