Merkez noktadan kare x mil koordinatlarını mı hesaplıyorsunuz?

Meydanın tüm taraflarının merkezden x mil olacağı merkezi bir noktanın etrafında balta mil karesi (veya daire) oluşturmaya çalışıyorum. 4 köşe koordinatına ihtiyacım var.

Kafamı döndürmeye çalışan beynimi karıştırıyor mu? Haversin formülünü kullanarak iki nokta arasındaki mesafeyi çalışabilirim ama matematik ciddi bir şekilde güçlü noktam değil ve günahı, çünkü vs'yi anlamıyorum.

Enlem / Boylam X milini noktadan hesaplamakla karşılaştım ? ama anlamıyorum!

Bunu elma ve armut açısından nasıl açıklayacağımı açıklayacak kadar nazik olacak mı?

Tam olarak ne yapmaya çalıştığımı açıklamak için;

Kullanıcıların belirli bir alandaki binaları arayabileceği bir web sitem var. Onlar bir kasaba veya yer (ki lat uzun bilecek) girecek ve belirli bir yarıçap içinde söylemek demek 10 mil yer.

Min / max lat ve 10 mil yarıçapı uzun bulmak gerekir, bu yüzden benzer bir nerede yan tümcesi kullanarak veritabanımı sorgulayabilirsiniz:

Where buildingLat <= maxLat 
  and buildingLat <= minLat 
  and buildingLong >= minLong 
   or buildingLong >= maxLong

Bir çeşit formüle ihtiyacım var!

Koordinatlarım ondalık derecelerde

Bu amaçla basit yaklaşımlar yeterince iyidir. Kuzey veya güney, bir derece yaklaşık 69 mil, doğu veya batı, sadece 69 * cos (enlem) mil. Enlemler on mil açıklıkta fazla değişmediği için, "kare" nin enleminin kosinüsünü güvenle kullanabilirsiniz. Bu nedenle, merkezi konumdan (f, l) r mil uzaklıktaki kare köşeler için istenen koordinatlar , lat-lon olarak verilir.

df = r/69        // North-south distance in degrees
dl = df / cos(f) // East-west distance in degrees
{(f-df,l-dl), (f+df,l-dl), (f+df,l+dl), (f-df,l+dl)} // List of vertices

Örneğin, r = 10 mil olduğunu ve merkezi konumun enlem 50 derece kuzey, boylam 1 derece batı olduğunu, böylece (f, l) = (50, -1) derece olduğunu varsayalım. Sonra

df = 10/69 = 0.145
dl = 0.145 / cos(50 degrees) = 0.145 / 0.6428 = 0.225
f - df = 50 - 0.145 = 49.855 (southernmost latitude)
f + df = 50 + 0.145 = 50.145 (northernmost latitude)
l - dl = -1 - 0.225 = -1.225 (western longitude)
l + dl = -1 + 0.225 = -0.775 (eastern longitude)

ve koordinatları güneybatı köşesinden başlayarak meydanın etrafında saat yönünde yürürken (49.855, -1.225), (50.145, -1.225), (50.145, -0.775) ve (49.855, -0.775) şeklindedir.

Bu yaklaşımı kutupların yakınında veya bir taraftaki birkaç dereceden daha büyük kareler için kullanmayın. Ayrıca, CBS'nin sınırlamalarına bağlı olarak, genellikle + -180 derecede alınan, boylamda küresel kesim çevresinde biraz özen gösterilmesi gerekebilir.

Merkezin X koordinatını alın ve bundan x mil çıkarın, bu karenizin sol tarafıdır. Sonra merkezin Y koordinatını alın ve X milini buradan çıkarın, bu karenizin alt kısmıdır. Bu adımları tekrarlayın, ancak sağ ve üst kenarları elde etmek için çıkartmak yerine ekleyin. Artık karenizin dört köşesini inşa edebilirsiniz.

