geotools 代码实例

最新推荐文章于 2024-08-13 17:33:49 发布

kone666

最新推荐文章于 2024-08-13 17:33:49 发布

阅读量717

点赞数

分类专栏： GeoTools

GeoTools 专栏收录该内容

6 篇文章 6 订阅

订阅专栏

前段时间的一个项目本来用ae完成了种种的查询，空间分析等等功能的代码，但是不幸的是这是一个web项目，无奈 ae各种错误，显然ae放在server端是不好使的无奈一咬牙一跺脚全部换换成geotools 看文档看api 从零开始算是把原来AE实现的东西统统改了过来用起来反而觉得既稳定效率还不错哈！
以下是部分功能总结：

1、连接数据库这里使用的postgis 链接代码如下：

private static void conn(String dbtype, String host, String port,
            String database, String userName, String password) {
        Map<String, Object> params = new HashMap<String, Object>();
        // params.put(PostgisNGDataStoreFactory.DBTYPE.key, "postgis");    // 两种代码方式
        // params.put(PostgisNGDataStoreFactory.HOST.key, "localhost");
        // params.put(PostgisNGDataStoreFactory.PORT.key, new Integer(5432));
        // params.put(PostgisNGDataStoreFactory.DATABASE.key, "postgis");
        // params.put(PostgisNGDataStoreFactory.SCHEMA.key, "public");
        // params.put(PostgisNGDataStoreFactory.USER.key, "postgres");
        // params.put(PostgisNGDataStoreFactory.PASSWD.key, "root");
        params.put(PostgisNGDataStoreFactory.DBTYPE.key, dbtype);
        params.put(PostgisNGDataStoreFactory.HOST.key, host);
        params.put(PostgisNGDataStoreFactory.PORT.key, new Integer(port));
        params.put(PostgisNGDataStoreFactory.DATABASE.key, database);
        params.put(PostgisNGDataStoreFactory.SCHEMA.key, "public");
        params.put(PostgisNGDataStoreFactory.USER.key, userName);
        params.put(PostgisNGDataStoreFactory.PASSWD.key, password);
        try {
            pgDatastore = DataStoreFinder.getDataStore(params);
            if (pgDatastore != null) {
                System.out.println("系统连接到位于：" + host + "的空间数据库" + database
                        + "成功！");
            } else {
                System.out.println("系统连接到位于：" + host + "的空间数据库" + database
                        + "失败！请检查相关参数");
            }
        } catch (IOException e) {
            e.printStackTrace();
            System.out.println("系统连接到位于：" + host + "的空间数据库" + database
                    + "失败！请检查相关参数");
        }
 
    }

调用方法为：conn("postgis", "localhost", 5432, "postgis", "postgres", "root");
2、图层的操作

2.1 查询

public static ArrayList<SimpleFeature> queryMethod(String filterStr,
            String layerName) {
        //pgDatastore为上文连接数据库获取相当于AE中的workspace
        //SimpleFeatureSource相当于AE中的featureClass
        SimpleFeatureSource featureSource =pgDatastore.getFeatureSource(layerName); 
        ArrayList<SimpleFeature> featureList = new ArrayList<SimpleFeature>();
        if(featureSource==null)
            return featureList;
        try {
            Filter filter;
            filter = CQL.toFilter(filterStr); // filterStr形式 如  name='武汉大学' or code like 'tt123%'
            SimpleFeatureCollection result = featureSource.getFeatures(filter);
 
            FeatureIterator<SimpleFeature> itertor = result.features();
            while (itertor.hasNext()) {
                SimpleFeature feature = itertor.next();
                featureList.add(feature);
            }
            itertor.close();
            return featureList;
        } catch (CQLException e) {
            // TODO Auto-generated catch block
            e.printStackTrace();
        } catch (IOException e) {
            // TODO Auto-generated catch block
            e.printStackTrace();
        }
        return null;
    }

2.2 要素操作对上面4.1中的 SimpleFeature操作

//获取feature的geometry
Geometry geo=(Geometry) feature.getDefaultGeometry();
//获取geometry中的坐标 这里用string的方式保存
int geoUnm = geo.getNumGeometries();  // 一个geometry可能含有n个geometry
for (int i = 0; i < geoUnm; i++) {
    Geometry singleGeo = geo.getGeometryN(i); //获取其中每一个geometry
    int pointCount = singleGeo.getNumPoints();
    Coordinate[] coords = singleGeo.getCoordinates();
    for (int j = 0; j < pointCount; j++) {
        if (j == pointCount - 1)
            sBuilder.append(coords[j].x + "," + coords[j].y);
        else {
            sBuilder.append(coords[j].x + "," + coords[j].y
                                    + ";");
        }
    }
    if (i != geoUnm - 1) {
        sBuilder.append("|");
    }
} 
//获取feature中的属性
feature.getAttribute(arg0);

