openCVを使ったハフ変換で角度を限定して検出したい
Posted: 2015年7月03日(金) 12:00
はじめまして。
現在、c++でopenCVを使用してハフ変換である範囲の角度の直線のみを検出しようとしています。
OSはwindows7、ソフトはvisual studio 2010を使用しています。
私がやった方法では一度すべての直線を検出した後、
描画する前に直線の角度を見て、指定した範囲の角度だった場合に描画します。
直線を検出することはできましたが、角度の限定がうまくいきません。
検出する直線の角度を限定できればよかったのですが、
方法がわからずこのような方法になってしまいました。
私が教えていただきたいのは
・ある範囲の角度の直線のみを検出する方法
・検出した直線の角度を判別し、描画するかしないかを判断する方法
です。
以下はソースコードです。
よろしくお願いいたします。
現在、c++でopenCVを使用してハフ変換である範囲の角度の直線のみを検出しようとしています。
OSはwindows7、ソフトはvisual studio 2010を使用しています。
私がやった方法では一度すべての直線を検出した後、
描画する前に直線の角度を見て、指定した範囲の角度だった場合に描画します。
直線を検出することはできましたが、角度の限定がうまくいきません。
検出する直線の角度を限定できればよかったのですが、
方法がわからずこのような方法になってしまいました。
私が教えていただきたいのは
・ある範囲の角度の直線のみを検出する方法
・検出した直線の角度を判別し、描画するかしないかを判断する方法
です。
以下はソースコードです。
よろしくお願いいたします。
#include <opencv2/core/core.hpp>
#include <opencv2/imgproc/imgproc.hpp>
#include <opencv2/highgui/highgui.hpp>
#include <opencv2/opencv.hpp>
#include <opencv2/opencv_lib.hpp>
#include <math.h>
#define PI 3.141592653589793
int main(int argc, char *argv[])
{
int i1;
char suuzi[249];
FILE *f;
int i=0;
char filename01[100];
char filename02[100];
for(i=1;i<240;i++){
sprintf(filename01,"00[%d].bmp",i);
sprintf(filename02,"out[%d].bmp",i);
//f=fopen(filename,"r");
//fclose(f);
cv::Mat src_img = cv::imread(filename01);//Canny1000.jpg
if(!src_img.data) return -1;
cv::Mat dst_img, work_img;
dst_img = src_img.clone();
cv::cvtColor(src_img, work_img, CV_RGB2GRAY);
cv::Canny(work_img, work_img, 50, 200, 3);
// （古典的）Hough変換
std::vector<cv::Vec2f> lines;
// 入力画像，出力，距離分解能，角度分解能，閾値，*,*
cv::HoughLines(work_img, lines, 1, CV_PI/180, 120, 0, 0);
std::vector<cv::Vec2f>::iterator it = lines.begin();
for(; it!=lines.end(); ++it) {
float rho = (*it)[0], theta = (*it)[1];
cv::Point pt1, pt2;
double x1,x2,y1,y2;
double a = cos(theta), b = sin(theta);
double x0 = a*rho, y0 = b*rho;
double x3,y3;
double tantan;
int amari,amo,c,d;
pt1.x= cv::saturate_cast<int>(x0 + 1000*(-b));
pt1.y = cv::saturate_cast<int>(y0 + 1000*(a));
pt2.x = cv::saturate_cast<int>(x0 - 1000*(-b));
pt2.y = cv::saturate_cast<int>(y0 - 1000*(a));
x1=pt1.x;
y1=pt1.y;
x2=pt2.x;
y2=pt2.y;
if(x1<x2){
x3=x1-x2;
}else{
x3=x2-x1;
}
if(y1>y2){
y3=y1-y2;
}else{
y3=y2-y1;
}
tantan=atan2(y3,x3);
amo=fabs(tantan);
d=amo*180/PI;
amari=d%360;
printf("//////////////////////\n");
printf("x1\n");
printf("%f\n",x1);
printf("y1\n");
printf("%f\n",y1);
printf("x2\n");
printf("%f\n",x2);
printf("y2\n");
printf("%f\n",y2);
printf("x3\n");
printf("%f\n",x3);
printf("y3\n");
printf("%f\n",y3);
printf("tantan\n");
printf("%f\n",tantan);
printf("amo\n");
printf("%d\n",amo);
printf("d\n");
printf("%d\n",d);
printf("amari\n");
printf("%d\n",amari);
if(0<amari&&amari<360){
cv::line(dst_img, pt1, pt2, cv::Scalar(0,0,255), 3, CV_AA);
}
}
cv::namedWindow("HoughLines", CV_WINDOW_AUTOSIZE|CV_WINDOW_FREERATIO);
//cv::imshow("HoughLines", dst_img);
cv::imwrite(filename02,dst_img);
//cv::waitKey(0);
}
cv::waitKey(0);
}