matlab怎么栅格化地图,用栅格化的思维进行地图散点数据绘制

原理上，这个课题利用point做粒子效果就可以解决，本不需要通过shader去处理。不过作为用gl的思维进行地图栅格数据绘制课题的后续，还是提出以下。本方案的核心思路还是数据栅格化，降纬处理数据减少碰撞测试次数提升性能。

Demo：http://lrdcq.com/test/mapwebgltiles/scatter/

思路

对于地图上的散点数据，我们还是将他们落到格子里：

1. 将地图区域划为MxN个格子，每个格子里会有0-x个点的经纬度数据，如果划分得相对合理，我们的x将极小(demo里总数据量约7w+，x为6)。

2. 假设我的栅格的大小为(x,y)，如果我的每个点的绘制区域的大小也正好是(m,n)=(x,y)，如下图

那我每一个栅格里点的绘制区域都有可能覆盖到3x3共计9个栅格的像素。同理，对于任意(m,n)的绘制区域，我们需要检查(ceil(m/x) * 2 + 1, ceil(n/y) * 2 + 1)的区域的像素。那对应的，对于我们的shader处理每一个像素时，我们需要检测周围这么大的区域的栅格里是否有数据。

3. 对每个像素处理，对拿到的可能覆盖范围内的所有数据点做矩形(即绘制区域)碰撞测试，如果到当前像素在某个数据点的绘制范围内，则取对应图片进行绘制。

相对于纯栅格绘制，这个散点的栅格绘制方案的成本包括多检索(ceil(m/x) * 2 + 1, ceil(n/y) * 2 + 1)的区域，并且每个区域需要做做多0-x次碰撞测试。以demo的数据来说，性能消耗最多是3x3x6攻击54倍——当然对于shader来说，这个也不算什么的样子。

Code

这里有一个思路转换。原本方案的数据图像(即一个像素表示一个栅格的图片)，每个像素里根本没存储啥东西，所以像素数据非常无所谓。但是现在这个方案，我们需要在数据图像中储存这个像素，就是这个栅格里包含了多少个点，还有每一个点的经纬度。

1. 储存经纬度，一个像素0xFFFFFFFF储存经纬度可以储存到小数点后7位了。而如果用rgb，不用a的六位0xFFFFFF来储存经纬度，可以储存到小数点后5位。5位的进度已经是米级的了，所以五位就够了。所以我们可以用两个像素0xAABBCCFF,0xAABBCCFF来足够准确的表述以个数据点的经纬度还可以储存256*256的其他信息。完全够了。

2. 因为我们的数据图像的像素要表示栅格，因此我们开了第二个图片，数据图像2，来储存所有的经纬度点数据。只需要两个两个一组把所有数据按即有顺序罗列在图片里面就行了。

生成数据图像2的数据写入片段是&#