本發(fā)明涉及衛(wèi)星超大幅寬擺掃成像方法,特別涉及一種衛(wèi)星快速旋轉(zhuǎn)超大幅寬擺掃成像方法。
背景技術(shù):
衛(wèi)星的推掃成像一般是通過將探測器安置在垂直于衛(wèi)星的飛行方向上,衛(wèi)星向前飛行的時候,一次收集一行圖像。圖像的范圍大小依據(jù)探測器中ccd的視場范圍而定,一般采用多個ccd拼接或者同時增大相機(jī)側(cè)擺能力可實(shí)現(xiàn)更大范圍的視場。這是目前較通用的光學(xué)遙感衛(wèi)星成像方式。衛(wèi)星的擺掃成像是使用反射鏡將光線反射到探測器內(nèi),利用反射鏡的往返擺動收集一個像素上的測量值。這種成像方式移動部件非常昂貴且容易損壞,視場范圍有限。
本發(fā)明在常規(guī)衛(wèi)星推掃和擺掃成像的基礎(chǔ)上,提出一種利用衛(wèi)星快速旋轉(zhuǎn)實(shí)現(xiàn)超大幅寬擺掃成像的方法。線陣ccd與衛(wèi)星飛行方向平行,通過探測器360°連續(xù)自旋,同時實(shí)現(xiàn)擺掃和推掃,從而實(shí)現(xiàn)衛(wèi)星對地成像覆蓋范圍的最大化。
現(xiàn)階段常規(guī)衛(wèi)星均為靜態(tài)推掃成像或大角度動態(tài)推掃成像,專利號cn201410151609.9“一種基于增量式偏流角的動中成像衛(wèi)星姿態(tài)控制方法”,涉及一種動中成像衛(wèi)星姿態(tài)控制方法,選取衛(wèi)星當(dāng)前目標(biāo)姿態(tài)為參考基準(zhǔn),求解偏流角的增量,并將更新后的參考姿態(tài)矩陣作為姿態(tài)控制時的目標(biāo)姿態(tài)矩陣,使得衛(wèi)星姿態(tài)能夠跟蹤上目標(biāo)姿態(tài)。專利號cn201410163903.1“一種新型敏捷衛(wèi)星機(jī)動中成像方法”,該方法在衛(wèi)星模型中建立ccd像平面,通過投影計算得到像移速度矢量和偏流角,控制衛(wèi)星的偏航角進(jìn)行偏流角的修正,滿足機(jī)動中的成像要求,實(shí)現(xiàn)衛(wèi)星在三軸姿態(tài)機(jī)動過程中開啟光學(xué)有效載荷進(jìn)行成像的動態(tài)成像。專利號cn201510466057.5“反射鏡擺動寬幅成像系統(tǒng)及成像方法”,提供了一種反射鏡擺動寬幅成像系統(tǒng),通過擺動反射鏡的轉(zhuǎn)動擴(kuò)大可見光面陣成像相機(jī)的成像視場,實(shí)現(xiàn)可見光面陣成像相機(jī)對星下點(diǎn)的寬度方向進(jìn)行多角度拍照。
上述技術(shù)方案只能針對衛(wèi)星均為靜態(tài)推掃成像、大角度動態(tài)推掃成像或衛(wèi)星靜態(tài)反射鏡擺動成像,無法滿足地面區(qū)域超大幅寬成像,不能解決單個衛(wèi)星相鄰兩軌成像區(qū)域之間的無縫拼接成像問題。
技術(shù)實(shí)現(xiàn)要素:
本發(fā)明的目的是為了解決只能針對衛(wèi)星均為靜態(tài)推掃成像、大角度動態(tài)推掃成像或衛(wèi)星靜態(tài)反射鏡擺動成像的方法,無法滿足地面區(qū)域超大幅寬成像,不能解決單個衛(wèi)星相鄰兩軌成像區(qū)域之間的無縫拼接成像問題,而提出的一種衛(wèi)星快速旋轉(zhuǎn)超大幅寬擺掃成像方法。
上述的發(fā)明目的是通過以下技術(shù)方案實(shí)現(xiàn)的:
步驟一、假設(shè)衛(wèi)星軌道高度h,取地球半徑r,計算圓心角2θ所對應(yīng)的衛(wèi)星軌道的曲線上的ab兩點(diǎn)之間的曲線距離
步驟二、假設(shè)探測器視場角為η,則根據(jù)視場角和軌道高度計算探測器的視場范圍l1即為飛行方向的幅寬;
步驟三、兩次成像時探測器光軸中心在地表的星下點(diǎn)軌跡上的距離l2必須小于等于飛行方向的幅寬l1,l2的臨界值是飛行方向的幅寬l1;
步驟四、取臨界值l2=l1帶入如下公式則計算出探測器沿軌方向的自旋速度:
探測器繞沿軌方向的自旋速度大于等于vz時,滿足探測器的兩次相鄰成像區(qū)域之間沒有縫隙,如果低于vz,則兩次成像區(qū)域之間會有縫隙;
步驟五、結(jié)合探測器固有的特征參數(shù)選擇不同的側(cè)擺角度計算各種軌道條件下探測器實(shí)現(xiàn)超大幅寬成像時對應(yīng)的分辨率ac;其中,探測器固有的特征參數(shù)包括像元尺寸和焦距;
步驟六、根據(jù)探測器自旋速度vz及分辨率ac計算ccd行頻fp:
步驟七、已知衛(wèi)星軌道高度為h的太陽同步軌道,衛(wèi)星繞地球一圈為t,則一天衛(wèi)星繞地球約q圈,赤道周長為l,則在赤道處,相鄰兩軌之間的距離為l3=l/q,已知軌道傾角為σ,則兩軌之間的直線距離為l4=l3×sinσ,由計算已知垂直于軌道的幅寬l5,則計算相鄰兩軌成像區(qū)域之間的疊加距離為l6=l5-l4。
發(fā)明效果
本發(fā)明在常規(guī)衛(wèi)星推掃和擺掃成像的基礎(chǔ)上,提出一種利用衛(wèi)星快速旋轉(zhuǎn)實(shí)現(xiàn)超大幅寬擺掃成像的方法。線陣ccd與衛(wèi)星飛行方向平行,通過衛(wèi)星360°連續(xù)自旋,衛(wèi)星快速旋轉(zhuǎn)速度1~10°/s;同時實(shí)現(xiàn)擺掃和推掃,實(shí)現(xiàn)衛(wèi)星對地成像覆蓋范圍的最大化從而實(shí)現(xiàn)中低軌道衛(wèi)星千公里級的超大幅寬擺掃成像。
傳統(tǒng)的推掃、擺掃衛(wèi)星成像,如果要對某一區(qū)域完全覆蓋,需要多次對該區(qū)域重訪,即多軌數(shù)據(jù)拼接;雖然加上側(cè)擺能力可以加大對區(qū)域的覆蓋能力,但是一次性采集地面數(shù)據(jù)的時間分辨率和覆蓋能力都不能滿足當(dāng)前衛(wèi)星快速響應(yīng)的任務(wù)需求。
本發(fā)明綜合考慮了衛(wèi)星的機(jī)動能力和成像能力,依據(jù)衛(wèi)星動態(tài)成像的姿態(tài)角和姿態(tài)角速度的變化及動態(tài)變化對成像的影響,設(shè)計了一種在衛(wèi)星飛行過程中,探測器光軸隨著衛(wèi)星的前進(jìn)方向旋轉(zhuǎn),在沿軌方向或非沿軌方向成像,從而實(shí)現(xiàn)一次性大范圍覆蓋成像的方法。本發(fā)明公開了一種衛(wèi)星快速旋轉(zhuǎn)超大幅寬擺掃成像的方法,該方法針對傳統(tǒng)遙感衛(wèi)星上采用的推掃、擺掃成像方法中幅寬有限的問題,在衛(wèi)星飛行過程中,使探測器光軸隨著衛(wèi)星的前進(jìn)緩慢的向前;隨著衛(wèi)星的360°勻速旋轉(zhuǎn),向左或右推移。使沿軌方向或非沿軌方向成像幅寬增大,以實(shí)現(xiàn)大范圍覆蓋成像的目的。本發(fā)明方法簡單可靠、易行,擴(kuò)大了推擺掃成像范圍,擴(kuò)展了推擺掃衛(wèi)星的應(yīng)用領(lǐng)域。本發(fā)明是國內(nèi)首個針對擴(kuò)大成像范圍提出的衛(wèi)星快速旋轉(zhuǎn)超大幅寬擺掃成像及大覆蓋面積的可靠方法。
探測器的旋轉(zhuǎn)可以通過探測器相對于衛(wèi)星平臺旋轉(zhuǎn)或者隨整星結(jié)構(gòu)整體旋轉(zhuǎn)實(shí)現(xiàn)。
為實(shí)現(xiàn)衛(wèi)星進(jìn)動和地球自轉(zhuǎn)情況下衛(wèi)星對地面大區(qū)域探測,衛(wèi)星姿態(tài)控制系統(tǒng)需要保證探測器360°連續(xù)勻速旋轉(zhuǎn),成像時衛(wèi)星姿態(tài)控制系統(tǒng)需要對衛(wèi)星本體坐標(biāo)系相對于軌道坐標(biāo)系的三軸姿態(tài)角(橫滾
附圖說明
圖1為具體實(shí)施方式一提出的衛(wèi)星自旋飛行示意圖;
圖2為具體實(shí)施方式一提出的探測器覆蓋區(qū)域示意圖。
圖3為具體實(shí)施方式一提出的垂軌幅寬計算示意圖。
圖4為具體實(shí)施方式一提出的沿軌幅寬計算示意圖。
圖5為具體實(shí)施方式一提出的一定側(cè)擺角成像對應(yīng)的分辨率下降示意圖。
圖6為具體實(shí)施方式一提出的衛(wèi)星側(cè)擺角與地面成像分辨率對應(yīng)關(guān)系。
圖7為實(shí)施例提出的兩軌拼接示意圖。
具體實(shí)施方式
具體實(shí)施方式一:本實(shí)施方式的一種衛(wèi)星快速旋轉(zhuǎn)超大幅寬推掃成像方法,具體是按照以下步驟制備的:
步驟一、假設(shè)衛(wèi)星軌道高度h,取地球半徑r,計算圓心角2θ所對應(yīng)的衛(wèi)星軌道的曲線上的ab兩點(diǎn)之間的曲線距離
步驟二、假設(shè)探測器視場角為η,則根據(jù)視場角和軌道高度計算探測器的視場范圍l1即為飛行方向的幅寬如圖4所示;
步驟三、只有探測器兩次相鄰的成像區(qū)域之間無縫隙才能實(shí)現(xiàn)超大幅寬,成像區(qū)域示意圖見圖2,即兩次成像時探測器光軸中心在地表的星下點(diǎn)軌跡上的距離l2必須小于等于飛行方向的幅寬l1,因此,可認(rèn)為l2的臨界值是飛行方向的幅寬l1;
步驟四、取臨界值l2=l1帶入如下公式則計算出探測器沿軌方向的自旋速度:
探測器繞沿軌方向的自旋速度大于等于vz時,可以滿足探測器的兩次相鄰成像區(qū)域之間沒有縫隙,如果低于vz,則兩次成像區(qū)域之間會有縫隙,不能實(shí)現(xiàn)超大幅寬成像;
步驟五、探測器的地表分辨率與其側(cè)擺角度大小密切相關(guān),一種衛(wèi)星快速旋轉(zhuǎn)超大幅寬擺掃成像的方法,通過探測器的360°旋轉(zhuǎn)能夠提供大范圍的側(cè)擺角,但并不是隨著側(cè)擺角的增大探測器的分辨率越高,只是在一定的側(cè)擺角范圍內(nèi),相同的軌道條件下,側(cè)擺角越大探測器的分辨率越高;為兼顧高分辨率和地面條帶無縫拼接兩方面的要求,需要成像幅寬與分辨率ac相互制約;
這是因為衛(wèi)星側(cè)擺角度很大時,衛(wèi)星與地球表面相切程度增大,除了衛(wèi)星分辨率隨著側(cè)擺物距的變化,衛(wèi)星在地面上的投影角度也變大,對應(yīng)成像的幾何分辨率變形嚴(yán)重,會出現(xiàn)嚴(yán)重的分辨率下降,如圖5所示;此時,結(jié)合探測器固有的特征參數(shù)選擇不同的側(cè)擺角度計算各種軌道條件下探測器實(shí)現(xiàn)超大幅寬成像時對應(yīng)的分辨率ac;衛(wèi)星側(cè)擺角與分辨率對應(yīng)關(guān)系如圖6所示;其中,探測器固有的特征參數(shù)包括像元尺寸和焦距;
步驟六、探測器的轉(zhuǎn)速vz及分辨率ac也決定了ccd行頻fp的大小,而ccd器件均具有自身的極限行頻,超出該極限行頻則探測器無法正常成像;根據(jù)探測器自旋速度vz及分辨率ac可以計算ccd行頻fp;
步驟七、如圖2所示,探測器自旋成像時,成像區(qū)域兩側(cè)變形較大,考慮兩軌之間即兩次成像區(qū)域之間的距離,可以減小成像區(qū)域在垂直沿軌方向的距離l6,即兩軌之間重疊區(qū)域,使用兩軌數(shù)據(jù)拼接的方式;如圖7所示;
已知衛(wèi)星軌道高度為h的太陽同步軌道,衛(wèi)星繞地球一圈為t,則一天衛(wèi)星繞地球約q圈,赤道周長為l,則在赤道處,相鄰兩軌之間的距離為l3=l/q,已知軌道傾角為σ,則兩軌之間的直線距離為l4=l3×sinσ,由計算已知垂直于軌道的幅寬l5,則可計算相鄰兩軌的成像區(qū)域之間的疊加距離為l6=l5-l4;
如果成像區(qū)域兩側(cè)邊緣處變形較大,可通過圖像處理的過程中減小l6的取值,即在保證兩軌之間無縫拼接的情況下,減小兩軌之間重疊區(qū)域,以保證實(shí)現(xiàn)超大幅寬的同時避免成像區(qū)域兩側(cè)邊緣處的變形帶來的圖像變形。
本實(shí)施方式效果:
本實(shí)施方式在常規(guī)衛(wèi)星推掃和擺掃成像的基礎(chǔ)上,提出一種利用衛(wèi)星快速旋轉(zhuǎn)實(shí)現(xiàn)超大幅寬擺掃成像的方法。線陣ccd與衛(wèi)星飛行方向平行,通過衛(wèi)星360°連續(xù)自旋,衛(wèi)星快速旋轉(zhuǎn)速度1~10°/s;同時實(shí)現(xiàn)擺掃和推掃,實(shí)現(xiàn)衛(wèi)星對地成像覆蓋范圍的最大化從而實(shí)現(xiàn)中低軌道衛(wèi)星千公里級的超大幅寬擺掃成像。
傳統(tǒng)的推掃、擺掃衛(wèi)星成像,如果要對某一區(qū)域完全覆蓋,需要多次對該區(qū)域重訪,即多軌數(shù)據(jù)拼接;雖然加上側(cè)擺能力可以加大對區(qū)域的覆蓋能力,但是一次性采集地面數(shù)據(jù)的時間分辨率和覆蓋能力都不能滿足當(dāng)前衛(wèi)星快速響應(yīng)的任務(wù)需求。
本實(shí)施方式綜合考慮了衛(wèi)星的機(jī)動能力和成像能力,依據(jù)衛(wèi)星動態(tài)成像的姿態(tài)角和姿態(tài)角速度的變化及動態(tài)變化對成像的影響,設(shè)計了一種在衛(wèi)星飛行過程中,探測器光軸隨著衛(wèi)星的前進(jìn)方向旋轉(zhuǎn),在沿軌方向或非沿軌方向成像,從而實(shí)現(xiàn)一次性大范圍覆蓋成像的方法。本實(shí)施方式公開了一種衛(wèi)星快速旋轉(zhuǎn)超大幅寬擺掃成像的方法,該方法針對傳統(tǒng)遙感衛(wèi)星上采用的推掃、擺掃成像方法中幅寬有限的問題,在衛(wèi)星飛行過程中,使探測器光軸隨著衛(wèi)星的前進(jìn)緩慢的向前;隨著衛(wèi)星的360°勻速旋轉(zhuǎn),向左或右推移。使沿軌方向或非沿軌方向成像幅寬增大,以實(shí)現(xiàn)大范圍覆蓋成像的目的。本實(shí)施方式方法簡單可靠、易行,擴(kuò)大了推擺掃成像范圍,擴(kuò)展了推擺掃衛(wèi)星的應(yīng)用領(lǐng)域。本實(shí)施方式是國內(nèi)首個針對擴(kuò)大成像范圍提出的衛(wèi)星快速旋轉(zhuǎn)超大幅寬擺掃成像及大覆蓋面積的可靠方法。
探測器的旋轉(zhuǎn)可以通過探測器相對于衛(wèi)星平臺旋轉(zhuǎn)或者隨整星結(jié)構(gòu)整體旋轉(zhuǎn)實(shí)現(xiàn)。
為實(shí)現(xiàn)衛(wèi)星進(jìn)動和地球自轉(zhuǎn)情況下衛(wèi)星對地面大區(qū)域探測,衛(wèi)星姿態(tài)控制系統(tǒng)需要保證探測器360°連續(xù)勻速旋轉(zhuǎn),成像時衛(wèi)星姿態(tài)控制系統(tǒng)需要對衛(wèi)星本體坐標(biāo)系相對于軌道坐標(biāo)系的三軸姿態(tài)角(橫滾
具體實(shí)施方式二:本實(shí)施方式與具體實(shí)施方式一不同的是:步驟一中所述垂直于軌道的幅寬l5計算如下:
圓心角
垂直于軌道的幅寬
具體實(shí)施方式三:本實(shí)施方式與具體實(shí)施方式一或二不同的是:步驟二中所述飛行方向的幅寬l1具體為:
具體實(shí)施方式四:本實(shí)施方式與具體實(shí)施方式一至三之一不同的是:步驟三中l(wèi)2計算如下:
其中,衛(wèi)星在軌道上的角速度為
具體實(shí)施方式五:本實(shí)施方式與具體實(shí)施方式一至四之一不同的是:步驟五中所述不同側(cè)擺角成像對應(yīng)的分辨率ac計算方法如下,如圖6所示:
α=π-δ
則分辨率為:
ac=lc×a/f/cosα
δ為視線與地心半徑夾角,lc為物距,f為焦距;a為像元尺寸;
采用以下實(shí)施例驗證本實(shí)施方式的有益效果:
實(shí)施例一:
本實(shí)施例一種衛(wèi)星快速旋轉(zhuǎn)超大幅寬擺掃成像方法,具體是按照以下步驟制備的:
步驟一、假設(shè)衛(wèi)星軌道高度h=500km,取地球半徑r=6371km,計算圓心角2θ所對應(yīng)的衛(wèi)星軌道的曲線上的ab兩點(diǎn)之間的曲線距離
圓心角
垂直于軌道的幅寬
步驟二、假設(shè)探測器視場角η=33.2°,則根據(jù)視場角和軌道高度計算探測器的視場范圍l1具體為:
步驟三、只有探測器兩次相鄰成像區(qū)域之間無縫隙才能實(shí)現(xiàn)超大幅寬,成像區(qū)域示意圖見圖2,即兩次成像時探測器光軸中心在地表的星下點(diǎn)軌跡上的距離l2必須小于等于飛行方向的幅寬l1,因此,可認(rèn)為l2的臨界值是飛行方向的幅寬l1;l2計算如下:
其中,衛(wèi)星在軌道上的角速度為
步驟四、取臨界值l2=l1=298km帶入如下公式則計算出探測器沿軌方向的自旋速度:
探測器繞沿軌方向的自旋速度大于等于vz時,可以滿足探測器的兩次相鄰成像區(qū)域之間沒有縫隙實(shí)現(xiàn)超大幅寬成像,如果低于vz,則兩次成像區(qū)域之間會有縫隙,不能實(shí)現(xiàn)超大幅寬成像;
探測器繞沿軌方向的自旋速度不能低于9.2°/s,如果低于9.2°/s兩次成像區(qū)域之間會有縫隙。
步驟五、探測器的地表分辨率與其側(cè)擺角度大小密切相關(guān),一種衛(wèi)星快速旋轉(zhuǎn)超大幅寬擺掃成像的方法,通過探測器的360°旋轉(zhuǎn)能夠提供大范圍的側(cè)擺角,但并不是隨著側(cè)擺角的增大探測器的分辨率越高,只是在一定的側(cè)擺角范圍內(nèi),相同的軌道條件下,側(cè)擺角越大探測器的分辨率越高;為兼顧高分辨率和地面條帶無縫拼接兩方面的要求,需要成像幅寬與分辨率ac相互制約;
這是因為衛(wèi)星側(cè)擺角度很大時,衛(wèi)星與地球表面相切程度增大,除了衛(wèi)星分辨率隨著側(cè)擺物距的變化,衛(wèi)星在地面上的投影角度也變大,對應(yīng)成像的幾何分辨率變形嚴(yán)重,會出現(xiàn)嚴(yán)重的分辨率下降,如圖5所示;此時,結(jié)合探測器固有的特征參數(shù)選擇不同的側(cè)擺角度計算各種軌道條件下探測器實(shí)現(xiàn)超大幅寬成像時對應(yīng)的分辨率ac;衛(wèi)星側(cè)擺角與分辨率對應(yīng)關(guān)系如圖6所示;其中,探測器固有的特征參數(shù)包括像元尺寸和焦距;
所述不同側(cè)擺角成像對應(yīng)的分辨率計算方法如下如圖6所示:
α=π-δ
則分辨率為:
ac=lc×a/f/cosα
δ為視線與地心半徑夾角,lc為物距,f為焦距;a為像元尺寸;
步驟六、探測器的轉(zhuǎn)速vz及分辨率ac也決定了ccd行頻fp的大小,而ccd器件均具有自身的極限行頻,超出該極限行頻則探測器無法正常成像;根據(jù)探測器自旋速度vz及分辨率ac計算ccd行頻fp:
假設(shè)探測器的焦距f=1.18m,像元尺寸a=3.5μm,此時,衛(wèi)星星下點(diǎn)對應(yīng)分辨率為ac=500km×3.5μm/1.18m=1.5m。當(dāng)探測器側(cè)擺角度取為
步驟七、如圖2所示,探測器自旋成像時,成像區(qū)域兩側(cè)變形較大,考慮兩軌之間即相鄰兩次成像區(qū)域之間的距離,可以減小成像區(qū)域在垂直沿軌方向的距離l6,即兩軌之間重疊區(qū)域,使用兩軌數(shù)據(jù)拼接的方式;如圖7所示;
已知衛(wèi)星軌道高度為h=500km的太陽同步軌道,衛(wèi)星繞地球一圈為t=5676.98s,則一天衛(wèi)星繞地球約q=15.2圈,赤道周長為l=40000km,則在赤道處,相鄰兩軌之間的距離為l3=l/q=40000km/15.2≈2631.6km,已知軌道傾角為σ=97.4°,則兩軌之間的直線距離為l4=l3×sinσ=2631.6km×sin97.4°≈2609.7km,由計算的該軌道下已知垂直于軌道的幅寬l5最大為4890km,則可計算相鄰兩軌的成像區(qū)域之間的疊加距離為l6=l5-l4=4890km-2609.7km=2280.3km;
如果成像區(qū)域兩側(cè)邊緣處變形較大,可通過圖像處理的過程中減小l6的取值,即在保證兩軌之間無縫拼接的情況下,減小兩軌之間重疊區(qū)域,以保證實(shí)現(xiàn)超大幅寬的同時避免成像區(qū)域兩側(cè)邊緣處的變形帶來的圖像變形。
顯然,上述實(shí)施例僅僅是為清楚地說明所作的舉例,而并非對實(shí)施方式的限定。對于所屬領(lǐng)域的普通技術(shù)人員來說,在上述說明的基礎(chǔ)上還可以做出其它不同形式的變化或變動。這里無需也無法對所有的實(shí)施方式予以窮舉。而由此所引伸出的顯而易見的變化或變動仍處于本發(fā)明創(chuàng)造的保護(hù)范圍之中。
本實(shí)施方式還可有其它多種實(shí)施例,在不背離本實(shí)施方式精神及其實(shí)質(zhì)的情況下,本領(lǐng)域技術(shù)人員當(dāng)可根據(jù)本實(shí)施方式作出各種相應(yīng)的改變和變形,但這些相應(yīng)的改變和變形都應(yīng)屬于本實(shí)施方式所附的權(quán)利要求的保護(hù)范圍。