变化监测NR算子 Opencv C++代码 附注释

变化检测NR算子：NR算子+Kmeans

#include <opencv2\highgui\highgui.hpp> // imread imshow   waitKey 函数的申明 图形用户功能

#include "opencv2/core/core.hpp"  // 矩阵计算、数据结构功能
#include"opencv2/imgproc/imgproc.hpp"  //图像处理模块 img process
#include<iostream>
using namespace std;
using namespace cv;


int main(int argc, char* argv[]) // argc 代表主函数在命令窗口的参数个数，argv代表参数数组，argv[0]是程序名
{
Scalar mean;
Scalar stddev;
Mat img1 = imread("1.bmp", 0);// 0 代表返回灰度图像 1代表返回彩色图像,可自行将3通道图像变成灰度图像
Mat img2 = imread("2.bmp", 0);
Mat img3 = imread("3.bmp", 0);// 参考用图像
//int c = img1.channels();
int k = img1.dims;
int m = img1.rows;
int n = img1.cols;
Mat img_save = img1;


// 参数设置 窗口大小
int w_mean_var = 5;
int w_half = (w_mean_var - 1) / 2;


//图像均值与方差求解
Mat img1_mean;
blur(img1, img1_mean, Size(w_mean_var, w_mean_var), Point(-1, -1), BORDER_REPLICATE);
Mat img2_mean;
blur(img2, img2_mean, Size(w_mean_var, w_mean_var), Point(-1, -1), BORDER_REPLICATE);

Mat image1_var = Mat_<double>(img1); //转换数据类型,方差阵
Mat image2_var = Mat_<double>(img2);




//Mat patch1 = img1(Range(5, 9), Range(5, 9));// 新对象与原始对象共享同一个数据区域 是左闭右开区间，不包括右max
//Mat patch2 = img1(Range(5,9), Range(5, 9));
//cout << "patch1_p = " << patch1 << endl;