;IDL编写NDVI计算函数，数值范围控制在[-1,1]，实现在ENVI中BAND MATH中调用

在ENVI (Environment for Visualizing Images) 中，要编写一个 IDL (Image Data Language) 函数来计算归一化差值植被指数 (Normalized Difference Vegetation Index, NDVI)，并确保结果数值范围在 [-1, 1] 内，可以按照以下步骤进行：

; 定义NDVI计算函数
function NDVI(band_red, band_nir)
    ; 检查输入数据类型
    if ~(band_red~type=real || band_nir~type=real) then begin
        print, "Error: Both bands must be of real type."
        return
    endif

    ; 计算NDVI
    ndvi = (band_nir - band_red) / (band_nir + band_red)

    ; 对结果应用边界条件
    ndvi = where(ndvi < -1, -1, ndvi)
    ndvi = where(ndvi > 1, 1, ndvi)

    ; 返回处理后的NDVI值
    return ndvi
end

; 调用示例
new_bnd = NDVI(red_band, nir_band)

在这个例子中，`red_band` 和 `nir_band` 分别代表红光和近红外波段的数据。函数首先检查输入是否为实数，然后通过公式 `(NIR - RED) / (NIR + RED)` 计算NDVI，接着应用边界条件将超出范围的值调整到 [-1, 1] 内。