Slide 1

Slide 1 text

空间索引在地理服务中实践 W E L C O M E

Slide 2

Slide 2 text

⾃我介绍 ID:⼀缕殇流化隐半边冰霜(花名:霜菜) 了解 iOS 开发、前端⼩⽩、后端新⼿ 略懂 JavaScript、Go、C、C++、Objective-C、Swift Github:@halfrost 微博:@halfrost

Slide 3

Slide 3 text

F O R E W O R D N E X T S L I D E

Slide 4

Slide 4 text

PRESENTATION AGENDA Geohash 01 03 Google S2 02 04 Application Space filling curve

Slide 5

Slide 5 text

#01 GEOHASH S P AT I A L I N D E X N E X T S L I D E

Slide 6

Slide 6 text

EXAMPLE B A S E - 3 2 B A S E - 3 6

Slide 7

Slide 7 text

EXAMPLE

Slide 8

Slide 8 text

4 5 % 4 5 % 4 5 % 4 5 % 4 5 % 4 5 % GEOHASH (LEVEL-6) (31.1932993, 121.43960190000007) 纬 经

Slide 9

Slide 9 text

GEOHASH (LEVEL-6) 101011000101110 L AT I T U D E 纬 度 110101100101101 L O N G I T U D E 经 度 偶数位放经度,奇数位放纬度 (第 0 位为第⼀位) 111001100111100000110011110110

Slide 10

Slide 10 text

N E X T 11100 11001 11100 00011 00111 10110 28 25 28 3 7 22 wtw37q

Slide 11

Slide 11 text

EXAMPLE

Slide 12

Slide 12 text

