GDAL庫本身就可以讀取Envisat的圖像數據,具體鏈接爲:http://www.gdal.org/frmt_various.html#Envisat。
但是對於ASAR傳感器的數據來說,GDAL在處理的時候將入射角數據沒要進行處理,也就是通過GDAL提供的接口無法獲取到這個信息。
ASAR數據進行輻射定標的時候,入射角信息是必須要的信息,如果通過GDAL庫來對ASAR數據進行輻射定標,那麼因爲沒法獲取到入射角信息從而導致沒法進行輻射定標。
ASAR數據裏面入射角信息位於GEOLOCATION GRID ADS中,GDAL庫在解析該部分內容時,只將GCP信息進行了解析,而對於入射角、時間等信息直接跳過,沒有進行處理,所以要使用GDAL庫能夠讀取入射角信息,最簡單的方式就是將入射角信息解析到GDAL的元數據中,這樣就可以在不改動任何接口的情況下來獲取入射角信息。
基於上面的分析,需要對GDAL庫中讀取ENVISAT數據的驅動進行修改,增加解析入射角的內容,並將入射角構造成元數據信息寫入Dataset中。通過分析ENVISAT驅動的源碼文件envisatdataset.cpp,可以在函數ScanForGCPs_ASAR()中增加相應的代碼。
增加的代碼片段爲:
- 解析First Line的數據
//解析入射角信息
char szKey[32];
char szValue[32];
//解析Sample
memcpy(&unValue, abyRecord + 25 + iRecord * 4, 4);
int nSample = CPL_MSBWORD32(unValue);
snprintf(szKey, sizeof(szKey), "FIRST_LINE_SAMPLE_%d", iRecord);
snprintf(szValue, sizeof(szValue), "%d", nSample);
SetMetadataItem(szKey, szValue, "GEOLOCATION_GRID_ADSR");
//解析Time
float fValue;
memcpy(&fValue, abyRecord + 25 + 44 + iRecord * 4, 4);
CPL_SWAP32PTR(&fValue);
snprintf(szKey, sizeof(szKey), "FIRST_LINE_TIME_%d", iRecord);
snprintf(szValue, sizeof(szValue), "%.6f", fValue);
SetMetadataItem(szKey, szValue, "GEOLOCATION_GRID_ADSR");
//解析Angles
memcpy(&fValue, abyRecord + 25 + 88 + iRecord * 4, 4);
CPL_SWAP32PTR(&fValue);
snprintf(szKey, sizeof(szKey), "FIRST_LINE_ANGLES_%d", iRecord);
snprintf(szValue, sizeof(szValue), "%.6f", fValue);
SetMetadataItem(szKey, szValue, "GEOLOCATION_GRID_ADSR");
- 解析Last Line的數據
//解析入射角信息
char szKey[32];
char szValue[32];
//解析Sample
snprintf(szKey, sizeof(szKey), "LAST_LINE_SAMPLE_%d", iGCP);
snprintf(szValue, sizeof(szValue), "%d", nSample);
SetMetadataItem(szKey, szValue, "GEOLOCATION_GRID_ADSR");
//解析Time
float fValue;
memcpy(&fValue, abyRecord + 279 + 44 + iGCP * 4, 4);
CPL_SWAP32PTR(&fValue);
snprintf(szKey, sizeof(szKey), "LAST_LINE_TIME_%d", iGCP);
snprintf(szValue, sizeof(szValue), "%.6f", fValue);
SetMetadataItem(szKey, szValue, "GEOLOCATION_GRID_ADSR");
//解析Angles
memcpy(&fValue, abyRecord + 279 + 88 + iGCP * 4, 4);
CPL_SWAP32PTR(&fValue);
snprintf(szKey, sizeof(szKey), "LAST_LINE_ANGLES_%d", iGCP);
snprintf(szValue, sizeof(szValue), "%.6f", fValue);
SetMetadataItem(szKey, szValue, "GEOLOCATION_GRID_ADSR");
增加上面兩部分代碼之後,重新編譯GDAL庫,生成DLL即可。下面使用gdalinfo工具來進行測試,輸出的入射角、時間信息如下:
Metadata (GEOLOCATION_GRID_ADSR):
FIRST_LINE_ANGLES_0=14.296353
FIRST_LINE_ANGLES_1=15.133030
FIRST_LINE_ANGLES_10=22.277876
FIRST_LINE_ANGLES_2=15.960590
FIRST_LINE_ANGLES_3=16.779432
FIRST_LINE_ANGLES_4=17.589853
FIRST_LINE_ANGLES_5=18.389204
FIRST_LINE_ANGLES_6=19.186239
FIRST_LINE_ANGLES_7=19.972475
FIRST_LINE_ANGLES_8=20.750854
FIRST_LINE_ANGLES_9=21.521433
FIRST_LINE_SAMPLE_0=1
FIRST_LINE_SAMPLE_1=832
FIRST_LINE_SAMPLE_10=8304
FIRST_LINE_SAMPLE_2=1663
FIRST_LINE_SAMPLE_3=2494
FIRST_LINE_SAMPLE_4=3325
FIRST_LINE_SAMPLE_5=4153
FIRST_LINE_SAMPLE_6=4987
FIRST_LINE_SAMPLE_7=5818
FIRST_LINE_SAMPLE_8=6649
FIRST_LINE_SAMPLE_9=7480
FIRST_LINE_TIME_0=5414506.500000
FIRST_LINE_TIME_1=5432814.000000
FIRST_LINE_TIME_10=5637850.500000
FIRST_LINE_TIME_2=5452048.500000
FIRST_LINE_TIME_3=5472201.000000
FIRST_LINE_TIME_4=5493261.000000
FIRST_LINE_TIME_5=5515137.000000
FIRST_LINE_TIME_6=5538061.500000
FIRST_LINE_TIME_7=5561780.000000
FIRST_LINE_TIME_8=5586363.000000
FIRST_LINE_TIME_9=5611798.500000
LAST_LINE_ANGLES_0=14.357172
LAST_LINE_ANGLES_1=15.185051
LAST_LINE_ANGLES_10=22.280605
LAST_LINE_ANGLES_2=16.004770
LAST_LINE_ANGLES_3=16.816580
LAST_LINE_ANGLES_4=17.620653
LAST_LINE_ANGLES_5=18.414248
LAST_LINE_ANGLES_6=19.206030
LAST_LINE_ANGLES_7=19.987459
LAST_LINE_ANGLES_8=20.761423
LAST_LINE_ANGLES_9=21.527927
LAST_LINE_SAMPLE_0=1
LAST_LINE_SAMPLE_1=832
LAST_LINE_SAMPLE_10=8304
LAST_LINE_SAMPLE_2=1663
LAST_LINE_SAMPLE_3=2494
LAST_LINE_SAMPLE_4=3325
LAST_LINE_SAMPLE_5=4153
LAST_LINE_SAMPLE_6=4987
LAST_LINE_SAMPLE_7=5818
LAST_LINE_SAMPLE_8=6649
LAST_LINE_SAMPLE_9=7480
LAST_LINE_TIME_0=5414506.500000
LAST_LINE_TIME_1=5432814.000000
LAST_LINE_TIME_10=5637850.500000
LAST_LINE_TIME_2=5452048.500000
LAST_LINE_TIME_3=5472201.000000
LAST_LINE_TIME_4=5493261.000000
LAST_LINE_TIME_5=5515137.000000
LAST_LINE_TIME_6=5538061.500000
LAST_LINE_TIME_7=5561780.000000
LAST_LINE_TIME_8=5586363.000000
LAST_LINE_TIME_9=5611798.500000
從上面的信息中可以看出,ASAR數據裏面的入射角等信息都位於GEOLOCATION_GRID_ADSR元數據之下,後續只需從該元數據中解析入射角、Sample等信息進行插值即可得到輻射定標的相關參數。