postgis mvt矢量切片 django drf mapboxgl

postgis mvt矢量切片 django drf mapboxgl

目录

0.前提

1.sql代码

2.django drf后端服务代码

3.具体的应用（整体代码）

4.参考

0.前提

        [1] 静态的矢量切片可以采用 tippecanoe 生成，nginx代理，这种数据是不更新的；

        [2] 动态的矢量切片，一般采用postgis生成。基本上矢量切片80%的厂商都采用postgis，确实好用！不谈商业的。

        [3] postgis矢量切片使用到的函数：ST_AsMVT、ST_AsMVTGeom、ST_TileEnvelope、ST_Transform、ST_Intersects、ST_SRID。（点击查看每个函数的介绍，都是官网文档，很详细）

        [4] postgis api参考文档官网：Official Manual | PostGIS，有postgis3.0版本以上的，点击html，点击8. PostGIS Reference，即可查看矢量处理的函数。        

        [5] 矢量切片mvt，需要坐标为EPSG:3857，如果想用mvt这种开源的标准，就使用ST_Transform统一转换成这个坐标。然后入库的数据都设置成EPSG:4326。

1.sql代码

        [1] 获取表的字段名称：（zzz替换成传入的表名）【PS：不建议动态查询将1-2结合，不好】

        select column_name from information_schema.columns where table_name='zzz';

        [2] 动态获取矢量切片： (1.0.0)替换成传入的z,x,y参数，zzz替换成传入的表名，geom替换成geom几何对应的字段名称。       

1.  
   with mvtgeom as (
2.  
   select ST_AsMVTGeom(ST_Transform(geom, 3857), ST_TileEnvelope(1,0,0)) as geom, gid 
3.  
       from 
4.  
       zzz, 
5.  
       (select ST_SRID(geom) as srid from zzz where geom is not null limit 1) a
6.  
       where
7.  
         ST_Intersects(geom, ST_Transform(ST_TileEnvelope(1,0,0), srid))
8.  
   )
9.  
   select ST_AsMVT(mvtgeom.*, 'zzz', 4096, 'geom') as mvt from mvtgeom;

2.django drf后端服务代码

        基于APIView重写get函数，再注册到urls.py中

        前端访问  ip/`table_name`/`z`/`x`/`y`，eg：http://127.0.0.1:8080/getmap/zzz/2/1/1 这样就可以接收矢量切片mvt了。

1.  
   from django.db import connection
2.  
   from rest_framework.views import APIView
3.  
   from rest_framework.response import Response
4.  
    
5.  
   class MapView(APIView):
6.  
    
7.  
   def get(self, request, table, z, x, y):
8.  
    
9.  
   print(table, z, x, y,)
10.  
    table_name = table #'zzz'
11.  
     
12.  
    # 将数据库的数据导出为geojson（建议不要太大，太大用 wms服务吧、或者mvt矢量切片）
13.  
    # sql = f"""select json_build_object('type', 'FeatureCollection', 'name', '{table_name}', 'features', json_agg(ST_ASGeoJSON(t.*)::json)) from {table_name} AS t """
14.  
     
15.  
    # 获取表的字段名称 列表
16.  
    # sql = f"""select column_name from information_schema.columns where table_name='{table_name}'""" # where后面是string 不应被转成对象变量
17.  
    # 获取表的字段名称 字符串
18.  
    # sql = f"""select array_to_string(array(select column_name from information_schema.columns where table_name='{table_name}' and column_name != 'geom'), ',');"""
19.  
    # 动态矢量切片mvt gid字段也可替换成前端传入的字段
20.  
    geom_name = 'geom'
21.  
    sql = f"""with mvtgeom as (
22.  
    select ST_AsMVTGeom(ST_Transform({geom_name}, 3857), ST_TileEnvelope({z},{x},{y})) as geom, gid
23.  
    from
24.  
    {table_name},
25.  
    (select ST_SRID({geom_name}) as srid from {table_name} where {geom_name} is not null limit 1) a
26.  
    where
27.  
    ST_Intersects({geom_name}, ST_Transform(ST_TileEnvelope({z},{x},{y}), srid))
28.  
    )
29.  
    select ST_AsMVT(mvtgeom.*, '{table_name}', 4096, 'geom') as mvt from mvtgeom;"""
30.  
     
31.  
    print(sql)
32.  
    with connection.cursor() as cursor:
33.  
    cursor.execute(sql)
34.  
    results = cursor.fetchall()
35.  
     
36.  
    return Response(results)

3.具体的应用（整体代码）

        TODO

        [1] shp、geojson前端上传，后端对数据进行校验（校验坐标系-强制4326、等），然后再采用后端代码将shp、geojson导入postgres数据库中。

        [2] jwt、permission等

        [3] 前端如何加jwt等信息 请求数据

4.参考

        [1] 参考了这个项目的sql代码（比国内很多博客写的强太多，如果不考虑权限等，就是部署个后端、数据库，手动导入shp数据，那么直接用这个开源项目即可，前端采用mapboxgl类似的开源库加载，api访问部署后的localhost:3000即可查看）：https://github.com/tobinbradley/dirt-simple-postgis-http-api/blob/master/routes/mvt.js

这篇好文章是转载于：学新通技术网