OpenCV-消除凸性缺陷(Python到C++)
OpenCV - Removing convexity defects (Python to C++)
我一直在将其中一个答案的代码从Python移植到C++,我希望有人能帮助我完成这项工作,并澄清其中的一部分过程。
我说的答案是——https://stackoverflow.com/a/11366549
我正在对质心进行排序(一个在最后一个之前)。有问题的代码完成了这项工作,但我不确定如何在不延长代码的情况下使用C++获得相同的结果。
centroids = np.array(centroids,dtype = np.float32)
c = centroids.reshape((100,2))
c2 = c[np.argsort(c[:,1])]
b = np.vstack([c2[i*10:(i+1)*10][np.argsort(c2[i*10:(i+1)*10,0])] for i in xrange(10)])
bm = b.reshape((10,10,2))
我想实现这一点的方法是:
初始质心数组已经以某种方式进行了排序(最高的y值Point的索引为0,最低的为99),所以我想反转它,使数组从上到下排序。
之后,只需要对行进行排序(按x轴对10行10列进行排序)
我相信这就是所要做的一切,但我找不到一种优雅的方式来编码它(不使用向量和手动排序)。
我也不太理解最后一步,图像被扭曲(使其难以移植),所以如果有人能提供一些关于这一部分正在做的事情的见解,以及可能与此相当的C++,我将不胜感激。
output = np.zeros((450,450,3),np.uint8)
for i,j in enumerate(b):
ri = i/10
ci = i%10
if ci != 9 and ri!=9:
src = bm[ri:ri+2, ci:ci+2 , :].reshape((4,2))
dst = np.array( [ [ci*50,ri*50],[(ci+1)*50-1,ri*50],[ci*50,(ri+1)*50-1],[(ci+1)*50-1,(ri+1)*50-1] ], np.float32)
retval = cv2.getPerspectiveTransform(src,dst)
warp = cv2.warpPerspective(res2,retval,(450,450))
output[ri*50:(ri+1)*50-1 , ci*50:(ci+1)*50-1] = warp[ri*50:(ri+1)*50-1 , ci*50:(ci+1)*50-1].copy()
我只是在学习OpenCV和C++,我知道这可能是微不足道的,所以我希望有人能抽出一些时间,提供一些基本的答案。
编辑
根据请求,以下是代码。我已经处理了第一部分,仍然不确定这是否是正确的做法。
#pragma mark Correcting the defects
findContours(res, contours, RETR_LIST, CHAIN_APPROX_SIMPLE);
vector<cv::Point> centroids;
centroids.reserve(100);
for (int i = 0; i < contours.size(); i++) {
vector<cv::Point> contour = contours.at(i);
Moments mom = moments(contour);
cv::Point center = cv::Point(int(mom.m10 / mom.m00), int(mom.m01 / mom.m00));
centroids.push_back(center);
}
std::reverse(centroids.begin(), centroids.end());
struct {
bool operator()(const cv::Point p1, const cv::Point p2) {
return p1.x < p2.x;
}
} pointXGreater;
for (int i = 0; i < 10; i++) {
std::sort(centroids.begin() + i * 10, centroids.begin() + (i + 1) * 10, pointXGreater);
}
Mat b(centroids);
Mat bm = b.reshape(10, 10);
在用索引在图像上绘制质心后,它们似乎是正确的。
现在,我被困在最后一部分,试图破译Python代码,却不知道它能做什么。到目前为止,我得到了这个:
Mat output = Mat::zeros(450, 450, CV_8U);
for (int i = 0; i < centroids.size(); i++) {
cv::Point j = centroids.at(i);
int ri = i / 10;
int ci = i % 10;
if (ci != 9 && ri != 9) {
Mat src = ??
Mat dst = ??
Mat retval = getPerspectiveTransform(src, dst);
Mat warp;
warpPerspective(res2, warp, retval, (450, 450));
Mat output = ??
}
}
我会继续努力理解它,但如果你真的不知道在做什么,我会非常感谢你的帮助,因为这种Python语法不是很友好。
经过一段时间的努力,这就是我想出的
#pragma mark Correcting the defects
findContours(res, contours, RETR_LIST, CHAIN_APPROX_SIMPLE);
vector<cv::Point> centroids;
centroids.reserve(100);
for (int i = 0; i < contours.size(); i++) {
vector<cv::Point> contour = contours.at(i);
Moments mom = moments(contour);
cv::Point center = cv::Point(int(mom.m10 / mom.m00), int(mom.m01 / mom.m00));
centroids.push_back(center);
}
std::reverse(centroids.begin(), centroids.end());
struct {
bool operator()(const cv::Point p1, const cv::Point p2) {
return p1.x < p2.x;
}
} pointXGreater;
for (int i = 0; i < 10; i++) {
std::sort(centroids.begin() + i * 10, centroids.begin() + (i + 1) * 10, pointXGreater);
}
Mat bm = Mat(centroids);
bm = bm.reshape(2, 10);
Mat output(450, 450, CV_8UC3, CV_RGB(1, 1, 1));
for (int i = 0; i < centroids.size(); i++) {
int ri = i / 10;
int ci = i % 10;
if (ci != 9 && ri != 9) {
cv::Point2f src[4];
src[0] = cv::Point2f(bm.at<cv::Point>(ri, ci).x, bm.at<cv::Point>(ri, ci).y);
src[1] = cv::Point2f(bm.at<cv::Point>(ri, ci + 1).x, bm.at<cv::Point>(ri, ci + 1).y);
src[2] = cv::Point2f(bm.at<cv::Point>(ri + 1, ci).x, bm.at<cv::Point>(ri + 1, ci).y);
src[3] = cv::Point2f(bm.at<cv::Point>(ri + 1, ci + 1).x, bm.at<cv::Point>(ri + 1, ci + 1).y);
cv::Point2f dst[4];
dst[0] = cv::Point2f(ci * 50, ri * 50);
dst[1] = cv::Point2f((ci + 1) * 50 - 1, ri * 50);
dst[2] = cv::Point2f(ci * 50, (ri + 1) * 50 - 1);
dst[3] = cv::Point2f((ci + 1) * 50 - 1, (ri + 1) * 50 - 1);
Mat retval = getPerspectiveTransform(src, dst);
Mat warp;
warpPerspective(res2, warp, retval, cv::Size(450, 450));
for (int j = ri * 50; j < (ri + 1) * 50 - 1; j++) {
for (int k = ci * 50; k < (ci + 1) * 50 - 1; k++) {
output.at<Vec3b>(j, k) = warp.at<Vec3b>(j, k);
}
}
}
}
相关文章:
- 如何运行位于boost/libs/python/example/tutorial目录中的hello.cpp和Jamfil
- Pybind11:将元组列表从Python传递到C++
- 如何在c++中使用引用实现类似python的行为
- 是否可以通过C++扩展强制多个python进程共享同一内存
- 递归列出所有目录中的C++与Python与Ruby的性能
- IPC使用多个管道和分支进程来运行Python程序
- 从python中调用C++函数并获取返回值
- Python 3.7 和 excess_args 的 SWIG 问题
- Python中的for循环与C++有何不同
- 使用Pybind11向Python公开Eigen::张量
- Python str to C++ to Python str
- 如何使用Python从C++中读取谷物序列化数据
- 如何在C++中使用pybind11加载一个pickle python列表
- 如何在c++中使用system()来运行包含空格的python脚本
- python集合的C++等价物是什么.计数器
- 如果C++对象的类在另一个boost模块中声明,如何使用boost将指向该对象的指针返回到python
- 从python调用openMP共享库时,未定义opnMP函数
- 使用JsonCpp将数据返回到带有pybind11的python会在python调用中产生Symbol not foun
- 如何将真正的字符串从python c-api转换为python脚本
- Python ctype 'c_char_p' Memory Leak