如何在两个其他点之间找到地理点

时间:2011-10-06 15:41:50

标签: java android maps latitude-longitude

对于我的应用程序,我必须在谷歌地图上找到一个点的位置,只知道它位于其他两个点和捕获坐标的时间(以毫秒为单位)之间。

在我的代码中,假设A和B为给定点,X为要查找的点,I:

  1. 计算A和B之间的距离

  2. 根据时间我发现从A到B的速度(微米/毫秒)

  3. 我找到了从A点到X点的距离(使用时间和速度)

  4. 使用类似的三角形规则,我从点A计算点X的纬度和经度

  5. 此工作流程会在地图上显示错误,因此,X标记通常不在A和B标记之间。

    如何让它更好用?它是地球球形的一个问题吗?

    谢谢大家。

    以下是代码:

        int ax = oldPoint.getLatitude();
        int ay = oldPoint.getLongitude();
    
        int bx = currentPoint.getLatitude();
        int by = currentPoint.getLongitude();
    
        long at = oldPoint.getDataRilevamento(); //get time first point
        long bt = currentPoint.getDataRilevamento(); // get time second point
        long xt = x.getDate(); // time of point to find
    
        int c1 = bx-ax;
        int c2 = by-ay;
        double hyp =  Math.sqrt(Math.pow(c1, 2) + Math.pow(c2, 2));
    
        double vel = hyp / (bt-at);
    
        double pos = vel*(xt - at);
    
        int posx = (int)((pos*c1)/hyp);
        int posy = (int)((pos*c2)/hyp);
    
        x.setLatitude(ax+posx); //set the latitude of X
        x.setLongitude(ay+posy); // set the longitude of X
    

1 个答案:

答案 0 :(得分:14)

您的问题可以通过以下步骤解决。

  • 计算距离A点和B点的距离(使用Haversine formula,这在这里已经足够好了,或者更复杂的Vincenty formula)。要使用 公式正确,Android的microdegrees,这就是getLatitude 和getLongitude返回,必须转换为弧度, 使用这样的公式:

     double radians = Math.toRadians((double)microdegrees/1000000);
    
  • 从A点和B点计算轴承(方向)(使用公式 在同一页上)。这与毕达哥拉斯公式不同,因为 地球是圆的,不是平的。

  • 然后你可以选择一个新的距离并计算给定的X点 A点和上一步中找到的轴承(参见同一页面上的“目标点给定距离和起点承受”或Vincenty's direct formula)。
  • 使用以下公式将生成点的弧度转换为微度:

     int microdegrees = (int)(Math.toDegrees(radians)*1000000);
    

总而言之,我们有以下功能,我将其置于公共领域:

    public static int[] getIntermediatePoint(
        int startLatMicroDeg,
        int startLonMicroDeg,
        int endLatMicroDeg,
        int endLonMicroDeg,
        double t // How much of the distance to use, from 0 through 1
    ){
        // Convert microdegrees to radians
        double alatRad=Math.toRadians((double)startLatMicroDeg/1000000);
        double alonRad=Math.toRadians((double)startLonMicroDeg/1000000);
        double blatRad=Math.toRadians((double)endLatMicroDeg/1000000);
        double blonRad=Math.toRadians((double)endLonMicroDeg/1000000);
        // Calculate distance in longitude
        double dlon=blonRad-alonRad;
        // Calculate common variables
        double alatRadSin=Math.sin(alatRad);
        double blatRadSin=Math.sin(blatRad);
        double alatRadCos=Math.cos(alatRad);
        double blatRadCos=Math.cos(blatRad);
        double dlonCos=Math.cos(dlon);
        // Find distance from A to B
        double distance=Math.acos(alatRadSin*blatRadSin +
                                  alatRadCos*blatRadCos *
                                  dlonCos);
        // Find bearing from A to B
        double bearing=Math.atan2(
            Math.sin(dlon) * blatRadCos,
            alatRadCos*blatRadSin -
            alatRadSin*blatRadCos*dlonCos);
        // Find new point
        double angularDistance=distance*t;
        double angDistSin=Math.sin(angularDistance);
        double angDistCos=Math.cos(angularDistance);
        double xlatRad = Math.asin( alatRadSin*angDistCos +
                                   alatRadCos*angDistSin*Math.cos(bearing) );
        double xlonRad = alonRad + Math.atan2(
            Math.sin(bearing)*angDistSin*alatRadCos,
            angDistCos-alatRadSin*Math.sin(xlatRad));
        // Convert radians to microdegrees
        int xlat=(int)Math.round(Math.toDegrees(xlatRad)*1000000);
        int xlon=(int)Math.round(Math.toDegrees(xlonRad)*1000000);
        if(xlat>90000000)xlat=90000000;
        if(xlat<-90000000)xlat=-90000000;
        while(xlon>180000000)xlon-=360000000;
        while(xlon<=-180000000)xlon+=360000000;
        return new int[]{xlat,xlon};
    }

以下是它的使用方法:

int ax = oldPoint.getLatitude();
int ay = oldPoint.getLongitude();

int bx = currentPoint.getLatitude();
int by = currentPoint.getLongitude();

long at = oldPoint.getDataRilevamento(); //get time first point
long bt = currentPoint.getDataRilevamento(); // get time second point
long xt = x.getDate(); // time of point to find

// Find relative time from point A to point B
double t=(bt==at) ? 0 : ((double)(xt-at))/((double)(bt-at));
// Find new point given the start and end points and the relative time
int[] xpos=getIntermediatePoint(ax,ay,bx,by,t);
x.setLatitude(xpos[0]); //set the latitude of X
x.setLongitude(xpos[1]); // set the longitude of X