删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于距离方程组的机载SAR图像立体定位方法

本站小编 Free考研考试/2021-12-25

闂傚倷娴囬褏鈧稈鏅犻、娆撳冀椤撶偟鐛ラ梺鍦劋椤ㄥ懐澹曟繝姘厵闁告挆鍛闂佹娊鏀遍崹鍫曞Φ閸曨垰绠涢柛鎾茬劍閸嬔冾渻閵堝繒鍒扮€殿喖澧庨幑銏犫攽鐎n亞鍔﹀銈嗗笒鐎氼剛绮婚妷锔轰簻闁哄啠鍋撻柛搴″暱閻g兘濡烽妷銏℃杸濡炪倖姊婚悺鏂库枔濡眹浜滈柨鏂垮⒔閵嗘姊婚崒姘偓鐑芥倿閿旈敮鍋撶粭娑樻噽閻瑩鏌熼悜姗嗘畷闁稿孩顨嗛妵鍕棘閸喒鎸冮梺鍛婎殕瀹€鎼佸箖濡も偓閳藉鈻庣€n剛绐楅梻浣哥-缁垰螞閸愵喖钃熸繛鎴欏灩鍞梺闈涚箚閸撴繈鎮甸敃鈧埞鎴︽倷閹绘帗鍊悗鍏夊亾闁归棿绀侀拑鐔兼煏閸繍妲哥紒鐙欏洦鐓曟い顓熷灥閺嬬喐绻涢崼婵堝煟婵﹨娅g槐鎺懳熼悡搴樻嫛闂備胶枪缁ㄦ椽宕愬Δ鍐ㄥ灊婵炲棙鍔曠欢鐐烘煙闁箑澧版い鏃€甯″娲嚃閳圭偓瀚涢梺鍛婃尰閻╊垶鐛繝鍌楁斀閻庯綆鍋嗛崢浠嬫⒑缂佹◤顏勵嚕閼搁潧绶為柛鏇ㄥ幐閸嬫挾鎲撮崟顒傤槰闂佹寧娲忛崹浠嬪极閹扮増鍊风痪鐗埫禍楣冩煥濠靛棝顎楀ù婊冨⒔缁辨帡骞夌€n剛袦闂佸搫鐬奸崰鎰缚韫囨柣鍋呴柛鎰ㄦ櫓閳ь剙绉撮—鍐Χ閸℃ê鏆楅梺纭呮珪閹瑰洦淇婇幘顔肩闁规惌鍘介崓鐢告⒑閹勭闁稿妫濇俊瀛樼節閸屾鏂€闂佺粯锕╅崑鍕妤e啯鈷戦柛娑橈功閳藉鏌f幊閸旀垵顕i弻銉晢闁告洦鍓欓埀顒€鐖奸弻锝夊箛椤撶偟绁烽梺鎶芥敱濮婅绌辨繝鍕勃闁稿本鑹鹃~鍥⒑閸濆嫮鐒跨紒缁樼箓閻i攱绺介崜鍙夋櫇闂侀潧绻掓慨瀵哥不閹殿喚纾介柛灞剧懅閸斿秵銇勯妸銉﹀殗閽樻繈姊婚崼鐔恒€掗柡鍡檮閹便劌顫滈崱妤€浼庣紓浣瑰敾缁蹭粙婀侀梺鎸庣箓鐎氼垶顢楅悢璁垮綊鎮℃惔銏犳灎濠殿喖锕ュ钘夌暦閵婏妇绡€闁稿本绮庨幊鍡樼節绾版ɑ顫婇柛瀣噽閹广垽宕奸妷褍绁﹂梺鍦濠㈡﹢鏌嬮崶顒佺厸闁搞儮鏅涢弸鎴炵箾閸涱喚澧紒缁樼⊕濞煎繘宕滆琚f繝鐢靛仜閹锋垹绱炴担鍝ユ殾闁炽儲鏋奸崼顏堟煕椤愩倕鏋庨柍褜鍓涢弫濠氬蓟閿濆顫呴柣妯哄悁缁敻姊洪幖鐐测偓鎰板磻閹剧粯鈷掑ù锝堫潐閸嬬娀鏌涢弬璺ㄐら柟骞垮灲瀹曠喖顢橀悙鑼喊闂佽崵濮村ú銈咁嚕椤掑嫬绫嶉柛灞绢殔娴滈箖鏌ㄥ┑鍡涱€楀褌鍗抽弻銊モ槈閾忣偄顏�
547闂傚倸鍊搁崐椋庣矆娴i潻鑰块梺顒€绉查埀顒€鍊圭粋鎺斺偓锝庝簽閿涙盯姊洪悷鏉库挃缂侇噮鍨堕崺娑㈠箳濡や胶鍘遍梺鍝勬处椤ㄥ棗鈻嶉崨瀛樼厽闊浄绲奸柇顖炴煛瀹€瀣埌閾绘牠鏌嶈閸撶喖寮绘繝鍥ㄦ櫜濠㈣泛锕﹂悿鍥⒑鐟欏嫬绀冩い鏇嗗懐鐭嗛柛鎰ㄦ杺娴滄粓鐓崶銊﹀鞍妞ゃ儲绮撻弻锝夊箻鐎靛憡鍒涘┑顔硷攻濡炶棄鐣峰Δ鍛闁兼祴鏅涢崵鎺楁⒒娴e憡鎲搁柛锝冨劦瀹曟垿宕熼娑樹患闂佺粯鍨兼慨銈夊疾閹间焦鐓ラ柣鏇炲€圭€氾拷1130缂傚倸鍊搁崐鎼佸磹妞嬪海鐭嗗〒姘e亾閽樻繃銇勯弽銊х煂闁活厽鎹囬弻锝夊閵忊晜姣岄梺绋款儐閹瑰洤鐣疯ぐ鎺濇晝闁挎繂娲﹂濠氭⒒娓氣偓閳ь剛鍋涢懟顖涙櫠閸欏浜滄い鎰╁焺濡叉椽鏌涢悩璇у伐妞ゆ挸鍚嬪鍕節閸愵厾鍙戦梻鍌欑窔閳ь剛鍋涢懟顖涙櫠閹绢喗鐓涢悘鐐登规晶鑼偓鍨緲鐎氼噣鍩€椤掑﹦绉靛ù婊勭矒閿濈偞鎯旈埦鈧弨浠嬫煟閹邦垰鐨哄褎鐩弻娑㈠Ω閵壯傝檸闂佷紮绲块崗姗€寮幘缁樺亹闁肩⒈鍓﹀Σ浼存煟閻斿摜鐭婄紒缁樺笧閸掓帒鈻庨幘宕囧€為梺鍐叉惈閸熶即鏁嶅⿰鍕瘈闁靛骏绲剧涵楣冩煥閺囶亪妾柡鍛劦濮婄粯鎷呴崨濠傛殘闁煎灕鍥ㄧ厱濠电姴鍟版晶杈╃磽閸屾稒宕岄柟绋匡攻缁旂喖鍩¢崒娑辨閻庤娲︽禍婵嬪箯閸涱垱鍠嗛柛鏇ㄥ幗琚欓梻鍌氬€风粈浣革耿闁秴鍌ㄧ憸鏃堝箖濞差亜惟闁宠桨鑳堕鍥⒑閸撴彃浜濇繛鍙夌墵閹偤宕归鐘辩盎闂佺懓顕崑娑㈩敋濠婂懐纾煎ù锝呮惈椤eジ鏌曢崶褍顏い銏℃礋婵偓闁宠桨绀佹竟澶愭⒒娴g懓顕滅紒瀣浮瀹曟繂鈻庨幘璺虹ウ闁诲函缍嗛崳顕€寮鍡欑瘈濠电姴鍊规刊鍏间繆閺屻儲鏁辩紒缁樼箞閹粙妫冨☉妤佸媰闂備焦鎮堕崝宀€绱炴繝鍌ゅ殨妞ゆ劑鍊楅惌娆愪繆椤愩倖鏆╅柛搴涘€楅幑銏犫攽鐎n亞鍊為梺闈浨归崕鏌ヮ敇濞差亝鈷戦柛婵嗗濡叉悂鏌eΔ浣虹煉鐎规洘鍨块獮鎺懳旈埀顒勫触瑜版帗鐓涢柛鎰╁妿婢ф盯鏌i幘宕囩闁哄本鐩崺鍕礃閳哄喚妲烽梻浣呵圭换鎰版儔閼测晜顫曢柟鐑橆殢閺佸﹪鏌涜箛鎿冩Ц濞存粓绠栭幃娲箳瀹ュ棛銈板銈庡亜椤︾敻鐛崱娑樻閹煎瓨鎸婚~宥夋⒑閸︻厼鍔嬮柛銊ㄦ珪缁旂喖寮撮悢铏诡啎闁哄鐗嗘晶浠嬪箖婵傚憡鐓涢柛婊€绀佹禍婊堝础闁秵鐓曟い鎰Т閸旀粓鏌i幘瀛樼闁哄瞼鍠栭幃婊兾熺拠鏌ョ€洪梻浣呵归鍥ㄧ箾閳ь剟鏌$仦鐣屝ユい褌绶氶弻娑滅疀閺冨倶鈧帗绻涢崱鎰仼妞ゎ偅绻勯幑鍕洪鍜冪船婵犲痉鏉库偓褏寰婃禒瀣柈妞ゆ牜鍋涚粻鐘虫叏濡顣抽柛瀣崌閻涱噣宕归鐓庮潛闂備礁鎽滈崰鎾寸箾閳ь剛鈧娲橀崹鍧楃嵁濡皷鍋撳☉娅亪顢撻幘缁樷拺缂備焦锚閻忥箓鏌ㄥ鑸电厓鐟滄粓宕滃☉銏犵;闁绘梻鍘ч悞鍨亜閹烘垵鏋ゆ繛鍏煎姍閺岀喖顢欓懖鈺佺厽閻庤娲樺ú鐔笺€佸☉銏″€烽柤纰卞墮婵附淇婇悙顏勨偓鏍垂婵傜ǹ纾垮┑鐘宠壘缁€鍌炴倶閻愭澘瀚庡ù婊勭矒閺岀喖骞嗚閹界娀鏌涙繝鍐ㄥ闁哄瞼鍠栭、娆撴嚃閳轰胶鍘介柣搴ゎ潐濞茬喐绂嶉崼鏇犲祦闁搞儺鍓欐儫闂侀潧顦崐鏇⑺夊顑芥斀闁绘劘鍩栬ぐ褏绱掗懠顒€浜剧紒鍌氱Ч閹崇偤濡疯濞村嫰姊洪幐搴㈢5闁稿鎹囧Λ浣瑰緞閹邦厾鍘遍棅顐㈡处濞叉牜鏁崼鏇熺厵闁稿繐鍚嬮崐鎰版煛鐏炵晫啸妞ぱ傜窔閺屾稖绠涢弮鍌楁闂傚洤顦甸弻娑㈠Ψ椤旂厧顫╃紒鐐劤閵堟悂寮婚弴鐔虹瘈闊洦娲滈弳鐘差渻閵堝棙绀夊瀛樻倐楠炲牓濡搁妷搴e枔缁瑩宕归纰辨綍闂傚倷鑳舵灙妞ゆ垵妫濋獮鎰節濮橆剛顔嗛梺鍛婁緱閸ㄩ亶宕伴崱娑欑厱闁哄洢鍔屾晶浼存煛閸℃ê鍝烘慨濠勭帛閹峰懘宕崟顐$帛闁诲孩顔栭崰妤呭磿婵傜ǹ桅闁圭増婢樼粈鍐┿亜韫囨挻顥犲璺哄娣囧﹪濡惰箛鏇炲煂闂佸摜鍣ラ崹璺虹暦閹达附鍋愮紓浣贯缚閸橀亶姊洪弬銉︽珔闁哥噥鍋呴幈銊╁焵椤掑嫭鈷戠紒瀣儥閸庢劙鏌熺粙娆剧吋妤犵偛绻樺畷銊р偓娑櫭禒鎯ь渻閵堝棛澧柤鐟板⒔缁骞嬮敂瑙f嫽婵炶揪绲介幉锟犲箚閸儲鐓曞┑鐘插閸︻厼寮查梻渚€娼х换鍫ュ磹閺囥垺鍊块柛顭戝亖娴滄粓鏌熺€电ǹ浠滄い鏇熺矌缁辨帗鎷呯憴鍕嚒濡炪値鍙€濞夋洟骞夐幘顔肩妞ゆ巻鍋撶痪鐐▕閹鈻撻崹顔界亾闂佽桨绀侀…鐑藉Υ娴g硶妲堟俊顖涚矌閸犲酣鎮鹃埄鍐跨矗濞达絽澹婂Λ婊勭節閻㈤潧浠╅柟娲讳簽缁辩偤鍩€椤掍降浜滄い鎰╁焺濡偓闂佽鍣换婵嬪春閳ь剚銇勯幒鎴濐仾闁抽攱甯¢弻娑氫沪閹规劕顥濋梺閫炲苯鍘哥紒顔界懇閵嗕礁鈻庨幇顔剧槇闂佸憡娲﹂崜锕€岣块悢鍏尖拺闁告挻褰冩禍婵囩箾閸欏澧辩紒顔垮吹缁辨帒螣闂€鎰泿闂備礁婀遍崑鎾翅缚濞嗘拲澶婎潩閼哥數鍘遍柣搴秵閸嬪懐浜告导瀛樼厵鐎瑰嫮澧楅崵鍥┾偓瑙勬礈閸忔﹢銆佸Ο琛℃敠闁诡垎鍌氼棜濠电姷鏁告慨鏉懨洪敃鍌氱9闁割煈鍋嗙粻楣冩煙鐎涙ḿ绠橀柡瀣暟缁辨帡鍩€椤掑倵鍋撻敐搴℃灍闁绘挸鍟伴幉绋库堪閸繄顦у┑鐐村灦濮樸劑鎯岄崱妞曞綊鏁愰崼鐔粹偓鍐煟閹烘埊韬柡宀€鍠庨埢鎾诲垂椤旂晫浜愰梻浣呵归鍡涘箰閹间礁鐓″璺哄閸嬫捇宕烽鐐愩儲銇勯敂鍨祮婵﹥妞介弻鍛存倷閼艰泛顏梺鍛娒幉锛勬崲濞戙垹绾ч柟瀵稿仜閺嬬姴顪冮妶鍐ㄧ仾闁挎洏鍨归悾鐑筋敃閿曗偓鍞悷婊冪灱缁厽寰勬繛鐐杸闁圭儤濞婂畷鎰板箻缂佹ê鈧潡鏌ㄩ弮鈧畷妯绘叏閾忣偅鍙忔俊顖氱仢閻撴劙鏌i幘宕囩闁哄本鐩崺鍕礃閳哄喚妲舵俊鐐€х拋锝嗕繆閸ヮ剙鐒垫い鎺嗗亾婵犫偓鏉堛劎浠氭俊鐐€ら崢濂稿床閺屻儲鍋╅柣鎴eГ閺呮煡鏌涢妷顖炴闁告洖鍟村铏圭矙閹稿孩鎷卞銈冨妼閹冲繒绮嬪澶婄畾妞ゎ兘鈧磭绉洪柡浣瑰姍瀹曘劑顢欓崗鍏肩暭闂傚倷绀侀幉鈥趁洪悢铏逛笉闁哄稁鍘奸拑鐔兼煥濠靛棭妲归柛濠勫厴閺屾稑鈻庤箛锝嗏枔濠碘槅鍋呴崹鍨潖濞差亝鐒婚柣鎰蔼鐎氫即鏌涘Ο缁樺€愰柡宀嬬秮楠炴帡鎮欓悽鍨闁诲孩顔栭崳顕€宕滈悢椋庢殾闁圭儤鍩堝ḿ鈺呮煥濠靛棙顥犻柛娆忓暞缁绘繂鈻撻崹顔界亾闂佺娅曢幐鍝ュ弲闂佺粯枪椤曆呭婵犳碍鐓欓柟顖嗗懏鎲兼繝娈垮灡閹告娊寮诲☉妯锋婵鐗婇弫楣冩⒑闂堚晝绋婚柟顔煎€垮濠氭晲閸℃ê鍔呴梺闈涚箳婵挳寮稿▎鎾寸厽闁绘ê鍟挎慨澶愭煕閻樺磭澧电€规洘妞介崺鈧い鎺嶉檷娴滄粓鏌熺€电ǹ浠滄い鏇熺矋閵囧嫰鏁冮崒銈嗩棖缂備浇椴搁幐鎼侇敇婵傜ǹ妞藉ù锝嚽规竟搴ㄦ⒒娴d警鏀版繛鍛礋閹囨偐鐠囪尙鐤囬梺缁樕戝鍧楀极閸℃稒鐓曢柟閭﹀枛娴滈箖鏌﹂幋婵愭Ш缂佽鲸鎹囧畷鎺戔枎閹存繂顬夐梻浣告啞閸旀洟鈥﹂悜鐣屽祦闊洦绋掗弲鎼佹煥閻曞倹瀚�28缂傚倸鍊搁崐鎼佸磹妞嬪孩顐介柨鐔哄Т绾捐顭块懜闈涘Е闁轰礁顑囬幉鎼佸籍閸稈鍋撴担鑲濇棃宕ㄩ闂寸盎闂備焦鍎崇换鎰耿闁秵鍋傞悗锝庡枟閳锋垿鎮峰▎蹇擃仾闁稿孩顨婇弻娑㈠Ω閵壯嶇礊婵犮垼顫夊ú鐔煎极閹剧粯鏅搁柨鐕傛嫹
潘志刚, 王超, 曹舸
中国科学院空天信息创新研究院, 北京 100190
2019年3月11日 收稿; 2019年11月27日 收修改稿
基金项目: 国家自然科学基金(61471340)资助
通信作者: 潘志刚,E-mail:zgpan@mail.ie.ac.cn

