首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >在不均匀的列表中查找接近的LON值

在不均匀的列表中查找接近的LON值
EN

Stack Overflow用户
提问于 2019-03-21 21:39:17
回答 1查看 29关注 0票数 0

我有两个经度值矩阵。我想找出两个列表中距离很近的点。我的问题是列表有不同的长度

其中一个名为Lon_pair,另一个名为Lon_epa。

Lon_pair的大小

代码语言:javascript
复制
size(Lon_pair) 
12635       1

Lon_epa的大小

代码语言:javascript
复制
size(Lon_epa)
20560       1

如何在Lon_pair中找到与Lon_epa中的点接近的点?还是找出它们之间的区别?

编辑:添加代码

代码语言:javascript
复制
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
EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2019-03-24 09:04:30

对于非常大的数据集,这不是最有效的,但它可能对您有效(或者至少是一个起点)。下面的双循环遍历两个数据集,找到点之间的距离,保持与最小距离相对应的伴生数据集的最小和行数。

代码语言:javascript
复制
% 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'};
票数 1
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/55281731

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档