PostGIS代码操作简介
1. 代码操作POSTGIS的可选方案
- jdbc
- postgis-java
- geotools
- gdal
2. JDBC
public void testJdbc() {
String sql = "select st_area(st_geomfromtext('MULTIPOLYGON (((39364656.2504190132021904 2701523.9713633288629353, 39364650.82893280684947968 2701491.44244607863947749, 39364683.77488745748996735 2701504.16208679834380746, 39364683.77488745748996735 2701504.16208679834380746, 39364656.2504190132021904 2701523.9713633288629353)))'))";
DataSource ds = new SimpleDataSource("jdbc:postgresql://localhost:5432/postgres",
"postgres", "123456");
Connection connection = ds.getConnection();
ResultSet rs = connection.createStatement().executeQuery(sql);
while (rs.next()) {
System.out.println(rs.getString(1));
}
connection.close();
}
3. PostGIS-Java
public void testPostgisJava(){
String wkt = "MULTIPOLYGON (((39364656.2504190132021904 2701523.9713633288629353, 39364650.82893280684947968 2701491.44244607863947749, 39364683.77488745748996735 2701504.16208679834380746, 39364683.77488745748996735 2701504.16208679834380746, 39364656.2504190132021904 2701523.9713633288629353)))";
Connection connection = DriverManager.getConnection("jdbc:postgresql://localhost:5432/postgres",
"postgres", "123456");
PGgeometry geom = new PGgeometry(wkt);
PreparedStatement ps = connection.prepareStatement("select st_area(?)");
ps.setObject(1, geom);
ResultSet rs = ps.executeQuery();
while (rs.next()) {
System.out.println(rs.getString(1));
}
connection.close();
}
4. GeoTools
public void testGeotools(){
Map<String, Object> params = new HashMap<>();
params.put("dbtype", "postgis");
params.put("host", "127.0.0.1");
params.put("port", "5432");
params.put("schema", "public");
params.put("database", "postgres");
params.put("user", "postgres");
params.put("passwd", "123456");
params.put("preparedStatements", true);
params.put("encode functions", true);
String sql = "select st_area(st_geomfromtext('MULTIPOLYGON (((39364656.2504190132021904 2701523.9713633288629353, 39364650.82893280684947968 2701491.44244607863947749, 39364683.77488745748996735 2701504.16208679834380746, 39364683.77488745748996735 2701504.16208679834380746, 39364656.2504190132021904 2701523.9713633288629353)))'))";
JDBCDataStore jdbcDataStore = (JDBCDataStore) DataStoreFinder.getDataStore(params);
Statement statement = jdbcDataStore.getConnection(Transaction.AUTO_COMMIT).createStatement();
ResultSet rs = statement.executeQuery(sql);
while (rs.next()) {
System.out.println(rs.getString(1));
}
statement.close();
jdbcDataStore.dispose();
}
5. GDAL
public void testGdal(){
String path = "PG: host=127.0.0.1 port=5432 dbname=postgres user=postgres password=123456 active_schema=public";
ogr.RegisterAll();
org.gdal.ogr.Driver driver = ogr.GetDriverByName("PostgreSQL");
org.gdal.ogr.DataSource dataSource = driver.Open(path, 1);
String sql = "select st_area(st_geomfromtext('MULTIPOLYGON (((39364656.2504190132021904 2701523.9713633288629353, 39364650.82893280684947968 2701491.44244607863947749, 39364683.77488745748996735 2701504.16208679834380746, 39364683.77488745748996735 2701504.16208679834380746, 39364656.2504190132021904 2701523.9713633288629353)))'))";
Layer layer = dataSource.ExecuteSQL(sql);
System.out.println(layer.GetFeature(0).GetFieldAsDouble(0));
dataSource.delete();
}
6. Geometry与Geography区别
PostGIS提供了两种主要的空间数据类型：GEOMETRY和GEOGRAPHY。这两种数据类型都可以用来存储地理空间数据，如点、线和多边形，但它们在处理这些数据时有一些重要的区别：
-
坐标系统：
GEOMETRY数据类型使用笛卡尔坐标系统（平面坐标系统），而GEOGRAPHY数据类型使用球面坐标系统（经纬度）。这意味着GEOMETRY类型更适合处理本地或区域数据，而GEOGRAPHY类型更适合处理全球数据。 -
计算：
GEOMETRY类型在计算距离、面积和长度时使用的是笛卡尔坐标系统，这可能会导致结果不准确，特别是在处理大范围的数据时。相反，GEOGRAPHY类型在计算这些值时会考虑地球的曲率，因此结果更准确。 -
性能：
GEOMETRY类型通常比GEOGRAPHY类型更快，因为笛卡尔计算通常比球面计算更简单。然而，这种性能差异通常只在处理大量数据时才显著。 -
空间函数支持：
GEOMETRY类型支持更多的空间函数，而GEOGRAPHY类型只支持一部分空间函数。
在选择使用GEOMETRY还是GEOGRAPHY时，你应该根据你的具体需求和数据来决定。如果你需要处理全球数据，并且需要准确的距离和面积计算，那么GEOGRAPHY可能是更好的选择。如果你需要处理本地或区域数据，或者需要使用更多的空间函数，那么GEOMETRY可能是更好的选择。