java - 三边测量法 Android Java

Question

我有一个小问题...有谁知道我该如何实现这个？

我需要一个像这样的java方法:

public /*COORDINATE OBJECT*/ getTrilaterationPointBetween(
/*COORDINATE OBJECT*/ coordinate1, double distance1,
/*COORDINATE OBJECT*/ coordinate2, double distance2,
/*COORDINATE OBJECT*/ coordinate3, double distance3){

    //MAGIC

    return /*COORDINATE OBJECT*/;
}

实际上，此方法将接收 3 个坐标以及“myLocation”与其中之一之间的距离...

我需要三边测量算法!!

我做了一些搜索并做了一些测试，结果让我进入了太平洋中部!!

关于OBJECT，很简单，有2个变量代表纬度和经度。

我尝试过这个方法:

static double[] MyTrilateration(Ponto ponto1, double dist1,
                                    Ponto ponto2, double dist2,
                                    Ponto ponto3, double dist3) {

        double[] tmpWAP1 = new double[3];
        double[] tmpWAP2 = new double[3];
        double[] tmpWAP3 = new double[3];

        double tmpLat2, tmpLong2, tmpLat3, tmpLong3;
        double tmpSlide, deg;
        double MyLat, MyLong;

        double[] MyLocation = new double[2];

        tmpLat2  = ponto2.getX() - ponto1.getX();
        tmpLong2 = ponto2.getY() - ponto1.getY();
        tmpLat3  = ponto3.getX() - ponto1.getX();
        tmpLong3 = ponto3.getY() - ponto1.getY();

        tmpSlide = Math.sqrt(Math.pow(tmpLat2,2)+Math.pow(tmpLong2,2));

        deg = (180/Math.PI)*Math.acos( Math.abs(tmpLat2)/Math.abs(tmpSlide));

        if( (tmpLat2>0 && tmpLong2>0) ) {
            deg = 360 - deg;
        }
        else if( (tmpLat2<0 && tmpLong2>0) ) {
            deg = 180 + deg;
        }
        else if( (tmpLat2<0 && tmpLong2<0)){
            deg = 180 - deg;
        }
        else if( (tmpLat2>0 && tmpLong2<0)) {
            deg = deg;
        }

        tmpWAP1[0] = 0.0;
        tmpWAP1[1] = 0.0;
        tmpWAP1[2] = dist1;
        tmpWAP2 = myRotation(tmpLat2, tmpLong2, dist2, deg);
        tmpWAP3 = myRotation(tmpLat3, tmpLong3, dist3, deg);


        MyLat = (Math.pow(tmpWAP1[2],2)-Math.pow(tmpWAP2[2],2)+Math.pow(tmpWAP2[0],2))/(2*tmpWAP2[0]);

        MyLong = (Math.pow(tmpWAP1[2],2)-Math.pow(tmpWAP3[2],2)-Math.pow(MyLat,2)
                +Math.pow(MyLat-tmpWAP3[0],2)+Math.pow(tmpWAP3[1], 2))/(2*tmpWAP3[1]);

        MyLocation = myRotation(MyLat, MyLong, 0, -deg);

        MyLocation[0] = MyLocation[0] + ponto1.getX();
        MyLocation[1] = MyLocation[1] + ponto1.getY();

        return MyLocation;
    }

这个方法:

public String getCoordinateWith(
            Ponto a, Ponto b, Ponto c,
            Float dA,
            Float dB,
            Float dC) {
        Float W, Z, x, y, y2;
        W = dA*dA - dB*dB - a.getX()*a.getX() - a.getY()*a.getY() + b.getX()*b.getX() + b.getY()*b.getY();
        Z = dB*dB - dC*dC - b.getX()*b.getX() - b.getY()*b.getY() + c.getX()*c.getX() + c.getY()*c.getY();

        x = (W*(c.getY()-b.getY()) - Z*(b.getY()-a.getY())) / (2 * ((b.getX()-a.getX())*(c.getY()-b.getY()) - (c.getX()-b.getX())*(b.getY()-a.getY())));
        y = (W - 2*x*(b.getX()-a.getX())) / (2*(b.getY()-a.getY()));
        //y2 is a second measure of y to mitigate errors
        y2 = (Z - 2*x*(c.getX()-b.getX())) / (2*(c.getY()-b.getY()));
        y = (y + y2) / 2;
        return "Position: " + x + " , " + y;
    }

