Geos判断点是否在多边形内

使用的geo版本是3.5.1

 #include <iostream>
 #include "geos.h"
 using namespace std;
 GeometryFactory factory;
 Point* createGeosPoint(double x,double y)
 {
     Coordinate pt(x,y);
     Point* p=factory.createPoint(pt);
     return p;
 }
 Polygon* createGeosPolygon(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+*offset));
     cas->add(Coordinate(x+*offset,y+*offset));
     cas->add(Coordinate(x+*offset,y));
     cas->add(Coordinate(x,y)); //与第一个点相等，形成闭合
     LinearRing *lr=factory.createLinearRing(cas);//构建闭合环线
     Polygon *poly=factory.createPolygon(lr,NULL); //把闭合线转成多边形；如果多边形中间没有孔洞，第二个参数设为NULL
     return poly;
 }
 int main()
 {
     cout<<"GEOS库版本为："<<GEOS_VERSION<<endl;
     double x=;
     double y=;
     Polygon *poly=createGeosPolygon(,,);
     Point* p = createGeosPoint(x,y);
     if(poly->contains(p))//判断点p是否在多边形poly内，与边缘重合不算在内
     {
         cout<<"p在多边形内"<<endl;
     }
 }

遇到的问题：提示GeometryFactory是protected型，无法调用
解决办法：  一开始版本是3.6.2，将版本降低到3.5.1即可。