摘要: 提出一种基于距离方程组的机载SAR图像立体定位方法,该方法使用同一成像区域的3幅SAR图像进行无控制点定位,首先确定目标在3幅图像中对应的同名像点,解析该像点对应的雷达平台瞬时经纬度坐标、海拔高度以及目标到雷达平台的斜距值,在地心坐标系下建立3幅图像对应同名像点的距离方程组,对非线性方程组进行求解得到同名像点基于地心坐标系的坐标值,并通过坐标转换得到目标点的经纬度坐标及海拔高度,实现三维定位。对影响定位精度的系统误差因素进行仿真分析,利用实飞的机载SAR图像进行目标定位实验,结果验证了该方法的有效性。
关键词: 合成孔径雷达立体定位距离多普勒模型距离方程组
Stereo localization method for air-borne SAR image based on range equations
PAN Zhigang, WANG Chao, CAO Ge
Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100190, China


Abstract: A method of stereo localization for airborne SAR image based on range equations is proposed. This method uses three SAR images of the same imaging area to implement stereo localization without control points. Firstly, the corresponding image points (CIP) in the three images are determined, and the platform parameters corresponding to the CIP such as instantaneous geographic coordinates, altitude, and slant range are analyzed. The range equations for the CIP of three images are established in the geocentric coordinate system. The geocentric coordinates of the CIP are obtained by solving range equations, the geographic coordinates and altitude are obtained by coordinate transformation, so as to realize the three-dimensional positioning. The system error factors affecting the accuracy of stereo localization are simulated and analyzed. The positioning experiment is carried out using the actual flight airborne SAR image. The results show that the proposed method is effective.
Keywords: synthetic aperture radar(SAR)stereo localizationrange-Doppler modelrange equations
合成孔径雷达(synthetic aperture radar, SAR)具有全天时、全天候成像等优点,是高分辨率对地观测的重要工具;对目标进行准确定位是SAR的重要用途之一[1],而SAR影像立体定位是获取目标点三维信息的有效手段,对于成像区域的DEM生成、正射影像图制作,地形图绘制以及军事目标的精准解译等方面具有重要的应用价值。
距离-多普勒(range-Doppler)模型[2],简称R-D模型,在SAR图像立体定位中得到了广泛的应用。R-D模型基于SAR侧视斜距成像和多普勒频移两个特点,利用距离方程和多普勒频移方程建立SAR斜距图像的像点坐标与雷达平台坐标之间的函数关系,具有明确的几何物理意义[3];如果SAR满足正侧视成像条件,则R-D模型可以简化为多普勒中心频移为零的F.Leberl公式[4]
针对机载SAR图像立体定位问题,张红敏等[5-6]基于R-D模型设计利用单个地面控制点的SAR图像立体定位方案。除R-D模型之外,韦北余等[7]研究在无地面控制点的情况下,利用和差波束SAR图像直接对地面目标进行定位的方法,可获取地面目标的三维坐标信息。程春泉等[8]通过对SAR影像成像几何机理的分析,以传感器位置和姿态作为定向参数,构建了考虑多普勒参数的SAR影像距离-共面方程及基于该方程的SAR影像空中三角测量模型。
星载SAR图像定位目前普遍采用R-D模型,例如加拿大的RADARSAT-2雷达系统[9]、德国的TerraSAR-X雷达系统[10-12]以及中国自主研制的GF-3雷达卫星[13-14]。使用R-D模型进行定位时,影响定位精度的主要因素在于天线相位中心的位置和速度误差,此外SAR系统时间误差、大气传播延迟误差以及成像处理的引入误差也是影响SAR图像定位精度的重要因素;对于星载SAR系统而言,采用高精度的系统定标方法可显著降低上述系统性误差,例如标定后的GF-3卫星[14]雷达天线相位中心的位置精度可达到5 cm,速度精度为0.05 mm/s。降低卫星的系统误差对于提高图像定位精度至关重要。
与星载SAR相比,机载SAR的平台稳定性要差很多,尤其是小型的飞机平台在飞行时容易受惯导漂移、风速、风向以及空中气流的影响,使得飞行轨迹不再满足SAR成像时所需的匀速直线运动要求,因此机载SAR成像时需要做运动补偿以保证图像质量,但是运动补偿精度依赖于飞机平台的传感器精度;在没有安装POS系统时,飞机平台的位置及速度误差较大,如果单纯依赖飞机平台提供的实时数据去求解定位参数,则由于数据源存在较大误差且各参数之间存在相关性,对定位模型的参数解算提出了很高要求。
把R-D模型应用于机载SAR图像定位,模型中的多普勒频移条件方程说明了飞机平台的速度分量与雷达斜距分量的乘积是关于多普勒中心频率的函数;多普勒中心频率在成像处理过程中可估计,雷达斜距误差(距离延迟)可标定,因此飞机平台的位置误差和速度误差成为影响定位精度的主要因素。与雷达距离延迟不同,平台位置误差和速度误差受限于机载定位导航设备的精度,具有随机性和难以标定等缺点。本文使用成像区域内的3幅SAR图像进行立体定位,针对3幅图像中的同名像点建立距离方程组,在定位模型上使定位的误差因素局限于平台位置误差和雷达斜距误差,尽可能消除平台速度和多普勒频移参数估计等误差因素对定位精度的影响,有利于提高SAR图像无控制点立体定位精度。
1 立体定位方法1.1 定位基本原理假设机载SAR从3条不同的航线(非平行航线、相互交叉、夹角大于30°,如图 1所示)分别对目标区域进行条带模式成像,则目标位置坐标可由目标成像时刻雷达天线相位中心的位置坐标及其与目标之间的斜距距离确定,在地心坐标系下,目标点P的坐标与3条航线上对应的雷达天线相位中心坐标的几何关系可通过式(1)所示的距离方程组表达,求解式(1)所示方程组,即得到目标的三维坐标值,这一点同GPS定位原理类似。
$\left\{ {\begin{array}{*{20}{l}}{{f_1} = {{\left( {X - {X_1}} \right)}^2} + {{\left( {Y - {Y_1}} \right)}^2} + {{\left( {Z - {Z_1}} \right)}^2} - R_1^2 = 0, }\\{{f_2} = {{\left( {X - {X_2}} \right)}^2} + {{\left( {Y - {Y_2}} \right)}^2} + {{\left( {Z - {Z_2}} \right)}^2} - R_2^2 = 0, }\\{{f_3} = {{\left( {X - {X_3}} \right)}^2} + {{\left( {Y - {Y_3}} \right)}^2} + {{\left( {Z - {Z_3}} \right)}^2} - R_3^2 = 0.}\end{array}} \right.$ (1)
Fig. 1
Download: JPG
larger image
图 1 机载SAR通过3条航线对目标成像示意图 Fig. 1 Airborne SAR imaging through three air lines
图 1 机载SAR通过3条航线对目标成像示意图

Fig. 1 Airborne SAR imaging through three air lines -->

式中:(X, Y, Z)为目标点的三维地心坐标,(X1, Y1, Z1)为飞行航线1的目标像点所对应的雷达天线相位中心坐标,(X2, Y2, Z2)为飞行航线2的目标像点所对应的雷达天线相位中心坐标,(X3, Y3, Z3)为飞行航线3的目标像点所对应的雷达天线相位中心坐标,R1R2R3为雷达天线相位中心至目标之间的斜距距离。雷达天线相位中心的经纬度坐标和海拔高度可通过机载卫星导航定位系统如GPS或者北斗获取,并可转换到地心坐标系,雷达天线相位中心至目标之间的斜距距离可通过雷达斜距测量确定。
由式(1)可见,当使用3条交叉航线所得SAR图像对目标进行定位处理时,目标的坐标值仅由雷达天线相位中心的坐标及目标与雷达之间的斜距距离决定,而与其他因素无关。如果雷达天线相位中心坐标、目标与雷达之间的斜距距离能够精确测量,则目标定位精度足够精确。该方法与R-D定位模型相比,至少使用3幅图像实现目标的立体定位,而R-D模型至少需2幅图像即可实现,但是R-D模型除距离方程外,还需多普勒频移方程;在多普勒频移方程中,飞机平台的速度误差和多普勒中心频率估计误差会对定位精度产生较大影响。
1.2 方程组求解式(1)所示方程组是一个非线性方程组,求解之前需要进行线性化,即把式(1)所示方程组进行泰勒展开,并保留0阶和1阶导数项,得到如式(2)所示的线性化方程组:
$\left\{ {\begin{array}{*{20}{l}}{{f_1} = f_1^0 + \frac{{\partial {f_1}}}{{\partial X}}\left( {X - {X_0}} \right) + \frac{{\partial {f_1}}}{{\partial Y}}\left( {Y - {Y_0}} \right) + \frac{{\partial {f_1}}}{{\partial Z}}\left( {Z - {Z_0}} \right) = 0, }\\{{f_2} = f_2^0 + \frac{{\partial {f_2}}}{{\partial X}}\left( {X - {X_0}} \right) + \frac{{\partial {f_2}}}{{\partial Y}}\left( {Y - {Y_0}} \right) + \frac{{\partial {f_2}}}{{\partial Z}}\left( {Z - {Z_0}} \right) = 0, }\\{{f_3} = f_3^0 + \frac{{\partial {f_3}}}{{\partial X}}\left( {X - {X_0}} \right) + \frac{{\partial {f_3}}}{{\partial Y}}\left( {Y - {Y_0}} \right) + \frac{{\partial {f_3}}}{{\partial Z}}\left( {Z - {Z_0}} \right) = 0.}\end{array}} \right.$ (2)
式中,令ΔX=X-X0,ΔY=Y-Y0,ΔZ=Z-Z0,改写成矩阵形式如下:
$\left[ {\begin{array}{*{20}{l}}{\frac{{\partial {f_1}}}{{\partial X}}}&{\frac{{\partial {f_1}}}{{\partial Y}}}&{\frac{{\partial {f_1}}}{{\partial Z}}}\\{\frac{{\partial {f_2}}}{{\partial X}}}&{\frac{{\partial {f_2}}}{{\partial Y}}}&{\frac{{\partial {f_2}}}{{\partial Z}}}\\{\frac{{\partial {f_3}}}{{\partial X}}}&{\frac{{\partial {f_3}}}{{\partial Y}}}&{\frac{{\partial {f_3}}}{{\partial Z}}}\end{array}} \right] \times \left[ {\begin{array}{*{20}{l}}{\Delta X}\\{\Delta Y}\\{\Delta Z}\end{array}} \right] + \left[ {\begin{array}{*{20}{l}}{f_1^0}\\{f_2^0}\\{f_3^0}\end{array}} \right] = \left[ {\begin{array}{*{20}{l}}{{f_1}}\\{{f_2}}\\{{f_3}}\end{array}} \right].$ (3)
采用牛顿迭代法求解式(3)所示方程组,需先设定(X, Y, Z)的初始值为(X0, Y0, Z0),由于存在平台位置误差和距离误差,无法得到线性方程组的精确解,因此采用线性方程组的最小二乘解法求得改正量(ΔX, ΔY, ΔZ),则(X, Y, Z)的修正值为X′=X0XY′=Y0YZ′=Z0Z,当(ΔX, ΔY, ΔZ)较大时,可令当前的(X′, Y′, Z′)值作为新的初值,重复求解式(3)所示方程组,求出新的(ΔX, ΔY, ΔZ)与(X′, Y′, Z′),上述求解过程反复迭代,直到(ΔX, ΔY, ΔZ)满足所要求的精度为止。迭代初值设为(0, 0, 0)即可得到优化的结果。
1.3 定位实现流程图 2给出利用交叉航线得到的3幅SAR图像进行立体定位的实现流程。
Fig. 2
Download: JPG
larger image
图 2 立体定位流程图 Fig. 2 Flow chart of stereo positioning
图 2 立体定位流程图

Fig. 2 Flow chart of stereo positioning -->

1) 在3幅图像中选择定位目标点对应的同名像点,得到目标点在3幅图像中对应的像素坐标;
2) 分别针对3幅图像,解析目标点像素所在方位线对应的辅助数据,得到雷达平台的经度、纬度、海拔高度以及图像近距对应的雷达斜距值R0,由同名像点的像坐标(x, y)和距离向采样间隔Mx,求得同名像点与雷达平台的斜距距离R=R0+x·Mx
3) 将得到的雷达平台经度、纬度和高度值等转换到地心坐标系,并构建如式(1)所示的距离方程组;
4) 对方程组进行线性化,解线性方程组,得到最小二乘法解,即为目标点在地心坐标系下的坐标值;
5) 将目标点在地心坐标系下的值转换到地理坐标系,完成三维立体定位。
2 误差因素分析由式(1)的定位模型可以看出,采用3幅SAR图像进行立体定位时,影响SAR图像定位精度的直接因素是雷达平台的位置误差和雷达斜距误差(距离延迟);其中平台位置误差取决于机载导航定位系统的精度水平,可进一步细分为平台位置误差(经纬度误差)和平台高度误差;雷达斜距误差取决于雷达系统数据采集的时间精度以及大气传输延迟。下面分别通过实例仿真来分析平台位置误差、高度误差、雷达距离延迟等3个参数的测量误差对定位精度的影响,仿真使用的参数见论文第3节第1段。
图 3给出平台位置误差与目标定位精度的关系曲面,由于目标定位误差是由多种误差因素共同作用的结果,为了模拟平台位置误差对定位精度的影响,本文假设平台高度和雷达斜距不存在误差,相当于对影响定位精度的耦合误差因素做简化处理,后面的处理方式均与此类似。图 3中平台位置误差的范围为[-20, 20]m,可以看出平台位置误差不仅影响目标的平面精度,对目标高程精度也有较大的影响。
Fig. 3
Download: JPG
larger image
图 3 平台位置误差对定位精度的影响 Fig. 3 Influence of position error of platform on positioning accuracy
图 3 平台位置误差对定位精度的影响

