Sınırlayıcı kutu problemi üzerinde, statik bir LAT, LONG noktanın SrcRad yarıçapı içindeki tüm noktaları bulmanın bir yan sorunu olarak çalışıyordum. Kullanan epeyce hesaplama var
maxLon = $lon + rad2deg($rad/$R/cos(deg2rad($lat)));
minLon = $lon - rad2deg($rad/$R/cos(deg2rad($lat)));
boylam sınırlarını hesaplamak için, ancak bunun gerekli olan tüm cevapları vermediğini buldum. Çünkü gerçekten yapmak istediğin şey
(SrcRad/RadEarth)/cos(deg2rad(lat))
Biliyorum, cevabın aynı olması gerektiğini biliyorum, ama olmadığını anladım. (SRCrad / RadEarth) yaptığımdan emin olmadığım için önce ve sonra Cos kısmına bölerek bazı konum noktalarını dışarıda bırakıyordum.
Tüm sınırlayıcı kutu noktalarınızı aldıktan sonra, enlem verilen Noktadan Noktaya Mesafeyi hesaplayan bir fonksiyonunuz varsa, sadece sabit noktadan belirli bir mesafe yarıçapı olan noktaları elde etmek kolaydır. İşte yaptığım şey. Birkaç adım daha attığını biliyorum ama bana yardımcı oldu
-- GLOBAL Constants
gc_pi CONSTANT REAL := 3.14159265359; -- Pi
-- Conversion Factor Constants
gc_rad_to_degs CONSTANT NUMBER := 180/gc_pi; -- Conversion for Radians to Degrees 180/pi
gc_deg_to_rads CONSTANT NUMBER := gc_pi/180; --Conversion of Degrees to Radians
lv_stat_lat -- The static latitude point that I am searching from
lv_stat_long -- The static longitude point that I am searching from
-- Angular radius ratio in radians
lv_ang_radius := lv_search_radius / lv_earth_radius;
lv_bb_maxlat := lv_stat_lat + (gc_rad_to_deg * lv_ang_radius);
lv_bb_minlat := lv_stat_lat - (gc_rad_to_deg * lv_ang_radius);
--Here's the tricky part, accounting for the Longitude getting smaller as we move up the latitiude scale
-- I seperated the parts of the equation to make it easier to debug and understand
-- I may not be a smart man but I know what the right answer is... :-)
lv_int_calc := gc_deg_to_rads * lv_stat_lat;
lv_int_calc := COS(lv_int_calc);
lv_int_calc := lv_ang_radius/lv_int_calc;
lv_int_calc := gc_rad_to_degs*lv_int_calc;
lv_bb_maxlong := lv_stat_long + lv_int_calc;
lv_bb_minlong := lv_stat_long - lv_int_calc;
-- Now select the values from your location datatable
SELECT * FROM (
SELECT cityaliasname, city, state, zipcode, latitude, longitude,
-- The actual distance in miles
spherecos_pnttopntdist(lv_stat_lat, lv_stat_long, latitude, longitude, 'M') as miles_dist
FROM Location_Table
WHERE latitude between lv_bb_minlat AND lv_bb_maxlat
AND longitude between lv_bb_minlong and lv_bb_maxlong)
WHERE miles_dist <= lv_limit_distance_miles
order by miles_dist
;