微波遙感集中編程實習(xí)報告二(專業(yè)版)
【正文】 然后就是計算對應(yīng)SAR影像的行列號處，每次運行完我編寫的程序后得到的影像都是一片黑，設(shè)斷點調(diào)試程序時發(fā)現(xiàn)每次的SAR影像上對應(yīng)的行號i計算出來的都是負(fù)值，不在SAR影像包括的范圍內(nèi)。利用公式求出其強度信息并讀出來存儲起來。注意為解求矩陣方便，可將第一點成像時間當(dāng)作0。掌握距離多普勒方程進(jìn)行校正的原理和難點，認(rèn)識校正所需參數(shù)意義與讀取方法。即利用DEM中各點的經(jīng)緯度坐標(biāo)BL結(jié)合DEM影像上的高度信息H，加上投影橢球參數(shù)，計算出各點的大地坐標(biāo)XYZ，再由該點的位置矢量結(jié)合當(dāng)時衛(wèi)星的位置矢量、速度矢量和加速度矢量代入距離多普勒方程，確定成像時刻，從而找到SAR影像上對應(yīng)的（i，j），取出此處的灰度值，賦給DEM對應(yīng)的大地坐標(biāo)處。可以通過設(shè)置一個閾值，當(dāng)時，迭代停止，此時求出的時間為真實方位向時間。當(dāng)然還是很感謝老師在編程過程中提供的幫助，老師十分耐心地予以指導(dǎo)，在我數(shù)組指針無法動態(tài)提供內(nèi)存時，讓我們不需要保存這些中間數(shù)據(jù)，只需在使用時再直接計算即可。代碼部分如下：圖7 結(jié)果的存儲與輸出圖8 運行結(jié)果四、 實習(xí)心得通過這次幾何校正編程實習(xí)，我弄清楚了利用RD方程間接幾何校正的具體流程，需要結(jié)合DEM數(shù)據(jù)，迭代求出SAR影像上的行列號。注意一些數(shù)據(jù)需要進(jìn)行轉(zhuǎn)換，如投影橢球的長短半軸長單位是千米，需要轉(zhuǎn)換到米為單位；經(jīng)緯度是以度為單位的，在由經(jīng)緯度坐標(biāo)BLH計算到大地坐標(biāo)XYZ時需要把經(jīng)緯度轉(zhuǎn)為弧度為單位；考慮到雙程傳輸，時間需要進(jìn)
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1