根据高度图计算体积等。

根据高度图计算体积等。

1、2017年8月2日14:31:18。

//计算高度图的面积(mm^2)、体积(mm^3)、平均高度(mm)。
//输入高度图：高度单位mm。
//每像素18.3um=0.0183mm。
int CalVolumeByHeightMat(const cv::Mat & HeightMat,double & dArea, double & dVolume, double & dAveHeight, double dMmPerPixel);

//计算高度图的面积(mm^2)、体积(mm^3)、平均高度(mm)。
//输入高度图：高度单位mm。
//输入掩码矩阵.
//每像素18.3um=0.0183mm。
int CalVolumeByHeightSrcMat(const cv::Mat & HeightSrcMat, const cv::Mat & MaskMat, double & dArea, double & dVolume, double & dAveHeight);

int TestCalVolumeByHeightSrcMat(const cv::Mat & HeightSrcMat, double & dArea, double & dVolume, double & dAveHeight);

int CalVolumeByHeightMat(const cv::Mat & HeightMat, double & dArea, double & dVolume, double & dAveHeight, double dMmPerPixel)
{
	int nCount = countNonZero(HeightMat);
	if (nCount <= 0)
	{
		dArea = 0;
		dVolume = 0;
		dAveHeight = 0;
		return 1;
	}
	cv::Scalar  vSum = sum(HeightMat);
	double  dTmpVolume = vSum[0];
	dAveHeight = dTmpVolume / nCount;
	dArea = nCount * dMmPerPixel * dMmPerPixel;
	dVolume = dTmpVolume * dMmPerPixel * dMmPerPixel;
	return 1;
}

int CalVolumeByHeightSrcMat(const cv::Mat & HeightSrcMat, const cv::Mat & MaskMat, double & dArea, double & dVolume, double & dAveHeight)
{
	int nHeightRows = HeightSrcMat.rows;
	int nHeightCols = HeightSrcMat.cols;

	int nMaskRows = MaskMat.rows;
	int nMaskCols = MaskMat.cols;

	cv::Mat HeightMat;
	if ((nHeightRows == nMaskRows) && (nHeightCols == nMaskCols))
	{
		HeightSrcMat.copyTo(HeightMat, MaskMat);
	}
	else
	{
		HeightSrcMat.copyTo(HeightMat);
	}

	//每像素18.3um=0.0183mm。
	double dMmPerPixel = 0.0183;

	CalVolumeByHeightMat( HeightMat, dArea, dVolume, dAveHeight, dMmPerPixel);
	return 1;
}

int TestCalVolumeByHeightSrcMat(const cv::Mat & HeightSrcMat, double & dArea, double & dVolume, double & dAveHeight)
{
	//cv::Mat MaskMat(HeightSrcMat.size(), CV_8UC1, cv::Scalar::all(255));
	cv::Mat MaskMat(HeightSrcMat.size(), CV_8UC1, cv::Scalar::all(0));

	cv::Rect ROIRect = cv::Rect(0, 0, 15,8);
	cv::Mat ROIMat = MaskMat(ROIRect);
	ROIMat.setTo(cv::Scalar::all(255));

	CalVolumeByHeightSrcMat(HeightSrcMat, MaskMat, dArea,dVolume, dAveHeight);
	return 1;
}

	cv::Mat mat = cv::imread(strFileName, cv::IMREAD_GRAYSCALE);

	double minVal = 0;
	double maxVal = 0;
	cv::minMaxLoc(mat, &minVal,&maxVal);

	const int channels[1] = { 0 };
	const int histSize[1] = { 256 };
	float hranges[2] = { 0,256 };
	const float* ranges[1] = { hranges };
	cv::MatND hist;
	cv::calcHist(&mat, 1, channels, cv::Mat(), hist, 1, histSize, ranges);
	int nNumOf0 = hist.at<float>(0);
	int nNumOfAll = mat.rows * mat.cols;
	hist.at<float>(0) = 0;
	cv::Scalar scalar = cv::sum(hist);
	int nNumOf1 = scalar[0];

	double dArea = 0;
	double dVolume = 0;
	double dAveHeight = 0;
	TestCalVolumeByHeightSrcMat(mat, dArea, dVolume, dAveHeight);

说明：

1、根据高度图计算体积、面积、平均高度。

2、熟悉cv::Mat的部分操作。

3、熟悉直方图使用。

4、熟悉ROI使用。