以 Landsat8 近红外波段,即 band5 为例。从头文件中获取: Radiance_max=64723 Radiance_min=0 DN_max=64516 DN_min=0 已知辐射定标公式:L=(a*DN+b)*0.1 将 Radinance 于 DN 值的最大值最小值带入定标公式中,结算获得 a 与 b 参数的值。 利用 a 与 b 参数值、定标公式、ENVI 中的 band_math 工具对近红 外波段整幅影像进行辐射定标。
时间: 2024-03-04 13:50:35 浏览: 84
根据辐射定标公式,可以得到:
a = (Radiance_max - Radiance_min) / (DN_max - DN_min) / 0.1
b = Radiance_min / 10 - a * DN_min / 10
将 Radiance_max=64723,Radiance_min=0,DN_max=64516,DN_min=0 带入上述公式,可以得到:
a = 0.099998
b = 0.000000
然后,利用 a 和 b 带入辐射定标公式 L=(a*DN+b)*0.1,即可对近红外波段整幅影像进行辐射定标。可以使用 ENVI 中的 band_math 工具进行操作,公式如下:
(Landsat8_B5 * 0.099998 + 0.000000) * 0.1
其中,Landsat8_B5 表示 Landsat8 的 band5 波段。
相关问题
以 Landsat8 近红外波段,即 band5 为例。从头文件中获取: Radiance_max=377.15604 Radiance_min=-31.14567 DN_max=65535 DN_min=1 已知辐射定标公式:L=(a*DN+b)*0.1 将 Radinance 于 DN 值的最大值最小值带入定标公式中,结算获得 a 与 b 参数的值。 利用 a 与 b 参数值、定标公式、ENVI 中的 band_math 工具对近红 外波段整幅影像进行辐射定标。
可以按照以下步骤进行:
1. 打开 ENVI 软件,选择 File -> Open Image File,选择 Landsat8 数据中的 "B5.TIF" 文件并打开。
2. 在 ENVI 中打开 "B5.TIF.hdr" 文件,找到 Radiance_max、Radiance_min、DN_max 和 DN_min 的值,记录下来。
3. 根据公式 L = (a * DN + b) * 0.1,将 Radiance_max 和 Radiance_min 代入得到以下两个方程:
Radiance_max = (a * DN_max + b) * 0.1
Radiance_min = (a * DN_min + b) * 0.1
将上述两个方程联立,可以解得 a 和 b 的值:
a = (Radiance_max - Radiance_min) * 10 / (DN_max - DN_min)
b = Radiance_max * 0.1 - a * DN_max
将已知数据代入计算,可以得到:
a = (377.15604 - (-31.14567)) * 10 / (65535 - 1) ≈ 0.005831
b = 377.15604 * 0.1 - 0.005831 * 65535 ≈ -193.465
4. 选择 Raster -> Spectral -> Band Math,打开 band_math 工具。
5. 在 band_math 工具中,输入下列表达式进行辐射定标:
(0.005831 * b5) - 193.465
其中,0.005831 和 -193.465 分别是上一步中计算得到的 a 和 b 参数值,b5 表示输入的 Landsat8 近红外波段数据。
6. 点击 OK 按钮,等待处理完成。
7. 处理完成后,保存输出结果。可以选择 File -> Save As,设置输出文件名和文件格式,并保存结果。
经过以上步骤,就可以利用 a 和 b 参数值、定标公式、ENVI 中的 band_math 工具对 Landsat8 近红外波段整幅影像进行辐射定标。
阅读全文