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