matlab 到 C++/openCV 归一化函数
matlab to C++/openCV normalization function
这是我的 matlab 代码:
imageData = imageData ./ toolbox.c3d.p.tprctile(imageData(xy),99.2);
imageData(imageData>1) = 1;
这里是我的openCV/c++代码,矩阵dst是一个openCV矩阵
cv::Mat dst
std::vector<float> result;
for (std::vector<int>::iterator it = index.begin() ; it != index.end(); ++it)
{
int ind = *it;
float temp = dst.at<float>(ind - 1);
result.push_back(temp);
}
float divider = tprctile(result,99.2);
dst = dst/ divider;
百分位数的效用函数
float Utils::tprctile(std::vector<float> channel, double pt)
{
std::sort(channel.begin(),channel.end());
int ptInd = Utilities::MatlabRound (pt/100 * channel.size() );
return channel[ptInd];
// Matlab code
// function val = tprctile(data, pt)
// data = sort(data);
// ptInd = round( pt/100 * length(data) );
// val = data(ptInd);
}
我的问题是关于 imageData(imageData>1) = 1
实现此功能的最有效方法是什么——我当然可以像以前一样遍历 dst。有没有更好的办法?
这就是我现在正在做的事情
int matrixSize = dst.rows *dst.cols;
cv::MatConstIterator_<float> it = dst.begin<float>(), it_end = dst.end<float>();
for(int i = 0 ; i < matrixSize ; ++i, ++it)
{
float value = *it;
if(value > 1.0)
{
dst.at<float>(i) = 1.0;
}
}
你想要的是用cv::threshold截断图像。
以下应该满足您的要求:
cv::threshold(dst, dst, 1, 1, CV_THRESH_TRUNC);
这将截断所有大于 1 的值并将结果存储在 dst 中。
这是我的 matlab 代码:
imageData = imageData ./ toolbox.c3d.p.tprctile(imageData(xy),99.2);
imageData(imageData>1) = 1;
这里是我的openCV/c++代码,矩阵dst是一个openCV矩阵
cv::Mat dst
std::vector<float> result;
for (std::vector<int>::iterator it = index.begin() ; it != index.end(); ++it)
{
int ind = *it;
float temp = dst.at<float>(ind - 1);
result.push_back(temp);
}
float divider = tprctile(result,99.2);
dst = dst/ divider;
百分位数的效用函数
float Utils::tprctile(std::vector<float> channel, double pt)
{
std::sort(channel.begin(),channel.end());
int ptInd = Utilities::MatlabRound (pt/100 * channel.size() );
return channel[ptInd];
// Matlab code
// function val = tprctile(data, pt)
// data = sort(data);
// ptInd = round( pt/100 * length(data) );
// val = data(ptInd);
}
我的问题是关于 imageData(imageData>1) = 1
实现此功能的最有效方法是什么——我当然可以像以前一样遍历 dst。有没有更好的办法?
这就是我现在正在做的事情
int matrixSize = dst.rows *dst.cols;
cv::MatConstIterator_<float> it = dst.begin<float>(), it_end = dst.end<float>();
for(int i = 0 ; i < matrixSize ; ++i, ++it)
{
float value = *it;
if(value > 1.0)
{
dst.at<float>(i) = 1.0;
}
}
你想要的是用cv::threshold截断图像。
以下应该满足您的要求:
cv::threshold(dst, dst, 1, 1, CV_THRESH_TRUNC);
这将截断所有大于 1 的值并将结果存储在 dst 中。