CODE package geohash import ( "bytes" ) const ( BASE32 = "0123456789bcdefghjkmnpqrstuvwxyz" MAX_LATITUDE float64 = 90 MIN_LATITUDE float64 = -90 MAX_LONGITUDE float64 = 180 MIN_LONGITUDE float64 = -180 ) var ( bits = []int{16, 8, 4, 2, 1} base32 = []byte(BASE32) ) type Box struct { MinLat, MaxLat float64 // 纬度 MinLng, MaxLng float64 // 经度 } func (this *Box) Width() float64 { return this.MaxLng - this.MinLng } func (this *Box) Height() float64 { return this.MaxLat - this.MinLat } // 输⼊值:纬度,经度,精度(geohash的⻓度) // 返回geohash, 以及该点所在的区域 func Encode(latitude, longitude float64, precision int) (string, *Box) { var geohash bytes.Buffer var minLat, maxLat float64 = MIN_LATITUDE, MAX_LATITUDE var minLng, maxLng float64 = MIN_LONGITUDE, MAX_LONGITUDE var mid float64 = 0 bit, ch, length, isEven := 0, 0, 0, true for length < precision { if isEven { if mid = (minLng + maxLng) / 2; mid < longitude { ch |= bits[bit] minLng = mid } else { maxLng = mid } } else { if mid = (minLat + maxLat) / 2; mid < latitude { ch |= bits[bit] minLat = mid } else { maxLat = mid } } isEven = !isEven if bit < 4 { bit++ } else { geohash.WriteByte(base32[ch]) length, bit, ch = length+1, 0, 0 } } b := &Box{ MinLat: minLat, MaxLat: maxLat, MinLng: minLng, MaxLng: maxLng, } return geohash.String(), b }

Slide 13

Slide 13 text

SCALE

Slide 14

Slide 14 text

ERROR

Slide 15

Slide 15 text

SHARDING RULE 根据经纬度计算 Shard ID 及 Cell ID 若上⼀次在同⼀个 Cell,直接更新 若上⼀次在不同 Cell,先删除上⼀ 个 Cell,再加⼊当前 Cell 的队列 中。 W R I T E ( U P D AT E L O C AT I O N ) 计算相交的 Shard 每个 Shard 做并⾏ Nearby 计算。 Geohash 只能先按照⼴度优先查 找相应的 Cell,过滤出 Cell 中符 合条件的点。 R E A D ( N E A R B Y S E A R C H )

Slide 16

Slide 16 text

WHY why? 偶 数 位 放 经 度 , 奇 数 位 放 纬 度 why? 理 论 基 础 why? 有 缺 点 么 ?

Slide 17

Slide 17 text

#02 SPACE FILLING CURVE S P AT I A L F I L L I N G C U R V E N E X T S L I D E

Slide 18

Slide 18 text

No content

Slide 19

Slide 19 text

Z-ORDER

Slide 20

Slide 20 text

Z-ORDER

Slide 21

Slide 21 text

DISADVANTAGE

Slide 22

Slide 22 text

DISADVANTAGE

Slide 23

Slide 23 text

DISADVANTAGE Z 阶曲线可以将⼆维或者多 维空间⾥的所有点都转换成 ⼀维曲线。在数学上成为分 形维。 降 维 局 部 保 序 性 突 变 性 搜 索 查 找 邻 近 点 ⽐ 较 快

Slide 24

Slide 24 text

EXAMPLE

Slide 25

Slide 25 text

PEANO CURVE ⽪亚诺曲线是⼀条连续的但处处不可导的曲线。

Slide 26

Slide 26 text

DRAGON CURVE

Slide 27

Slide 27 text

GOSPER CURVE

Slide 28

Slide 28 text

KOCH CURVE

Slide 29

Slide 29 text

MOORE CURVE

Slide 30

Slide 30 text

SIERPIŃSKI CURVE

Slide 31

Slide 31 text

HILBERT CURVE

Slide 32

Slide 32 text

HILBERT CURVE

Slide 33

Slide 33 text

HAUSDORFF FRACTALS DIMENSION

Slide 34

Slide 34 text

HILBERT CURVE

Slide 35

Slide 35 text

HILBERT CURVE

Slide 36

Slide 36 text

HILBERT CURVE

Slide 37

Slide 37 text

HILBERT CURVE

Slide 38

Slide 38 text

HILBERT CURVE

Slide 39

Slide 39 text

HILBERT CURVE

Slide 40

Slide 40 text

HILBERT CURVE

Slide 41

Slide 41 text

HILBERT CURVE

Slide 42

Slide 42 text

HILBERT CURVE

Slide 43

Slide 43 text

CHARACTERISTIC 降 维 稳 定 连 续 希尔伯特曲线是连续的,所以能保证⼀定 可以填满空间。连续性是需要数学证明 的。具体证明⽅法这⾥就不细说了,感兴 趣的可以点⽂章末尾⼀篇关于希尔伯特曲 线的论⽂,那⾥有连续性的证明。

Slide 44

Slide 44 text

CONTINUITY

Slide 45

Slide 45 text

ANY QUESTIONS? T H A N K S F O R Y O U R AT T E N T I O N ! N E X T S L I D E

Slide 46

Slide 46 text

#03 GOOGLE S2 S P AT I A L I N D E X N E X T S L I D E

Slide 47

Slide 47 text

LAT / LNG x = r * sin θ * cos φ y = r * sin θ * sin φ z = r * cos θ

Slide 48

Slide 48 text

LAT / LNG s(lat,lng) -> f(x,y,z)

Slide 49

Slide 49 text

PROJECTION

Slide 50

Slide 50 text

PROJECTION

Slide 51

Slide 51 text

FRACTAL

Slide 52

Slide 52 text

FRACTAL

Slide 53

Slide 53 text

FRACTAL 3 0 % 4 5 % 5 0 % f(x,y,z) -> g(face,u,v)

Slide 54

Slide 54 text

FIXED g(face,u,v) -> h(face,s,t)

Slide 55

Slide 55 text

PROGRAM 线性变换 u = 0.5 * ( u + 1) tan() 三⻆变换 u = 2 / pi * (atan(u) + pi / 4) = 2 * atan(u) / pi + 0.5 ⼆次变换 u >= 0,u = 0.5 * sqrt(1 + 3*u) u < 0,u = 1 - 0.5 * sqrt(1 - 3*u) 1 2 3

Slide 56

Slide 56 text

TRANSFORM h(face,s,t) -> H(face,i,j)

Slide 57

Slide 57 text

HILBERT CURVE

Slide 58

Slide 58 text

HILBERT CURVE

Slide 59

Slide 59 text

HILBERT CURVE

Slide 60

Slide 60 text

HILBERT CURVE LEVEL H(face,i,j) -> CellID

Slide 61

Slide 61 text

CELL ID

Slide 62

Slide 62 text

N E X T S L I D E S(lat,lng) -> f(x,y,z) -> g(face,u,v) -> h(face,s,t) -> H(face,i,j) -> CellID GOOGLE S2

Slide 63

Slide 63 text

CLACULATE

Slide 64

Slide 64 text

GEOHASH VS GOOGLE S2 各种向量计算,⾯ 积计算,多边形覆 盖,距离问题,球 ⾯球体上的问题 ⼏ 何 计 算 S2 还能解决多 边形覆盖的问题 多 边 形 覆 盖 S2 有30级,从 0.7cm² 到 85,000,000km² 。S2 的存 储只需要⼀个 uint64 即可 存下 突 变 性 Geohash 有12级,从5000km 到 3.7cm。中间每⼀ 级的变化⽐较⼤。有时候可能选择上⼀级会⼤很 多,选择下⼀级⼜会⼩⼀些。⽐如选择字符串⻓度 为4,它对应的 cell 宽度是39.1km,需求可能是 50km,那么选择字符串⻓度为5,对应的 cell 宽度 就变成了156km,瞬间⼜⼤了3倍了。Geohash 需 要 12 bytes 存储 L E V E L 精 细 度

Slide 65

Slide 65 text

GOOGLE S2 1.涉及到⻆度,间隔,纬度经度点,单位⽮量等的表示,以及对这些 类型的各种操作。 2.单位球体上的⼏何形状,如球冠(“圆盘”),纬度 - 经度矩形,折 线和多边形。 3.⽀持点,折线和多边形的任意集合的强⼤的构造操作(例如联合) 和布尔谓词(例如,包含)。 4.对点,折线和多边形的集合进⾏快速的内存索引。 5.针对测量距离和查找附近物体的算法。 6.⽤于捕捉和简化⼏何的稳健算法(该算法具有精度和拓扑保证)。 7.⽤于测试⼏何对象之间关系的有效且精确的数学谓词的集合。 8.⽀持空间索引,包括将区域近似为离散“S2单元”的集合。此功能可 以轻松构建⼤型分布式空间索引。

Slide 66

Slide 66 text

#04 APPLICATION F I N A L N E X T S L I D E

Slide 67

Slide 67 text

APPLICATION

Slide 68

Slide 68 text

APPLICATION

Slide 69

Slide 69 text

APPLICATION

Slide 70

Slide 70 text

APPLICATION

Slide 71

Slide 71 text

APPLICATION

Slide 72

Slide 72 text

APPLICATION

Slide 73

Slide 73 text

APPLICATION

Slide 74

Slide 74 text

APPLICATION 流量是每秒钟⼤概数万条消息,⼀天⼤概是⼏亿,并且每条消息包含⼏⼗个字段 1.⽀持时序和地理空间的切⽚ 2.⽀持⼤流量数据 3.⽀持秒级(毫秒级?)查询 4.⽀持原始数据查询 ElasticSearch + Kafka

Slide 75

Slide 75 text

ONE MORE THING https://github.com/halfrost/Halfrost-Field

Slide 76

Slide 76 text

THANKS B Y E B Y E