Yukarıdakilerin, merkez noktanızın mil cinsinden olduğunu varsayalım. İlk önce değilse yeniden projekte edin. Diğer bilge tüm bahisler kapalı ve kare kare olmayacak.

Sonunda cevabım: (c #)

Muhtemelen 4 koordinatlara ihtiyacım yok ama bence oldukça doğrular.

 public static void GetBoundingCoords(double centerLat, double centerLong,  double distance)
    {
     Coordinate top=   MaxLatLongOnBearing(centerLat, centerLong,45,10);
     Coordinate right = MaxLatLongOnBearing(centerLat, centerLong, 135, 10);
     Coordinate bottom = MaxLatLongOnBearing(centerLat, centerLong, 225, 10);
     Coordinate left = MaxLatLongOnBearing(centerLat, centerLong, 315, 10);
    }

    public static Coordinate MaxLatLongOnBearing(double centerLat, double centerLong, double bearing, double distance)
    {

        var lonRads = ToRadian(centerLong);
        var latRads = ToRadian(centerLat);
        var bearingRads = ToRadian(bearing);
        var maxLatRads = Math.Asin(Math.Sin(latRads) * Math.Cos(distance / 6371) + Math.Cos(latRads) * Math.Sin(distance / 6371) * Math.Cos(bearingRads));
        var maxLonRads = lonRads + Math.Atan2((Math.Sin(bearingRads) * Math.Sin(distance / 6371) * Math.Cos(latRads)), (Math.Cos(distance / 6371) - Math.Sin(latRads) * Math.Sin(maxLatRads)));

        var maxLat = RadiansToDegrees(maxLatRads);
        var maxLong = RadiansToDegrees(maxLonRads);

        return new Coordinate(){Latitude=maxLat, Longitude=maxLong};
    }

DÜZENLE

Kare x milimin köşelerini merkez noktadan belirlediğimde fark ettim, karemin kenarları aynı x mil olmayacak. (matematik benim güçlü nokta değil dedi) Yani kareler kenarları x mil istiyorum köşe noktalarını mesafe elde etmek için Pythagoras Teorem çapraz mesafe çalışmak için kullandım. (dik açılı üçgende, hipotenüs üzerindeki kare (diyagonal) diğer iki tarafın karesine eşittir)

Uzamsal olarak farkında bir veritabanı kullanıyorsanız, ilgi alanınızı verilerinizin depolandığı aynı koordinat sistemine dönüştürebilir ve daha sonra bir elma ile elma karşılaştırmasını yapabilirsiniz.

Örneğin:

1. Kullanıcı bir konum seçerek lat / lon ile sonuçlanır.
2. Uzamsal veri tabanından bu noktayı alana uygun olarak öngörülen bir koordinat sistemine dönüştürmesini isteyin (fit veya metre birimi vb.).
3. İlgi alanınızı öngörülen nokta etrafında oluşturun.
4. Uzamsal veritabanından bu ilgi alanını tekrar lat / lon'a dönüştürmesini isteyin.
5. Yapmanız gereken her türlü karşılaştırmayı yapın.

Ben içinde ne kullanılmış bu sayfada

Başlangıç ​​noktasından uzaklık ve yön verilen hedef noktası

Formül:
lat2 = asin (sin (lat1) * cos (d / R) + cos (lat1) * sin (d / R) * cos (θ))
lon2 = lon1 + atan2 (sin (θ) * sin (d / R) * cos (lat1), cos (d / R) −sin (lat1) * sin (lat2))

θ yataktır (radyan cinsinden, kuzeyden saat yönünde); d / R açısal mesafedir (radyan cinsinden), burada d gidilen mesafedir ve R dünyanın yarıçapıdır

Θ "Sol üst nokta" için -45 derece (radyan cinsinden) ve "sağ alt" için 135 derece kullandım

(Geçenlerde aynı soruyu matematik sitesinde sordum )