postgis计算矢量切片

    没写错,是使用postgis计算出来矢量切片。在这之前先准备一个数据:一个GIS数据表(本例中数据为一百万的点数据,坐标:4326),并在表中添加x,y字段,方便后面的数据筛选。sql中用到了
ST_AsMVTST_AsMVTGeom

    本文中创建矢量切片很简单,就是使用下方的一个sql,运行结果如下图。接着写一个矢量切片的http服务(参考go-vtile-example,这个例子中矢量切片压缩率更高),并且使用mapbox进行前端展示(小贴士:sql中‘points’的字符串与渲染中mapbox里的source-layer一致).代码见最下方

SELECT ST_AsMVT(tile,'points') tile  FROM(
SELECT ST_AsMVTGeom(geom,ST_MakeEnvelope(100,10,125,22, 4326),4096, 0, true)    
 AS geom FROM grid20180322 ) 
AS  tile where  tile.geom is not null
准备的数据
sql查询

前端展示

切片请求控制台截图
package main

import 
(   
    _ "github.com/lib/pq"
    "database/sql"
    "time"
    "log"
    "math"
    "errors"
    "fmt"
    "net/http"
    "regexp"
    "strconv"
    "strings"
)

func tilePathToXYZ(path string) (TileID, error) {
    xyzReg := regexp.MustCompile("(?P<z>[0-9]+)/(?P<x>[0-9]+)/(?P<y>[0-9]+)")
    matches := xyzReg.FindStringSubmatch(path)
    if len(matches) == 0 {
        return TileID{}, errors.New("Unable to parse path as tile")
    }
    x, err := strconv.ParseUint(matches[2], 10, 32)
    if err != nil {
        return TileID{}, err
    }
    y, err := strconv.ParseUint(matches[3], 10, 32)
    if err != nil {
        return TileID{}, err
    }
    z, err := strconv.ParseUint(matches[1], 10, 32)
    if err != nil {
        return TileID{}, err
    }
    return TileID{x: uint32(x), y: uint32(y), z: uint32(z)}, nil
}
type LngLat struct {
    lng float64
    lat float64
}
type TileID struct {
    x uint32
    y uint32
    z uint32
}

func tile2lon( x int,  z int)(a float64) {
    return float64(x) /math.Pow(2, float64(z)) * 360.0 - 180;
 }

 func tile2lat( y int,  z int)(a float64) {
   n := math.Pi - (2.0 * math.Pi * float64(y)) / math.Pow(2, float64(z));
   return math.Atan(math.Sinh(n))*180/math.Pi;
 }
 

func FloatToString(input_num float64) string {
    // to convert a float number to a string
    return strconv.FormatFloat(input_num, 'f', 6, 64)
}

 func loadData(xyz TileID)(a []byte){
    ymax :=FloatToString(tile2lat(int(xyz.y), int(xyz.z)));
    ymin := FloatToString(tile2lat(int(xyz.y+1), int(xyz.z)));
    xmin := FloatToString(tile2lon(int(xyz.x), int(xyz.z)));
    xmax := FloatToString(tile2lon(int(xyz.x+1), int(xyz.z)));
    fmt.Println("ymax: ", ymax)
    fmt.Println("ymin: ",ymin)
    fmt.Println("xmin : ",xmin )
    fmt.Println("xmax : ",xmax )
    connStr := "dbname=xx user=xx password=xx host=localhost port=5433  sslmode=disable"
    db, err := sql.Open("postgres", connStr)
    if err != nil {
        panic(err)
    }
    defer db.Close()

    err = db.Ping()
    if err != nil {
        panic(err)
    }

    fmt.Println("Successfully connected!")
    var tile []byte
     s := []string{xmin,ymin,xmax,ymax}
     maxmin:=strings.Join(s, ",") 
     s2 := []string{" where (x between", xmin,"and",xmax,") and ( y between",ymin,"and",ymax,")"}
    wmaxmin:=strings.Join(s2, " ") 
    sql:="SELECT ST_AsMVT(tile,'points') tile  FROM(SELECT ST_AsMVTGeom(w.geom,ST_MakeEnvelope("+maxmin+", 4326),4096, 0, true)  AS geom FROM (select geom from grid20180322"+wmaxmin+") w) AS  tile where  tile.geom is not null"
    rows1:= db.QueryRow(sql)
    err1 := rows1.Scan(&tile)
    if err1 != nil {
        log.Fatal(err1)
    }
    fmt.Println(sql)
    //defer rows1.Close()
    return tile
 }
