python怎么判断两条线段相交

Python012

python怎么判断两条线段相交,第1张

先判断两条线段是否不平行(最好同时判断是否有交点并且不平行,因为浮点运算不精确),然后计算两条线段的交点。以下是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

}

1、在源码 python2.7.3/Modules/Setup.dist 搜索zlib ,打开该模块

2、继续下去,又提示缺少zlib.h,也就是我们的交叉编译库中缺少zlib头文件和库文件

3、在http://www.zlib.net/ 下载zlib源码,交叉编译zlib

由于zlib 编译不支持--host选项,所以在zlib目录下执行如下命令

~/zlib1.2.8/# mkdir _install

~/zlib1.2.8/#CC=arm-xxxxx(你的交叉编译器) ./configure --host=./_install(最好是绝对路径,我偷懒了,这个是安装目录)

~/zlib1.2.8/# make &&make install

你会在_install目录下发现include 和 lib 目录

下一步,我们回到第1步

为了防止交叉编译链污染,我们在python目录下创建otherinclude otherlib

把zlib下的include和lib下的内容分别按照对应名字拷贝过来

4.在Modules/setup 找到我们前面setup.dist那句话,添加 -L../otherlib -I../otherinclude 这个应该都懂的

直接make

5,如果你之前make distclean过了,没有setup文件,请修改setup.dist.文件,重新configure!

1. 简介。

图像处理是一门应用非常广的技术,而拥有非常丰富第三方扩展库的 Python 当然不会错过这一门盛宴。PIL (Python Imaging Library)是 Python 中最常用的图像处理库,目前版本为 1.1.7,我们可以 在这里 下载学习和查找资料。

Image 类是 PIL 库中一个非常重要的类,通过这个类来创建实例可以有直接载入图像文件,读取处理过的图像和通过抓取的方法得到的图像这三种方法。

2. 使用。

导入 Image 模块。然后通过 Image 类中的 open 方法即可载入一个图像文件。如果载入文件失败,则会引起一个 IOError ;若无返回错误,则 open 函数返回一个 Image 对象。现在,我们可以通过一些对象属性来检查文件内容,即:

1 >>>import Image

2  >>>im = Image.open("j.jpg")

3  >>>print im.format, im.size, im.mode

4 JPEG (440, 330) RGB

这里有三个属性,我们逐一了解。

format : 识别图像的源格式,如果该文件不是从文件中读取的,则被置为 None 值。

size : 返回的一个元组,有两个元素,其值为象素意义上的宽和高。

mode : RGB(true color image),此外还有,L(luminance),CMTK(pre-press image)。

现在,我们可以使用一些在 Image 类中定义的方法来操作已读取的图像实例。比如,显示最新载入的图像:

1 >>>im.show()

2  >>>

输出原图:

3. 函数概貌。

3.1    Reading and Writing Images : open( infilename ) , save( outfilename )

3.2    Cutting and Pasting and Merging Images :

crop() : 从图像中提取出某个矩形大小的图像。它接收一个四元素的元组作为参数,各元素为(left, upper, right, lower),坐标系统的原点(0, 0)是左上角。

paste() :

merge() :

1 >>>box = (100, 100, 200, 200)

2  >>>region = im.crop(box)

3  >>>region.show()

4  >>>region = region.transpose(Image.ROTATE_180)

5  >>>region.show()

6  >>>im.paste(region, box)

7  >>>im.show()

其效果图为:

旋转一幅图片:

1 def roll(image, delta):

2     "Roll an image sideways"

3

4     xsize, ysize = image.size

5

6     delta = delta % xsize

7     if delta == 0: return image

8

9     part1 = image.crop((0, 0, delta, ysize))

10     part2 = image.crop((delta, 0, xsize, ysize))

11     image.paste(part2, (0, 0, xsize-delta, ysize))

12     image.paste(part1, (xsize-delta, 0, xsize, ysize))

13

14     return image

3.3    几何变换。

3.3.1    简单的几何变换。

1 >>>out = im.resize((128, 128))                     #

2  >>>out = im.rotate(45)                             #逆时针旋转 45 度角。

3  >>>out = im.transpose(Image.FLIP_LEFT_RIGHT)       #左右对换。

4  >>>out = im.transpose(Image.FLIP_TOP_BOTTOM)       #上下对换。

5  >>>out = im.transpose(Image.ROTATE_90)             #旋转 90 度角。

6  >>>out = im.transpose(Image.ROTATE_180)            #旋转 180 度角。

7 >>>out = im.transpose(Image.ROTATE_270)            #旋转 270 度角。

各个调整之后的图像为:

图片1:

图片2:

图片3:

图片4:

3.3.2    色彩空间变换。

convert() : 该函数可以用来将图像转换为不同色彩模式。

3.3.3    图像增强。

Filters : 在 ImageFilter 模块中可以使用 filter 函数来使用模块中一系列预定义的增强滤镜。

1 >>>import ImageFilter

2 >>>imfilter = im.filter(ImageFilter.DETAIL)

3 >>>imfilter.show()

3.4    序列图像。

即我们常见到的动态图,最常见的后缀为 .gif ,另外还有 FLI / FLC 。PIL 库对这种动画格式图也提供了一些基本的支持。当我们打开这类图像文件时,PIL 自动载入图像的第一帧。我们可以使用 seek 和 tell 方法在各帧之间移动。

1 import Image

2 im.seek(1)        # skip to the second frame

3

4 try:

5     while 1:

6         im.seek( im.tell() + 1)

7         # do something to im

8 except EOFError:

9     pass

3.5    更多关于图像文件的读取。

最基本的方式:im = Image.open("filename")

类文件读取:fp = open("filename", "rb")im = Image.open(fp)

字符串数据读取:import StringIOim = Image.open(StringIO.StringIO(buffer))

从归档文件读取:import TarIOfp = TarIo.TarIO("Image.tar", "Image/test/lena.ppm")im = Image.open(fp)

基本的 PIL 目前就练习到这里。其他函数的功能可点击 这里 进一步阅读。