Fig. 3 Influence of position error of platform on positioning accuracy -->

图 4给出平台高度误差与目标定位精度的关系曲线,这是假定平台位置误差和雷达斜距误差为0时所得到的结果。由图 4可见,平台高度误差对目标的平面精度和高程精度均产生影响,并且对目标平面精度的影响较大。
Fig. 4
Download: JPG
larger image
图 4 平台高度误差对定位精度的影响 Fig. 4 Influence of altitude error of platform on positioning accuracy
图 4 平台高度误差对定位精度的影响

Fig. 4 Influence of altitude error of platform on positioning accuracy -->

图 5给出雷达距离延迟与目标定位精度的关系曲线,这里假设平台位置误差和平台高度误差均为0。由图中可见,雷达距离延迟对目标平面精度的影响较小,当距离延迟在[-20, 20]m范围内变化时,目标平面误差最大值为1.1 m;与之相反的是,距离延迟对目标高程精度的影响甚为显著,[-20, 20]m范围内的斜距误差会导致最大72 m的目标高程误差。
Fig. 5
Download: JPG
larger image
图 5 雷达距离延迟对定位精度的影响 Fig. 5 Influence of range delay error of radar on positioning accuracy
图 5 雷达距离延迟对定位精度的影响

Fig. 5 Influence of range delay error of radar on positioning accuracy -->

