地图经纬度 区域平分

业务需求：

• 需要在百度地图 某一块区域上 例如 整个四川的区域
• 取四川区域的东南角 跟 西北角的 经纬度
• 在整个四川省的区域进行平分，平分成10*10的网格
• 获取每个网格的中心经纬度坐标
• 目的是为了算出每个网格里大概人数有多少
• 如图所示：
  

函数实现： 其实就是等差数列

type Aera struct {
	NorthEast *Point `json:"north_east"` //地图东北角经纬度
	SouthWest *Point `json:"south_west"` //地图西南角经纬度
}

type Point struct {
	Lng   float64 `json:"lng"`
	Lat   float64 `json:"lat"`
}

// 区域平分成 step x step 份
func heatmapSplit(aera *Aera, step float64) [][]*Point {
	lng0, lng1 := aera.SouthWest.Lng, aera.NorthEast.Lng
	lat0, lat1 := aera.SouthWest.Lat, aera.NorthEast.Lat
	lngD := (lng1 - lng0) / step
	latD := (lat1 - lat0) / step
	dis := math.Min(lngD, latD)
	maxi := int(math.Ceil((lng1 - lng0) / dis))
	maxj := int(math.Ceil((lat1 - lat0) / dis))
	points := make([][]*Point, 0)
	for i := 0; i < maxi; i++ {
		lng10 := lng0 + dis*float64(i)
		lng11 := lng0 + dis*float64(i+1)
		for j := 0; j < maxj; j++ {
			lat10 := lat0 + dis*float64(j)
			lat11 := lat0 + dis*float64(j+1)
			points = append(points, []*Point{
				{
					Lng: lng11,
					Lat: lat11,
				},
				{
					Lng: lng11,
					Lat: lat10,
				},
				{
					Lng: lng10,
					Lat: lat10,
				},
				{
					Lng: lng10,
					Lat: lat11,
				},
				{
					Lng: (lng10 + lng11) / 2,
					Lat: (lat10 + lat11) / 2,
				}, // 中心点经纬度
			})
		}
	}
	return points
}