雷達(dá)和AIS是水上交通監(jiān)控里最常用也最容易“打架”的兩類信息來(lái)源。雷達(dá)每個(gè)掃描周期都能給出目標(biāo)的位置但不知道這條船是誰(shuí)AIS每隔幾秒到幾分鐘報(bào)一次船名、航速、航向但目標(biāo)一旦關(guān)機(jī)或者報(bào)文延遲監(jiān)控畫面就會(huì)出現(xiàn)一條“隱身船”。把兩者的航跡融合成一條準(zhǔn)確、連續(xù)、帶身份的目標(biāo)軌跡是這個(gè)領(lǐng)域繞不開(kāi)的核心工作。我最近用一個(gè)基于MATLAB的雷達(dá)與AIS航跡融合處理項(xiàng)目把從數(shù)據(jù)預(yù)處理、時(shí)空對(duì)齊、航跡關(guān)聯(lián)到融合濾波的完整鏈路跑通了一遍這里把思路、代碼和踩過(guò)的坑一起整理出來(lái)。內(nèi)容主要面向海事雷達(dá)工程師、做水上交通信息處理的學(xué)生以及正在寫多傳感器融合代碼的朋友如果你只想快速跑通一個(gè)Demo也可以直接照著核心邏輯改。1. 融合需求剖析與整體策略1.1 雷達(dá)與AIS的互補(bǔ)關(guān)系以及為什么非融合不可雷達(dá)是主動(dòng)傳感器自己發(fā)射電磁波、接收回波好處是不依賴目標(biāo)配合任何有雷達(dá)反射面的目標(biāo)都能被探測(cè)到壞處也很明顯它只給“這里有個(gè)東西”不給“這是誰(shuí)”。在近岸、彎曲航道、雨雪天氣下雷達(dá)點(diǎn)跡里還混著大量海雜波和假回波距離誤差隨距離增大而增大方位誤差在遠(yuǎn)區(qū)會(huì)被放大成好幾倍的橫向偏差。AIS則完全是另一套邏輯。它通過(guò)船載設(shè)備主動(dòng)廣播本船的靜態(tài)信息船名、MMSI、船長(zhǎng)船寬和動(dòng)態(tài)信息經(jīng)緯度、對(duì)地航速、對(duì)地航向定位用的是GNSS位置精度通常在幾米到十幾米遠(yuǎn)遠(yuǎn)好于雷達(dá)。但它是“合作式”系統(tǒng)船上沒(méi)裝設(shè)備、設(shè)備關(guān)機(jī)或報(bào)文被干擾時(shí)就一點(diǎn)數(shù)據(jù)都沒(méi)有。而且AIS更新率隨航行狀態(tài)變化錨泊船可能幾分鐘才發(fā)一條船多了報(bào)文還可能排隊(duì)延遲。我實(shí)際調(diào)試時(shí)遇到過(guò)這種情況監(jiān)控畫面里一個(gè)雷達(dá)回波和一個(gè)AIS目標(biāo)相距約200米看起來(lái)是兩條船實(shí)際上雷達(dá)回波的位置噪聲偏大、AIS定位又很準(zhǔn)兩者指向的其實(shí)是同一艘船。如果這時(shí)候直接拿兩個(gè)位置做平均身份和位置都會(huì)一起漂掉。所以融合系統(tǒng)的第一件事不是算平均而是先回答“這兩個(gè)測(cè)量到底是不是同一個(gè)目標(biāo)”這個(gè)問(wèn)題答錯(cuò)了后面所有的濾波和加權(quán)都沒(méi)有意義。1.2 融合質(zhì)量的評(píng)價(jià)指標(biāo)融合做得好不好不能靠肉眼和感覺(jué)。我自己的習(xí)慣是先把評(píng)價(jià)指標(biāo)定下來(lái)再動(dòng)手寫算法。做過(guò)這個(gè)方向的同行應(yīng)該都有同感沒(méi)有量化指標(biāo)調(diào)試時(shí)很難說(shuō)清楚一個(gè)改動(dòng)是變好了還是變壞了。指標(biāo)工程含義怎么衡量位置誤差融合航跡與真值的偏差與高精度定位或人工標(biāo)定真值比較取均方根誤差身份保持率航跡號(hào)和MMSI綁定的穩(wěn)定程度統(tǒng)計(jì)一段時(shí)間內(nèi)身份未跳變的軌跡占比軌跡完整率目標(biāo)可連續(xù)跟蹤的程度目標(biāo)全程未斷裂、未丟失的軌跡占比航向航速跳變率融合輸出的平滑程度相鄰兩次更新中航向或航速突變超過(guò)設(shè)定閾值的次數(shù)誤關(guān)聯(lián)率關(guān)聯(lián)算法張冠李戴的比例錯(cuò)誤配對(duì)點(diǎn)跡數(shù)占總關(guān)聯(lián)次數(shù)的比例這些指標(biāo)不需要做得很重一個(gè)簡(jiǎn)單的離線統(tǒng)計(jì)腳本就夠了。每次改動(dòng)算法或參數(shù)后跑同一段回放數(shù)據(jù)對(duì)比這些數(shù)字的變化比憑感覺(jué)判斷可靠得多。2. 航跡關(guān)聯(lián)前的時(shí)空對(duì)齊處理2.1 坐標(biāo)統(tǒng)一把經(jīng)緯度與極坐標(biāo)放到同一平面上雷達(dá)原始輸出通常是極坐標(biāo)距離和方位角AIS輸出是WGS84經(jīng)緯度加對(duì)地速度、對(duì)地航向。要把兩者做關(guān)聯(lián)第一步是放到同一個(gè)平面坐標(biāo)系里。工程上最常用的公共參考系是以雷達(dá)站為原點(diǎn)的本地切平面直角坐標(biāo)也就是東向x、北向y的ENU坐標(biāo)系。小范圍水域幾十公里以內(nèi)可以用等距近似轉(zhuǎn)換公式很簡(jiǎn)單function [x, y] ll2xy(lat, lon, lat0, lon0) % 簡(jiǎn)化等距轉(zhuǎn)換適合幾十公里內(nèi)的監(jiān)控場(chǎng)景 R 6371000; % 地球近似半徑單位m x (lon - lon0) * pi/180 * R * cos(lat0*pi/180); y (lat - lat0) * pi/180 * R; end這段代碼夠用但要注意兩點(diǎn)。第一參考點(diǎn)lat0、lon0一定要選在監(jiān)控區(qū)域中心偏離太遠(yuǎn)時(shí)等距近似的誤差會(huì)明顯變大。第二更嚴(yán)謹(jǐn)?shù)膱?chǎng)景建議直接用投影轉(zhuǎn)換把經(jīng)緯度投影到UTM坐標(biāo)MATLAB里用自帶函數(shù)創(chuàng)建投影對(duì)象后做正向轉(zhuǎn)換這樣在大范圍、高緯度區(qū)域也不會(huì)出問(wèn)題。雷達(dá)極坐標(biāo)轉(zhuǎn)直角坐標(biāo)也要統(tǒng)一口徑x rho * sin(theta); % theta從正北順時(shí)針單位弧度 y rho * cos(theta);不同雷達(dá)設(shè)備輸出的方位角定義可能不一樣有的是從正北順時(shí)針有的是從正東逆時(shí)針轉(zhuǎn)坐標(biāo)前必須確認(rèn)設(shè)備手冊(cè)里的角度定義否則整個(gè)關(guān)聯(lián)都會(huì)錯(cuò)位。2.2 時(shí)間同步用運(yùn)動(dòng)模型把航跡外推到同一時(shí)刻AIS和雷達(dá)的采樣時(shí)刻幾乎不可能完全一致。雷達(dá)掃描周期通常是2到6秒AIS更新周期從幾秒到幾分鐘不等二者在同一時(shí)刻點(diǎn)上的概率很低。時(shí)間差小于一個(gè)閾值時(shí)可以直接忽略大于閾值就必須做外推。外推用最樸素的運(yùn)動(dòng)模型就夠了假設(shè)目標(biāo)在短時(shí)間內(nèi)外推窗口內(nèi)保持勻速直線運(yùn)動(dòng)function [x_est, y_est] extrapolate(x0, y0, v, cog, dt) % v為對(duì)地航速單位m/scog為對(duì)地航向單位弧度dt為外推時(shí)間差單位s x_est x0 v * sin(cog) * dt; y_est y0 v * cos(cog) * dt; end這里的方向定義要跟坐標(biāo)轉(zhuǎn)換保持一致。如果目標(biāo)航速接近零比如錨泊船直接沿用舊位置不要硬外推。外推時(shí)間越長(zhǎng)位置誤差越大所以一般建議外推窗口不超過(guò)10到30秒。AIS的SOG和COG本身就是GNSS解算出來(lái)的平滑值短時(shí)間外推的精度還可以雷達(dá)點(diǎn)跡如果也要做時(shí)間對(duì)齊它的速度估計(jì)噪聲大外推時(shí)要加約束防止位置被推到離譜的地方。時(shí)間對(duì)齊做完后每個(gè)待關(guān)聯(lián)的點(diǎn)都必須打上“測(cè)量時(shí)間戳”和“外推后的位置”后續(xù)關(guān)聯(lián)用的是外推后的位置但原始測(cè)量值要保留下來(lái)方便排查問(wèn)題。3. 航跡關(guān)聯(lián)判定與融合濾波實(shí)現(xiàn)3.1 關(guān)聯(lián)門限與最近鄰分配關(guān)聯(lián)判定的核心邏輯是在公共坐標(biāo)系下計(jì)算雷達(dá)航跡點(diǎn)和AIS航跡點(diǎn)之間的距離、速度差、航向差然后判斷它們是否屬于同一個(gè)目標(biāo)。位置門限可以用誤差傳播的思路來(lái)定。設(shè)雷達(dá)測(cè)距誤差標(biāo)準(zhǔn)差為σ_RAIS定位誤差標(biāo)準(zhǔn)差為σ_A那么兩者差的方差是兩者方差之和取3倍標(biāo)準(zhǔn)差作為門限覆蓋大約99.7%的隨機(jī)誤差同時(shí)給目標(biāo)機(jī)動(dòng)留一點(diǎn)余量G_p 3 * sqrt(σ_R^2 σ_A^2)舉個(gè)例子σ_R取50米σ_A取10米則G_p約為153米。實(shí)際工程里為了應(yīng)對(duì)目標(biāo)機(jī)動(dòng)和航跡外推誤差我一般會(huì)把門限再放寬到300米左右再配合速度和航向門限做細(xì)篩。速度門限取Δv ≤ 5節(jié)航向門限取Δφ ≤ 30度這兩個(gè)約束能擋掉大量近距平行航行的錯(cuò)誤配對(duì)。匹配策略上簡(jiǎn)單場(chǎng)景用最近鄰就夠了多目標(biāo)交叉、編隊(duì)航行的場(chǎng)景強(qiáng)烈建議用全局最優(yōu)分配避免出現(xiàn)一個(gè)AIS目標(biāo)同時(shí)匹配兩個(gè)雷達(dá)目標(biāo)這種沖突。MATLAB里可以用匹配函數(shù)做最小代價(jià)指派代價(jià)矩陣就是距離矩陣超過(guò)門限的項(xiàng)設(shè)置成大數(shù)或者直接排除function [pairs] associate_gnn(Xr, Xa, gate) % Xr: 雷達(dá)點(diǎn)跡坐標(biāo) [Nx2]Xa: AIS點(diǎn)跡坐標(biāo) [Mx2] N size(Xr,1); M size(Xa,1); D sqrt((Xr(:,1) - Xa(:,1)).^2 (Xr(:,2) - Xa(:,2)).^2); % 超過(guò)門限的代價(jià)置為很大的數(shù)讓求解器不要選它 D(D gate) 1e9; [assignment, ~] matchpairs(D, min); % assignment列出的是有效配對(duì)的索引 pairs assignment; end沒(méi)裝優(yōu)化工具箱的讀者也別急這個(gè)最小距離分配用兩層循環(huán)加一個(gè)占位標(biāo)記也能寫性能差別在點(diǎn)跡數(shù)量不大的場(chǎng)景下完全體現(xiàn)不出來(lái)。提醒一句不要每幀都對(duì)全部點(diǎn)跡做一次全量關(guān)聯(lián)那樣容易造成身份抖動(dòng)。推薦以AIS航跡為“錨”雷達(dá)點(diǎn)跡來(lái)配對(duì)更新AIS沒(méi)有覆蓋到的目標(biāo)雷達(dá)點(diǎn)跡可以獨(dú)立建航跡以后再嘗試關(guān)聯(lián)。3.2 位置融合從加權(quán)平均到卡爾曼濾波關(guān)聯(lián)配對(duì)成功后才輪到真正的“融合”。需要說(shuō)清楚的是融合不是簡(jiǎn)單的坐標(biāo)平均而是按誤差逆加權(quán)。雷達(dá)位置誤差σ_R遠(yuǎn)大于AIS的σ_A但AIS更新慢所以融合結(jié)果應(yīng)當(dāng)偏向“當(dāng)前時(shí)刻更可信”的信息源。加權(quán)融合公式如下x_f (σ_A^2 * x_r σ_R^2 * x_a) / (σ_R^2 σ_A^2) y_f (σ_A^2 * y_r σ_R^2 * y_a) / (σ_R^2 σ_A^2)這個(gè)公式本質(zhì)是最小方差估計(jì)。MATLAB實(shí)現(xiàn)起來(lái)就是幾行function [xf, yf] fuse_weighted(xr, yr, xa, ya, sigmaR, sigmaA) wR sigmaA^2 / (sigmaR^2 sigmaA^2); wA sigmaR^2 / (sigmaR^2 sigmaA^2); xf wR * xr wA * xa; yf wR * yr wA * ya; end想更進(jìn)一步可以用常速度模型做卡爾曼濾波狀態(tài)取[x, y, vx, vy]雷達(dá)和AIS作為兩類量測(cè)輸入。MATLAB自帶跟蹤濾波函數(shù)可以直接調(diào)用但我更建議自己把濾波核心寫一遍狀態(tài)轉(zhuǎn)移矩陣和量測(cè)更新邏輯才不會(huì)變成黑盒。量測(cè)更新部分的核心代碼大概是function [x_est, P_est] kf_update(x_pred, P_pred, z, H, R) K P_pred * H / (H * P_pred * H R); x_est x_pred K * (z - H * x_pred); P_est (eye(size(P_pred)) - K * H) * P_pred; end預(yù)測(cè)部分就是狀態(tài)轉(zhuǎn)移矩陣乘狀態(tài)、狀態(tài)協(xié)方差乘轉(zhuǎn)移矩陣轉(zhuǎn)置再加過(guò)程噪聲協(xié)方差。工程上Q取多大、R取多大跟場(chǎng)景強(qiáng)相關(guān)不要迷信某個(gè)固定數(shù)值要根據(jù)回放數(shù)據(jù)里的跳動(dòng)幅度來(lái)調(diào)。3.3 航向航速融合與航跡身份管理位置融合之外SOG和COG的融合也有講究。AIS的SOG、COG來(lái)自GNSS解算統(tǒng)計(jì)特性比雷達(dá)點(diǎn)跡差分出來(lái)的速度穩(wěn)得多所以在AIS報(bào)文“新鮮”的時(shí)候直接用AIS的SOG和COG作為融合輸出不做額外計(jì)算。AIS報(bào)文超過(guò)一定時(shí)間后就不新鮮了比如超過(guò)30秒這時(shí)需要用雷達(dá)點(diǎn)跡的位置差分估算速度航向作為輔助。雷達(dá)差分速度噪聲大我一般會(huì)再做一次中值濾波把跳變點(diǎn)壓掉。這里的原則是優(yōu)先用可信源備用源做兜底。身份綁定是融合系統(tǒng)里容易被忽略卻很重要的環(huán)節(jié)。融合航跡號(hào)與AIS的MMSI綁定后即使AIS信號(hào)中斷航跡身份也要保持住不能等AIS重新出現(xiàn)時(shí)又當(dāng)成新目標(biāo)。身份管理的核心就是“綁定一次保持到底”直到確認(rèn)目標(biāo)離開(kāi)或長(zhǎng)時(shí)間失聯(lián)。4. MATLAB工程實(shí)現(xiàn)與調(diào)試心得4.1 數(shù)據(jù)輸入與模擬數(shù)據(jù)生成調(diào)試融合算法最怕沒(méi)有數(shù)據(jù)。項(xiàng)目早期如果拿不到真實(shí)雷達(dá)和AIS記錄可以先寫一個(gè)模擬數(shù)據(jù)生成器把傳感器的時(shí)序特征造出來(lái)。模擬器要還原三個(gè)關(guān)鍵特征雷達(dá)掃描周期短但噪聲大AIS更新周期長(zhǎng)且不一定規(guī)律兩個(gè)傳感器的時(shí)間戳不同步。function gen_demo() t 0:1:300; % 仿真時(shí)間軸單位秒 % 目標(biāo)勻速直線運(yùn)動(dòng)位置單位為米 x_true 1000 5 * t; y_true 2000 3 * t; % AIS每10秒上報(bào)一次定位噪聲小 aisT t(1:10:end); aisX x_true(1:10:end) randn(size(aisT)) * 10; aisY y_true(1:10:end) randn(size(aisT)) * 10; % 雷達(dá)每3秒掃描一次位置噪聲大 radT t(1:3:end); radX x_true(1:3:end) randn(size(radT)) * 30; radY y_true(1:3:end) randn(size(radT)) * 30; save(demo.mat, aisT,aisX,aisY,radT,radX,radY, x_true,y_true); end這個(gè)模擬器生成的數(shù)據(jù)足夠把關(guān)聯(lián)、融合、航跡管理整條鏈路跑通。建議模擬場(chǎng)景里再放一條近距離平行航行的干擾船專門用來(lái)驗(yàn)證關(guān)聯(lián)算法在多目標(biāo)場(chǎng)景下的表現(xiàn)。4.2 關(guān)鍵參數(shù)配置與調(diào)試順序我把這套融合系統(tǒng)里最關(guān)鍵的參數(shù)整理成了一張表初值是基于普通VTS監(jiān)控場(chǎng)景的經(jīng)驗(yàn)值實(shí)際使用時(shí)需要按具體數(shù)據(jù)調(diào)整。參數(shù)推薦初值說(shuō)明坐標(biāo)參考點(diǎn)監(jiān)控區(qū)域中心影響坐標(biāo)轉(zhuǎn)換精度時(shí)間對(duì)齊窗口2秒小于窗口直接用大于窗口外推位置關(guān)聯(lián)門限300至800米按σ組合計(jì)算再留機(jī)動(dòng)余量速度門限5節(jié)防止誤關(guān)聯(lián)航向門限30度防止大角度錯(cuò)誤配對(duì)AIS新鮮度閾值30秒超過(guò)閾值的AIS信息不參與融合航跡保持周期5個(gè)雷達(dá)周期防止軌跡頻繁斷裂航跡起始條件連續(xù)3次關(guān)聯(lián)成功防止雜波起批外推時(shí)長(zhǎng)上限30秒防止外推發(fā)散調(diào)試順序也很重要。我建議先從坐標(biāo)轉(zhuǎn)換入手畫一張俯視圖把雷達(dá)點(diǎn)、AIS點(diǎn)、轉(zhuǎn)換后的點(diǎn)疊在一起看確認(rèn)沒(méi)有系統(tǒng)性偏移然后調(diào)時(shí)間對(duì)齊觀察軌跡是否出現(xiàn)鋸齒接著調(diào)關(guān)聯(lián)門限確認(rèn)配對(duì)關(guān)系符合直覺(jué)最后才動(dòng)濾波器參數(shù)。不要一上來(lái)就調(diào)卡爾曼的Q和R問(wèn)題定位會(huì)很痛苦。4.3 可視化調(diào)試三板斧融合算法的調(diào)試可視化比看日志高效得多。我常用的有三招。第一招把所有原始點(diǎn)和融合軌跡畫在一張圖上。雷達(dá)點(diǎn)用圓圈、AIS點(diǎn)用三角、融合軌跡用連續(xù)的線不同目標(biāo)用不同顏色。一眼就能看出關(guān)聯(lián)錯(cuò)誤和位置偏移比對(duì)著數(shù)字猜快太多了。第二招打印關(guān)聯(lián)矩陣。每次關(guān)聯(lián)計(jì)算后把配對(duì)雙方的航跡編號(hào)、距離、速度差、航向差輸出到命令行人工檢查有沒(méi)有一對(duì)多、跳號(hào)、距離異常等狀況。第三招畫軌跡生命周期視圖。橫軸是時(shí)間縱軸是航跡號(hào)一條目標(biāo)的生命過(guò)程就是一條水平線斷裂和跳變?cè)谶@個(gè)視圖里十分直觀。目標(biāo)分裂、丟失、身份跳變的時(shí)刻都能一眼定位。這三招配合起來(lái)基本能把算法跑飛的原因鎖定到具體環(huán)節(jié)。5. 典型問(wèn)題與排查實(shí)錄5.1 五個(gè)高頻問(wèn)題與處理思路問(wèn)題之一多目標(biāo)交叉時(shí)身份亂跳。這是關(guān)聯(lián)門限設(shè)得過(guò)大造成的。我把位置門限從3σ收緊到2.5σ同時(shí)強(qiáng)制航向門限誤關(guān)聯(lián)率明顯下降。如果場(chǎng)景里經(jīng)常出現(xiàn)近距離并行船可以再加一個(gè)船長(zhǎng)門限做輔助判斷。問(wèn)題之二融合軌跡來(lái)回?cái)[動(dòng)。雷達(dá)點(diǎn)噪聲偏大加權(quán)系數(shù)沒(méi)有反映誤差變化。處理辦法是對(duì)雷達(dá)點(diǎn)先做一次α-β平滑再把平滑后的位置參與融合。α取0.7、β取0.1左右效果比較溫和不會(huì)把真實(shí)機(jī)動(dòng)削掉。問(wèn)題之三融合軌跡整體滯后。AIS報(bào)文的傳輸延遲會(huì)導(dǎo)致軌跡整體偏向歷史位置。短外推解決不了系統(tǒng)性的滯后我在項(xiàng)目里記錄AIS報(bào)文“發(fā)送時(shí)刻”與“接收時(shí)刻”的時(shí)間差估算固定延遲后做補(bǔ)償滯后問(wèn)題基本消除。問(wèn)題之四目標(biāo)短暫丟失后重新出現(xiàn)被當(dāng)成新目標(biāo)。這是航跡生命周期參數(shù)設(shè)置不合理。把航跡保持時(shí)間拉長(zhǎng)在丟失階段用AIS信息或運(yùn)動(dòng)外推維持航跡再用連續(xù)計(jì)數(shù)邏輯把“短暫消失”和“真正離開(kāi)”區(qū)分開(kāi)。問(wèn)題之五近岸雜波引起的假航跡。雷達(dá)在近岸會(huì)產(chǎn)生大量假回波單點(diǎn)起批很容易把雜波當(dāng)成船。對(duì)策是航跡起始必須滿足連續(xù)3次關(guān)聯(lián)成功并且加最小航速約束靜止雜波一般過(guò)不了這個(gè)條件。5.2 問(wèn)題排查速查表現(xiàn)象可能原因處理動(dòng)作身份跳變位置門限過(guò)大、特征約束弱收縮門限增加航速航向門限軌跡擺動(dòng)雷達(dá)噪聲大、加權(quán)失衡先平滑雷達(dá)點(diǎn)按誤差更新權(quán)重整體滯后AIS鏈路傳輸延遲估算延遲并補(bǔ)償優(yōu)先采用新鮮源軌跡分裂航跡保持周期太短加長(zhǎng)“丟失保持”時(shí)間放寬重新關(guān)聯(lián)條件誤匹配局部最近鄰沖突改用全局最優(yōu)分配排除一對(duì)多6. 寫在最后的一些體會(huì)融合算法做出來(lái)不難難的是在真實(shí)數(shù)據(jù)環(huán)境下一直不壞?；乜催@個(gè)項(xiàng)目我大部分時(shí)間其實(shí)花在數(shù)據(jù)清洗和參數(shù)調(diào)試上——AIS字段缺失、雷達(dá)野值、時(shí)間戳不齊這些問(wèn)題不處理干凈再漂亮的濾波公式輸出也是垃圾。另一個(gè)體會(huì)是調(diào)試時(shí)一定要保留可追溯的中間結(jié)果每次改動(dòng)參數(shù)都能復(fù)現(xiàn)問(wèn)題否則調(diào)參只能靠猜。這套基于MATLAB的實(shí)現(xiàn)思路覆蓋了普通VTS監(jiān)控場(chǎng)景的典型需求如果后續(xù)要做高動(dòng)態(tài)編隊(duì)或者復(fù)雜航道多目標(biāo)還需要引入更強(qiáng)的數(shù)據(jù)關(guān)聯(lián)方法比如概率數(shù)據(jù)關(guān)聯(lián)或帶航跡管理能力的多假設(shè)跟蹤。希望這份記錄能幫正在做同類工作的朋友少踩幾個(gè)坑。