综合图 3~图 5的仿真结果可以看出,对目标平面精度影响最大的误差因素是平台位置误差,对目标高程精度影响最大的因素是雷达距离延迟。值得一提的是,虽然经过标定后,雷达距离延迟可以很小,但是如果飞行过程中飞机姿态不稳定并且机载雷达位置坐标(由机载导航设备测量得到)存在较大误差,依然会造成SAR图像上的目标到雷达的斜距距离有不可忽略的误差值,从而造成目标高程精度降低。
3 试验结果与分析定位试验采用的SAR图像是由中国科学院电子学研究所研制的某试验型机载SAR系统获取的详查模式图像。雷达从3条不同的航线分别对同一目标区域成像,3条航线对应的航迹角分别为65.96°、128.83°和10.61°,雷达作用距离为近距10 km,飞机飞行高度3 100 m,飞行地速40 m/s,正侧视条带成像,所得3幅图像的有效重叠率为30%。雷达平台的实时位置坐标由GPS测量得到,该GPS没有采用差分定位,定位误差在2~10 m;经过标定后的雷达距离延迟小于1 m。
在3幅SAR图像中分别选择10个同名目标点用于定位测试,选择含有楼房、电线杆、平地角点等标志物的像素作为定位目标点,目标点的经纬度坐标和高程事先由GPS测量得到。
按照第1.3节所述的3幅图像联合立体定位的实现流程进行同名点选择、辅助数据解析和三维坐标解算,得到目标点的经纬度坐标和高程值,结果如表 1所示;表 2同时给出采用R-D模型[5]的立体定位结果,R-D定位模型采用2幅SAR图像,由于是正侧视成像,R-D模型中的多普勒中心频率参数设置为0。应保密要求,经纬度坐标的中间两位数字以“**”代替。
Table 1
表 1 本文所提方法的立体定位结果Table 1 Stereo positioning results of the proposed method
序号 测量经度/(°) 测量纬度/(°) 测量高程/m 解算经度/(°) 解算纬度/(°) 解算高程/m
1 112.2** 956 30.3** 992 30 112.2** 950 30.3** 967 29.87
2 112.2** 859 30.3** 036 60 112.2** 904 30.3** 203 75.61
3 112.2** 757 30.3** 193 40 112.2** 752 30.3** 212 43.57
4 112.2** 285 30.3** 744 30 112.2** 299 30.3** 794 42.47
5 112.2** 253 30.3** 420 48 112.2** 268 30.3** 414 58.30
6 112.2** 468 30.3** 615 28 112.2** 462 30.3** 635 37.75
7 112.2** 966 30.3** 831 28 112.2** 002 30.3** 846 35.89
8 112.2** 299 30.3** 669 29 112.2** 272 30.3** 704 45.15
9 112.2** 711 30.3** 330 30 112.2** 728 30.3** 302 22.57
10 112.2** 068 30.3** 712 30 112.2** 066 30.3** 719 39.65

