如圖所示，我提取的邊緣有很多鋸齒，想通過曲線擬合來解決，現在遇到幾個問題
1、我的曲線擬合法選擇的是貝塞爾方程曲線擬合，但是這個方法需要利用前一個點和后一個點來進行控制點的計算，也就是要求資料點要有順序，我用opencv中的findcontours函式獲得的該影像輪廓點集，并不知道有沒有順序。。。
2、在做曲線擬合前需要先提取角點么？
3、有其他對資料點的順序沒有要求的曲線擬合演算法么？
uj5u.com熱心網友回復：
可以試試5點三次平滑uj5u.com熱心網友回復：
僅供參考：#include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include "opencv2/imgproc/imgproc_c.h"
using namespace std;
using namespace cv;
Mat img,smallImg,gray,bw;
vector<Vec4i> hierarchy;
vector<vector<Point> > contours;
int threshval=128;
Rect r;
Rect maxrect,brect;
int idx,n;
const static Scalar colors[15]={
CV_RGB( 0, 0,128),
CV_RGB( 0,128, 0),
CV_RGB( 0,128,128),
CV_RGB(128, 0, 0),
CV_RGB(128, 0,128),
CV_RGB(128,128, 0),
CV_RGB(128,128,128),
CV_RGB(160,160,160),
CV_RGB( 0, 0,255),
CV_RGB( 0,255, 0),
CV_RGB( 0,255,255),
CV_RGB(255, 0, 0),
CV_RGB(255, 0,255),
CV_RGB(255,255, 0),
CV_RGB(255,255,255),
};
Scalar color;
void gamma_correct(Mat& img, Mat& dst, double gamma) {
Mat temp;
CvMat tmp;
img.convertTo(temp, CV_32FC1, 1.0/255.0, 0.0);
tmp=temp;
cvPow(&tmp, &tmp, gamma);
temp.convertTo(dst , CV_8UC1 , 255.0 , 0.0);
}
int main() {
cvNamedWindow("display",1);
img=imread("image.jpg",1);
r.x=img.cols/10;
r.y=img.rows/3;
r.width=img.cols*8/10;
r.height=img.rows*2/3;
smallImg=img(r);
cvtColor(smallImg,gray,CV_BGR2GRAY);
// medianBlur(gray,gray,5);
equalizeHist(gray,gray);
gamma_correct(gray,gray,4.0);
imshow("display",gray);
waitKey(0);
bw=(gray>threshval);
imshow("display",bw);
waitKey(0);
Mat Structure0=getStructuringElement(MORPH_RECT,Size(3,3));
erode(bw,bw,Structure0,Point(-1,-1));
Mat Structure1=getStructuringElement(MORPH_RECT,Size(6,6));
dilate(bw,bw,Structure1, Point(-1,-1));
imshow("display",bw);
waitKey(0);
findContours(bw,contours,hierarchy,RETR_EXTERNAL,CHAIN_APPROX_SIMPLE);
if (!contours.empty()&&!hierarchy.empty()) {
idx=0;
n=0;
vector<Point> approx;
for (;idx>=0;idx=hierarchy[idx][0]) {
color=colors[idx%15];
// drawContours(smallImg,contours,idx,color,1,8,hierarchy);
approxPolyDP(Mat(contours[idx]), approx, arcLength(Mat(contours[idx]), true)*0.005, true);//0.005為將毛邊拉直的系數
const Point* p = &approx[0];
int m=(int)approx.size();
polylines(smallImg, &p, &m, 1, true, color);
circle(smallImg,Point(p[0].x,p[0].y),3,color);
circle(smallImg,Point(p[1].x,p[1].y),2,color);
for (int i=2;i<m;i++) circle(smallImg,Point(p[i].x,p[i].y),1,color);
n++;
if (1==n) {
maxrect=boundingRect(Mat(contours[idx]));
} else {
brect=boundingRect(Mat(contours[idx]));
CvRect mr(maxrect),br(brect);
maxrect=cvMaxRect(&mr,&br);
}
}
circle(smallImg,Point(maxrect.x+maxrect.width/2,maxrect.y+maxrect.height/2),2,CV_RGB(255,0,0));
}
imshow("display",smallImg);
waitKey(0);
cvDestroyWindow("display");
return 0;
}
uj5u.com熱心網友回復：
五點三次平滑不是對等距點的平滑方式么？我這是一張圖，point的x，y都是亂的呀，不等距uj5u.com熱心網友回復：
貝塞爾方程與點的順序無關，你的影像呈現不光滑是因為你傳遞到貝塞爾方程中的變數少于三個，曲線必須高于二階貝塞爾方程。uj5u.com熱心網友回復：
我是做貝賽爾插值，怎么會和點的順序無關呢？插值程序中的控制點是需要根據每個點前后的資料點來求取的呀。傳遞到貝賽爾方程中的變數？
ABCD都是根據資料點求的呀
uj5u.com熱心網友回復：
你需要理解貝塞爾方程，你的錯誤是因為你不理解方程意義導致的。
uj5u.com熱心網友回復：
那你能給我講一下么？ 我的理解是：這是一個三階貝塞爾曲線公式，我只需要提供我要做插值的起始點（A）和終止點（D）以及三界貝塞爾曲線中所需的兩個控制點（B、C），控制點求取方式我是在網上查的。每兩個點之間的三階貝賽爾曲線只要將上面的ABCD代入公式即可，x用ABCD四點的X值，y用ABCD四點的Y值。這樣理論上應該是沒問題的呀。我剛又查了一下貝賽爾方程的推導程序，這幾個點的概念還是沒問題的呀？難道是我控制點求的不對？
uj5u.com熱心網友回復：
ABCD是4個點，大于3了呀
uj5u.com熱心網友回復：
opencv中的findcontours函式獲得的該影像輪廓點集是按照順序的uj5u.com熱心網友回復：
參考：拋物樣條曲線的原理說明及畫拋物曲線的一個類uj5u.com熱心網友回復：
宏觀的看你的圖形，感覺細節資料權重太大，導致大的區域不夠光滑。建議你，走走特殊路線。比如將每3個點的資料取平均，作為一個點資料使用。 具體操作幾個點效果好，你可以測驗。
uj5u.com熱心網友回復：
是這個道理
uj5u.com熱心網友回復：
...插值？uj5u.com熱心網友回復：
考慮下樣條插值(spline)?你把數值輸入excel看看效果如何uj5u.com熱心網友回復：
5點三次平滑 這樣可以嗎？uj5u.com熱心網友回復：
我覺得這個不用擬合吧，簡單做個濾波就夠了uj5u.com熱心網友回復：
說實話還是感覺你說的靠譜些，我之前是對插值點進行了篩選，大概是每8個點選一個進行顯示，然后確實光滑了不少，就是影像細節方面丟失的更多一點。看樣子暫時是找不到又有細節又沒有毛刺的方法了
uj5u.com熱心網友回復：
嗯嗯，樣條插值看起來比較復雜就沒選用，我都畢業了，以后看興趣再研究研究吧。uj5u.com熱心網友回復：
學習了,謝謝轉載請註明出處，本文鏈接：https://www.uj5u.com/gongcheng/65505.html
標籤：圖形處理/算法
下一篇：MFC上EXCEL操作問題