程序运行出现vector subscript out of range
#include"opencv2/core/core.hpp"
#include"opencv2/highgui/highgui.hpp"
#include"opencv2/imgproc/imgproc.hpp"
#include"iostream"
#include"vector"
using namespace std;
using namespace cv;
int getA(Mat& src, Mat &dark)
{
int row=src.rows;
int col=src.cols;
vector<uchar>list;
for(int i=0;i<row;++i)
{
for(int j=0;j<col;++j)
{
list.push_back(dark.at<uchar>(i,j));
}
}
sort(list.begin(),list.end());
int thresh=list[0.999*row*col];
vector<uchar>list1;
int k=0;
for(int i=0;i<row;++i)
{
for(int j=0;j<col;++j)
{
if(dark.at<uchar>(i,j)>thresh)
list1.push_back(src.at<uchar>(i,j));
++k;
}
}
int Sum=0;
for(int i=0;i<k;++i)
{
Sum=Sum+list1[i];
}
int A;
A=Sum/k;
return A;
}
Mat getMax(Mat& src)
{
Mat image(src.size(),CV_8UC1);
int row=src.rows;
int col=src.cols;
for(int i=0;i<row;i++ )
{
for(int j=0;j<col;j++)
{
Vec3b intensity=src.at<Vec3b>(i,j);
image.at<uchar>(i,j)=intensity.val[0];
if(image.at<uchar>(i,j)<intensity.val[1])
image.at<uchar>(i,j)=intensity.val[1];
else if(image.at<uchar>(i,j)<intensity.val[2])
image.at<uchar>(i,j)<intensity.val[2];
}
}
return image;
}
int main()
{
Mat srcImage=imread("G:\\study resource\\underwater_image\\haze18.jpg");
if(srcImage.empty())
return -1;
imshow("srcGray",srcImage);
Mat dstImage;
std::vector<Mat>vInts(3);
split(srcImage,vInts);
imshow("dstImageB",vInts[0]);
imshow("dstImageG",vInts[1]);
imshow("dstImageR",vInts[2]);
int row,col;
row=srcImage.rows;
col=srcImage.cols; //get dark image
Mat dark(srcImage.size(),CV_8UC1);
for(int i=0;i<row;i++ )
{
for(int j=0;j<col;j++)
{
Vec3b intensity=srcImage.at<Vec3b>(i,j);
dark.at<uchar>(i,j)=intensity.val[0];
if(dark.at<uchar>(i,j)>intensity.val[1])
dark.at<uchar>(i,j)=intensity.val[1];
else if(dark.at<uchar>(i,j) > intensity.val[2])
dark.at<uchar>(i,j) < intensity.val[2];
}
}
imshow("dark",dark); // end
Mat dst(dark.size(),CV_8UC1); //dark image经过最小值滤波后图像begin
int radius=7; //窗口半径
vector<uchar>vlist;
for(int i=radius;i<row-radius;++i)
{
for(int j=radius;j<col-radius;++j)
{
for(int pi=i-radius;pi<=i+radius;pi++)
{
for(int pj=j-radius;pj<=j+radius;pj++)
{
vlist.push_back(dark.at<uchar>(pi,pj));
}
}
sort(vlist.begin(),vlist.end());
dst.at<uchar>(i,j)=vlist[0];
vlist.clear();
}
} //end
//求环境光A
Mat Aimage=getMax(srcImage);
int A=getA(Aimage,dark);
cout<<"a"<<A<<endl;
imshow("dst",dst);
waitKey(0);
return 0;
}