微波遙感集中編程實(shí)習(xí)報(bào)告二(文件)
【正文】 程實(shí)習(xí)SAR影像幾何校正。注意一些數(shù)據(jù)需要進(jìn)行轉(zhuǎn)換，如投影橢球的長短半軸長單位是千米，需要轉(zhuǎn)換到米為單位；經(jīng)緯度是以度為單位的，在由經(jīng)緯度坐標(biāo)BLH計(jì)算到大地坐標(biāo)XYZ時(shí)需要把經(jīng)緯度轉(zhuǎn)為弧度為單位；考慮到雙程傳輸，時(shí)間需要進(jìn)行改正。+虛部178。對(duì)多普勒方程求微分，有式中VS、RS表示衛(wèi)星速度矢量和位置矢量，RP表示DEM大地坐標(biāo)，對(duì)DEM的任一點(diǎn)來說，成像時(shí)間未知，多普勒方程不為0，可針對(duì)方位向時(shí)間進(jìn)行迭代。再利用第一行成像時(shí)間和方位向成像時(shí)間間隔，求出行號(hào)。代碼部分如下：圖7 結(jié)果的存儲(chǔ)與輸出圖8 運(yùn)行結(jié)果四、 實(shí)習(xí)心得通過這次幾何校正編程實(shí)習(xí)，我弄清楚了利用RD方程間接幾何校正的具體流程，需要結(jié)合DEM數(shù)據(jù)，迭代求出SAR影像上的行列號(hào)。我原本希望將計(jì)算出的經(jīng)緯度值都以矩陣的方式存儲(chǔ)起來，但程序運(yùn)行到這里時(shí)總是提示錯(cuò)誤，錯(cuò)誤理由檢查后發(fā)現(xiàn)是內(nèi)存沒分配到，在詢問老師后，老師說每個(gè)程序運(yùn)行時(shí)內(nèi)存都是一定的，我這樣寫內(nèi)存不足沒法實(shí)現(xiàn)分配存儲(chǔ)，建議我在每次需要利用BLH求解大地坐標(biāo)XYZ時(shí)再計(jì)算經(jīng)緯度BL，也不要把它存儲(chǔ)起來。總的來說，這次編程實(shí)習(xí)，我收獲很大。當(dāng)然還是很感謝老師在編程過程中提供的幫助，老師十分耐心地予以指導(dǎo)，在我數(shù)組指針無法動(dòng)態(tài)提供內(nèi)存時(shí)，讓我們不需要保存這些中間數(shù)據(jù)，只需在使用時(shí)再直接計(jì)算即可。經(jīng)檢查后發(fā)現(xiàn)是自己的經(jīng)緯度BL弄反了，并且由經(jīng)緯度坐標(biāo)BLH計(jì)算地理坐標(biāo)XYZ的公式就寫錯(cuò)了，導(dǎo)
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1