基于Solr的空间搜索学习笔记
在Solr中基于空间地址查询主要围绕2个概念实现：
（1） Cartesian Tiers 笛卡尔层
Cartesian Tiers是通过将一个平面地图的根据设定的层次数，将每层的分解成若干个网格，如下图所示:


 
    每层以2的评方递增，所以第一层为4个网格，第二层为16 个，所以整个地图的经纬度将在每层的网格中体现：

笛卡尔层在Lucene中对空间地理位置查询最大的用处在查找周边地址的时候有效的减少查询量，即将查询量可以控制在分层后最小的网格中的若干docId。
那么如何构建这样的索引结构呢，其实很简单，只需要对应笛卡尔层的层数来构建域即可。

也即是tiers0->field_0，tiers1->field_1,tiers2-field_2,……，tiers19->field_19。（一般20层即可）。每个对应笛卡尔层次的域将根据当前这条记录的经纬度通过笛卡尔算法计算出归属于当前层的网格，然后将gridId（网格唯一标示）以term的方式存入索引。这样每条记录关于笛卡尔0-19的域将都会有一个gridId对应起来。但是查询的时候一般是需要查周边的地址，那么可能周边的范围超过一个网格的范围，那么实际操作过程是根据经纬度和一个距离确定出需要涉及查询的从19-0（从高往低查，留给读者思考）若干层对应的若干网格的数据（关于代码实现在后面的文章内容阐述）。那么一个经纬度周边地址的查询只需要如下图圆圈内的数据：

所以通过这样的数据过滤，将极大的减少计算量。
（2） GeoHash算法
在Lucene索引中将经纬度的二维坐标通过geohash，变成一个一维的字符串base32的坐标，例如，经纬度对应一个base32的坐标为DRT2Y，那这个base32的字符串什么意思呢，
其实编码中每个字符都是代表一个区域，并且前面的字符是后面字符的父区域，即R是D区域内的子区域，T又为D区域的子区域，大家可以从如下图片获得base32的层级关系（以下图片均来自互联网）：

进入D区域，则看到又分为若干区域，而R为其子区域：

继续进入R区域，可以继续看到有子区域T区域：

而2Y也是基于以上的关系类推，所以一个base32的编码是标示一个区域，而编码过程中会根据经纬度的精度来确定这个区域大小。从上面的解释大家肯定会想到编码的前缀是表示更大的区域。例如wx4g0ec1，它的前缀wx4g0e表示包含编码wx4g0ec1在内的更大区域。所以根据这个特点，利用模糊查询是可以达到一种附近地点的查询。
geohash算法实现其实非常简单，网上有很多例子，在这里借用下这些例子再加上比较详细的说明。基本算法流程是基于多轮的收敛，以达到满足精度要求为止。具体流程以(39.92324 纬度, 116.3906 经度)为例，首先将纬度的范围(-90, 90)平分成两个区间(-90, 0)、(0, 90)，如果目标纬度位在（-90,0），则编码为0，在（0,90）则编码为1。由于上面的例子中维度39.92324是属于(0, 90)，所以第一轮获得的编码位取1。接下来再将(0, 90)分成 (0, 45), (45, 90)两个区间，而39.92324位于(0, 45)，所以编码为0。以此类推，直到精度符合要求为止，如下图所示：

所以通过16轮的计算后得到经度39.92324的编码为：1011 1000 1100 0111 1001
经度也用同样的算法，对(-180, 180)多轮的依次细分计算：

得到经度116.3906的编码为1101 0010 1100 0100 0100
经纬度的编码都计算完毕后，接下来就需要合并经纬度的编码，规则是以经度开始，依次每次取一位合并成5位的新编码，如上图红色字标示顺序所示:

完成合并编码后就需要将该编码和base32编码表对应起来，做法是每5位为一个十进制数，以11100为例，它的十进制数是28，所以对应的base32编码表示W，如下图所示：

其他的五位编码依次从表中找到对应位置后，(39.92324 纬度, 116.3906 经度)的base32编码为：wx4g0ec1
解码算法与编码算法相反，先进行base32解码，然后分离出经纬度，最后根据二进制编码对经纬度范围进行细分即可，这里不再赘述。不过由于geohash表示的是区间，编码越长越精确，但不可能解码出完全一致的地址

