Matlab下地理热度图的绘制

    Matlab自带有强大的Mapping工具箱，但网上的相关信息较少，本文通过分享Matlab下地理热度图的绘制方法，介绍Matlab地理工具箱的部分使用方法。

    本文使用的数据是我国三大运营商的部分移动基站位置数据，格式为（lngt,lat）。我们将地图分割分等大小的地理栅格（Grid），Matlab里称为Raster，通过统计每个栅格内的基站数量，得到基站的热度分布，如下图所示。

    

    算法的实现代码如下：

clear
close all
clc

load BSposition.mat %基站数据， 存储在result中，格式为（lngt,lat）

% 将每度分割成的网格数量，取20相当于3min一个网格
Scale=20; 

China=shaperead('maps/bou1_4p.shp','UseGeoCoords',true); %读入中国的国界文件shp文件
lat=[guo.Lat];
lon=[guo.Lon];
[Z,refvec]=vec2mtx(lat,lon,Scale,'filled');%得到中国的国土栅格图Z，refvec 是一个三元素的参考向量，表示为[cells/degree north-latitude west-longitude]，用于将经纬度值转换为网格的行列号

% 构建初始热度图heat，热度值均为0