Tính toán kinh độ mới, vĩ độ từ + n mét cũ

Tôi muốn tạo 2 kinh độ mới và 2 vĩ độ mới dựa trên tọa độ và khoảng cách tính bằng mét, tôi muốn tạo một hộp giới hạn đẹp xung quanh một điểm nhất định. Nó dành cho một phần của thành phố và tối đa là ± 1500 mét. Do đó, tôi không nghĩ rằng phải tính đến độ cong của trái đất.

Vì vậy, tôi có 50.0452345(x) và 4.3242234(y) và tôi muốn biết x + 500 mét, x - 500 mét, y - 500 mét, y + 500 mét

Tôi đã tìm thấy nhiều thuật toán nhưng hầu như tất cả đều giải quyết được khoảng cách giữa các điểm.

Số km trên mỗi độ kinh là khoảng

(2*pi/360) * r_earth * cos(theta)

nơi thetalà vĩ độ và r_earthkhoảng 6378 km.

Số km trên mỗi độ vĩ độ là gần như nhau ở tất cả các vị trí, khoảng

(2*pi/360) * r_earth = 111 km / degree 

Vì vậy, bạn có thể làm:

new_latitude  = latitude  + (dy / r_earth) * (180 / pi);
new_longitude = longitude + (dx / r_earth) * (180 / pi) / cos(latitude * pi/180);

Miễn là dxvà dynhỏ so với bán kính của trái đất và bạn không đến quá gần các cực.

Câu trả lời được chấp nhận là hoàn toàn đúng và hoạt động. Tôi đã thực hiện một số chỉnh sửa và biến thành điều này:

double meters = 50;

// number of km per degree = ~111km (111.32 in google maps, but range varies
   between 110.567km at the equator and 111.699km at the poles)
// 1km in degree = 1 / 111.32km = 0.0089
// 1m in degree = 0.0089 / 1000 = 0.0000089
double coef = meters * 0.0000089;

double new_lat = my_lat + coef;

// pi / 180 = 0.018
double new_long = my_long + coef / Math.cos(my_lat * 0.018);

Hy vọng điều này cũng sẽ giúp.

Đối với vĩ độ làm:

var earth = 6378.137,  //radius of the earth in kilometer
    pi = Math.PI,
    m = (1 / ((2 * pi / 360) * earth)) / 1000;  //1 meter in degree

var new_latitude = latitude + (your_meters * m);

Đối với kinh độ làm:

var earth = 6378.137,  //radius of the earth in kilometer
    pi = Math.PI,
    cos = Math.cos,
    m = (1 / ((2 * pi / 360) * earth)) / 1000;  //1 meter in degree

var new_longitude = longitude + (your_meters * m) / cos(latitude * (pi / 180));

Biến your_meterscó thể chứa một giá trị dương hoặc một giá trị âm.

Bạn đã kiểm tra chưa: Làm cách nào để tìm vĩ độ / kinh độ là x km về phía bắc của một vĩ độ / kinh độ nhất định ?

Những tính toán này rất khó chịu, tôi đã thực hiện rất nhiều. Công thức hasrsine sẽ là người bạn của bạn.

Một số tài liệu tham khảo: http://www.movable-type.co.uk/scripts/latlong.html

Tôi đã phải dành khoảng hai giờ để tìm ra giải pháp bằng @nibot, tôi chỉ cần một phương pháp để tạo một hộp ranh giới với điểm trung tâm và chiều rộng / chiều cao (hoặc bán kính) tính bằng km:

Tôi không hiểu giải pháp về mặt toán học / địa lý. Tôi đã điều chỉnh giải pháp (bằng cách thử và sai) để có được bốn tọa độ:

Bắc:

private static Position FromKmToNPosition(Position p, double km)
{
    double r_earth = 6378;
    var pi = Math.PI;
    var new_latitude = p.Lat + (km / r_earth) * (180 / pi);
    return new Position(new_latitude, p.Long);
}

Phía đông:

private static Position FromKmToEPosition(Position p, double km)
{
    double r_earth = 6378;
    var pi = Math.PI;
    var new_longitude = p.Long + (km / r_earth) * (180 / pi) / Math.Cos(p.Lat * pi / 180);
    return new Position(p.Lat, new_longitude);
}

Miền Nam:

private static Position FromKmToSPosition(Position p, double km)
{
    double r_earth = 6378;
    var pi = Math.PI;
    var new_latitude = p.Lat - (km / r_earth) * (180 / pi);
    return new Position(new_latitude, p.Long);
}

Hướng Tây:

private static Position FromKmToWPosition(Position p, double km)
{
    double r_earth = 6378;
    var pi = Math.PI;
    var new_longitude = p.Long - (km / r_earth) * (180 / pi) / Math.Cos(p.Lat * pi / 180);
    return new Position(p.Lat, new_longitude);
}

nếu bạn không cần phải chính xác lắm thì: mỗi 10000 mét là khoảng 0,1 cho vĩ độ và kinh độ. ví dụ: tôi muốn tải các vị trí 3000 mét xung quanh point_A từ cơ sở dữ liệu của mình:

double newMeter =  3000 * 0.1 / 10000;
double lat1 = point_A.latitude - newMeter;
double lat2 = point_A.latitude + newMeter;
double lon1 = point_A.longitude - newMeter;
double lon1 = point_A.longitude + newMeter;
Cursor c = mDb.rawQuery("select * from TABLE1  where lat >= " + lat1 + " and lat <= " + lat2 + " and lon >= " + lon1 + " and lon <= " + lon2 + " order by id", null);

public double MeterToDegree(double meters, double latitude)
{
    return meters / (111.32 * 1000 * Math.Cos(latitude * (Math.PI / 180)));
}

var meters = 50;
var coef = meters * 0.0000089;
var new_lat = map.getCenter().lat.apply() + coef;
var new_long = map.getCenter().lng.apply() + coef / Math.cos(new_lat * 0.018);
map.setCenter({lat:new_lat, lng:new_long});