关于Solr+Lucene使用Cartesian Tiers 笛卡尔层和GeoHash的构建索引和查询的细节介绍将在新的Blog中阐述。

在Solr中其实支持很多默认距离函数，但是基于坐标构建索引和查询的主要会基于2种方案：

（1）GeoHash

（2）Cartesian Tiers+GeoHash

而这块的源码实现都在lucene-spatial.jar中可以找到。接下来我将根据这2种方案展开关于构建索引和查询细节进行阐述，都是代码分析，感兴趣的看官可以继续往下看。

GeoHash

 构建索引阶段

定义geohash域，在schema.xml中定义：

<fieldtype name="geohash" class="solr.GeoHashField"/>

接下来再构建索引的时候使用到lucene-spatial.jar的GeoHashUtils类：

String geoHash = GeoHashUtils.encode(latitude, longitude);//通过geoHash算法将经纬度变成base32的编码

document.addField("geohash", geoHash); //将经纬度对应的bash32编码存入索引。

查询阶段

在solrconfig.xml中配置好QP，该QP将对用户的请求Query进行QParser，查询语法规范是

{!spatial  sfield=geofield pt= latitude, longitude d=xx, sphere_radius=xx }

sfield:geohash对应的域名

pt:经纬度字符串

d=球面距离

sphere_radius：圆周半径

接下来看看QP是如何解析上述查询语句，然后生成基于GeoHash的Query的，见如下代码，代码来源SpatialFilterQParser的parse()方法：

Java代码  
  1. //GeohashType一定是继承SpatialQueryable的
  2. if (type instanceof SpatialQueryable) {
  3. double radius = localParams.getDouble(SpatialParams.SPHERE_RADIUS, DistanceUtils.EARTH_MEAN_RADIUS_KM); //圆周半径
  4. //pointStr=经纬度串，dist=距离，DistanceUnits.KILOMETERS 距离单位
  5. SpatialOptions opts = new SpatialOptions(pointStr, dist, sf, measStr, radius, DistanceUnits.KILOMETERS);
  6. opts.bbox = bbox;
  7. //通过GeoHashField 创建查询Query
  8. result = ((SpatialQueryable)type).createSpatialQuery(this, opts);
  9. }

其中最核心的方法便是GeoHashField的createSpatialQuery(),该方法负责生成基于geoHash的查询Query，展开看该方法：

Java代码  
  1. public Query createSpatialQuery(QParser parser, SpatialOptions options) {
  2. double [] point = new double[0];
  3. try {
  4. //解析经纬度
  5. point = DistanceUtils.parsePointDouble(null, options.pointStr, 2);
  6. } catch (InvalidGeoException e) {
  7. throw new SolrException(SolrException.ErrorCode.BAD_REQUEST, e);
  8. }
  9. //将经纬度编码成bash32，对如何编码请看本文geohash算法解析篇幅
  10. String geohash = GeoHashUtils.encode(point[0], point[1]);
  11. //TODO: optimize this
  12. return new SolrConstantScoreQuery(new ValueSourceRangeFilter(new GeohashHaversineFunction(getValueSource(options.field, parser),
  13. new LiteralValueSource(geohash), options.radius), "0", String.valueOf(options.distance), true, true));
  14. }

从源码中可以看到代码作者有标示TODO：optimize this，笔者从源码中看到这块的实现，也觉得确实有疑惑，整个大体实现流程是基于Lucene的Filter的方式来过滤命中docId,但是其过滤的范围让笔者看起来觉得性能会出现问题，可能也是源码中有TODO：optimize this的缘故吧。

接下来继续讲下核心处理流程，Lucene的查询规则是Query->Weight->Scorer,而主要负责查询遍历结果集合的就是Scorer，该例子也不例外，同样是：

SolrConstantScoreQueryà ConstantWeightà ConstantScorer，通过Query生成Weight，Weight生成Scorer，熟悉Lucene的读者应该很清楚了，这里不再累述，

其中ConstantScorer的通过docIdSetIterator遍历获取满足条件的docId

