一种逆向渲染绘制热力图的方法技术

本发明专利技术公开一种逆向渲染绘制热力图的方法，包括以若干具有时空特性的热力分析数据点为样本，训练核密度分析函数概率模型；由图像对应坐标范围及分辨率计算各个像素点对应的地理坐标；将其输入核密度分析函数概率模型，拟合模型的核密度估计函数计算获得像素区域的影响力；将像素区域的影响力映射到范围[0，255]的灰度区间内，渲染出像素点的影响力灰度；根据像素点的影响力灰度获取色彩梯度的分割阈值；依据所得的色彩梯度的分割阈值对像素点的影响力灰度赋色生成热力图。本发明专利技术将像素点替代数据点作为分析元，极大减小了地理空间分析的计算量，可在适应亿级海量数据热力图渲染的需求的同时保证热力图渲染的质量。

【技术实现步骤摘要】
一种逆向渲染绘制热力图的方法
本专利技术涉及热力图绘制领域，特别涉及一种逆向渲染绘制热力图的方法。
技术介绍
热力图是对热力分析数据的一种直观表示方法，在空间大数据挖掘和知识发现的研究中具有良好的展示效果。一般的热力图，将数据进行投影聚类，确定数据点影响半径，然后绘制出热力灰度图，并用渐变色为热力图赋色。很多开源网络地图如：国外的Openlayer、Leaflet以及国内的百度地图、天地图和超图地图等在线地图都为开发者提供了API接口或插件，可以在这些平台上便捷地实现一般热力图的绘制。目前的热力图可视化研究和一般热力图绘制也存在一些不足，主要有以下2点：1自适应能力差，不同尺度下变形失真，在不同的地图比例下，热力图所展示的数据特征差异很大；2大数据支持粒度不足，一般只能支持10w以下的点热力图渲染；3热力图在绘制时由于将数据投影到屏幕视窗，未充分考虑数据的地理空间特性。
技术实现思路
本专利技术的主要目的是提出一种逆向渲染绘制热力图的方法，旨在克服以上问题。为实现上述目的，本专利技术提出的一种逆向渲染绘制热力图的方法，包括如下步骤：S10以若干具有时空特性的热力分析数据点为样本，训练核密度分析函数概率模型；S20由图像对应坐标范围及分辨率计算各个像素点对应的地理坐标；S30将各个像素点对应的地理坐标输入核密度分析函数概率模型，拟合模型的核密度估计函数计算获得各个像素中心点的概率密度，即该像素区域的影响力；S40将该像素区域的影响力映射到范围[0，255]的灰度区间内，渲染出像素点的影响力灰度；S50根据像素点的影响力灰度获取色彩梯度的分割阈值；S60依据所得的色彩梯度的分割阈值对像素点的影响力灰度赋色生成热力图。优选地，所述S50根据像素点的影响力灰度，采用最大熵阈值分割算法获取色彩梯度的分割阈值，最大熵阈值分割算法具体为：S501将图像分割为若干子图像；S502根据像素点的影响力灰度，计算各子图像的熵；S503遍历穷举分割子图像的各个阈值，得到使得各个子图像熵的和最大的子图像分割阈值，以此为色彩梯度的分割阈值。优选地，所述方法采用Openlayer交互式地图开源JavaScript库。优选地，所述渲染是采用HtmlCanvasl图形容器的脚本程序绘制。优选地，所述脚本程序采用JavaScript语言编写。优选地，所述像素中心点的概率密度是与该点临近的像素点的概率密度值按照地理空间坐标距离衰减权重加权而得。本专利技术以像素替代数据点为分析元，极大减小了地理空间分析的粒度，可在适应亿级海量数据热力图渲染的需求的同时保证热力图渲染的质量。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图示出的结构获得其他的附图。图1为本专利技术一实施例的方法流程图；图2为所述S50采用最大熵阈值分割算法的方法流程图；图3为一般的HeatMap热力图绘制方法和本专利的逆向绘制地理热力图方法，在(1:20000)比例尺下做对比测试图；图4为一般的HeatMap热力图绘制方法和本专利的逆向绘制地理热力图方法，在(1:2000)比例尺下做对比测试图；图5为一般的HeatMap热力图绘制方法和本专利的逆向绘制地理热力图方法，在(1:500)比例尺下做对比测试图，本专利技术目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。需要说明，若本专利技术实施例中有涉及方向性指示(诸如上、下、左、右、前、后……)，则该方向性指示仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。另外，若本专利技术实施例中有涉及“第一”、“第二”等的描述，则该“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本专利技术要求的保护范围之内。本专利技术提出的一种逆向渲染绘制热力图的方法，包括如下步骤：S10以若干具有时空特性的热力分析数据点为样本，训练核密度分析函数概率模型；S20由图像对应坐标范围及分辨率计算各个像素点对应的地理坐标；S30将各个像素点对应的地理坐标输入核密度分析函数概率模型，拟合模型的核密度估计函数计算获得各个像素中心点的概率密度，即该像素区域的影响力；S40将该像素区域的影响力映射到范围[0，255]的灰度区间内，渲染出像素点的影响力灰度；S50根据像素点的影响力灰度获取色彩梯度的分割阈值；S60依据所得的色彩梯度的分割阈值对像素点的影响力灰度赋色生成热力图。在本专利技术实施例中，本专利以视窗像素为对象，分析屏幕像素点所表示的地理区域和周围邻近数据点之间的相互联系(即像素区域的影响力)计算获得像素点的影响力灰度(并非灰度值为0-255的传统意义灰度图)，然后按照灰度图的整体分布规律划分得到渐变色分级阈值(即色彩梯度的分割阈值)，以这种自适应的方式获得热力图的色彩梯度,然后赋色得到热力图。本专利技术通过色彩梯度自适应并自定义影响力的叠加，实现了一种逆向渲染热力图的绘制方法。本专利技术以视窗单个像素为分析元，以渲染器像素映射的地理空间作为计算分析的空间粒度，极大减小了地理空间分析的计算量，解决了大数据支持粒度不足的问题，也解决了热力图影响力叠加规则依赖于渲染器机制的缺点，可在适应亿级海量数据热力图渲染的需求的同时保证热力图渲染的质量。优选地，所述S50根据像素点的影响力灰度，采用最大熵阈值分割算法获取色彩梯度的分割阈值，最大熵阈值分割算法具体为：S501将图像分割为若干子图像；S502根据像素点的影响力灰度，计算各子图像的熵；S503遍历穷举分割子图像的各个阈值，得到使得各个子图像熵的和最大的子图像分割阈值，以此为色彩梯度的分割阈值。优选地，所述方法采用Openlayer交互式地图开源JavaScript库。优选地，所述渲染是采用HtmlCanvasl图形容器的脚本程序绘制。优选地，所述脚本程序采用JavaScript语言编写。优选地，所述像素中心点的概率密度是与该点临近的像素点的概率密度值按照地理空间坐标距离衰减权重加权而得。在本专利技术实施例中，针对自适应能力差的问题，逆向渲染热力图方法使用地理距离与绘制像素结合计算得到分析点缓冲区半径系数与影响力参数，每个像素点的影响力不仅仅是该点的核密度概率值，而是与该点临近的像素点的概率密度值按照距离衰减权重加权而得，以此来减弱在不同的地图尺度下热力图的形变程度。实操实例：实验截取了工作日20本文档来自技高网...

