मैं एक स्थिर LAT, LONG बिंदु के SrcRad त्रिज्या के भीतर सभी बिंदुओं को खोजने के लिए एक साइड इश्यू के रूप में बाउंडिंग बॉक्स समस्या पर काम कर रहा था। काफी कुछ गणनाओं का उपयोग किया गया है
maxLon = $lon + rad2deg($rad/$R/cos(deg2rad($lat)));
minLon = $lon - rad2deg($rad/$R/cos(deg2rad($lat)));
देशांतर सीमा की गणना करने के लिए, लेकिन मुझे यह सभी उत्तरों की आवश्यकता नहीं होने के लिए मिला। क्योंकि आप वास्तव में क्या करना चाहते हैं
(SrcRad/RadEarth)/cos(deg2rad(lat))
मुझे पता है, मुझे पता है कि उत्तर समान होना चाहिए, लेकिन मैंने पाया कि यह नहीं था। ऐसा प्रतीत होता है कि यह सुनिश्चित न करके कि मैं (SRCrad / RadEarth) पहले कर रहा था और फिर कॉस भाग से विभाजित करके मैं कुछ स्थान बिंदुओं को छोड़ रहा था।
अपने सभी बाउंडिंग बॉक्स पॉइंट्स प्राप्त करने के बाद, यदि आपके पास एक ऐसा फंक्शन है जो पॉइंट टू पॉइंट डिस्टेंस को दिए गए लेट को कैलकुलेट करता है, तो केवल उन पॉइंट्स को प्राप्त करना आसान है, जो निश्चित पॉइंट से एक निश्चित दूरी के दायरे हैं। मैंने जो किया था यह रहा। मुझे पता है कि इसने कुछ अतिरिक्त कदम उठाए लेकिन इससे मुझे मदद मिली
-- 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
;