表 1 本文所提方法的立体定位结果Table 1 Stereo positioning results of the proposed method


Table 2
表 2 采用R-D模型的立体定位结果Table 2 Stereo positioning results of the R-D model
序号 测量经度/(°) 测量纬度/(°) 测量高度/m 解算经度/(°) 解算纬度/(°) 解算高度/m
1 112.2** 956 30.3** 992 30 112.2** 959 30.3** 042 54.23
2 112.2** 859 30.3** 036 60 112.2** 902 30.3** 188 69.96
3 112.2** 757 30.3** 193 40 112.2** 755 30.3** 277 66.0
4 112.2** 285 30.3** 744 30 112.2** 294 30.3** 817 49.81
5 112.2** 253 30.3** 420 48 112.2** 264 30.3** 453 70.83
6 112.2** 468 30.3** 615 28 112.2** 444 30.3** 613 29.21
7 112.2** 966 30.3** 831 28 112.2** 989 30.3** 773 11.39
8 112.2** 299 30.3** 669 29 112.2** 26 30.3** 675 34.8
9 112.2** 711 30.3** 330 30 112.2** 728 30.3** 320 28.74
10 112.2** 068 30.3** 712 30 112.2** 058 30.3** 662 21.05

表 2 采用R-D模型的立体定位结果Table 2 Stereo positioning results of the R-D model

利用表 1表 2中的数值,可得到本文方法与R-D定位模型的立体定位误差曲线,分别如图 6(a)6 (b)所示。图 6(a)给出两种方法的目标平面精度比较,6(b)给出两种方法的目标高程精度比较。依据表 1表 2中数值,计算得到本文方法与R-D模型所得定位误差的均方根值(RMS)如表 3所示。由图 6表 3可以看出,本文方法的平面精度和高程精度总体上优于R-D模型。对于R-D模型来说,除雷达平台的位置误差、斜距误差外,还引入了飞机的速度误差以及多普勒中心频率估计误差。由于机载SAR成像运动补偿的影响,导致在图像中的某些点距离误差较大,影响了定位精度,使得图 6中的少数点采用R-D模型定位时精度占优。
Fig. 6
Download: JPG
larger image
图 6 两种模型定位精度的比较 Fig. 6 Comparison of positioning accuracy between two models
图 6 两种模型定位精度的比较

Fig. 6 Comparison of positioning accuracy between two models -->


Table 3
表 3 两种模型定位精度的RMS值比较Table 3 Comparison of RMS value of positioning accuracy between two models
定位方法 平面精度RMS/m 高程精度RMS/m
本文方法 3.74 10.41
R-D定位模型 5.96 16.33

表 3 两种模型定位精度的RMS值比较Table 3 Comparison of RMS value of positioning accuracy between two models

4 结束语R-D定位模型在星载SAR图像定位中得到了广泛应用并取得了很好的结果,机载SAR的平台稳定性远不如星载SAR,表现在速度、航迹的实时变化比较大。对于R-D模型来说,多普勒频移条件方程中的速度和多普勒中心频率参数在机载SAR图像定位时会成为影响定位精度的重要因素,在速度误差难以标定的前提下,如果只使用SAR辅助数据中的平台轨道参数而不采用地面控制点解算定向参数,则R-D模型针对机载SAR的无控制点定位误差较大,尤其是高程误差。
与R-D模型采用2幅SAR图像进行立体定位不同,本文使用成像区域不同航线所对应的3幅SAR图像进行联合立体定位,针对3幅图像中同一目标点所对应的同名像点建立距离方程组,多普勒频移条件方程则舍弃不用,在定位模型中仅需要知道平台的位置、高度和雷达的斜距距离即可完成定位解算。与R-D模型相比,本文方法可以减少定向参数并降低参数误差对定位精度的影响,代价为多使用1幅SAR图像参与定位解算。
本文采用3幅机载SAR条带模式图像进行立体定位的方法探索,并针对实飞的SAR条带图像进行定位试验,结果验证了方法的有效性。由于缺少聚束模式的试验数据,无法验证本文方法在聚束模式或者聚束与条带模式混合情况下的定位精度,这也是后续研究的方向。

