指南：將任意點(diǎn)吸附到線上的最近點(diǎn)計算與球面幾何實現(xiàn)原理)
數(shù)據(jù)分析【免費(fèi)下載鏈接】turfA modular geospatial engine written in JavaScript and TypeScript項目地址https://gitcode.com/gh_mirrors/tu/turf點(diǎn)擊查看免費(fèi)下載turf/nearest-point-on-line是 Turf 模塊化地理引擎中專門用于點(diǎn)吸附snap to line的核心工具給定一條LineString或MultiLineString與一個參考點(diǎn)它返回線上距離該點(diǎn)最近的那個位置點(diǎn)并附帶完整的距離與索引信息。本指南將帶你完整掌握它的 API、返回值語義含 v7.4 新增屬性與舊屬性對照、球面大圓幾何的底層實現(xiàn)以及如何通過倉庫源碼與測試用例驗證其行為邊界——讀完后你可以在路徑糾偏、地圖吸附、導(dǎo)航匹配等場景中直接落地使用。功能概述與典型應(yīng)用場景nearestPointOnLine解決的問題非常明確求一條線上離某個給定點(diǎn)最近的點(diǎn)。它接收任意形式的線段數(shù)據(jù)LineString/MultiLineString的 Geometry 或 Feature和任意形式的點(diǎn)數(shù)據(jù)Point的 Geometry、Feature甚至直接是[lng, lat]坐標(biāo)數(shù)組返回一個帶豐富屬性的PointFeature。典型的實戰(zhàn)場景包括地圖交互吸附用戶在地圖上點(diǎn)擊時把點(diǎn)擊位置吸附到道路、河流等線性要素上實現(xiàn)自動貼線。路徑糾偏 / 軌跡匹配GPS 采樣點(diǎn)偏移到實際道路之外時將軌跡點(diǎn)拉回到路網(wǎng)上。路線輔助計算結(jié)合返回的totalDistance、segmentIndex等屬性定位點(diǎn)在線路上的里程位置進(jìn)而配合turf/along等模塊做分段分析。該模塊屬于 Turf 的獨(dú)立子包版本為 7.4.0被 turf 主包統(tǒng)一導(dǎo)出同時依賴turf/distance測距、turf/helpers要素構(gòu)造、turf/invariant坐標(biāo)校驗與turf/meta要素遍歷四個基礎(chǔ)模塊見 package.json。安裝與引入單獨(dú)安裝單個模塊$ npm install turf/nearest-point-on-line該包聲明為 ESM 模塊type: module運(yùn)行環(huán)境要求 Node.js 22見 package.json。安裝全量 turf/turf如果你希望使用 Turf 全家桶一次性引入所有功能函數(shù)$ npm install turf/turf之后即可通過turf.nearestPointOnLine(...)調(diào)用全量包的導(dǎo)出定義在 packages/turf/index.ts。在代碼中引入import { nearestPointOnLine } from turf/nearest-point-on-line; // 或從全量包引入 // import { nearestPointOnLine } from turf/turf;API 簽名與參數(shù)詳解函數(shù)簽名function nearestPointOnLineG extends LineString | MultiLineString( lines: FeatureG | G, // 要吸附的線 inputPoint: Coord, // 參考點(diǎn)Geometry / FeaturePoint / number[] options: { units?: Units } {} // 可選參數(shù) ): FeaturePoint, {...屬性...}完整的類型定義見 index.ts。參數(shù)說明參數(shù)類型說明linesGeometry \| FeatureLineString \| MultiLineString要吸附到的線要素可以是LineString或MultiLineString的 Geometry 或 FeatureinputPointGeometry \| FeaturePoint \| number[]參考點(diǎn)支持三種形態(tài)Point的 Geometry、PointFeature或直接傳[lng, lat]坐標(biāo)數(shù)組options.unitsUnits距離單位默認(rèn)kilometersoptions.units 支持的取值options.units支持 Turf 的全部合法單位定義于 packages/turf-helpers/index.tsmeters | metres | m | millimeters | millimetres | mm | centimeters | centimetres | cm | kilometers | kilometres | km | miles | mi | nauticalmiles | nmi | inches | in | yards | yd | feet | ft | radians | rad | degrees | deg注意radians/degrees等角度單位在距離語義上具有緯度依賴性——經(jīng)度方向上每一度所代表的實際距離隨緯度變化在地球上不同位置度量差異很大請僅在明確理解其含義時使用。返回值屬性v7.4 新屬性與舊屬性對照函數(shù)返回線上離inputPoint最近的點(diǎn)FeaturePoint。自 v7.4 起Turf 開始遷移返回值屬性命名新的推薦屬性與舊屬性一一對應(yīng)如下新屬性v7.4 起推薦舊屬性仍兼容已標(biāo)記 deprecated含義lineStringIndexmultiFeatureIndex最近點(diǎn)落在輸入MultiLineString的第幾條 LineString 上輸入為 LineString 時恒為 0segmentIndexindex最近點(diǎn)落在該 LineString 的第幾個線段上totalDistancelocation從整條 MultiLineString 絕對起點(diǎn)算起的沿線距離lineDistance—舊版無對應(yīng)從最近點(diǎn)所在 LineString 起點(diǎn)算起的沿線距離segmentDistance—舊版無對應(yīng)從最近點(diǎn)所在線段起點(diǎn)算起的沿線距離pointDistancedist參考點(diǎn)與最近點(diǎn)之間的直線距離官方保證multiFeatureIndex、index、location、dist至少在下一個大版本發(fā)布之前繼續(xù)可用見 README.md 與 index.ts 的 JSDoc 說明。源碼在生成新屬性的同時會同步將multiFeatureIndex/index/location/dist指向?qū)?yīng)的新屬性值見 index.ts因此新舊寫法可以并存。類型層面新屬性與舊屬性都被聲明在返回類型的 properties 中其中舊屬性帶有deprecated標(biāo)記見 index.ts。建議新代碼直接使用新屬性以便在未來的大版本升級中平滑遷移?；A(chǔ)用法示例官方示例來源 README.md演示一條折線與一個參考點(diǎn)以英里為單位求吸附點(diǎn)var line turf.lineString([ [-77.031669, 38.878605], [-77.029609, 38.881946], [-77.020339, 38.884084], [-77.025661, 38.885821], [-77.021884, 38.889563], [-77.019824, 38.892368] ]); var inputPoint turf.point([-77.037076, 38.884017]); var snapped turf.nearestPointOnLine(line, inputPoint, {units: miles}); // addToMap var addToMap [line, inputPoint, snapped]; snapped.properties[marker-color] #00f;上面的例子中snapped即為吸附后的最近點(diǎn)其geometry.coordinates是線上的經(jīng)緯度坐標(biāo)properties中攜帶上文表格里的全部索引與距離字段。需要展示時可將原始線、參考點(diǎn)和吸附點(diǎn)一并加入地圖圖層并給吸附點(diǎn)設(shè)置marker-color以突出顯示。更完整的可運(yùn)行示例import { nearestPointOnLine } from turf/nearest-point-on-line; import { point, lineString, multiLineString } from turf/helpers; const pt point([1.5, 1.5]); const line lineString([ [0, 0], [1, 1], ]); const multiLine multiLineString([ [ [0, 0], [1, 1], [2, 2], [0, 0], ], ]); // 各種合法的參數(shù)組合 nearestPointOnLine(line, pt); // 默認(rèn)單位 kilometers nearestPointOnLine(multiLine, pt); // 支持 MultiLineString nearestPointOnLine(line.geometry, pt); // 直接傳 Geometry nearestPointOnLine(multiLine.geometry, pt); // 直接傳 MultiLineString Geometry nearestPointOnLine(line, pt, { units: miles }); // 指定單位 nearestPointOnLine(line, [1.5, 1.5]); // 傳坐標(biāo)數(shù)組 // 輸出可以作為下一次輸入二次吸附 const output nearestPointOnLine(line, pt); nearestPointOnLine(line, output); // 新舊屬性均可讀取 output.properties.dist; // 舊屬性 output.properties.multiFeatureIndex; // 舊屬性 output.properties.index; // 舊屬性 output.properties.location; // 舊屬性 output.properties.pointDistance; // 新屬性推薦 output.properties.totalDistance; // 新屬性推薦上述參數(shù)組合與屬性訪問方式均由 types.ts 中的類型測試覆蓋保證編譯期可驗證。源碼實現(xiàn)原理逐段求最近點(diǎn)與球面幾何深入 index.ts可以看清整個算法的執(zhí)行流程1. 參數(shù)校驗與逐要素遍歷函數(shù)首先校驗lines與inputPoint是否都提供了否則拋出lines and inputPoint are required arguments見 index.ts。隨后通過getCoord提取參考點(diǎn)坐標(biāo)并用turf/meta的flattenEach將輸入包括 MultiLineString扁平化為逐條 LineString 遍歷同時記錄lineStringIndex見 index.ts。2. 逐線段計算并維護(hù)累計距離對每條 LineString 的每一段coords[i]到coords[i 1]算法先計算該段長度segmentLength distance(start, stop, options)再求參考點(diǎn)在該段上的垂足最近點(diǎn)短路優(yōu)化如果參考點(diǎn)恰好等于某線段的端點(diǎn)則直接采用該端點(diǎn)不再做投影計算見 index.ts。這里特意先比較終點(diǎn)再比較起點(diǎn)以保證起點(diǎn)與終點(diǎn)重合時行為一致。計算候選點(diǎn)與參考點(diǎn)的距離pointDistance若小于當(dāng)前最優(yōu)解則更新closestPt并同步累計出segmentDistance、lineDistance、totalDistance。每處理完一段就把totalDistance與lineDistance累加該段長度當(dāng)flattenEach切換 LineString 時lineDistance重置為 0見 index.ts。這樣totalDistance恒為從整條 MultiLineString 絕對起點(diǎn)到最近點(diǎn)的沿線里程lineDistance為該最近點(diǎn)所在 LineString 內(nèi)的里程segmentDistance為該最近點(diǎn)所在線段內(nèi)的里程。3. 核心幾何nearestPointOnSegment 的大圓投影真正的投影邏輯在nearestPointOnSegment(posA, posB, posC)中見 index.ts其核心思想是把經(jīng)緯度坐標(biāo)轉(zhuǎn)換到三維單位球面的直角坐標(biāo)系中再通過球面大圓幾何求最近點(diǎn)坐標(biāo)球面化lngLatToVector將[lng, lat]轉(zhuǎn)為單位向量(x, y, z)見 index.ts。構(gòu)造大圓平面法向量線段兩端點(diǎn)向量A、B的叉積cross(A, B)即該線段所在大圓平面的法向量segmentAxis。求過參考點(diǎn)的大圓交點(diǎn)依次計算targetAxis cross(segmentAxis, C)、intersectionAxis cross(targetAxis, segmentAxis)歸一化后得到大圓的兩個交點(diǎn)I1/I2用點(diǎn)積余弦相似度挑選離參考點(diǎn)C更近的那個作為候選交點(diǎn)I。判斷交點(diǎn)是否落在線段弧上通過兩次叉積點(diǎn)積比較cross(A, I)與cross(I, B)相對segmentAxis的方向一致性判斷旋轉(zhuǎn)方向若候選交點(diǎn)確實位于 A→B 的弧段內(nèi)則直接返回否則退化為選擇距離參考點(diǎn)更近的端點(diǎn) A 或 B見 index.ts。整個算法對每條線段只做常數(shù)次向量運(yùn)算配合前面先比端點(diǎn)、再算投影的短路邏輯復(fù)雜度為 O(線段數(shù))。4. 退化情況與反足點(diǎn)處理源碼對兩類數(shù)學(xué)退化情況做了顯式處理線段兩端向量平行aligned或反足antipodal此時叉積結(jié)果為[0, 0, 0]無法構(gòu)造大圓平面。平行dot(A, B) 0意味著 A、B 重合算法直接返回終點(diǎn)保持一致行為反足意味著 A、B 之間存在無數(shù)個大圓弧未定義此時選擇返回參考點(diǎn)本身語義上等同于該點(diǎn)位于線上見 index.ts。這正是 README 開頭所說明的若輸入 LineString 中存在反足線段弧未定義函數(shù)將返回點(diǎn)位于線上。候選交點(diǎn)與參考點(diǎn)重合/反足targetAxis為 0 向量參考點(diǎn)到整個線段等距為保持一致返回終點(diǎn)見 index.ts。此外vectorToLngLat在反解經(jīng)緯度時對 z 分量做了[-1, 1]鉗制避免浮點(diǎn)誤差導(dǎo)致asin越界見 index.ts。邊界行為與測試驗證模塊的邊界行為在 test.ts 中有大量回歸測試覆蓋可直接作為行為文檔閱讀參考點(diǎn)恰為線起點(diǎn)/終點(diǎn)點(diǎn)不移動totalDistance分別為 0 或全線長。參考點(diǎn)在第一個頂點(diǎn)之前吸附到第一個頂點(diǎn)totalDistance為 0見points behind first point測試。參考點(diǎn)在最后一個頂點(diǎn)之后吸附到最后一個頂點(diǎn)見points in front of last point測試。參考點(diǎn)恰在折點(diǎn)joint上原地不動見points on joints測試。參考點(diǎn)位于線上位移不超過1e-6量級見points on top of line、point along line測試其中用turf/along生成線上點(diǎn)做交叉驗證。線段中點(diǎn)與索引segmentIndex與pointDistance的精確數(shù)值被斷言如segmentIndex 8、pointDistance 0.823802。Geometry / 坐標(biāo)數(shù)組輸入四種輸入組合均可正常調(diào)用見Geometry Support測試。MultiLineString 索引驗證lineStringIndex正確指到第 1 條 LineString見lineStringIndex測試。歷史 issue 回歸包括 #691totalDistance不得為 NaN、#2753MultiLineString 的lineDistance精度、#1514經(jīng)線吸附精度、#965長線吸附、#2808 與重復(fù)點(diǎn)支持含重復(fù)頂點(diǎn)不得拋錯、#2934跨半球時須選中正確的南/北端點(diǎn)、#2939微小線段在瀏覽器與 Node 的浮點(diǎn)差異下仍返回正確端點(diǎn)等。測試數(shù)據(jù)位于 test/in 目錄每個 fixture 是一個 FeatureCollection第一個要素為線第二個要素為參考點(diǎn)例如 line1.geojson、multiLine1.geojson、route1.geojson近 5000 行的真實路線數(shù)據(jù)、line-northern-latitude-#344.geojson高緯地區(qū)反足退化場景。如需復(fù)現(xiàn)可在包目錄下運(yùn)行pnpm run test:tape基于 tape 的斷言測試與pnpm run test:types類型檢查見 package.json。性能基準(zhǔn)模塊自帶 Benchmark 基準(zhǔn)腳本 bench.ts用test/in下的 fixture 逐一測速。從腳本中保留的歷史結(jié)果看短線段數(shù)據(jù)集如line1可達(dá)數(shù)十萬 ops/sec而真實長路線route1、route2約為百余 ops/sec——這符合算法按線段數(shù)線性遍歷的特性。你可以通過pnpm run bench在當(dāng)前硬件上復(fù)現(xiàn)自己的基準(zhǔn)數(shù)據(jù)。與其他 Turf 模塊的協(xié)同turf/turf全量包nearestPointOnLine與nearestPoint、nearestPointToLine、nearestNeighborAnalysis一并導(dǎo)出見 packages/turf/index.ts三者適用場景不同nearestPointOnLine求線上最近點(diǎn)nearestPointToLine求點(diǎn)到線的垂足/最近點(diǎn)距離nearestPoint求點(diǎn)集內(nèi)最近鄰。turf/distance本模塊的pointDistance、segmentDistance等距離字段全部由turf/distance按指定units計算單位語義與 Turf 全局一致。turf/along/turf/length測試中常用于交叉驗證——先along取線上點(diǎn)再nearestPointOnLine吸附斷言位移趨近于零也可用于組合實現(xiàn)沿線的里程定位類業(yè)務(wù)。turf/meta的flattenEach是支撐MultiLineString輸入的關(guān)鍵基礎(chǔ)能力。總結(jié)nearestPointOnLine是一個實現(xiàn)簡潔、行為嚴(yán)謹(jǐn)?shù)狞c(diǎn)吸附工具對外提供三種點(diǎn)輸入形態(tài)、完整的新舊屬性語義與全部 Turf 單位支持對內(nèi)以逐段短路比對 三維單位球大圓投影完成最近點(diǎn)求解并顯式處理了反足、重合、微小線段等退化場景。無論是做地圖交互吸附、軌跡糾偏還是路線里程定位它都能直接作為可靠的底層函數(shù)使用其行為邊界也可以通過 test.ts 與 test/in 中的真實數(shù)據(jù) fixtures 充分驗證。贊分享數(shù)據(jù)分析【免費(fèi)下載鏈接】turfA modular geospatial engine written in JavaScript and TypeScript項目地址https://gitcode.com/gh_mirrors/tu/turf點(diǎn)擊查看免費(fèi)下載相關(guān)推薦搞定幾何算法從凸包到最近點(diǎn)對的Java實現(xiàn)指南搞定幾何算法從凸包到最近點(diǎn)對的Java實現(xiàn)指南 你是否還在為計算平面點(diǎn)集的凸包而煩惱是否想知道如何高效找到海量點(diǎn)中的最近點(diǎn)對本文將帶你一步掌握這兩個核心幾示例工程gh_mirrors/cs/C-Sharp中的幾何算法點(diǎn)、線、面計算實現(xiàn)gh_mirrors/cs/C Sharp中的幾何算法點(diǎn)、線、面計算實現(xiàn) 幾何算法是計算機(jī)科學(xué)中的重要基礎(chǔ)廣泛應(yīng)用于地圖服務(wù)、空間分析、游戲開發(fā)等領(lǐng)域。在示例工程科學(xué)計算TDengine 發(fā)布?xì)v史全覽版本號規(guī)則、升級兼容性與 3.x Release Notes 解讀TDengine 發(fā)布?xì)v史全覽版本號規(guī)則、升級兼容性與 3.x Release Notes 解讀 本文以 docs/en/17 release history數(shù)據(jù)分析上一篇在 Zed 中接入 daisyUI MCP ServerBlueprint、Context7 與 GitMCP 三種方案配置實戰(zhàn)下一篇Deepagents在線教育在線教育平臺的AI代理終極指南創(chuàng)作聲明：本文部分內(nèi)容由AI輔助生成（AIGC），僅供參考 Turf nearestPointOnLine 實戰(zhàn)指南：將任意點(diǎn)吸附到線上的最近點(diǎn)計算與球面幾何實現(xiàn)原理
數(shù)據(jù)分析【免費(fèi)下載鏈接】turfA modular geospatial engine written in JavaScript and TypeScript項目地址https://gitcode.com/gh_mirrors/tu/turf點(diǎn)擊查看免費(fèi)下載turf/nearest-point-on-line是 Turf 模塊化地理引擎中專門用于點(diǎn)吸附snap to line的核心工具給定一條LineString或MultiLineString與一個參考點(diǎn)它返回線上距離該點(diǎn)最近的那個位置點(diǎn)并附帶完整的距離與索引信息。本指南將帶你完整掌握它的 API、返回值語義含 v7.4 新增屬性與舊屬性對照、球面大圓幾何的底層實現(xiàn)以及如何通過倉庫源碼與測試用例驗證其行為邊界——讀完后你可以在路徑糾偏、地圖吸附、導(dǎo)航匹配等場景中直接落地使用。功能概述與典型應(yīng)用場景nearestPointOnLine解決的問題非常明確求一條線上離某個給定點(diǎn)最近的點(diǎn)。它接收任意形式的線段數(shù)據(jù)LineString/MultiLineString的 Geometry 或 Feature和任意形式的點(diǎn)數(shù)據(jù)Point的 Geometry、Feature甚至直接是[lng, lat]坐標(biāo)數(shù)組返回一個帶豐富屬性的PointFeature。典型的實戰(zhàn)場景包括地圖交互吸附用戶在地圖上點(diǎn)擊時把點(diǎn)擊位置吸附到道路、河流等線性要素上實現(xiàn)自動貼線。路徑糾偏 / 軌跡匹配GPS 采樣點(diǎn)偏移到實際道路之外時將軌跡點(diǎn)拉回到路網(wǎng)上。路線輔助計算結(jié)合返回的totalDistance、segmentIndex等屬性定位點(diǎn)在線路上的里程位置進(jìn)而配合turf/along等模塊做分段分析。該模塊屬于 Turf 的獨(dú)立子包版本為 7.4.0被 turf 主包統(tǒng)一導(dǎo)出同時依賴turf/distance測距、turf/helpers要素構(gòu)造、turf/invariant坐標(biāo)校驗與turf/meta要素遍歷四個基礎(chǔ)模塊見 package.json。安裝與引入單獨(dú)安裝單個模塊$ npm install turf/nearest-point-on-line該包聲明為 ESM 模塊type: module運(yùn)行環(huán)境要求 Node.js 22見 package.json。安裝全量 turf/turf如果你希望使用 Turf 全家桶一次性引入所有功能函數(shù)$ npm install turf/turf之后即可通過turf.nearestPointOnLine(...)調(diào)用全量包的導(dǎo)出定義在 packages/turf/index.ts。在代碼中引入import { nearestPointOnLine } from turf/nearest-point-on-line; // 或從全量包引入 // import { nearestPointOnLine } from turf/turf;API 簽名與參數(shù)詳解函數(shù)簽名function nearestPointOnLineG extends LineString | MultiLineString( lines: FeatureG | G, // 要吸附的線 inputPoint: Coord, // 參考點(diǎn)Geometry / FeaturePoint / number[] options: { units?: Units } {} // 可選參數(shù) ): FeaturePoint, {...屬性...}完整的類型定義見 index.ts。參數(shù)說明參數(shù)類型說明linesGeometry \| FeatureLineString \| MultiLineString要吸附到的線要素可以是LineString或MultiLineString的 Geometry 或 FeatureinputPointGeometry \| FeaturePoint \| number[]參考點(diǎn)支持三種形態(tài)Point的 Geometry、PointFeature或直接傳[lng, lat]坐標(biāo)數(shù)組options.unitsUnits距離單位默認(rèn)kilometersoptions.units 支持的取值options.units支持 Turf 的全部合法單位定義于 packages/turf-helpers/index.tsmeters | metres | m | millimeters | millimetres | mm | centimeters | centimetres | cm | kilometers | kilometres | km | miles | mi | nauticalmiles | nmi | inches | in | yards | yd | feet | ft | radians | rad | degrees | deg注意radians/degrees等角度單位在距離語義上具有緯度依賴性——經(jīng)度方向上每一度所代表的實際距離隨緯度變化在地球上不同位置度量差異很大請僅在明確理解其含義時使用。返回值屬性v7.4 新屬性與舊屬性對照函數(shù)返回線上離inputPoint最近的點(diǎn)FeaturePoint。自 v7.4 起Turf 開始遷移返回值屬性命名新的推薦屬性與舊屬性一一對應(yīng)如下新屬性v7.4 起推薦舊屬性仍兼容已標(biāo)記 deprecated含義lineStringIndexmultiFeatureIndex最近點(diǎn)落在輸入MultiLineString的第幾條 LineString 上輸入為 LineString 時恒為 0segmentIndexindex最近點(diǎn)落在該 LineString 的第幾個線段上totalDistancelocation從整條 MultiLineString 絕對起點(diǎn)算起的沿線距離lineDistance—舊版無對應(yīng)從最近點(diǎn)所在 LineString 起點(diǎn)算起的沿線距離segmentDistance—舊版無對應(yīng)從最近點(diǎn)所在線段起點(diǎn)算起的沿線距離pointDistancedist參考點(diǎn)與最近點(diǎn)之間的直線距離官方保證multiFeatureIndex、index、location、dist至少在下一個大版本發(fā)布之前繼續(xù)可用見 README.md 與 index.ts 的 JSDoc 說明。源碼在生成新屬性的同時會同步將multiFeatureIndex/index/location/dist指向?qū)?yīng)的新屬性值見 index.ts因此新舊寫法可以并存。類型層面新屬性與舊屬性都被聲明在返回類型的 properties 中其中舊屬性帶有deprecated標(biāo)記見 index.ts。建議新代碼直接使用新屬性以便在未來的大版本升級中平滑遷移?；A(chǔ)用法示例官方示例來源 README.md演示一條折線與一個參考點(diǎn)以英里為單位求吸附點(diǎn)var line turf.lineString([ [-77.031669, 38.878605], [-77.029609, 38.881946], [-77.020339, 38.884084], [-77.025661, 38.885821], [-77.021884, 38.889563], [-77.019824, 38.892368] ]); var inputPoint turf.point([-77.037076, 38.884017]); var snapped turf.nearestPointOnLine(line, inputPoint, {units: miles}); // addToMap var addToMap [line, inputPoint, snapped]; snapped.properties[marker-color] #00f;上面的例子中snapped即為吸附后的最近點(diǎn)其geometry.coordinates是線上的經(jīng)緯度坐標(biāo)properties中攜帶上文表格里的全部索引與距離字段。需要展示時可將原始線、參考點(diǎn)和吸附點(diǎn)一并加入地圖圖層并給吸附點(diǎn)設(shè)置marker-color以突出顯示。更完整的可運(yùn)行示例import { nearestPointOnLine } from turf/nearest-point-on-line; import { point, lineString, multiLineString } from turf/helpers; const pt point([1.5, 1.5]); const line lineString([ [0, 0], [1, 1], ]); const multiLine multiLineString([ [ [0, 0], [1, 1], [2, 2], [0, 0], ], ]); // 各種合法的參數(shù)組合 nearestPointOnLine(line, pt); // 默認(rèn)單位 kilometers nearestPointOnLine(multiLine, pt); // 支持 MultiLineString nearestPointOnLine(line.geometry, pt); // 直接傳 Geometry nearestPointOnLine(multiLine.geometry, pt); // 直接傳 MultiLineString Geometry nearestPointOnLine(line, pt, { units: miles }); // 指定單位 nearestPointOnLine(line, [1.5, 1.5]); // 傳坐標(biāo)數(shù)組 // 輸出可以作為下一次輸入二次吸附 const output nearestPointOnLine(line, pt); nearestPointOnLine(line, output); // 新舊屬性均可讀取 output.properties.dist; // 舊屬性 output.properties.multiFeatureIndex; // 舊屬性 output.properties.index; // 舊屬性 output.properties.location; // 舊屬性 output.properties.pointDistance; // 新屬性推薦 output.properties.totalDistance; // 新屬性推薦上述參數(shù)組合與屬性訪問方式均由 types.ts 中的類型測試覆蓋保證編譯期可驗證。源碼實現(xiàn)原理逐段求最近點(diǎn)與球面幾何深入 index.ts可以看清整個算法的執(zhí)行流程1. 參數(shù)校驗與逐要素遍歷函數(shù)首先校驗lines與inputPoint是否都提供了否則拋出lines and inputPoint are required arguments見 index.ts。隨后通過getCoord提取參考點(diǎn)坐標(biāo)并用turf/meta的flattenEach將輸入包括 MultiLineString扁平化為逐條 LineString 遍歷同時記錄lineStringIndex見 index.ts。2. 逐線段計算并維護(hù)累計距離對每條 LineString 的每一段coords[i]到coords[i 1]算法先計算該段長度segmentLength distance(start, stop, options)再求參考點(diǎn)在該段上的垂足最近點(diǎn)短路優(yōu)化如果參考點(diǎn)恰好等于某線段的端點(diǎn)則直接采用該端點(diǎn)不再做投影計算見 index.ts。這里特意先比較終點(diǎn)再比較起點(diǎn)以保證起點(diǎn)與終點(diǎn)重合時行為一致。計算候選點(diǎn)與參考點(diǎn)的距離pointDistance若小于當(dāng)前最優(yōu)解則更新closestPt并同步累計出segmentDistance、lineDistance、totalDistance。每處理完一段就把totalDistance與lineDistance累加該段長度當(dāng)flattenEach切換 LineString 時lineDistance重置為 0見 index.ts。這樣totalDistance恒為從整條 MultiLineString 絕對起點(diǎn)到最近點(diǎn)的沿線里程lineDistance為該最近點(diǎn)所在 LineString 內(nèi)的里程segmentDistance為該最近點(diǎn)所在線段內(nèi)的里程。3. 核心幾何nearestPointOnSegment 的大圓投影真正的投影邏輯在nearestPointOnSegment(posA, posB, posC)中見 index.ts其核心思想是把經(jīng)緯度坐標(biāo)轉(zhuǎn)換到三維單位球面的直角坐標(biāo)系中再通過球面大圓幾何求最近點(diǎn)坐標(biāo)球面化lngLatToVector將[lng, lat]轉(zhuǎn)為單位向量(x, y, z)見 index.ts。構(gòu)造大圓平面法向量線段兩端點(diǎn)向量A、B的叉積cross(A, B)即該線段所在大圓平面的法向量segmentAxis。求過參考點(diǎn)的大圓交點(diǎn)依次計算targetAxis cross(segmentAxis, C)、intersectionAxis cross(targetAxis, segmentAxis)歸一化后得到大圓的兩個交點(diǎn)I1/I2用點(diǎn)積余弦相似度挑選離參考點(diǎn)C更近的那個作為候選交點(diǎn)I。判斷交點(diǎn)是否落在線段弧上通過兩次叉積點(diǎn)積比較cross(A, I)與cross(I, B)相對segmentAxis的方向一致性判斷旋轉(zhuǎn)方向若候選交點(diǎn)確實位于 A→B 的弧段內(nèi)則直接返回否則退化為選擇距離參考點(diǎn)更近的端點(diǎn) A 或 B見 index.ts。整個算法對每條線段只做常數(shù)次向量運(yùn)算配合前面先比端點(diǎn)、再算投影的短路邏輯復(fù)雜度為 O(線段數(shù))。4. 退化情況與反足點(diǎn)處理源碼對兩類數(shù)學(xué)退化情況做了顯式處理線段兩端向量平行aligned或反足antipodal此時叉積結(jié)果為[0, 0, 0]無法構(gòu)造大圓平面。平行dot(A, B) 0意味著 A、B 重合算法直接返回終點(diǎn)保持一致行為反足意味著 A、B 之間存在無數(shù)個大圓弧未定義此時選擇返回參考點(diǎn)本身語義上等同于該點(diǎn)位于線上見 index.ts。這正是 README 開頭所說明的若輸入 LineString 中存在反足線段弧未定義函數(shù)將返回點(diǎn)位于線上。候選交點(diǎn)與參考點(diǎn)重合/反足targetAxis為 0 向量參考點(diǎn)到整個線段等距為保持一致返回終點(diǎn)見 index.ts。此外vectorToLngLat在反解經(jīng)緯度時對 z 分量做了[-1, 1]鉗制避免浮點(diǎn)誤差導(dǎo)致asin越界見 index.ts。邊界行為與測試驗證模塊的邊界行為在 test.ts 中有大量回歸測試覆蓋可直接作為行為文檔閱讀參考點(diǎn)恰為線起點(diǎn)/終點(diǎn)點(diǎn)不移動totalDistance分別為 0 或全線長。參考點(diǎn)在第一個頂點(diǎn)之前吸附到第一個頂點(diǎn)totalDistance為 0見points behind first point測試。參考點(diǎn)在最后一個頂點(diǎn)之后吸附到最后一個頂點(diǎn)見points in front of last point測試。參考點(diǎn)恰在折點(diǎn)joint上原地不動見points on joints測試。參考點(diǎn)位于線上位移不超過1e-6量級見points on top of line、point along line測試其中用turf/along生成線上點(diǎn)做交叉驗證。線段中點(diǎn)與索引segmentIndex與pointDistance的精確數(shù)值被斷言如segmentIndex 8、pointDistance 0.823802。Geometry / 坐標(biāo)數(shù)組輸入四種輸入組合均可正常調(diào)用見Geometry Support測試。MultiLineString 索引驗證lineStringIndex正確指到第 1 條 LineString見lineStringIndex測試。歷史 issue 回歸包括 #691totalDistance不得為 NaN、#2753MultiLineString 的lineDistance精度、#1514經(jīng)線吸附精度、#965長線吸附、#2808 與重復(fù)點(diǎn)支持含重復(fù)頂點(diǎn)不得拋錯、#2934跨半球時須選中正確的南/北端點(diǎn)、#2939微小線段在瀏覽器與 Node 的浮點(diǎn)差異下仍返回正確端點(diǎn)等。測試數(shù)據(jù)位于 test/in 目錄每個 fixture 是一個 FeatureCollection第一個要素為線第二個要素為參考點(diǎn)例如 line1.geojson、multiLine1.geojson、route1.geojson近 5000 行的真實路線數(shù)據(jù)、line-northern-latitude-#344.geojson高緯地區(qū)反足退化場景。如需復(fù)現(xiàn)可在包目錄下運(yùn)行pnpm run test:tape基于 tape 的斷言測試與pnpm run test:types類型檢查見 package.json。性能基準(zhǔn)模塊自帶 Benchmark 基準(zhǔn)腳本 bench.ts用test/in下的 fixture 逐一測速。從腳本中保留的歷史結(jié)果看短線段數(shù)據(jù)集如line1可達(dá)數(shù)十萬 ops/sec而真實長路線route1、route2約為百余 ops/sec——這符合算法按線段數(shù)線性遍歷的特性。你可以通過pnpm run bench在當(dāng)前硬件上復(fù)現(xiàn)自己的基準(zhǔn)數(shù)據(jù)。與其他 Turf 模塊的協(xié)同turf/turf全量包nearestPointOnLine與nearestPoint、nearestPointToLine、nearestNeighborAnalysis一并導(dǎo)出見 packages/turf/index.ts三者適用場景不同nearestPointOnLine求線上最近點(diǎn)nearestPointToLine求點(diǎn)到線的垂足/最近點(diǎn)距離nearestPoint求點(diǎn)集內(nèi)最近鄰。turf/distance本模塊的pointDistance、segmentDistance等距離字段全部由turf/distance按指定units計算單位語義與 Turf 全局一致。turf/along/turf/length測試中常用于交叉驗證——先along取線上點(diǎn)再nearestPointOnLine吸附斷言位移趨近于零也可用于組合實現(xiàn)沿線的里程定位類業(yè)務(wù)。turf/meta的flattenEach是支撐MultiLineString輸入的關(guān)鍵基礎(chǔ)能力。總結(jié)nearestPointOnLine是一個實現(xiàn)簡潔、行為嚴(yán)謹(jǐn)?shù)狞c(diǎn)吸附工具對外提供三種點(diǎn)輸入形態(tài)、完整的新舊屬性語義與全部 Turf 單位支持對內(nèi)以逐段短路比對 三維單位球大圓投影完成最近點(diǎn)求解并顯式處理了反足、重合、微小線段等退化場景。無論是做地圖交互吸附、軌跡糾偏還是路線里程定位它都能直接作為可靠的底層函數(shù)使用其行為邊界也可以通過 test.ts 與 test/in 中的真實數(shù)據(jù) fixtures 充分驗證。贊分享數(shù)據(jù)分析【免費(fèi)下載鏈接】turfA modular geospatial engine written in JavaScript and TypeScript項目地址https://gitcode.com/gh_mirrors/tu/turf點(diǎn)擊查看免費(fèi)下載相關(guān)推薦搞定幾何算法從凸包到最近點(diǎn)對的Java實現(xiàn)指南搞定幾何算法從凸包到最近點(diǎn)對的Java實現(xiàn)指南 你是否還在為計算平面點(diǎn)集的凸包而煩惱是否想知道如何高效找到海量點(diǎn)中的最近點(diǎn)對本文將帶你一步掌握這兩個核心幾示例工程gh_mirrors/cs/C-Sharp中的幾何算法點(diǎn)、線、面計算實現(xiàn)gh_mirrors/cs/C Sharp中的幾何算法點(diǎn)、線、面計算實現(xiàn) 幾何算法是計算機(jī)科學(xué)中的重要基礎(chǔ)廣泛應(yīng)用于地圖服務(wù)、空間分析、游戲開發(fā)等領(lǐng)域。在示例工程科學(xué)計算TDengine 發(fā)布?xì)v史全覽版本號規(guī)則、升級兼容性與 3.x Release Notes 解讀TDengine 發(fā)布?xì)v史全覽版本號規(guī)則、升級兼容性與 3.x Release Notes 解讀 本文以 docs/en/17 release history數(shù)據(jù)分析上一篇在 Zed 中接入 daisyUI MCP ServerBlueprint、Context7 與 GitMCP 三種方案配置實戰(zhàn)下一篇Deepagents在線教育在線教育平臺的AI代理終極指南創(chuàng)作聲明：本文部分內(nèi)容由AI輔助生成（AIGC），僅供參考