我有两个数组代表两个不同的 GPS 路径。每个数组包含偶数索引中的纬度(从 0 开始)和奇数索引中的经度,如下所示:
48.855002219371706,2.3472976684570312,48.855050000000006,2.34735,48.85508,2.3473200000000003,48.85584,2.3477300000000003,48.8562,2.3465000000000003,...
我想计算这两条路径的最近点之间的距离之和。我使用 Haversine 方法计算纬度和经度对之间的距离。
这就是我所做的:我在当前数组中搜索距离最小的纬度和经度对。然后我将最小距离添加到一个汇总所有最小距离的变量中。然后我从数组中删除处理的纬度和经度。
这是代码:
public class Distance {
private String filename1;
private String filename2;
public Distance(String filename1, String filename2) {
this.filename1 = filename1;
this.filename2 = filename2;
}
public double distance() {
if(this.filename1.equals(this.filename2))
return 0;
double distance = 0;
BufferedReader br1 = null;
BufferedReader br2 = null;
try {
br1 = new BufferedReader(new FileReader("files/" + this.filename1));
br2 = new BufferedReader(new FileReader("files/" + this.filename2));
String []data1 = br1.readLine().split(",");
String []data2 = br2.readLine().split(",");
System.out.println("Number of points in " + this.filename1 + " = " + data1.length/2);
System.out.println("Number of points in " + this.filename2 + " = " + data2.length/2);
while((isPresent(data1)) && (isPresent(data2))) {
int posi = -1;
int posj = -1;
double minD = Double.MAX_VALUE;
for(int i = 0; i < data1.length; i+=2) {
if(!data1[i].equals("-9999")) {
double lat1 = Double.parseDouble(data1[i]);
double lon1 = Double.parseDouble(data1[i+1]);
for(int j = 0; j < data2.length; j+=2) {
if(!data2[j].equals("-9999")) {
double lat2 = Double.parseDouble(data2[j]);
double lon2 = Double.parseDouble(data2[j+1]);
double d = getDistance(lat1, lon1, lat2, lon2);
if(minD > d) {
minD = d;
posi = i;
posj = j;
}
}
}
}
}
if(posi != -1){
data1[posi] = data1[posi+1] = data2[posj] = data2[posj+1] = "-9999";
distance += minD;
}
}
} catch (FileNotFoundException e) {
e.printStackTrace();
} catch (IOException e) {
e.printStackTrace();
} finally {
if(br1 != null) {
try {
br1.close();
} catch (IOException e) {
e.printStackTrace();
}
}
if(br2 != null) {
try {
br2.close();
} catch (IOException e) {
e.printStackTrace();
}
}
}
return distance;
}
private boolean isPresent(String arr[]) {
for(int i = 0; i < arr.length; i+=2) {
if(!arr[i].equals("-9999"))
return true;
}
return false;
}
private double getDistance(double lat1, double lon1, double lat2,
double lon2) {
double R = 6378.1; // Radius of the earth in km
double dLat = deg2rad(lat2-lat1); // deg2rad below
double dLon = deg2rad(lon2-lon1);
double a =
Math.sin(dLat/2) * Math.sin(dLat/2) +
Math.cos(deg2rad(lat1)) * Math.cos(deg2rad(lat2)) *
Math.sin(dLon/2) * Math.sin(dLon/2)
;
double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
double d = R * c; // Distance in km
return d*1000;
}
private double deg2rad(double d) {
return d * (Math.PI/180);
}
}
你认为这是一种笨拙的方法吗?我怎样才能快速做到这一点?
注意我已经看到最近的点对问题,但不知道如何使用纬度和经度来实现它,以便比较两个不同数组之间的点。