参考文献
[1] 潘志刚, 潘卓, 曹舸. 机载SAR无控制点直接定位方法[J]. 中国科学院大学学报, 2015, 32(4): 536-541.
[2] Curlander J C. Location of pixels in space-borne SAR imagery[J]. IEEE Transactions on Geosciences and Remote Sensing, 1982, 20(3): 359-364.
[3] Johnsen H, Lauknes L, Guneriussen T. Geocoding of fast-delivery ERS-1 SAR image mode product using DEM data[J]. International Journal of Remote Sensing, 1995, 16(11): 1957-1968. Doi:10.1080/01431169508954532
[4] Leberl F W. Radargrammetric image processing[M]. Massachusetts: Artech House, 1990.
[5] 张红敏. SAR图像高精度定位技术研究[D].郑州: 信息工程大学, 2013.
[6] 张红敏, 靳国旺, 徐青, 等. 利用单个地面控制点的SAR图像高精度立体定位[J]. 雷达学报, 2014, 3(1): 85-91.
[7] 韦北余, 朱岱寅, 吴迪, 等. 一种基于和差波束的机载SAR定位方法[J]. 电子与信息学报, 2013, 35(6): 1464-1470.
[8] 程春泉, 张继贤, 黄国满, 等. 考虑多普勒参数的SAR影像距离-共面方程及其定位[J]. 遥感学报, 2013, 17(6): 1444-1450.
[9] Schubert A, Small D, Jehle M. COSMO-Skymed, TerraSAR-X, and RADARSAT-2 geolocation accuracy after compensation for earth-system effects[C]//2012 IEEE International Geoscience and Remote Sensing Symposium (IGARSS2012), Germany, 2012: 3301-3304.
[10] Gisinger C, Balss U, Pail R, et al. Precise three-dimensional stereo localization of corner reflectors and persistent scatters with TerraSAR-X[J]. IEEE Transactions on Geosciences and Remote Sensing, 2015, 53(4): 1782-1802. Doi:10.1109/TGRS.2014.2348859
[11] Eineder M, Minet C, Steigenberger P, et al. Imaging geodesy:toward centimeter-level ranging accuracy with TerraSAR-X[J]. IEEE Transactions on Geosciences and Remote Sensing, 2011, 49(2): 661-671. Doi:10.1109/TGRS.2010.2060264
[12] Cong X Y, Balss U, Eineder M. Imaging geodesy:centimeter-level ranging accuracy with TerraSAR-X:an update[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(5): 948-952. Doi:10.1109/LGRS.2012.2187042
[13] 云菲. 高分三号卫星[J]. 卫星应用, 2016(8): F0004.
[14] 丁赤飚, 刘佳音, 雷斌, 等. 高分三号SAR卫星系统级几何定位精度处探[J]. 雷达学报, 2017, 6(1): 11-16.


相关话题/图像 测量 系统 控制 图片

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于拉普拉斯人工势场的无人机避障控制
    顾育津1,2,3,宋孝成3,刘晓培3,陆疌31.中国科学院上海微系统与信息技术研究所,上海200050;2.中国科学院大学,北京100049;3.上海科技大学,上海2012102019年1月29日收稿;2019年4月15日收修改稿基金项目:国家自然科学基金(61603254)资助通信作者:顾育津摘要 ...
    本站小编 Free考研考试 2021-12-25
  • 一种基于差分算法的共形聚焦式微波热疗系统
    徐丽凡1,2,3,王雄11.上海科技大学,上海201210;2.中国科学院大学,北京100049;3.中国科学院上海微系统与信息技术研究所,上海2000502019年3月1日收稿;2019年5月6日收修改稿基金项目:国家自然科学基金(61701305),上海市浦江人才计划(17PJ1406600)和 ...
    本站小编 Free考研考试 2021-12-25
  • 基于净空技术的ISM频段可靠无线通信系统
    吉行健1,2,3,梁广2,3,孙思月2,3,姜泉江2,3,余金培2,31.中国科学院上海微系统与信息技术研究所,上海200050;2.中国科学院大学,北京100049;3.中国科学院微小卫星创新研究院,上海2012032018年12月27日收稿;2019年5月8日收修改稿基金项目:国家自然科学基金( ...
    本站小编 Free考研考试 2021-12-25
  • 基于仿真SAR图像深度迁移学习的自动目标识别
    王泽隆1,2,徐向辉1,张雷31.中国科学院电子学研究所,北京100190;2.中国科学院大学,北京100049;3.中国人民解放军95899部队,北京1000762019年1月29日收稿;2019年3月20日收修改稿基金项目:国家重点研发计划(2017YFB0503001)资助通信作者:王泽隆,E ...
    本站小编 Free考研考试 2021-12-25
  • 基于生成对抗网络的半监督遥感图像飞机检测
    陈国炜1,2,刘磊1,郭嘉逸1,2,潘宗序1,胡文龙11.中国科学院电子学研究所中国科学院空间信息处理与应用系统技术重点实验室,北京100190;2.中国科学院大学,北京1000492018年12月18日收稿;2019年3月27日收修改稿基金项目:国家自然科学基金(61701478)资助通信作者:刘 ...
    本站小编 Free考研考试 2021-12-25
  • 基于卷积神经网络与主动学习的高光谱图像分类
    宋晗,杨炜暾,耿修瑞,赵永超中国科学院电子学研究所,北京100190;中国科学院空间信息处理与应用系统技术重点实验室,北京100190;中国科学院大学,北京1000492018年12月20日收稿;2019年3月6日收修改稿基金项目:国家自然科学基金委重大科研仪器研制项目(41427805)和国防科工 ...
    本站小编 Free考研考试 2021-12-25
  • 自适应阈值收缩算子的稀疏正则化图像重建算法
    张胜男,许燕斌,董峰天津大学电气自动化与信息工程学院天津市过程检测与控制重点实验室,天津3000722019年3月12日收稿;2019年5月8日收修改稿基金项目:国家自然科学基金(61671322,61571321)和天津市自然科学基金(16JCYBJC18600)资助通信作者:许燕斌,E-mail ...
    本站小编 Free考研考试 2021-12-25
  • 基于傅里叶叠层成像的光学图像加密
    许文慧1,李拓1,史祎诗1,21.中国科学院大学,北京100049;2.中国科学院光电研究院,北京1000942015年09月06日收稿;2016年03月26日收修改稿基金项目:国家自然科学基金(61575197)资助通信作者:史祎诗?E-mail:shiyishi@ucas.ac.cn摘要:提出一 ...
    本站小编 Free考研考试 2021-12-25
  • 双随机相位加密系统的无约束最优化攻击
    王国华1,2,李拓1,3,张三国1,2,史祎诗1,31.中国科学院大学,北京100049;2.中国科学院大数据挖掘与知识管理重点实验室,北京100049;3.中国科学院光电研究院,北京1000942016年02月24日收稿;2016年03月31日收修改稿基金项目:国家自然科学基金(61575197) ...
    本站小编 Free考研考试 2021-12-25
  • 基于动态K均值聚类算法的SAR图像分割
    邢涛1,黄友红2,胡庆荣1,李军1,王冠勇11.中国航天二院二十三所,北京100854;2.中国人民解放军驻航天二院二十三所军代表室,北京1008542016年01月11日收稿;2016年04月05日收修改稿基金项目:国家自然科学基金(61271417)和高分专项青年创新基金(GFZX0406010 ...
    本站小编 Free考研考试 2021-12-25
婵犵數濮烽弫鎼佸磻閻愬搫鍨傞柛顐f礀缁犱即鏌熺紒銏犳灈缁炬儳顭烽弻鐔煎礈瑜忕敮娑㈡煟閹惧鈽夋い顓炴健閹虫粌顕ュΔ濠侀偗闁诡喗锕㈤幃鈺冪磼濡厧甯鹃梻浣稿閸嬪懐鎹㈤崟顖氭槬闁挎繂顦伴悡娆戔偓瑙勬礀濞层倝鍩㈤崼鈶╁亾鐟欏嫭绀冪紒顔肩Ч楠炲繘宕ㄩ弶鎴炲祶濡炪倖鎸鹃崰鎰邦敊韫囨稒鈷掗柛灞捐壘閳ь剙鍢查湁闁搞儺鐏涘☉銏犵妞ゆ劑鍊栧▓鎯ь渻閵堝棗鍧婇柛瀣尰閵囧嫰顢曢敐鍥╃杽婵犵鍓濋幃鍌炲春閳╁啯濯撮柧蹇曟嚀楠炩偓婵犵绱曢崑鎴﹀磹閺嶎厽鍋嬫俊銈呮噺閸嬶繝鏌曢崼婵囩┛濠殿喗濞婇弻鈩冨緞婵犲嫭鐨戝┑鈩冨絻閻楁捇寮婚敓鐘茬闁挎繂鎳嶆竟鏇熺節閻㈤潧袨闁搞劍妞介弫鍐閻樺灚娈鹃梺鍛婄箓鐎氼噣寮抽崱娑欑厱闁哄洢鍔屾晶顔界箾閸繄鐒告慨濠冩そ瀹曘劍绻濋崒姣挎洘绻涚€涙ḿ鐭岄柛瀣ㄥ€曢悾宄懊洪鍕紜闂佸搫鍊堕崕鏌ワ綖瀹ュ鈷戦悷娆忓閸斻倝鏌f幊閸斿孩绂嶉幖渚囨晝闁靛牆娲ㄩ敍婊冣攽鎺抽崐鏇㈠疮椤愶箑鍑犻柡鍐ㄧ墛閻撴瑥顪冪€n亪顎楅柍璇茬墛椤ㄣ儵鎮欓弶鎴犱紝濡ょ姷鍋涘ú顓€€佸▎鎾充紶闁告洦浜i崺鍛存⒒閸屾艾鈧绮堟笟鈧獮鏍敃閿曗偓绾惧湱鎲搁悧鍫濈瑲闁稿绻濆鍫曞醇濮橆厽鐝曞銈庡亝濞茬喖寮婚妸鈺傚亞闁稿本绋戦锟�
2婵犵數濮烽弫鎼佸磻閻愬搫鍨傞柛顐f礀缁犳壆绱掔€n偓绱╂繛宸簻鎯熼梺鍐叉惈椤戝洨绮欒箛娑欌拺闁革富鍘奸崝瀣亜閵娿儲顥㈢€规洜鏁婚崺鈧い鎺戝閳锋垿鏌涘☉姗堝伐濠殿噯绠戦湁婵犲﹤鎳庢禒杈┾偓瑙勬礃濡炰粙寮幘缁樺亹鐎规洖娲ら獮妤呮⒒娓氣偓濞佳呮崲閸儱纾归柡宓偓濡插牏鎲搁弮鍫濊摕闁挎繂顦悞娲煕閹板吀绨奸柛锝庡幘缁辨挻鎷呴崜鎻掑壈闂佹寧娲︽禍顏勵嚕椤愶箑纾奸柣鎰綑濞堟劙姊洪崘鍙夋儓闁哥姵鑹惧嵄闁告鍋愰弨浠嬫煃閽樺顥滃ù婊呭仜椤儻顦虫い銊ワ躬瀵偆鈧綆鍓涚壕钘壝归敐澶嬫锭濠殿喖鍊搁湁婵犲﹤妫楅悡鎰庨崶褝鍔熼柍褜鍓氱粙鎺曟懌婵犳鍨伴顓犳閹烘垟妲堟慨妤€妫楅崜杈╃磽閸屾氨孝闁挎洏鍎茬粚杈ㄧ節閸ヨ埖鏅濋梺闈涚墕閹峰寮抽銏♀拺闁告捁灏欓崢娑㈡煕閵娿儳鍩g€规洘妞介崺鈧い鎺嶉檷娴滄粓鏌熸潏鍓хɑ缁绢叀鍩栭妵鍕晜閼测晝鏆ら梺鍝勬湰缁嬫垿鍩㈡惔銈囩杸闁哄洨濯崬鍦磽閸屾瑧绐旂紓鍌涜壘铻為柛鏇ㄥ枤娴滄瑩姊绘担鍛婂暈婵炶绠撳畷銏c亹閹烘垹锛涢梺鍦劋椤ㄥ棝鍩涢幋锔界厱婵犻潧妫楅鈺呮煃瑜滈崜娆戠礊婵犲洤绠栭梺鍨儐缂嶅洭鏌嶉崫鍕簽婵炶偐鍠庨埞鎴︻敊鐟欐帞鎳撻埢鏂库槈閵忊€冲壒濠德板€愰崑鎾绘煃鐟欏嫬鐏撮柟顔规櫊楠炴捇骞掗崱妞惧闂佸綊妫跨粈渚€鏌ㄩ妶鍛斀闁绘ɑ褰冮弸銈嗙箾閸粎鐭欓柡宀嬬秮楠炲洭顢楁担鍙夌亞闂備焦鎮堕崐妤呭窗閹邦喗宕叉繝闈涱儏閻掑灚銇勯幒鎴濐仼闁绘帗妞介弻娑㈠箛椤栨稓銆婇梺娲诲幗椤ㄥ懘鍩為幋锔绘晩缂佹稑顑嗛悾鍫曟⒑缂佹﹩娈旂紒缁樺笧閸掓帡宕奸悢椋庣獮闁诲函缍嗛崜娑㈩敊閺囥垺鈷戦柣鐔煎亰閸ょ喎鈹戦鐐毈鐎殿喗濞婇崺锟犲磼濠婂拋鍟庨梺鑽ゅТ濞壯囧礋椤愵偂绱�547闂傚倸鍊搁崐椋庣矆娴i潻鑰块梺顒€绉查埀顒€鍊圭粋鎺斺偓锝庝簽閿涙盯姊洪悷鏉库挃缂侇噮鍨堕崺娑㈠箳濡や胶鍘遍梺鍝勬处椤ㄥ棗鈻嶉崨瀛樼厽闊浄绲奸柇顖炴煛瀹€瀣埌閾绘牠鎮楅敐搴′簻妞ゅ骏鎷�4婵犵數濮烽弫鎼佸磻閻愬搫鍨傞柛顐f礀缁犳壆绱掔€n偓绱╂繛宸簼閺呮煡鏌涢妷銏℃珖妞わ富鍨跺娲偡闁箑娈堕梺绋款儑閸犳牠宕洪姀銈呯睄闁逞屽墴婵$敻宕熼鍓ф澑闂佽鍎抽顓⑺囬柆宥嗏拺缂佸顑欓崕鎰版煙閻熺増鎼愰柣锝呭槻椤粓鍩€椤掑嫨鈧線寮崼婵嗚€垮┑掳鍊曢崯顐︾嵁閹扮増鈷掗柛灞剧懅椤︼箓鏌涘顒夊剰妞ゎ厼鐏濋~婊堝焵椤掆偓閻g兘顢涢悜鍡樻櫇闂侀潧绻堥崹鍝勨枔妤e啯鈷戦梻鍫熶緱濡狙冣攽閳ヨ櫕鍠橀柛鈹垮灲瀵噣宕奸悢鍝勫箥闂備胶顢婇~澶愬礉閺囥垺鍎嶆繛宸簼閻撶喖鏌i弮鍫熸暠閻犳劧绱曠槐鎺撴綇閵娿儳鐟查悗鍨緲鐎氼噣鍩€椤掑﹦绉靛ù婊呭仦缁傛帡鎮℃惔妯绘杸闂佺粯鍔樺▔娑氭閿曞倹鐓曟俊銈呭閻濐亜菐閸パ嶅姛闁逞屽墯缁嬫帟鎽繝娈垮灡閹告娊骞冨畡鎵虫瀻婵炲棙鍨甸崺灞剧箾鐎涙ḿ鐭掔紒鐘崇墵瀵鈽夐姀鐘电杸闂佺ǹ绻愰幗婊堝极閺嶎厽鈷戠紒顖涙礃濞呮梻绱掔紒妯肩疄鐎殿喛顕ч埥澶娾堪閸涱垱婢戦梻浣瑰缁诲倿骞婃惔顭掔稏闁冲搫鎳忛埛鎴︽煕濞戞﹫鍔熼柟铏礈缁辨帗娼忛妸锔绢槹濡ょ姷鍋涚换姗€骞冮埡鍐╁珰闁肩⒈鍓﹂崯瀣⒒娴e憡鍟炲〒姘殜瀹曞綊骞庨崜鍨喘閸╋繝宕ㄩ瑙勫闂佽崵鍋炵粙鍫ュ焵椤掆偓閸樻牗绔熼弴銏♀拻濞达絽鎲$拹锟犲几椤忓棛纾奸柕濞垮妼娴滃湱绱掗鍛箺鐎垫澘瀚伴獮鍥敇閻樻彃绠婚梻鍌欑閹碱偆鈧凹鍓涢幑銏ゅ箳閺冨洤小闂佸湱枪缁ㄧ儤绂嶅⿰鍫熺厸闁搞儺鐓侀鍫熷€堕柤纰卞厴閸嬫挸鈻撻崹顔界彯闂佺ǹ顑呴敃銈夘敋閿濆洦宕夐悶娑掑墲閻庡姊虹拠鈥崇€婚柛蹇庡嫎閸婃繂顫忕紒妯诲闁荤喖鍋婇崵瀣磽娴e壊鍎愰柛銊ㄥ劵濡喎顪冮妶鍡樺蔼闁搞劌缍婇幃鐐哄垂椤愮姳绨婚梺鍦劋閸╁﹪寮ㄦ繝姘€垫慨妯煎亾鐎氾拷40缂傚倸鍊搁崐鎼佸磹妞嬪海鐭嗗〒姘e亾閽樻繃銇勯弽銊х煂闁活厽鎸鹃埀顒冾潐濞叉牕煤閵娧呬笉闁哄啫鐗婇悡娆撴煙椤栧棗鑻▓鍫曟⒑瀹曞洨甯涙慨濠傜秺楠炲牓濡搁妷顔藉缓闂侀€炲苯澧版繛鎴犳暬楠炴牗鎷呴崨濠勨偓顒勬煟鎼搭垳绉靛ù婊冪埣閹垽宕卞☉娆忎化闂佹悶鍎荤徊娲磻閹捐绀傞柛娑卞弾濡粎绱撻崒姘偓宄懊归崶銊d粓闁归棿鐒﹂崑锟犳煃閸濆嫭鍣归柦鍐枔閳ь剙鍘滈崑鎾绘煕閺囥劌浜炴い鎾存そ濮婃椽骞愭惔锝囩暤濠电偠灏欐繛鈧€规洘鍨块獮妯肩磼濡鍔掗梺鑽ゅ枑閻熴儳鈧凹鍓熷畷銏c亹閹烘挴鎷洪梺鍛婄箓鐎氼厼顔忓┑瀣厱閹兼番鍨归悘鈺備繆閸欏濮囨顏冨嵆瀹曞ジ鎮㈤崫鍕闂傚倷鑳剁涵鍫曞礈濠靛枹娲冀椤愩儱小缂備緡鍋勭€殿剟姊婚崒姘偓椋庢濮橆兗缂氱憸宥堢亱闂佸搫鍟崐濠氭儗閸℃褰掓晲閸偄娈欓梺鑽ゅ枑鐎氬牓寮崼婵嗙獩濡炪倖妫侀~澶屸偓鍨墵濮婄粯鎷呴崨濠傛殘婵炴挻纰嶉〃濠傜暦閵忋倖瀵犲璺烘閻庢椽鎮楅崗澶婁壕闂佸憡娲﹂崜娑㈠储闁秵鈷戦柛婵嗗閺嗙偤鏌熺粙鍨挃濠㈣娲熼獮鎰償濞戞鐩庨梻渚€娼ф蹇曟閺団偓鈧倿鎳犻鍌滐紲闂佸搫鍟崐鎼佸几濞戞瑣浜滈柕蹇婂墲缁€瀣煙椤旇娅婃い銏℃礋閿濈偤顢橀悜鍡橆棥濠电姷鏁搁崑鐘诲箵椤忓棛绀婇柍褜鍓氶妵鍕敃閵忊晜鈻堥梺璇″櫙缁绘繈宕洪埀顒併亜閹烘垵顏柍閿嬪浮閺屾稓浠﹂幑鎰棟闂侀€炲苯鍘哥紒顔界懇閵嗕礁鈻庨幇顔剧槇闂佸憡娲﹂崜锕€岣块悢鍏尖拺闁告挻褰冩禍婵囩箾閸欏澧辩紒顔垮吹缁辨帒螣闂€鎰泿闂備浇顫夊畷妯衡枖濞戙埄鏁佺€光偓閸曨剛鍘告繛杈剧到婢瑰﹪宕曡箛鏂讳簻妞ゆ挴鍓濈涵鍫曟煙妞嬪骸鈻堥柛銊╃畺瀹曟宕ㄩ娑樼樆闂傚倸鍊风欢姘跺焵椤掑倸浠滈柤娲诲灦瀹曘垽骞栨担鍦幘闂佸憡鍔樼亸娆撳春閿濆應鏀介柨娑樺閺嗩剟鏌熼鐣屾噰鐎殿喖鐖奸獮瀣敇閻愭惌鍟屾繝鐢靛У椤旀牠宕板Δ鍛櫇闁冲搫鎳庣粈鍌涚箾閹寸偟顣叉い顐f礋閺屻劌鈹戦崱妯轰痪閻熸粎澧楃敮妤呭疾閺屻儲鐓曢柍鈺佸暟閹冲懘鏌i幘鍐测偓鎼佲€旈崘顔嘉ч柛鎰╁妿娴犲墽绱掗悙顒佺凡缂佸澧庨崚鎺楀煛閸涱喖浜滅紒鐐妞存悂寮插┑瀣拺闂傚牊绋撴晶鏇熺箾鐠囇呯暤妤犵偛妫濋弫鎰緞鐎Q勫闂備礁婀辨灙婵炲鍏橀崺銉﹀緞鐎c劋绨婚梺鎸庢椤曆冾嚕椤曗偓閺屾盯鍩為幆褌澹曞┑锛勫亼閸婃牜鏁幒妤佹櫇闁靛/鈧崑鎾愁潩閻愵剙顏�28缂傚倸鍊搁崐鎼佸磹妞嬪孩顐介柨鐔哄Т绾捐顭块懜闈涘Е闁轰礁顑囬幉鎼佸籍閸垹绁﹂梺鍛婂姦閸犳牜绮绘繝姘厱闁规崘灏欑粣鏃堟煃閻熸壆绠茬紒缁樼箞婵偓闁挎繂妫涢妴鎰斿Δ濠佺凹闁圭ǹ鍟块悾宄扳攽鐎n亜绐涢柣搴㈢⊕宀e潡宕㈤柆宥嗏拺闁告繂瀚弳濠囨煕鐎n偅灏电紒杈ㄥ笧閳ь剨缍嗛崑鍛暦瀹€鈧埀顒侇問閸n噣宕戞繝鍥х畺濞寸姴顑呴崹鍌涖亜閹扳晛鐏╂鐐村灴濮婄粯鎷呴崨濠冨創濠电偠顕滅粻鎴︼綖濠靛惟闁冲搫鍊告禒顓㈡⒑鐎圭姵銆冮悹浣瑰絻鍗遍柛顐犲劜閻撴瑩鏌i幇闈涘缂傚秵鍨块弻鐔煎礂閸忕厧鈧劙鏌$仦鐣屝ユい褌绶氶弻娑㈠箻閸楃偛顫囧Δ鐘靛仜缁绘﹢寮幘缁樻櫢闁跨噦鎷�1130缂傚倸鍊搁崐鎼佸磹妞嬪海鐭嗗〒姘e亾閽樻繃銇勯弽銊х煂闁活厽鎹囬弻娑㈠箻閼碱剦妲梺鎼炲妽缁诲牓寮婚妸鈺傚亜闁告繂瀚呴姀銏㈢<闁逞屽墴瀹曟帡鎮欑€电ǹ骞堟繝鐢靛仦閸ㄥ爼鏁冮锕€缁╃紓浣贯缚缁犻箖鏌涢锝囩畼闁绘帗鎮傞弻锛勪沪缁嬪灝鈷夐悗鍨緲鐎氼噣鍩€椤掑﹦绉靛ù婊勭矒閿濈偤宕堕浣叉嫼闂備緡鍋嗛崑娑㈡嚐椤栨稒娅犲Δ锝呭暞閻撴瑩鏌涢幋娆忊偓鏍偓姘炬嫹