2.3 拓扑查询

public static Filter getGeoFilter(FilterFactory2 ff,                //构建拓扑查询的filter
            String geometryAttributeName, Geometry refGeo,
            SpatialReltionType.TopoRelTypeEnum relType) {   //这个SpatialReltionType是我自己定义的。。。
 
        switch (relType) {
        case intersect:
            return ff.intersects(ff.property(geometryAttributeName), ff
                    .literal(refGeo));
        case contains:
            return ff.contains(ff.property(geometryAttributeName), ff
                    .literal(refGeo));
        case within:
            return ff.within(ff.property(geometryAttributeName), ff
                    .literal(refGeo));
        case cross:
            return ff.crosses(ff.property(geometryAttributeName), ff
                    .literal(refGeo));
        case overlaps:
            return ff.overlaps(ff.property(geometryAttributeName), ff
                    .literal(refGeo));
        case touches:
            return ff.touches(ff.property(geometryAttributeName), ff
                    .literal(refGeo));
        case equals:
            return ff.equals(ff.property(geometryAttributeName), ff
                    .literal(refGeo));
        case disjoint:
            return ff.disjoint(ff.property(geometryAttributeName), ff
                    .literal(refGeo));
        default:
            return null;
        }
    }
// 普通的拓扑查询
public static ArrayList<Geometry> topoQueryMethod(Geometry refGeo,
            String layerName, SpatialReltionType.TopoRelTypeEnum relType) {
        ArrayList<SimpleFeature> featurelist=new ArrayList<SimpleFeature>();
        FilterFactory2 ff = CommonFactoryFinder.getFilterFactory2(null);
        SimpleFeatureSource featureSource=pgDatastore.getFeatureSource(layerName); 
    
        SimpleFeatureType schema = featureSource.getSchema();
        String geometryAttributeName = schema.getGeometryDescriptor().getLocalName();
        Filter filter1= getGeoFilter(ff,geometryAttributeName, refGeo, relType);    //上面的方法
        SimpleFeatureCollection result=null;
        try {
            result = featureSource.getFeatures(filter1);
        } catch (IOException e) {
            // TODO Auto-generated catch block
            e.printStackTrace();
        }
        if(result==null)
            return null;
        FeatureIterator<SimpleFeature> itertor = result.features();
        while (itertor.hasNext()) {
            SimpleFeature feature = itertor.next();
            featurelist.add(feature);
        }
        //这个方法是将feature转为geometry 自己定义的
        return SpatialUtil.ConverToGeoList(featurelist);  
    }
//联合属性的拓扑查询
public static ArrayList<Geometry> topoQueryMethod(Geometry refGeo,
            String queryName, String layerName,
            SpatialReltionType.TopoRelTypeEnum relType) {
        FilterFactory2 ff = CommonFactoryFinder.getFilterFactory2(null);
        ArrayList<SimpleFeature> featurelist=new ArrayList<SimpleFeature>();
        SimpleFeatureSource featureSource=pgDatastore.getFeatureSource(layerName); 
    
        SimpleFeatureType schema = featureSource.getSchema();
        String geometryAttributeName = schema.getGeometryDescriptor().getLocalName();
        Filter filter1= SpatialUtil.getGeoFilter(ff,geometryAttributeName, refGeo, relType);    
        Filter filter2=null;
        try {
            filter2=CQL.toFilter("StandName = '"+queryName+"'");
        } catch (CQLException e1) {
            // TODO Auto-generated catch block
            e1.printStackTrace();
        }
        List<Filter> match = new ArrayList<Filter>();
        match.add(filter1);
        match.add(filter2);
        Filter filter = ff.and(match);
 
        SimpleFeatureCollection result=null;
        try {
            result = featureSource.getFeatures(filter);
        } catch (IOException e) {
            // TODO Auto-generated catch block
            e.printStackTrace();
        }
        if(result==null)
            return null;
        FeatureIterator<SimpleFeature> itertor = result.features();
        while (itertor.hasNext()) {
            SimpleFeature feature = itertor.next();
            featurelist.add(feature);
        }
        return SpatialUtil.ConverToGeoList(featurelist);
        
    }