【技术保护点】
1.一种逆向渲染绘制热力图的方法，其特征在于，包括如下步骤：S10以若干具有时空特性的热力分析数据点为样本，训练核密度分析函数概率模型；S20由图像对应坐标范围及分辨率计算各个像素点对应的地理坐标；S30将各个像素点对应的地理坐标输入核密度分析函数概率模型，拟合模型的核密度估计函数计算获得各个像素中心点的概率密度，即该像素区域的影响力；S40将该像素区域的影响力映射到范围[0，255]的灰度区间内，渲染出像素点的影响力灰度；S50根据像素点的影响力灰度获取色彩梯度的分割阈值；S60依据所得的色彩梯度的分割阈值对像素点的影响力灰度赋色生成热力图。

【技术特征摘要】
1.一种逆向渲染绘制热力图的方法，其特征在于，包括如下步骤：S10以若干具有时空特性的热力分析数据点为样本，训练核密度分析函数概率模型；S20由图像对应坐标范围及分辨率计算各个像素点对应的地理坐标；S30将各个像素点对应的地理坐标输入核密度分析函数概率模型，拟合模型的核密度估计函数计算获得各个像素中心点的概率密度，即该像素区域的影响力；S40将该像素区域的影响力映射到范围[0，255]的灰度区间内，渲染出像素点的影响力灰度；S50根据像素点的影响力灰度获取色彩梯度的分割阈值；S60依据所得的色彩梯度的分割阈值对像素点的影响力灰度赋色生成热力图。2.如权利要求1所述的逆向渲染绘制热力图的方法，其特征在于，所述S50根据像素点的影响力灰度，采用最大熵阈值分割算法获取色彩梯度的分割阈值，最大熵阈值分割算法...

【专利技术属性】
技术研发人员：曾耀国，童昊昕，
申请(专利权)人：航天精一广东信息科技有限公司，
类型：发明
国别省市：广东,44