這個是2012年左右放在百度空間的，誰知百度空間關閉。。。轉移到博客園。
最近項目用到3.1.2版本的LEDAPS，新版本的使用情況會在后續文章中慢慢豐富。
HuSr是將LEDAPS項目中的TM/ETM+大氣校正流程系列算法中的大氣校正模塊由linux系統移植到windows下的產物，代碼本身改動不大，使用接口不變。
包含文件：
HuSr.exe
6sv1.0b.exe
命令行調用示例：HuSr.exe空格D:\LEDAPS\lndsr.L5014033_03320060801.txt
Vs2005下編譯，包含的靜態庫有：
xdr.lib
libjpeg.lib
zlib.lib
szlib.lib
HD426.lib
HM426.lib
ws2_32.lib
gctp.lib
hdfeos.lib
為結果重復定義，忽略的庫包括：libc.lib;libcmtd.lib
原始工程中包含CSALBR.f與CHAND.f兩個fotran代碼，這兩個文件原屬於輻射傳輸包6s中的代碼，為了在vs2005的c＋＋工程中加載它們，於是先在linux下使用f2c工具將其代碼轉換成標准c語言，然后再加載到vs2005的項目中，同時還需將f2c的頭文件f2c.h加入。
原始代碼中有幾個linux系統的函數使用標准c庫中的函數代替。Windows下缺少unistd.h函數定義的問題采用的解決方案與編譯前兩個模塊的相同。
程序需要外部調用6s，而代用外部程序的方法linux與windows不同，並且6s在linux與windows中的調用方式也不同，所以在sixs_runs.c文件中對調用6s的部分進行了修改。6s的調用版本是6sv1.0b.exe，原始的fotran代碼我在windows下面用f77沒有編譯成功，只好從網上下載了一個編譯好的exe用，結果表明計算結果與linux的相同，應該沒問題吧。
對原有的hdf庫做了更新，使用hdfeos庫為2.18b版本,HDF4庫使用4.26版本。沒有發現有啥問題。
校正一景TM數據在15分鍾左右，要比linux版本的稍快些。
0 (user defined)
29.28 129.03 0.00 0.00 8 1 (geometrical conditions sza saz vza vaz month day)
8 (option for water vapor and ozone)
4.20 0.30 (water vapor and ozone)
1 (continental model)
0 (option for optical thickness at 550 nm)
0.010 (value of aot550
0.000000 (target level)
-1000 (sensor level : -1000=satellite level)
25 (predefined band)
0 (homogeneous surface)
0 (no directional effects)
0 (constant value for rho)
0.140 (value of rho)
-1 (no atmospheric correction)
0
程序中需要的參數用黃色標記：
自己注釋用紅色字體
******************************* 6sV version 1.0B ******************************
* *
* geometrical conditions identity *
* ------------------------------- *
* user defined conditions *
* *
* month: 8 day : 1 需要的參數與具體日期無關 *
* solar zenith angle: 29.28 deg solar azimuthal angle: 129.03 deg *
* view zenith angle: 0.00 deg view azimuthal angle: 0.00 deg *
* scattering angle: 150.72 deg azimuthal angle difference: 129.03 deg *
* *
* atmospheric model description *
* ----------------------------- *
* atmospheric model identity : *
* user defined water content : uh2o= 4.200 g/cm2 *
* user defined ozone content : uo3 = 0.300 cm-atm *
* aerosols type identity : *
* Continental aerosol model *
* optical condition identity : *
* visibility :****** km opt. thick. 550 nm : 0.0100 *
* *
* spectral condition *
* ------------------ *
* tm 1 *
* value of filter function : *
* wl inf= 0.430 mic wl sup= 0.560 mic *
* *
* Surface polarization parameters *
* ---------------------------------- *
* *
* *
* Surface Polarization Q,U,Rop,Chi 0.00000 0.00000 0.00000 -1.#J *
* *
* *
* target type *
* ----------- *
* homogeneous ground *
* constant reflectance over the spectra 0.140 *
* *
* target elevation description *
* ---------------------------- *
* ground pressure [mb] 1013.00 *
* ground altitude [km] 0.000 高度一直為0 *
* *
*******************************************************************************
*******************************************************************************
* *
* integrated values of : *
* -------------------- *
* *
* apparent reflectance 0.1820603 appar. rad.(w/m2/sr/mic) 96.005 *
* total gaseous transmittance 0.987 *
* *
*******************************************************************************
* *
* coupling aerosol -wv : *
* -------------------- *
* wv above aerosol : 0.182 wv mixed with aerosol : 0.182 *
* wv under aerosol : 0.182 *
*******************************************************************************
* *
* integrated values of : *
* -------------------- *
* *
* app. polarized refl. 0.0075 app. pol. rad. (w/m2/sr/mic) 0.240 *
* direction of the plane of polarization -39.03 *
* total polarization ratio 0.041 *
* *
*******************************************************************************
* *
* int. normalized values of : *
* --------------------------- *
* % of irradiance at ground level *
* % of direct irr. % of diffuse irr. % of enviro. irr *
* 0.881 0.101 0.018 *
* reflectance at satellite level *
* atm. intrin. ref. background ref. pixel reflectance *
* 0.064 0.011 0.107 *
* *
* int. absolute values of *
* ----------------------- *
* irr. at ground level (w/m2/mic) *
* direct solar irr. atm. diffuse irr. environment irr *
* 1344.346 153.382 27.867 *
* rad at satel. level (w/m2/sr/mic) *
* atm. intrin. rad. background rad. pixel radiance *
* 33.758 5.576 56.671 *
* *
* *
* int. funct filter (in mic) int. sol. spect (in w/m2) *
* 0.0604850 114.879 *
* *
*******************************************************************************
*******************************************************************************
* *
* integrated values of : *
* -------------------- *
* *
* downward upward total *
* global gas. trans. : 0.99294 0.99384 0.98684 *
* water " " : 1.00000 1.00000 1.00000T_g_wv[i]*
* ozone " " : 0.99294 0.99384 0.98684tgoz*
* co2 " " : 1.00000 1.00000 1.00000tgco2*
* oxyg " " : 1.00000 1.00000 1.00000tgo2*
* no2 " " : 1.00000 1.00000 1.00000tgno2*
* ch4 " " : 1.00000 1.00000 1.00000tgch4*
* co " " : 1.00000 1.00000 1.00000tgco*
* T_g_og[i]=tgoz*tgco2*tgo2*tgno2*tgno2*tgch4*tgco*
* *
* rayl. sca. trans. : 0.91343T_r_down[i]0.92360T_r_up[i]0.84364T_r[i] *aeros. sca. " : 0.99733T_a_down[ij]0.99783T_a_up[ij]0.99517T_a[ij]
* total sca. " : 0.91057T_ra_down[ij]0.92124T_ra_up[ij]0.83885T_ra[ij]*
* *
* *
* *
* rayleigh aerosols total *
* *
* spherical albedo : 0.12771S_r[i]0.00454 0.13062S_ra[i][j]*
* optical depth total: 0.16504 0.01133 aot_wavelength[i][j] 0.17637*
* optical depth plane: 0.16504 0.01133 0.17637 *
* reflectance I : 0.06415rho_r[i]0.00059rho_a[ij]0.06478rho_ra[ij]*
* reflectance Q : 0.00160 -0.00003 0.00157 *
* reflectance U : -0.00756 0.00000 -0.00744 *
* polarized reflect. : 0.00773 0.00003 0.00760 *
* degree of polar. : 12.04 5.80 11.74 *
* dir. plane polar. : -39.03 90.00 -39.03 *
* phase function I : 1.30737 0.19812 1.23612 *
* phase function Q : -0.17199 0.05671 -0.15730 *
* phase function U : -1.25436 -0.06031 -1.17766 *
* primary deg. of pol: -0.13156 0.28625 -0.12726 *
* sing. scat. albedo : 1.00000 0.89911 0.99352 *
* *
* *
*******************************************************************************