3,编辑图层

3.1 添加要素

  //添加一个feature到图层中 在添加前要确定构造featureType
    public static SimpleFeatureType createFeatureType(String typeName,Class type) {
        SimpleFeatureTypeBuilder builder = new SimpleFeatureTypeBuilder();
        builder.setName(typeName);
        builder.setCRS(DefaultGeographicCRS.WGS84); // <- Coordinate reference
                                                    // system
        builder.add("the_geom", type);  //这个为地理属性字段 postgis中为 the——geom
 
        builder.add("StandName", String.class); // 这是其他属性字段 自己定义的....                            
        // build the type
        final SimpleFeatureType TYPE = builder.buildFeatureType();
 
        return TYPE;
    }
    
    //添加到图层的图层名，添加的要素空间属性和要素的某属性名
    public static boolean addFeature(String layerName,Geometry geo,String featureName){ 
        String type=geo.getGeometryType();
        Class TypeClass=null;
        if(type.toLowerCase().equals("point")){
            TypeClass=Point.class;
        }else if(type.toLowerCase().equals("polygon")){
            TypeClass=Polygon.class;
        }else if(type.toLowerCase().equals("polyline")){
            TypeClass=Polyline.class;
        }else if(type.toLowerCase().equals("multipolygon")){
            TypeClass=MultiPolygon.class;
        }
        SimpleFeatureType featureType=createFeatureType(layerName,TypeClass);
         SimpleFeatureCollection collection = FeatureCollections.newCollection();
 
        SimpleFeatureBuilder featureBuilder = new SimpleFeatureBuilder(featureType);
         /* Longitude (= x coord) first ! */
        GeometryFactory geometryFactory = JTSFactoryFinder.getGeometryFactory(null);
     
        featureBuilder.add(geo);
        featureBuilder.add(featureName);
       
        SimpleFeature feature = featureBuilder.buildFeature(null);
        collection.add(feature);
 
        FeatureSource featureSource=pgDatastore.getFeatureSource(layerName); 
        if (featureSource instanceof SimpleFeatureStore) {
            SimpleFeatureStore featureStore = (SimpleFeatureStore) featureSource;
            Transaction transaction = new DefaultTransaction("create");
            featureStore.setTransaction(transaction);
            try {
                featureStore.addFeatures(collection);
                transaction.commit();
                return true;
 
            } catch (Exception problem) {
                problem.printStackTrace();
                try {
                    transaction.rollback();
                } catch (IOException e) {
                    // TODO Auto-generated catch block
                    e.printStackTrace();
                }
 
            } finally {
                try {
                    transaction.close();
                } catch (IOException e) {
                    // TODO Auto-generated catch block
                    e.printStackTrace();
                }
            }
         
        } else {
            System.out.println(layerName + " does not support read/write access");      
        }
        return false;
    }

3.2 修改要素

   // 修改feacode为XX的要素的名字为featureName 地理方位为geo  (feacode StandName为你的属性字段自定义)
     public static boolean modifyFeature(String layerName,Geometry geo,String featureName,String FeaCode){
          FeatureSource featureSource=pgDatastore.getFeatureSource(layerName); 
            if (featureSource instanceof SimpleFeatureStore) {
                SimpleFeatureStore featureStore = (SimpleFeatureStore) featureSource;
                Transaction transaction = new DefaultTransaction("create");
                featureStore.setTransaction(transaction);
                try {
                    String filterStr="FeaCode= '"+FeaCode+"'";
                    String[] names=new String[2];
                    names[0]="StandName";
                    names[1]="the_geom";
                    Object[] values=new Object[2];
                    values[0]=featureName;
                    values[1]=geo;
                    featureStore.modifyFeatures(names, values, CQL.toFilter(filterStr));
                
                   
                    transaction.commit();
                    return true;
                } catch (Exception problem) {
                    problem.printStackTrace();
                    try {
                        transaction.rollback();
                    } catch (IOException e) {
                        // TODO Auto-generated catch block
                        e.printStackTrace();
                    }
 
                } finally {
                    try {
                        transaction.close();
                    } catch (IOException e) {
                        // TODO Auto-generated catch block
                        e.printStackTrace();
                    }
                }
             
            } else {
                System.out.println(layerName + " does not support read/write access");      
            }
            return false;
    }

