{"id":204521,"date":"2025-05-29T14:04:36","date_gmt":"2025-05-29T06:04:36","guid":{"rendered":"https:\/\/server.hk\/cnblog\/204521\/"},"modified":"2025-05-29T14:04:36","modified_gmt":"2025-05-29T06:04:36","slug":"%e5%a6%82%e4%bd%95%e4%bc%98%e5%8c%96%e6%89%b9%e9%87%8f%e7%bb%8f%e7%ba%ac%e5%ba%a6%e8%b7%9d%e7%a6%bb%e8%ae%a1%e7%ae%97%ef%bc%9f","status":"publish","type":"post","link":"https:\/\/server.hk\/cnblog\/204521\/","title":{"rendered":"\u5982\u4f55\u4f18\u5316\u6279\u91cf\u7ecf\u7eac\u5ea6\u8ddd\u79bb\u8ba1\u7b97\uff1f"},"content":{"rendered":"<p><b><\/b>     <\/p>\n<h1>\u5982\u4f55\u4f18\u5316\u6279\u91cf\u7ecf\u7eac\u5ea6\u8ddd\u79bb\u8ba1\u7b97\uff1f<\/h1>\n<p>\u4ece\u73b0\u5728\u5f00\u59cb\uff0c\u52aa\u529b\u5b66\u4e60\u5427\uff01\u672c\u6587\u4e3b\u8981\u8bb2\u89e3\u4e86<span style=\"color: #FF6600;, Helvetica, Arial, sans-serif;font-size: 14px;background-color: #FFFFFF\"><\/span>\u7b49\u7b49\u76f8\u5173\u77e5\u8bc6\u70b9\uff0c\u6211\u4f1a\u5728\u4e2d\u6301\u7eed\u66f4\u65b0\u76f8\u5173\u7684\u7cfb\u5217\u6587\u7ae0\uff0c\u6b22\u8fce\u5927\u5bb6\u5173\u6ce8\u5e76\u79ef\u6781\u7559\u8a00\u5efa\u8bae\u3002\u4e0b\u9762\u5c31\u5148\u4e00\u8d77\u6765\u770b\u4e00\u4e0b\u672c\u7bc7\u6b63\u6587\u5185\u5bb9\u5427\uff0c\u5e0c\u671b\u80fd\u5e2e\u5230\u4f60\uff01<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/www.17golang.com\/uploads\/20241109\/1731155971672f5803b3085.jpg\" class=\"aligncenter\"><\/p>\n<p><strong>\u6279\u91cf\u7ecf\u7eac\u5ea6\u8ddd\u79bb\u8ba1\u7b97\u4f18\u5316<\/strong><\/p>\n<p>\u5728\u6279\u91cf\u7ecf\u7eac\u5ea6\u8ddd\u79bb\u8ba1\u7b97\u95ee\u9898\u4e2d\uff0c\u9762\u5bf9\u6d77\u91cf\u6570\u636e\u65f6\uff0c\u5904\u7406\u65f6\u95f4\u5f80\u5f80\u6210\u4e3a\u4e00\u4e2a\u96be\u9898\u3002\u672c\u6587\u5c06\u4ecb\u7ecd\u5982\u4f55\u901a\u8fc7\u6539\u8fdb\u73b0\u6709\u4ee3\u7801\u6765\u7f29\u77ed\u5904\u7406\u65f6\u95f4\u3002<\/p>\n<p>\u4ee3\u7801\u4e2d\uff0c\u8017\u65f6\u8f83\u591a\u7684\u90e8\u5206\u662f\uff1a<\/p>\n<pre>df41 = df4[(df4['longitude']&amp;gt;df4_lon_min) &amp;amp; (df4['longitude']&amp;lt;df4_lon_max) &amp;amp; (df4['latitude']&amp;lt;df4_lat_max) &amp;amp; (df4['latitude']&amp;gt;df4_lat_min)]<\/pre>\n<p>\u8be5\u90e8\u5206\u901a\u8fc7\u7b5b\u9009\u7ecf\u7eac\u5ea6\u8303\u56f4\u6765\u7f29\u5c0f\u8ba1\u7b97\u89c4\u6a21\u3002\u7136\u800c\uff0c\u968f\u7740\u6570\u636e\u91cf\u7684\u589e\u52a0\uff0c\u7b5b\u9009\u8fc7\u7a0b\u8017\u65f6\u4f1a\u53d8\u5f97\u66f4\u957f\u3002<\/p>\n<p>\u4e3a\u4e86\u4f18\u5316\u8fd9\u4e00\u70b9\uff0c\u6211\u4eec\u53ef\u4ee5\u91c7\u7528\u66f4\u6709\u6548\u7684\u65b9\u6cd5\uff0c\u4f8b\u5982\u4f7f\u7528\u7a7a\u95f4\u7d22\u5f15\u6216\u7f51\u683c\u67e5\u8be2\u3002\u8fd9\u4e9b\u65b9\u6cd5\u53ef\u4ee5\u5feb\u901f\u5b9a\u4f4d\u5904\u4e8e\u76ee\u6807\u8303\u56f4\u5185\u7684\u70b9\uff0c\u4ece\u800c\u51cf\u5c11\u7b5b\u9009\u65f6\u95f4\u3002<\/p>\n<p>\u4ee5\u4e0b\u662f\u4e00\u6bb5\u6539\u8fdb\u540e\u7684\u4ee3\u7801\uff1a<\/p>\n<pre>import geopandas as gpd\n\n# \u5c06\u7ecf\u7eac\u5ea6\u6570\u636e\u8f6c\u6362\u4e3a GeoDataFrame\ngdf1 = gpd.GeoDataFrame(df1, geometry=gpd.points_from_xy(df1['Longitude'], df1['Latitude']))\ngdf2 = gpd.GeoDataFrame(df2, geometry=gpd.points_from_xy(df2['Longitude'], df2['Latitude']))\n\n# \u521b\u5efa\u7a7a\u95f4\u7d22\u5f15\ngdf1.sindex\ngdf2.sindex\n\n# \u67e5\u627e\u8ddd\u79bb\u76ee\u6807\u70b9 2km \u4ee5\u5185\u7684\u70b9\nbuffer = gdf1.buffer(2000)\nresult = gpd.sjoin(gdf2, buffer, op='within')<\/pre>\n<p>\u901a\u8fc7\u4f7f\u7528\u7a7a\u95f4\u7d22\u5f15\uff0c\u8ddd\u79bb\u8ba1\u7b97\u7684\u65f6\u95f4\u590d\u6742\u5ea6\u53ef\u4ee5\u4ece o(mn) \u964d\u4f4e\u5230 o(log(mn))\uff0c\u5176\u4e2d m \u548c n \u5206\u522b\u662f a \u8868\u548c b \u8868\u4e2d\u7684\u70b9\u6570\u3002<\/p>\n<p>\u6b64\u5916\uff0c\u4ee5\u4e0b\u5efa\u8bae\u4e5f\u53ef\u8fdb\u4e00\u6b65\u4f18\u5316\u4ee3\u7801\uff1a<\/p>\n<ul>\n<li>\u4f7f\u7528\u591a\u7ebf\u7a0b\u6216\u5206\u5e03\u5f0f\u8ba1\u7b97\u6765\u5e76\u884c\u5904\u7406\u4e0d\u540c\u5206\u7ec4\u3002<\/li>\n<li>\u8c03\u6574\u5206\u7ec4\u5927\u5c0f\u4ee5\u5e73\u8861\u901f\u5ea6\u548c\u5185\u5b58\u4f7f\u7528\u3002<\/li>\n<li>\u4f18\u5316 haversine \u8ddd\u79bb\u8ba1\u7b97\u51fd\u6570\u3002<\/li>\n<\/ul>\n<p>\u4eca\u5929\u5e26\u5927\u5bb6\u4e86\u89e3\u4e86\u7684\u76f8\u5173\u77e5\u8bc6\uff0c\u5e0c\u671b\u5bf9\u4f60\u6709\u6240\u5e2e\u52a9\uff1b\u5173\u4e8e\u6587\u7ae0\u7684\u6280\u672f\u77e5\u8bc6\u6211\u4eec\u4f1a\u4e00\u70b9\u70b9\u6df1\u5165\u4ecb\u7ecd\uff0c\u6b22\u8fce\u5927\u5bb6\u5173\u6ce8\u516c\u4f17\u53f7\uff0c\u4e00\u8d77\u5b66\u4e60\u7f16\u7a0b~<\/p>\n","protected":false},"excerpt":{"rendered":"<p>\u5982\u4f55\u4f18\u5316\u6279\u91cf\u7ecf\u7eac\u5ea6\u8ddd\u79bb\u8ba1\u7b97\uff1f \u4ece&#46;&#46;&#46;<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[4925],"tags":[],"class_list":["post-204521","post","type-post","status-publish","format-standard","hentry","category-4925"],"_links":{"self":[{"href":"https:\/\/server.hk\/cnblog\/wp-json\/wp\/v2\/posts\/204521","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/server.hk\/cnblog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/server.hk\/cnblog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/server.hk\/cnblog\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/server.hk\/cnblog\/wp-json\/wp\/v2\/comments?post=204521"}],"version-history":[{"count":0,"href":"https:\/\/server.hk\/cnblog\/wp-json\/wp\/v2\/posts\/204521\/revisions"}],"wp:attachment":[{"href":"https:\/\/server.hk\/cnblog\/wp-json\/wp\/v2\/media?parent=204521"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/server.hk\/cnblog\/wp-json\/wp\/v2\/categories?post=204521"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/server.hk\/cnblog\/wp-json\/wp\/v2\/tags?post=204521"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}