func main(){
    //t1 := time.Now() 
    mux := http.NewServeMux()
    tileBase := "/tiles/"
    mux.HandleFunc(tileBase, func(w http.ResponseWriter, r *http.Request) {
        t2 := time.Now() 
        log.Printf("url: %s", r.URL.Path)
        tilePart := r.URL.Path[len(tileBase):]
        fmt.Println("tilePart: ", tilePart)
        xyz, err := tilePathToXYZ(tilePart)
        fmt.Println("xyz: ", xyz)
        
        if err != nil {
            http.Error(w, "Invalid tile url", 400)
            return
        }
        tile:=loadData(xyz)
        
        // All this APi to be requests from other domains.
        w.Header().Set("Content-Type", "application/x-protobuf")
        w.Header().Set("Access-Control-Allow-Origin", "*")
        w.Header().Set("Access-Control-Allow-Methods", "GET, POST, OPTIONS")
        w.Write(tile)
        elapsed2 := time.Since(t2)
    
    fmt.Println("耗时: ", elapsed2)
    })
    log.Fatal(http.ListenAndServe(":8081", mux))
}

<!DOCTYPE html>
<html>

<head>
    <meta charset='utf-8' />
    <title>Add a third party vector tile source</title>
    <meta name='viewport' content='initial-scale=1,maximum-scale=1,user-scalable=no' />
    <script src='mapbox-gl.js'></script>
    <link href='mapbox-gl.css' rel='stylesheet' />
    <style>
        body {
            margin: 0;
            padding: 0;
        }
        
        #map {
            position: absolute;
            top: 0;
            bottom: 0;
            width: 100%;
        }
    </style>
</head>

<body>

    <div id='map'></div>
    <script>
        mapboxgl.accessToken = undenfined;
        var tileset = 'mapbox.streets';
        var map = new mapboxgl.Map({
            container: 'map',
            zoom: 12,
            center: [109.898625809072612, 19.106708155556731],
            style: 'mapbox://styles/mapbox/light-v9',
            hash: false
        });

        map.on('load', function loaded() {
            map.addSource('custom-go-vector-tile-source', {
                type: 'vector',
                tiles: ['http://localhost:8081/tiles/{z}/{x}/{y}']
            });
            map.addLayer({
                id: 'background',
                type: 'background',
                paint: {
                    'background-color': 'white'
                }
            });
            map.addLayer({
                "id": "custom-go-vector-tile-layer",
                "type": "circle",
                "source": "custom-go-vector-tile-source",
                "source-layer": "points",
                paint: {
                    'circle-radius': {
                        stops: [
                            [8, 0.1],
                            [11, 0.5],
                            [15, 3],
                            [20, 20]
                        ]
                    },
                    'circle-color': '#e74c3c',
                    'circle-opacity': 1
                }
            });
        });
    </script>

</body>

</html>
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 222,252评论 6 516
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 94,886评论 3 399
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 168,814评论 0 361
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,869评论 1 299
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,888评论 6 398
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 52,475评论 1 312
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 41,010评论 3 422
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,924评论 0 277
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 46,469评论 1 319
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,552评论 3 342
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,680评论 1 353
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 36,362评论 5 351
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 42,037评论 3 335
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,519评论 0 25
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,621评论 1 274
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 49,099评论 3 378
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,691评论 2 361

推荐阅读更多精彩内容

  • 方案背景     今年三月份写了一篇postgis计算矢量切片,参考了网上资料给出了一份很粗糙的相关方案(文章写的...
    polong阅读 1,891评论 3 3
  • from PostgreSQL+PostGIS 的使用 一、PostGIS中的几何类型 PostGIS支持所有O...
    山华水清阅读 6,977评论 0 5
  • 本文转载自PostGIS中的常用函数并稍加修改 说明:这里许多函数是以ST_[X]yyy形式命名的,事实上很多函数...
    AllanHao阅读 7,009评论 0 4
  • 若想发布自己矢量切片,除了使用mapbox推荐的工具与geoserver外还能够直接使用postgis+后端程序进...
    上岸躲雨阅读 9,036评论 2 5
  • 冒着隆冬的寒气,他细心地整理着材料,想着可能遇到的问题。心里稍微有点忐忑地向办事大厅走去。 “小伙子,...
    龙二二六阅读 378评论 0 0