geotools 使用 部分代码总结
2012-09-05 22:39
423 查看
geotools 使用 部分代码总结
分类:gis2012-04-08 20:19160人阅读评论(0)收藏举报
前段时间的一个项目 本来用ae完成了种种的 查询,空间分析等等功能的代码,但是不幸的是 这是一个web项目,无奈 ae各种错误,显然ae放在server端是不好使的 无奈 一咬牙一跺脚 全部换 换成geotools 看文档 看api 从零 开始算是把 原来AE实现的东西 统统改了过来 用起来 反而觉得既稳定 效率还不错哈!
以下是部分功能总结:
1、连接数据库 这里使用的postgis 链接代码如下:
[java]
view plaincopyprint?
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");
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、图层的操作
[java]
view plaincopyprint?
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.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; }
[java]
view plaincopyprint?
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.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);
[java]
view plaincopyprint?
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;
}
}
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; } }
[java]
view plaincopyprint?
// 普通的拓扑查询
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 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); }
[java]
view plaincopyprint?
//联合属性的拓扑查询
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);
}
//联合属性的拓扑查询 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,编辑图层
[java]
view plaincopyprint?
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.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; }
[java]
view plaincopyprint?
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;
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方法:这里转载一个别人写的比较好的
[java]
view plaincopyprint?
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.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; }
[java]
view plaincopyprint?
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.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; }
[java]
view plaincopyprint?
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.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; }
[java]
view plaincopyprint?
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.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; }
[java]
view plaincopyprint?
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;
}
相关文章推荐
- geotools 使用 部分代码总结
- GIS的学习(二十六)geotools 使用 部分代码总结
- geotools 使用 部分代码总结
- geotools 使用 部分代码总结
- Delphi listview使用部分总结代码
- 第二章后续部分代码,总结
- i2c-tools,i2cdump和i2cset等的使用总结
- 使用IdHttp获取网页代码出现的问题,只能得到部分代码,为什么,请大家帮忙
- asp.net部分控件使用和开发技巧总结
- android使用ndk-stack调试JNI部分的C/C++代码
- Android中使用代码截图的各种方法总结
- IOS常用代码总结 - 第三方库部分
- 使用cairo绘制图形的部分函数总结
- jQuery使用代码总结
- Java代码审查工具findbugs的使用总结
- 写高质量OC代码52建议总结:36.不要使用retainCount
- afinal 的部分使用总结
- MQTT(使用mosquitto做broker)做Android推送部分总结
- 【dedecms】DEDE代码使用总结
- Android中资源文件(非代码部分)的使用概览