android - 如何仅通过 GPS 读数确定方向并以 360 度显示？

Question

我正在开发 AS3 AIR for Android 应用程序。“我的问题是:给定一个 GPS 位置，我想检索运动方向(北、西北、西、西南、南、东南、东、东北)”

我想要有关如何使用 2 个不同的 GPS 读数来推断行进方向的代码。然后我想 360 度使用该数据。

我想我了解如何获得大致方向，如北、南、东和西。我会得到 2 个独立的 GPS 读数并比较类似这样的东西......

查看第一个或第二个 GPS 读数中的 X 是否更大。查看第一个或第二个 GPS 读数中的 Y 是否更大。从这里我应该能看到大方向。将 X1 和 X2 的差异与 Y1 和 Y2 的差异进行比较。然后你可以看到用户最常去的方向。然后显示方向。

我不确定如何获取数据并将其用于像指南针一样的 360 度解释...有人可以帮忙吗？

工作代码-------------------------------------------- ----------------------------------

package 

{
import flash.display.Sprite;
// Geolocation sensor stuff
import flash.sensors.Geolocation; 
import flash.events.GeolocationEvent;
//Timer setup stuff for setting frequesncy of GPS Update
import flash.utils.Timer;
import flash.events.TimerEvent;

// Sprite and Textfield display
import flash.display.Sprite;
import flash.text.TextField;
import flash.text.TextFieldAutoSize;
import flash.text.TextFormat;


    public class FlashTest extends Sprite 

    {
        //Variable to check if current or prior Geolocation Event fired
        private var gpsCheck:int = 1;
        public var dLon:Number


        //latitude and longitude in degrees (RAW info from Geolocation)
        public var gpsLat1:Number
        public var gpsLon1:Number
        private var gpsLat2:Number
        private var gpsLon2:Number
        public var bearing:Number


        // Latitude and longitude in radians converted from Degrees
        public var gpsLat1R:Number
        public var gpsLon1R:Number
        private var gpsLat2R:Number
        private var gpsLon2R:Number     
        private var yy:Number  
        private var xx:Number          


        public function FlashTest() 

        {
            // Text box for displaying results
            var my_txt:TextField = new TextField();
                my_txt.wordWrap=true; 
                my_txt.width = 300;
                my_txt.height = 300;
                 addChild(my_txt);

            /*
            //If GPS is on device create Geolocation object named "my_geo"
            //Request updates from my_geo every 2000 milliseconds. Run onGeoUpdate function
            //when Event is triggered. After that create the text box for displaying data.
            if (Geolocation.isSupported)
            {
                var my_geo:Geolocation = new Geolocation();
                my_geo.setRequestedUpdateInterval(2000);                
                my_geo.addEventListener(GeolocationEvent.UPDATE, onGeoUpdate);
                var my_txt:TextField = new TextField();
                my_txt.wordWrap=true; 
                my_txt.width = 300;
                my_txt.height = 300;
                addChild(my_txt);
            } 
            // If GPS is not supported on device display "No GPS Signal" 
            else
            {
                addChild(my_txt);
                my_txt.wordWrap=true; 
                my_txt.width = 300;
                my_txt.height = 300;
                addChild(my_txt);
                my_txt.text = "No GPS Signal ";
            }
            */


            // False GPS reading being passed for testing 
            //COMMENT THESE LINES OUT STARTING HERE---
                gpsLat1 =  42.1234584;
                gpsLon1 = -83.1234577;
                gpsLat2 =  42.1234583;
                gpsLon2 = -83.1234577;
           // END OF WHAT SHOULD BE COMMENTED OUT---

           // Equations to convert all RAW Geolocation information over to radians 
                gpsLat1R = gpsLat1 * Math.PI / 180;
                gpsLon1R = gpsLon1 * Math.PI / 180;
                gpsLat2R = gpsLat2 * Math.PI / 180;
                gpsLon2R = gpsLon2 * Math.PI / 180;
          // The rest of the math     
                dLon = gpsLon1 - gpsLon2;
                yy = Math.sin(dLon) * Math.cos(gpsLat2R);
                xx = Math.cos(gpsLat1R) * Math.sin(gpsLat2R) - Math.sin(gpsLat1R) * Math.cos(gpsLat2R) * Math.cos(dLon);
                bearing = Math.atan2(yy, xx) * 180 / Math.PI;
         // Run The Geoupdate function                
            onGeoUpdate();


        // onGeoUpdate basically displays the information that was collected and converted.
        // This is where you will put what you want the code to do with the results
            function onGeoUpdate():void
            {
            my_txt.text = "My Latitude is "+gpsLat1+ " and my Longitude is "+gpsLon1+ 
            "My 2nd Latitude is "+gpsLat2+" and my 2nd Longitude is "+gpsLon2+ 
            " Bearing is " +bearing;

            }            
        }
    }
}

Answer 1

你需要的公式(对于lon1，lon2经度，和lat1，lat2纬度的两个points) 在许多地方给出 - 例如，在 http://www.movable-type.co.uk/scripts/latlong.html .数学是这样的(将其转换为您最喜欢的语言...):

dLon = lon2 - lon1;
y = sin(dLon) * cos(lat2);
x = cos(lat1) * sin(lat2) - sin(lat1) * cos(lat2) *cos(dLon);
bearing = atan2(y, x) * 180 / Pi;

注意:atan2 将结果缩减为 [-pi, pi] 之间的值。乘以 180/pi 后，数字将介于 [-180,180] 之间。为了获得 0 到 360 之间的值，您可以这样做，例如:

if (bearing < 0) bearing += 360;

我希望您可以填写您需要的任何其他详细信息，以便为您完成这项工作。

编辑 我已经使用您提供的数字编写了一些代码:这是 Matlab，但它应该非常易读:

% bearings
lat1 = 42.1234584;
lat2 = 42.1234583;
lon1 = -83.1234577;
lon2 = -83.1234510;

% convert to radians:
g2r = pi/180;
lat1r = lat1 * g2r;
lat2r = lat2 * g2r;
lon1r = lon1 * g2r;
lon2r = lon2 * g2r;
dlonr = lon2r - lon1r;
y = sin(dlonr) * cos(lat2r);
x = cos(lat1r)*sin(lat2r) - sin(lat1r) * cos(lat2r)*cos(dlonr);

% compute bearning and convert back to degrees:
bearing = atan2(y, x) / g2r;

fprintf(1,'x: %+.3e\ny: %+.3e\nbearing: %.3f\n', x, y, bearing);

这导致输出:

x: -1.745e-09
y: +8.673e-08
bearing: 91.153

如您所见，x 和 y 很小 - 它们代表您移动的“地球周长的一部分”(大约正东 3.5 米，这似乎...)。您应该能够使用这些数字调试您的实现。

请注意，GPS 的“绝对”精度可能较差(在您的情况下，不确定性 > 100 m)，但“相对”精度仍然不错(它测量两个位置之间的差异远优于 100米)。