而docIdSetIterator便是前面源码中的ValueSourceRangeFilter，该Filter将会过滤掉不在一个指定球面距离范围内的数据，而ValueSourceRangeFilter并不是实际工作的类，它又将过滤交给了GeohashHaversineFunction，见ValueSourceRangeFilter如下代码：

Java代码  
  1. public DocIdSet getDocIdSet(final Map context, final IndexReader reader) throws IOException {
  2. return new DocIdSet() {
  3. lowerVal=0,upperVal=distance,includeLower=true,includeupper=true
  4. @Override
  5. public DocIdSetIterator iterator() throws IOException {
  6. valueSource= GeohashHaversineFunction,也是实际进行DocList过滤的类
  7. return valueSource.getValues(context, reader).getRangeScorer(reader, lowerVal, upperVal, includeLower, includeUpper);
  8. }
  9. };
  10. }

那么继续看GeohashHaversineFunction，首先看其 getRangeScorer（）方法，最核心的部分为：

Java代码  
  1. if (includeLower && includeUpper) {
  2. return new ValueSourceScorer(reader, this) {
  3. @Override
  4. public boolean matchesValue(int doc) {
  5. //计算docId对应的经纬度和查询传入的经纬度的距离
  6. float docVal = floatVal(doc);
  7. //如果返回的docVal(目标坐标和查询坐标的球面距离)在给定的distance之内则返回true
  8. //也就是说目标地址为待查询的周边范围内
  9. return docVal >= l && docVal <= u;
  10. }
  11. };
  12. }

所以再看看计算球面距离的GeohashHaversineFunction.floatVal()方法，可以从该方法最终调用的是distance()方法,如下所示：

Java代码  
  1. protected double distance(int doc, DocValues gh1DV, DocValues gh2DV) {
  2. double result = 0;
  3. String h1 = gh1DV.strVal(doc); //docId对应的经纬度的base32编码
  4. String h2 = gh2DV.strVal(doc); //查询的经纬度的base32编码
  5. if (h1 != null && h2 != null && h1.equals(h2) == false){
  6. //TODO: If one of the hashes is a literal value source, seems like we could cache it
  7. //and avoid decoding every time
  8. double[] h1Pair = GeoHashUtils.decode(h1); //base32解码
  9. double[] h2Pair = GeoHashUtils.decode(h2);
  10. //计算2个经度纬度之间的球面距离
  11. result = DistanceUtils.haversine(Math.toRadians(h1Pair[0]), Math.toRadians(h1Pair[1]),
  12. Math.toRadians(h2Pair[0]), Math.toRadians(h2Pair[1]), radius);
  13. } else if (h1 == null || h2 == null){
  14. result = Double.MAX_VALUE;
  15. }
  16. //返回2个经纬度之间球面距离
  17. return result;
  18. }

所以整个查询流程是将索引中的所有docId从第一个docId 0开始,对应的经度纬度和查询经纬度的球面距离是否在查询给定的distance之内，满足着将该docId返回，不满足则过滤。

大家可能看到是所有docId,这也是笔者觉得该过滤范围实现不靠谱的地方，也许是作者说需要进一步优化的地方。大家如果对怎么是所有docId进行过滤有疑惑，可以查看ValueSourceScorer的nextDoc() advance()方法，相信看过之后就明白了。到此Solr基于GeoHash的查询实现介绍完毕了。