4 、Geometry 与 JTS

geotools 构建 geometry方法：这里转载一个别人写的比较好的

4.1构建点

public Point createPoint(){  
        Coordinate coord = new Coordinate(109.013388, 32.715519);  
        Point point = geometryFactory.createPoint( coord );  
        return point;  
    } 
public Point createPointByWKT() throws ParseException{  
        WKTReader reader = new WKTReader( geometryFactory );  
        Point point = (Point) reader.read("POINT (109.013388 32.715519)");  
        return point;  
    } 
public MultiPoint createMulPointByWKT()throws ParseException{  
        WKTReader reader = new WKTReader( geometryFactory );  
        MultiPoint mpoint = (MultiPoint) reader.read("MULTIPOINT(109.013388 32.715519,119.32488 31.435678)");  
        return mpoint;  
    }

4.2 构建线

public LineString createLine(){  
        Coordinate[] coords  = new Coordinate[] {new Coordinate(2, 2), new Coordinate(2, 2)};  
        LineString line = geometryFactory.createLineString(coords);  
        return line;  
    }  
 public LineString createLineByWKT() throws ParseException{  
        WKTReader reader = new WKTReader( geometryFactory );  
        LineString line = (LineString) reader.read("LINESTRING(0 0, 2 0)");  
        return line;  
    }  
public MultiLineString createMLine(){  
        Coordinate[] coords1  = new Coordinate[] {new Coordinate(2, 2), new Coordinate(2, 2)};  
        LineString line1 = geometryFactory.createLineString(coords1);  
        Coordinate[] coords2  = new Coordinate[] {new Coordinate(2, 2), new Coordinate(2, 2)};  
        LineString line2 = geometryFactory.createLineString(coords2);  
        LineString[] lineStrings = new LineString[2];  
        lineStrings[0]= line1;  
        lineStrings[1] = line2;  
        MultiLineString ms = geometryFactory.createMultiLineString(lineStrings);  
        return ms;  
    }  
public MultiLineString createMLineByWKT()throws ParseException{  
        WKTReader reader = new WKTReader( geometryFactory );  
        MultiLineString line = (MultiLineString) reader.read("MULTILINESTRING((0 0, 2 0),(1 1,2 2))");  
        return line;  
    }

4.3　构建多边形

public Polygon createPolygonByWKT() throws ParseException{  
        WKTReader reader = new WKTReader( geometryFactory );  
        Polygon polygon = (Polygon) reader.read("POLYGON((20 10, 30 0, 40 10, 30 20, 20 10))");  
        return polygon;  
    }  
public MultiPolygon createMulPolygonByWKT() throws ParseException{  
        WKTReader reader = new WKTReader( geometryFactory );  
        MultiPolygon mpolygon = (MultiPolygon) reader.read("MULTIPOLYGON(((40 10, 30 0, 40 10, 30 20, 40 10),(30 10, 30 0, 40 10, 30 20, 30 10)))");  
        return mpolygon;  
    }

4.4 构建geo集合

public GeometryCollection createGeoCollect() throws ParseException{  
        LineString line = createLine();  
        Polygon poly =  createPolygonByWKT();  
        Geometry g1 = geometryFactory.createGeometry(line);  
        Geometry g2 = geometryFactory.createGeometry(poly);  
        Geometry[] garray = new Geometry[]{g1,g2};  
        GeometryCollection gc = geometryFactory.createGeometryCollection(garray);  
        return gc;  
    }

4.5 构建圆

public Polygon createCircle(double x, double y, final double RADIUS){  
        final int SIDES = 32;//圆上面的点个数  
        Coordinate coords[] = new Coordinate[SIDES+1];  
        for( int i = 0; i < SIDES; i++){  
            double angle = ((double) i / (double) SIDES) * Math.PI * 2.0;  
            double dx = Math.cos( angle ) * RADIUS;  
            double dy = Math.sin( angle ) * RADIUS;  
            coords[i] = new Coordinate( (double) x + dx, (double) y + dy );  
        }  
        coords[SIDES] = coords[0];  
        LinearRing ring = geometryFactory.createLinearRing( coords );  
        Polygon polygon = geometryFactory.createPolygon( ring, null );  
        return polygon;  
    }

PostGIS删除表：

 SELECT DropGeometryTable ('my_schema','my_spatial_table');

如： SELECT DropGeometryTable ('public','river');

kone666

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

专栏目录