GEOS库学习之四:几何关系判断
2024-10-19 05:22:14
原理上一篇已经介绍过了,这篇就直接进行程序练习
#include "geos.h" GeometryFactory factory; //创建一条环线,与线的区别就是环线是闭合的。即第一个点和最后一点重合
LinearRing* createGeosRing(double x,double y,double offset)
{
CoordinateArraySequenceFactory csf;
CoordinateSequence* cs = csf.create(,);
cs->setAt(Coordinate(x,y),);
cs->setAt(Coordinate(x,y+offset),);
cs->setAt(Coordinate(x+offset,y+offset),);
cs->setAt(Coordinate(x+*offset,y+*offset),);
cs->setAt(Coordinate(x+*offset,y),);
cs->setAt(Coordinate(x,y),); //与第一个点相等
LinearRing *lr=factory.createLinearRing(cs);
return lr;
} //创建一个多边形,如果多边形内部没有孔洞实际上与环线是一样的
Polygon* createGeosPolygon(double x,double y,double offset)
{
LinearRing *lr=createGeosRing(x,y,offset);
Polygon *poly=factory.createPolygon(lr,NULL); //如果多边形中间没有孔洞,第二个参数设为NULL
return poly;
} int main()
{
Polygon *p1=createGeosPolygon(,,); //创建第一个多边形 for(int i=;i<=;i++)
{
cout<<i<<": ";
Polygon *p2=createGeosPolygon(,,i); //创建第二个多边形
IntersectionMatrix *im=p2->relate(p1);
cout<<*im<<" "; //返回DE-9IM交叉矩阵
if(p2->disjoint(p1))
cout<<"不相交"<<endl;
else
{
if(p2->touches(p1))
cout<<"接触"<<endl;
else if(p2->overlaps(p1))
cout<<"部分重叠"<<endl;
else if(p2->covers(p1))
cout<<"覆盖"<<endl;
else
cout<<*im<<endl;
}
}
system("pause");
return ;
}
结果如下:
最新文章
- 转:MYSQL连接字符串参数解析(解释)
- 对Java中正则表达式的一些理解
- JS实现验证码倒计时效果
- ahjesus 捕获entity framework生成的sql语句
- SU Demos-03T-F Analysis-01Sugabor
- Maven仓库的理解及配置
- iOS - UITextView
- php运行步骤解析
- UVA 11078 Open Credit System(扫描 维护最大值)
- BZOJ 1600: [Usaco2008 Oct]建造栅栏
- Fatal error in launcher:Unable to create process using &#39;";&#39;
- 【docker简易笔记】docker基础信息的分享
- C语言实现快速翻转数组的顺序
- [c/c++] programming之路(30)、位运算(一)
- C++入门程序作业2
- HTML常用标签及属性
- 附9 elasticsearch-curator + Linux定时任务
- HP管理工具System Management Homepage安装配置
- 第五篇:Hadoop流
- css - 文字元素等的美化效果代码汇总(更新中...)