opencv编程中求近似两条平行线的中线，但老是出现内存出错的情况，特求助于大侠了

Posted 2023-03-31

#include "stdafx.h"
#include <stdio.h>
#include <math.h>
#include <cv.h>
#include <stdlib.h>
#include <highgui.h>
#include <cxcore.h>
int main(int argc, char*argv[])

IplImage* pImg=0; //声明IplImage指针
IplImage* image=0;
long int i,j,m=0,n=0;
uchar* imgArray;
uchar* array;
int col[387];//定义数组
for(i=0;i<388;i++)
col[i]=0;//数组的初始化
int row[507];
for(j=0;j<508;j++)
row[j]=0;
image = cvLoadImage("1.tif",0); //载入图像
cvNamedWindow( "Image", CV_WINDOW_AUTOSIZE);//创建窗口
IplImage * img = cvCreateImage( cvSize(194,254), IPL_DEPTH_8U, 3);//创建一个IplImage对象,尺寸100*100，字节型，3通道
cvZero(img);//将像素值全部置为0
imgArray=(uchar*)image->imageData;//图像数据的指针
printf("begin ");
for (i=0;i<image->height;i++)//i<image->height

array=imgArray+image->widthStep*i;//当前行的指针
for (j=0;j<image->width;j++)//j<image->width

int k=0;
if(*array>10)
row[m]=i;//将目标像素的坐标储存在数组里面
m++;
col[n]=j;
n++;
printf(" %d*%d ",i,j);
//物体近似为两条平行的线段，要求出这两条线的中轴线，只要提取的坐标的纵坐标相同，
//让横坐标相加除以2即可,然后在创建的背景中描出这些点即可
if((m>=1)&&(row[m]==row[m-1])&&(m<507))
k=(col[n-1]+col[n]);
//printf("%d*%d ",m,k);
cvCircle(img,cvPoint(m,k),1 , cvScalar(255,255,255,0.5), 1, 8, 0);//标出计算的点

//
array++;

printf("end ");
cvShowImage( "Image",image ); //显示图像
cvNamedWindow("show an image", CV_WINDOW_AUTOSIZE);
cvShowImage("show an image", img);
cvWaitKey(0); //等待按键

cvReleaseImage(&img);
cvDestroyWindow("show an image");
cvReleaseImage( &image ); //释放图像
cvDestroyWindow( "Image" );//销毁窗口
return -1;

1 hough变换获取直线，可以得到两条直线首尾两点坐标
2 根据坐标计算直线斜率，然后推算中轴线斜率
3 根据直线坐标计算它们的交点，然后有了交点和中轴线斜率，就出来结果了。追问

其实是两条曲线，不过这个方法挺好的，上面的程序已经调试好了。谢谢。

参考技术A IplImage * img = cvCreateImage( cvSize(194,254), IPL_DEPTH_8U, 3);//创建一个IplImage对象,尺寸100*100，字节型，3通道

这句话追问

运行的时候说内存冲突的是从
if(*array>10)
row[m]=i;//将目标像素的坐标储存在数组里面
m++;
col[n]=j;
n++;
printf(" %d*%d ",i,j);
开始的

追答

开始没越界，运行着运行发现越界了。

2013-11-16 23:41:58追加
不过你求平行线为什么不用hough检测出来直线然后将截距求均值呢?

追问

你是说求出直线的斜率和截距，然后求这两条线的截距，有斜率和截距就可以求出中线对吗？问题是这两条线并非是平行的，所以我想求出中轴线，从而得到中轴线倾斜角，最后求出这幅倾斜图像变为竖直图像要旋转的角度。也就是说最终是求旋转角的。不知道在这方面你有没有更好的想法。

追答

只要是直线就可以用hough做 很简单

两条线段间的位置 (叉积)

两条线段共有三种相对位置，平行，相交和重合，题目如果给出两条线段端点坐标，然后判断这两条线段的相对位置，如果是相交还要算出交点坐标。 
可以用叉积来做，将两条线段看做两条向量 
平行就是这两条向量的叉积为0 
重合就是他们首尾连接的向量与其中某条向量的叉积也为0 

 

//先判断两条线段是否不平行（最好同时判断是否有交点并且不平行，因为浮点运算不精确），然后计算两条线段的交点。以下是C语言代码：
#include<stdio.h>
#include<math.h>
#define eps 1e-8　　　　　　//精度，因为浮点数运算不精确
#define zero(x) (((x)>0?(x):-(x))<eps)
struct point
{
    double x,y;
};
 
//计算交叉乘积(P1-P0)x(P2-P0)
double xmult(point p1,point p2,point p0)
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
 
//判点是否在线段上,包括端点
int dot_online_in(point p,point l1,point l2)
{  
    return zero(xmult(p,l1,l2))&&(l1.x-p.x)*(l2.x-p.x)<eps&&(l1.y-p.y)*(l2.y-p.y)<eps;
}
 
//判两点在线段同侧,点在线段上返回0
int same_side(point p1,point p2,point l1,point l2)
{
    return xmult(l1,p1,l2)*xmult(l1,p2,l2)>eps;
}
 
//判两直线平行
int parallel(point u1,point u2,point v1,point v2)
{
    return zero((u1.x-u2.x)*(v1.y-v2.y)-(v1.x-v2.x)*(u1.y-u2.y));
}
 
//判三点共线
int dots_inline(point p1,point p2,point p3)
{
    return zero(xmult(p1,p2,p3));
}
 
//判两线段相交,包括端点和部分重合
int intersect_in(point u1,point u2,point v1,point v2)
{
    if (!dots_inline(u1,u2,v1)||!dots_inline(u1,u2,v2))
        return !same_side(u1,u2,v1,v2)&&!same_side(v1,v2,u1,u2);
    return dot_online_in(u1,v1,v2)||dot_online_in(u2,v1,v2)||dot_online_in(v1,u1,u2)||dot_online_in(v2,u1,u2);
}
 
//计算两线段交点,请判线段是否相交(同时还是要判断是否平行!)
point intersection(point u1,point u2,point v1,point v2)
{
    point ret=u1;
    double t=((u1.x-v1.x)*(v1.y-v2.y)-(u1.y-v1.y)*(v1.x-v2.x))
   /((u1.x-u2.x)*(v1.y-v2.y)-(u1.y-u2.y)*(v1.x-v2.x));
    ret.x+=(u2.x-u1.x)*t;
    ret.y+=(u2.y-u1.y)*t;
    return ret;
}
 
int main(void)
{
    point u1,u2,v1,v2,ans;
    printf("请输入线段1的两个端点:\n");
    scanf("%lf%lf%lf%lf",&u1.x,&u1.y,&u2.x,&u2.y);
    printf("请输入线段2的两个端点:\n");
    scanf("%lf%lf%lf%lf",&v1.x,&v1.y,&v2.x,&v2.y);
    if (parallel(u1,u2,v1,v2)||!intersect_in(u1,u2,v1,v2))
    {
        printf("无交点!\n");
    }
    else
    {
        ans=intersection(u1,u2,v1,v2);
        printf("交点为:(%lf,%lf)",ans.x,ans.y);
    }
    return 0;
}