c++ - 在 opencv 中创建 NDWI 矩阵

Question

首先，如果你不知道，我应该告诉你什么是ndwi。 Ndwi 代表归一化差异水指数。水的图形指示器，取值范围为[-1 1]。 Ndwi定义如下:

(Green - NIR) / (Green + NIR)

我正在开发一个基于 opencv 的简单海岸线提取工具。我在 MATLAB 中完成了它，结果如下所示:

然而，opencv 版本的结果看起来像二值化的:

当我调试程序时，我看到 ndwi 矩阵中的最小值为零，这是错误的，因为它应该是 -0.8057。负责ndwi计算的代码(opencv版)如下:

Mat ndwi = (greenRoi - nirRoi) / (greenRoi + nirRoi);
double min;
double max;
minMaxIdx(ndwi, &min, &max);
Mat adjNDWI;
convertScaleAbs(ndwi, adjNDWI, 255 / max);

这里有什么问题，我怎样才能计算出正确的 ndwi 值？

注意:

greenRoi 和 nirRoi 是这样创建的:

Rect rectangle = boundingRect(Mat(testCorners)); //vector<Point2f> testCorners(4);

Mat testImgGreen = imread((LPCSTR)testImgGreenPath, 0);
Mat testImgNir = imread((LPCSTR)testImgNirPath, 0);

Mat greenRoi(testImgGreen, rectangle);
Mat nirRoi(testImgNir, rectangle);

Answer 1

您需要显式创建一个浮点 cv::Mat

cv::Mat image(cols,rows,CV_32FC1) 或者 CV_64FC1 如果你需要 double