基于Solr的空间搜索学习笔记相关推荐

  1. 基于MVS的三维重建算法学习笔记（五）— 立体匹配经典算法PatchMatch论文翻译及要点解读

    基于MVS的三维重建算法学习笔记(五)- 立体匹配经典算法PatchMatch论文翻译及要点解读 声明 问题提出 问题建模 通过PatchMatch获取平面参数--Inference via Patc ...

  2. 基于MVS的三维重建算法学习笔记（一）— MVS三维重建概述与OpenMVS开源框架配置

    基于MVS的三维重建算法学习笔记(一)- MVS三维重建概述与OpenMVS开源框架配置 声明 1. MVS(Multi-view stereo)概述 稀疏重建与稠密重建的区别 稀疏重建--SFM(S ...

  3. 基于MVS的三维重建算法学习笔记（四）— 立体匹配经典算法Semi-Global Matching（SGM）论文翻译及要点解读

    基于MVS的三维重建算法学习笔记(四)- 立体匹配经典算法Semi-Global Matching(SGM)论文翻译及要点解读 声明 SGM概述 Cost Calculation(像素代价计算)--M ...

  4. ELK搜索学习笔记--Day1

    ELK搜索学习笔记–Day1 1． 课程简介 1.1 课程内容 ELK是包含但不限于Elasticsearch(简称es).Logstash.Kibana 三个开源软件的组成的一个整体.这三个软件合成 ...

  5. python小车行驶路线图_基于opencv-Python小车循线学习笔记

    基于opencv-Python小车循线学习笔记 加入摄像头模块,让小车实现自动循迹行驶 思路为:摄像头读取图像,进行二值化,将白色的赛道凸显出来 选择下方的一行像素,黑色为0,白色为255 找到白色值 ...

  6. 基于MVS的三维重建算法学习笔记（二）— 立体视觉的几何基础总结

    基于MVS的三维重建算法学习笔记(二)- 立体视觉的几何基础总结 声明 概述 1. 常见三维数据类型 2. 三维形状的几种表达形式 3. 三维空间刚体运动 4. 李群和李代数 5. 相机标定 6. 非 ...

  7. 基于python的数字图像处理--学习笔记（三）

    基于python的数字图像处理--学习笔记(三) 前言 一.灰度拉伸 二.幂律(伽马)变换 三.对数变换 前言 进入冈萨雷斯的第三章内容,并用python实现功能.我更改了代码源,之前找到太烂了,代码 ...

  8. 空间计量经济学学习笔记（一）

    空间计量经济学学习笔记(一) Edit by Linhao Cui 参考 James LeSage(2014) , Paul Elhost (2018), Luc Anselin (1988), Qi ...

  9. 基于tiva的匿名飞控学习笔记（1）

    基于tiva的匿名飞控学习笔记(1) 开关状态任务 遥控器数据处理任务 数传数据交换 延时存储任务 开关状态任务 匿名飞控的开关状态任务为函数Swtich_State_Task(u8 dT_ms),定 ...

最新文章

  1. 汽车车灯灯具系统（上）
  2. 一个“蝇量级” C 语言协程库
  3. Fragment详解之五——Fragment间参数传递
  4. debian部署mysql和tomcat_【超级详细】在Debian/kali/linux服务器下搭建Tomcat7+mysql+jdk8环境...
  5. 01. elastcsearch-monitor (es 监控)
  6. nc 模拟服务器_NC集群服务器使用详解
  7. bbb mmc_blk_probe 分析
  8. 基于Hadoop 2.6.0运行数字排序的计算
  9. 《应用软件安全编程指南》国标发布 奇安信代码卫士已全面支持
  10. PHP实现常用的数据结构之链表
  11. Zircon DDK入门指南
  12. 2014年上半年计算机试题,2014年上半年程序员考试上午试题及答案
  13. SQL service 数据库 某工厂的物料管理系统数据库设计与实现
  14. 分布式 | 如何与 DBLE 进行“秘密通话”
  15. MySQL -update语句流程总结
  16. Android 天气APP（九）细节优化、必应每日一图
  17. 2021年计算机专业工作规划,2021最新计算机课教学计划
  18. HTML的滚动文字标签marquee
  19. Java国际化ResourceBundle详解
  20. MATLAB （信号与系统）实验一：信号的MATLAB表示及信号运算

热门文章

  1. html5调用手机摄像头，实现拍照上传功能
  2. DPDK uio驱动实现（二十）
  3. leetcode算法题--两数之和
  4. Spring Boot系列——7步集成RabbitMQ
  5. iOS 录音转MP3格式
  6. FastDFS蛋疼的集群和负载均衡(十五)之lvs四层+Nginx七层负载均衡
  7. 迈向成功的关键在于执行（摘自李开复博士的《做最好的自己》）
  8. Maven多模块项目使用MyBatis Generator
  9. Oracle procedure调用实例
  10. Debian读写ntfs分区格式！