GEOS库的学习之二:简单几何图形的创建
几何图形(Geometry)是geos里面基本的操作对象,因此Geometry类就是最重要的一个类
几何图形中主要有三个要素:点,线,面。横纵坐标构成点,多个点构成线,环线构成面,点线面混合构成几何集合。对应的几个类为
坐标:Coordinate
点:Point、MultiPoint
线:LineString、MultiLineString(多条线)、LinearRing(环线)
面:Polygon、MultiPolygon
集合:GeometryCollection
在geos中,最小的组成单位是坐标,由Coordinate类来创建,如:Coordinate(2,3),表示一个点,横坐标是2,纵坐标是3.
所有的几何图形的创建,都是由类GeometryFactory来完成。因此,在创建几何图形前,可以先创建一个全局的GeometryFactory对象;
GeometryFactory factory; //全局对象,所有的图形都由此对象创建
1、点的创建
Point* createGeosPoint(double x,double y) { Coordinate pt(x,y); Point* p=factory.createPoint(pt); return p; }
2、非闭合线条的创建
LineString* createGeosLine(double x,double y, double offset) { CoordinateArraySequence *cas=new CoordinateArraySequence(); //构建点序列 cas->add(Coordinate(x,y)); cas->add(Coordinate(x,y+offset)); cas->add(Coordinate(x+offset,y+offset)); cas->add(Coordinate(x+offset,y+2*offset)); cas->add(Coordinate(x+2*offset,y+2*offset)); LineString *ls=factory.createLineString(cas); return ls; }
3、闭合线条的创建
//创建一条环线,与线的区别就是环线是闭合的。即第一个点和最后一点重合 LinearRing* createGeosRing(double x,double y,double offset) { CoordinateArraySequence *cas=new CoordinateArraySequence(); //构建点序列 cas->add(Coordinate(x,y)); cas->add(Coordinate(x,y+offset)); cas->add(Coordinate(x+offset,y+offset)); cas->add(Coordinate(x+offset,y+2*offset)); cas->add(Coordinate(x+2*offset,y+2*offset)); cas->add(Coordinate(x+2*offset,y)); cas->add(Coordinate(x,y)); //与第一个点相等 LinearRing *lr=factory.createLinearRing(cas); return lr; }
除了用add的方法来构建点序列,也可以用另外一种方法setAt
LinearRing* createGeosRing(double x,double y,double offset) { CoordinateArraySequenceFactory csf; CoordinateSequence* cs = csf.create(7,2); cs->setAt(Coordinate(x,y),0); cs->setAt(Coordinate(x,y+offset),1); cs->setAt(Coordinate(x+offset,y+offset),2); cs->setAt(Coordinate(x+offset,y+2*offset),3); cs->setAt(Coordinate(x+2*offset,y+2*offset),4); cs->setAt(Coordinate(x+2*offset,y),5); cs->setAt(Coordinate(x,y),6); //与第一个点相等 LinearRing *lr=factory.createLinearRing(cs); return lr; }
4、多边形的创建
//创建一个多边形,如果多边形内部没有孔洞实际上与环线是一样的 Polygon* createGeosPolygon(double x,double y,double offset) { LinearRing *lr=createGeosRing(x,y,offset); Polygon *poly=factory.createPolygon(lr,NULL); //如果多边形中间没有孔洞,第二个参数设为NULL return poly; }
测试:
#include "geos.h" int main() { LineString *ls=createGeosRing(10,10,5); cout<<"线条点数:"<<ls->getNumPoints()<<" 线条长度:"<<ls->getLength()<<endl; Polygon *poly=createGeosPolygon(10,10,5); cout<<"多边形面积:"<<poly->getArea()<<endl; system("pause"); return 1; }