postgis的地图切割方案

点线面图层的sh脚本：

点图层：

#!/usr/bin/sh<br /> export PGCLIENTENCODING=GBK<br /> LAYER=$1<br /> POINT_SQL=”SELECT user_label,b.the_geom FROM ${LAYER} b ,ST_SetSRID(ST_MakeBox2D(ST_SetSRID(ST_MakePoint($2,$3),4326),ST_SetSRID(ST_MakePoint($4,$5),4326)),4326) a_geom where st_intersects(b.the_geom,a_geom) ;”<br /> pgsql2shp -f “/home/postgres/${LAYER}.shp” postgis_test “${POINT_SQL}”<br />

线图层：

#!/usr/bin/sh<br /> export PGCLIENTENCODING=GBK<br /> LAYER=$1<br /> LINE_SQL=”SELECT b.user_label,(st_intersection(a_geom,b.the_geom)) as the_geom FROM ${LAYER} b ,ST_SetSRID(ST_MakeBox2D(ST_SetSRID(ST_MakePoint(120.0,20.0),4326),ST_SetSRID(ST_MakePoint(128.0,40.0),4326)),4326) a_geom where b.the_geom && a_geom and not ST_ISclosed(st_intersection(a_geom,b.the_geom)) ;”<br /> pgsql2shp -f “/home/postgres/${LAYER}.shp” postgis_test “${LINE_SQL}”

面图层：

#!/usr/bin/sh<br /> export PGCLIENTENCODING=GBK<br /> LAYER=$1<br /> PLOYGON_SQL=”SELECT b.user_label,(st_intersection(a_geom,b.the_geom)) as the_geom FROM ${LAYER} b ,ST_SetSRID(ST_MakeBox2D(ST_SetSRID(ST_MakePoint($2,$3),4326),ST_SetSRID(ST_MakePoint($4,$5),4326)),4326) a_geom where b.the_geom && a_geom and ST_ISclosed(st_intersection(a_geom,b.the_geom)) group by b.user_label,a_geom ,b.the_geom ;”<br /> pgsql2shp -f “/home/postgres/weih/Excursion/shp/${LAYER}.shp” postgis_test “${PLOYGON_SQL}”<br />

然后可以通过Java的telnet与Ftp的实现对这些图层数据的控制。

转载自：https://blog.csdn.net/no_cross_no_crown/article/details/6213735