遥感图像包络线去除

包络线去除是光谱识别中一个重要的步骤，算法如下：
    （1）通过求导得到光谱曲线上所有极大值点，即“峰”值点，然后比较大小，得到极大值点中的最大值点；
    （2）以最大值点作为包络线的一个端点，计算该点与长波方向（波长增长的方向）各个极大值连线的斜率，以斜率最大点作为包络线下一个端点，再以此点为起点循环，直至最后一点；
    （3）以最大值点作为包选线的一个端点，向短波（波长减少的方向）进行类似计算，以斜率最小点为下个端点，再以此点为起点循环，直至曲线上的开始点；

    （4）沿长波方向连接所有端点，可形成曲线的包络线，用实际光谱反射率去除包络线上相应波段的反射率值，可得到包络线消除法归一化后的值。

/**
    求取包络线，去除包络线(为了保存所有去除包络线数据)--->为了提高速度
    r:读取的光谱数据
    w:读取的波长
    h:输出的去包络线数据
*/

template <typename T> void RsEnvelopingLine(T* r, float* w, float* h, int bands)
{
    int i = 0, j = 1, m;
    int k = 0;
    float tempHH = 0;
    while (i < bands - 1)
    {
        if (j == bands - 1)
        {
            //从i-j循环
            float temp1 = ((float)r[j] - (float)r[i]) / (w[j] - w[i]);
            float temp2 = (float)r[i] - temp1 * w[i];
            for (k = i; k <= j; k++)
            {
                tempHH = temp1 * w[k] + temp2;
                if (tempHH == 0)
                {
                    h[k] = 0.0f;
                }
                else
                {
                    h[k] = (float)r[k] / tempHH;
                }
            }
            break;
        }
        else
        {
            for (m = j + 1; m <bands; m++)
            {
                float temp1 = ((float)r[j] - (float)r[i]) / (w[j] - w[i]);
                float temp2 = (float)r[i] - temp1*w[i];
                float y = temp1 * w[m] + temp2;
                if (y >= r[m])
                {
                    if (m == bands - 1)
                    {
                        //从i-j循环
                        for (k = i; k <= j; k++)
                        {
                            tempHH = temp1*w[k] + temp2;
                            if (tempHH == 0)
                            {
                                h[k] = 0.0f;
                            }
                            else
                            {
                                h[k] = (float)r[k] / tempHH;
                            }
                        }
                        i = j;
                        break;
                    }
                }
                else
                {
                    break;
                }
            }
        }
        j = j + 1;
    }
}