2634
FAQ: Radiometric Calibration of Worldview data
- Where to locate Gain and Offset values in Worldview metadata ?
Worldview data are not provided with Gain and Offset values. Only Absolute Scale Factor (ASF) and Effective Bandwith (EBW) are provided. Those values are located in the .mdi or .xml files provided with the dataset. ENVI reads those values them from the .xml file
For example: Wordlview-3 data SWIR Band S1
<ABSCALFACTOR>2.671600000000000e-04</ABSCALFACTOR>
<EFFECTIVEBANDWIDTH>3.300000000000000e-02</EFFECTIVEBANDWIDTH>
The absolute radiometric calibration values for Gain and Offset are not included in WorldView dataset. They are provided by Digital Globe and are revised regularly.
- Which absolute Gain and Offset values are used in ENVI for Worldview data calibration
In ENVI the absolute radiometric calibration values for Gain and Offset are hard coded in the following file
<INSTALL_DIR>\Harris\ENVI56\resource\filterfuncs\sensor_attributes.json
The file is updated when Digital Globe provides revised data
For example: Wordlview-3 data SWIR Band
- Offset is equal to 0 for all SWIR bands
- Gain values are:
"S1": 1.030,
"S2": 1.052,
"S3": 0.992,
"S4": 1.014,
"S5": 1.012,
"S6": 1.082,
"S7": 1.056,
"S8": 1.101
- Which version of absolute Gain and Offset values are hardcoded in ENVI 5.6.2
The data hardcoded in the sensor_attributes.json file for ENVI 5.6.2 are compliant with the slide 10 of MAXAR/DigitalGlobe reference file below (page 9 -10)
https://calval.cr.usgs.gov/apps/sites/default/files/jacie/MicheleKuester.pdf
They correspond to the Gain and Offset provided with calibration campaign 2018v0 for VNIR bands and 2019v0 for SWIR bands
- Why Gain values displayed in the ENVI metadata viewer differ from those hardcoded in the sensor_attributes.json file?
ENVI metadata viewer displays Gain and Offset values for Worldview data:
- The Offset is the one hardcoded in the sensor_attributes.json file
- However, Gain values as viewed in the ENVI Metadata viewer are different from the one hardcoded in the sensor_attributes.json file.
Gain values included in sensor_attributes.json file are absolute values. They need to be adjusted with a factor corresponding to the ratio ASF/EBW for each dataset
For example:
Gain values of WV-3 SWIR bands from sensor_attributes.json are listed below
G= [1.03,1.052,0.992,1.014,1.012,1.082,1.056,1.101]
ASF values from a given WV-3 data for SWIR bands, extracted from the .xml file associated to this dataset, are listed below
ASF=[2.671600e-04, 1.728110e-04,1.576480e-04,1.466560e-04,6.666700e-05,6.862700e-05,6.862700e-05,7.254900e-05]
EBW values from a given WV-3 data for SWIR bands, extracted from the .xml file associated to this dataset, are listed below
EBW=[3.300000e-02,3.970000e-02,3.730000e-02,4.160000e-02,3.890000e-02,4.090000e-02,4.760000e-02,6.790000e-02]
The following computation (ASF/EBW)*G will result in the below output:
G_ADJ=[0.00833863,0.00457927,0.00419268,0.00357474,0.00173437,0.00181551,0.00152248,0.00117638]
Those values are the Gain values displayed in ENVI metadata viewer for the corresponding SWIR bands of this dataset
- Why the equation provided by MAXAR/DigitalGlobe for radiance calibration differs from the one in ENVI help
ENVI helps shows the below equation to calibrate data in radiance
𝐿 = 𝐺𝐴𝐼𝑁 ∗ 𝐷𝑁 + 𝑂𝐹𝐹𝑆𝐸𝑇
MAXAR/DigitalGlobe (p.5 of the link below) provides a different equation
https://dg-cms-uploads-production.s3.amazonaws.com/uploads/document/file/207/Radiometric_Use_of_WorldView-3_v2.pdf?
𝐿 = 𝐺𝐴𝐼𝑁 ∗ (ASF/EBW) ∗ 𝐷𝑁 + 𝑂𝐹𝐹𝑆𝐸𝑇
As mentioned at point 4, Gain values displayed in ENVI metadata viewer are adjusted as follow
G_ADJ=(ASF/EBW)*G
The adjusted gain is then used by ENVI for the radiance calibration in the equation below
𝐿 = G_ADJ ∗ 𝐷𝑁 + 𝑂𝐹𝐹𝑆𝐸𝑇 = 𝐺 ∗ (ASF/EBW) ∗ 𝐷𝑁 + 𝑂𝐹𝐹𝑆𝐸𝑇
So the 2 equations for radiance calibration provided by ENVI and MAXAR/DigitalGlobe are equivalent
- Why the equation provided by MAXAR/DigitalGlobe for reflectance calibration differs from the one in ENVI help
ENVI helps shows the below equation to calibrate data in radiance
𝜌𝜆 = 𝐿𝜆𝑑²𝜋 / 𝐸SUN sin 𝜃
With 𝜃 the sun elevation
MAXAR/DigitalGlobe (p.7 of the link below) provides a different equation
https://dg-cms-uploads-production.s3.amazonaws.com/uploads/document/file/207/Radiometric_Use_of_WorldView-3_v2.pdf?
𝜌𝜆 = 𝐿𝜆𝑑²𝜋 / 𝐸𝜆 cos 𝜃𝑆
With 𝜃𝑆 the solar zenith angle
In ENVI 𝜃 angle is defined as the sun elevation while in for MAXAR/DigitalGlobe 𝜃𝑆 is the solar zenith angle.
Still from MAXAR/DigitalGlobe reference (p. 11), the solar zenith angle is equal to 90°- the sun elevation angle so
𝜃𝑆 =90- 𝜃= 𝜋 /2 - 𝜃
From trigonometry
sin(𝜃)=cos (𝜋 /2 - 𝜃)
So
cos(𝜃𝑆)=cos(𝜋 /2 - 𝜃)=sin(𝜃)
The solar irradiance has a similar definition in both cases:
Eλ =Esun
So the 2 equations for radiance calibration provided by ENVI and MAXAR/DigitalGlobe are equivalent
-------------------------------
created by BC on 7/2/2022
reviewed by MM on 7/18/22