微波遙感集中編程實習(xí)報告二(編輯修改稿)
2024-08-30 12:41 本頁面
【文章內(nèi)容簡介】 求解任意時刻的衛(wèi)星位置矢量、速度矢量和加速度矢量，這里對擬合參數(shù)的求解是利用矩陣最小二乘法求解的。在矩陣計算時，由于時間t較大，高次方運算可能會超過數(shù)據(jù)范圍，所以把初始成像時間設(shè)為0。注意一些數(shù)據(jù)需要進行轉(zhuǎn)換，如投影橢球的長短半軸長單位是千米，需要轉(zhuǎn)換到米為單位；經(jīng)緯度是以度為單位的，在由經(jīng)緯度坐標(biāo)BLH計算到大地坐標(biāo)XYZ時需要把經(jīng)緯度轉(zhuǎn)為弧度為單位；考慮到雙程傳輸，時間需要進行改正。需要利用到DEM數(shù)據(jù)和SAR影像，利用GDAL庫的RasterIO函數(shù)將影像數(shù)據(jù)讀入。，分為實部和虛部，數(shù)據(jù)格式是CInt16，而我們需要的是它的灰度信息，灰度信息是影像的強度，符合公式：強度178。=實部178。+虛部178。，利用公式求出其強度信息并讀出來存儲起來。編程實現(xiàn)DEM影像讀取和SAR影像讀取并存儲。圖4 DEM的讀取圖5 SAR的讀取2. 由RD方程計算行列號利用RD方程間接校正，即利用DEM的大地坐標(biāo)，結(jié)合衛(wèi)星位置矢量、速度矢量、加速度矢量代入RD方程，求解出該點對應(yīng)到SAR影像上的行列號。對多普勒方程求微分，有式中VS、RS表示衛(wèi)星速度矢量和位置矢量，RP表示DEM大地坐標(biāo)，對DEM的任一點來說，成像時間未知，多普勒方程不為0，可針對方位向時間進行迭代。以中間行成像時間為初始迭代時間，求出當(dāng)前時刻的VS、RS、aS，結(jié)合DEM上該點的RP，代入上式可以求出，令，再進行下一次迭代，隨著迭代次數(shù)的增加，會減小，成像時間會
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1