pg数据库初探-河北省-县级行政区划-CFANZ编程社区

从这两张图可以看出河北主要的县级行政区划主要集中在太行山以西的华北平原上。
综合考虑大概有这几方面的原因：

平原地形是基础——水土条件好，农业基础好；
山前大道是关键——河北省太行山前这片区域，就成了早期华北平原上最安全的地带；
历史文化太厚重——华北平原是中华文化发源地。

现在抛开历史和地理因素的影响，单纯的从GIS角度考虑。这么多的县级行政区划点做聚类的展示，使用PostGIS数据库中的Voronoi多边形分析。在多个点的影响下，如何对行政区进行等面积的划分。

（感谢遥想公瑾大神的搬运）

这一番操作就可以将河北省的面进行Voronoi多边形的切割。

CREATE TABLE hebei_points AS SELECT (ST_Dump(ST_GeneratePoints(geom, 1000))).geom AS geom FROM hebei;

CREATE TABLE hebei_pts_clustered AS SELECT geom, ST_ClusterKMeans(geom, 10) over () AS cluster FROM hebei_points;

CREATE TABLE hebei_centers AS SELECT cluster, ST_Centroid(ST_collect(geom)) AS geom FROM hebei_pts_clustered GROUP BY cluster;

CREATE TABLE hebei_voronoi AS SELECT (ST_Dump(ST_VoronoiPolygons(ST_collect(geom)))).geom AS geom FROM hebei_centers;

CREATE TABLE hebei_divided AS SELECT ST_Intersection(a.geom, b.geom) AS geom FROM hebei a CROSS JOIN hebei_voronoi b;

在地图展示部分选用的是Leaflet的聚类展示。大量点数据展示，使用聚类展示更加简洁明了。
展示效果：

同时加载了切分后的面数据。

代码如下：

<!DOCTYPE html>
<html>
<head>
    <meta charset="UTF-8">
    <meta name="viewport" content="width=device-width, initial-scale=1.0">
    <title>河北县级行政</title>
    <link rel="stylesheet" href="https://unpkg.com/leaflet@1.0.3/dist/leaflet.css" integrity="sha512-07I2e+7D8p6he1SIM+1twR5TIrhUQn9+I6yjqD53JQjFiMf8EtC93ty0/5vJTZGF8aAocvHYNEDJajGdNx1IsQ==" crossorigin="" />
    <script src="https://unpkg.com/leaflet@1.0.3/dist/leaflet-src.js" integrity="sha512-WXoSHqw/t26DszhdMhOXOkI7qCiv5QWXhH9R7CgvgZMHz1ImlkVQ3uNsiQKu5wwbbxtPzFXd1hK4tzno2VqhpA==" crossorigin=""></script>
    <link rel="stylesheet" href="css/screen.css" />
    <link rel="stylesheet" href="css/MarkerCluster.css" />
    <link rel="stylesheet" href="css/MarkerCluster.Default.css" />
    <script src="js-src/leaflet.markercluster-src.js"></script>
    <script src="JSON/pointJson.js"></script>
    <script src="JSON/hebei_divided.js"></script>
    <style>
        #map,html,body{
            height: 100%;
            width: 100%;
            padding: 0;
            margin: 0;
        }
    </style>
</head>
<body>
<div id="map"></div>
</body>
<script type="text/javascript">
    var tiles = L.tileLayer('https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png', {
        maxZoom: 18,
        //attribution: '&copy; <a href="https://www.openstreetmap.org/copyright">OpenStreetMap</a> contributors'
    });

    var heBei_xian = L.tileLayer.wms("http://localhost:8088/geoserver/hebei/wms",{
        layers: 'hebei:hebei_xian',
        format: 'image/png',
        transparent: true,
        zoom: 15
    });

    var heBei_xian = L.tileLayer.wms("http://localhost:8088/geoserver/hebei/wms",{
        layers: 'hebei:hebei_xian',
        format: 'image/png',
        transparent: true
    });

    var heBei_V = L.tileLayer.wms("http://localhost:8088/geoserver/hebei_N/wms",{
        layers: 'hebei_N:hebei_divided',
        format: 'image/png',
        transparent: true
    });

    var heBei = L.tileLayer.wms("http://localhost:8088/geoserver/hebei/wms",{
        layers: 'hebei:hebei_divided',
        format: 'image/png',
        transparent: true
    });

    var map = L.map('map');
    map.addLayer(tiles);
    map.addLayer(heBei_V);
    map.addLayer(heBei_xian);
    map.addLayer(heBei);

    var markers = L.markerClusterGroup();

    var geoJsonLayer = L.geoJson(geoJsonSample, {
        onEachFeature: function (feature, layer) {
            layer.bindPopup(feature.properties.NAME);
        }
    });
    markers.addLayer(geoJsonLayer);

    //设置加载GeoJSON面的格式
    var myStyle = {
        "color": "#b3b4ff",
        "weight": 3,
        "opacity": 0.8
    };

    L.geoJSON(divided, {
        style: myStyle
    }).addTo(map);

    map.addLayer(markers);
    map.fitBounds(markers.getBounds());

    var baseLayers = {
        "河北Vorino155切分": heBei_V,
        "河北Vorino10切分": heBei,
        "河北县级行政区划": heBei_xian
    };


    L.control.layers(baseLayers).addTo(map);
    map.setView(new L.LatLng(38.827739,115.448818),6);
</script>
</html>