поле широты и долготы для C #? - PullRequest
Новая
       51

поле широты и долготы для C #?

6 голосов
/ 17 июля 2010

Я просто хочу найти простой класс C #, который принимает начальную широту и долготу и расстояние и находит ограничивающий прямоугольник (max lat, min lat, max lon, min lon). Есть и другие, похожие вопросы, здесь, на SO, но ни один из них действительно не отвечает на этот вопрос, и те, которые не находятся в C #.

Помощь.

Ответы [ 3 ]

10 голосов
/ 14 января 2013

Вот что вы просите. Престижность Федерико А. Рампони, который написал оригинал на Python здесь . 

public class MapPoint
{
    public double Longitude { get; set; } // In Degrees
    public double Latitude { get; set; } // In Degrees
}

public class BoundingBox
{
    public MapPoint MinPoint { get; set; }
    public MapPoint MaxPoint { get; set; }
}        

// Semi-axes of WGS-84 geoidal reference
private const double WGS84_a = 6378137.0; // Major semiaxis [m]
private const double WGS84_b = 6356752.3; // Minor semiaxis [m]

// 'halfSideInKm' is the half length of the bounding box you want in kilometers.
public static BoundingBox GetBoundingBox(MapPoint point, double halfSideInKm)
{            
    // Bounding box surrounding the point at given coordinates,
    // assuming local approximation of Earth surface as a sphere
    // of radius given by WGS84
    var lat = Deg2rad(point.Latitude);
    var lon = Deg2rad(point.Longitude);
    var halfSide = 1000 * halfSideInKm;

    // Radius of Earth at given latitude
    var radius = WGS84EarthRadius(lat);
    // Radius of the parallel at given latitude
    var pradius = radius * Math.Cos(lat);

    var latMin = lat - halfSide / radius;
    var latMax = lat + halfSide / radius;
    var lonMin = lon - halfSide / pradius;
    var lonMax = lon + halfSide / pradius;

    return new BoundingBox { 
        MinPoint = new MapPoint { Latitude = Rad2deg(latMin), Longitude = Rad2deg(lonMin) },
        MaxPoint = new MapPoint { Latitude = Rad2deg(latMax), Longitude = Rad2deg(lonMax) }
    };            
}

// degrees to radians
private static double Deg2rad(double degrees)
{
    return Math.PI * degrees / 180.0;
}

// radians to degrees
private static double Rad2deg(double radians)
{
    return 180.0 * radians / Math.PI;
}

// Earth radius at a given latitude, according to the WGS-84 ellipsoid [m]
private static double WGS84EarthRadius(double lat)
{
    // http://en.wikipedia.org/wiki/Earth_radius
    var An = WGS84_a * WGS84_a * Math.Cos(lat);
    var Bn = WGS84_b * WGS84_b * Math.Sin(lat);
    var Ad = WGS84_a * Math.Cos(lat);
    var Bd = WGS84_b * Math.Sin(lat);
    return Math.Sqrt((An*An + Bn*Bn) / (Ad*Ad + Bd*Bd));
}
4 голосов
/ 31 мая 2012

Я знаю, что я воскрешаю двухлетний вопрос, но в случае, если кто-нибудь попадет сюда из Google, я написал библиотеку классов C #, которая делает именно это и немного больше. Я создал Github-репозиторий для размещения кода и вскоре напишу в блоге.

Хранилище геолокации

0 голосов
/ 19 июля 2010

Если у вас есть только начальная широта, долгота и высота, тогда ограничивающего прямоугольника нет, поскольку вы определили только одну точку в пространстве. Любой прямоугольник / ограничивающий прямоугольник потребует как минимум две точки для его определения!

...