谢谢!

Answer 1

我明白了!这是我的方法，对我有用，但不是100%准确!适用于 z = 0 并且仅适用于 3 个不同的点(2D 坐标)。

“Ponto”类是一个简单的类，具有两个值(X 和 Y/Lat 和 Lon)，即“double”。

public static Ponto getLocationByTrilateration(
                Ponto ponto1, double distance1,
                Ponto ponto2, double distance2,
                Ponto ponto3, double distance3){

        //DECLARACAO DE VARIAVEIS
        Ponto retorno = new Ponto();
        double[] P1   = new double[2];
        double[] P2   = new double[2];
        double[] P3   = new double[2];
        double[] ex   = new double[2];
        double[] ey   = new double[2];
        double[] p3p1 = new double[2];
        double jval  = 0;
        double temp  = 0;
        double ival  = 0;
        double p3p1i = 0;
        double triptx;
        double xval;
        double yval;
        double t1;
        double t2;
        double t3;
        double t;
        double exx;
        double d;
        double eyy;

        //TRANSFORMA OS PONTOS EM VETORES
        //PONTO 1
        P1[0] = ponto1.getX();
        P1[1] = ponto1.getY();
        //PONTO 2
        P2[0] = ponto2.getX();
        P2[1] = ponto2.getY();
        //PONTO 3
        P3[0] = ponto3.getX();
        P3[1] = ponto3.getY();

        //TRANSFORMA O VALOR DE METROS PARA A UNIDADE DO MAPA
        //DISTANCIA ENTRE O PONTO 1 E A MINHA LOCALIZACAO
        distance1 = (distance1 / 100000);
        //DISTANCIA ENTRE O PONTO 2 E A MINHA LOCALIZACAO
        distance2 = (distance2 / 100000);
        //DISTANCIA ENTRE O PONTO 3 E A MINHA LOCALIZACAO
        distance3 = (distance3 / 100000);

        for (int i = 0; i < P1.length; i++) {
            t1   = P2[i];
            t2   = P1[i];
            t    = t1 - t2;
            temp += (t*t);
        }
        d = Math.sqrt(temp);
        for (int i = 0; i < P1.length; i++) {
            t1    = P2[i];
            t2    = P1[i];
            exx   = (t1 - t2)/(Math.sqrt(temp));
            ex[i] = exx;
        }
        for (int i = 0; i < P3.length; i++) {
            t1      = P3[i];
            t2      = P1[i];
            t3      = t1 - t2;
            p3p1[i] = t3;
        }
        for (int i = 0; i < ex.length; i++) {
            t1 = ex[i];
            t2 = p3p1[i];
            ival += (t1*t2);
        }
        for (int  i = 0; i < P3.length; i++) {
            t1 = P3[i];
            t2 = P1[i];
            t3 = ex[i] * ival;
            t  = t1 - t2 -t3;
            p3p1i += (t*t);
        }
        for (int i = 0; i < P3.length; i++) {
            t1 = P3[i];
            t2 = P1[i];
            t3 = ex[i] * ival;
            eyy = (t1 - t2 - t3)/Math.sqrt(p3p1i);
            ey[i] = eyy;
        }
        for (int i = 0; i < ey.length; i++) {
            t1 = ey[i];
            t2 = p3p1[i];
            jval += (t1*t2);
        }
        xval = (Math.pow(distance1, 2) - Math.pow(distance2, 2) + Math.pow(d, 2))/(2*d);
        yval = ((Math.pow(distance1, 2) - Math.pow(distance3, 2) + Math.pow(ival, 2) + Math.pow(jval, 2))/(2*jval)) - ((ival/jval)*xval);

        t1 = ponto1.getX();
        t2 = ex[0] * xval;
        t3 = ey[0] * yval;
        triptx = t1 + t2 + t3;
        retorno.setX(triptx);
        t1 = ponto1.getY();
        t2 = ex[1] * xval;
        t3 = ey[1] * yval;
        triptx = t1 + t2 + t3;
        retorno.setY(triptx);

        return retorno;
    }