python 绘制中国地图并利用经纬度标注散点

python 绘制中国地图并利用经纬度标注散点

所需要的包：GeoPandas，安装教程有很多，自行百度即可。

用到的中国地级市shp文件：链接：https://pan.baidu.com/s/18aaxczrz4tIRMeCusOrDQA 
                                              提取码：rav1 

一、GeoPandas类简单介绍

  GeoPandas实现了两个主要的数据结构，GeoSeries和GeoDataFrame。它们分别是pandas中Series和DataFrame的子类。实际上，如果你了解DataFrame结构，这个就很好理解了。

  一个GeoSeries就是包含一个几何图形的序列，比如多边形、点。具体形式就在下图，如果想要绘图，很简单，只要用plot方法就可以，多边形GeoSeries绘制出来的就是多边图形， 点GeoSeries当然就是绘制点了。

多边形GeoSeries：

 点GeoSeries：

绘图实例：

 读取shp文件绘制中国地图

china_map=gp.GeoDataFrame.from_file("C:/Users/Administrator/Desktop/map_shp/shipnew
                               /中国行政区划2004_shp/dijishi_2004.shp",\
                                    encoding = 'gb18030')
china_map.plot(figsize=(20,12),color="white", edgecolor="black")

 

利用经纬度绘制散点图：

lng = data2002com['lng']
lat= data2002com['lat']
pts = gp.GeoSeries([Point(x, y) for x, y in zip(lng, lat)])
 
pts.plot(figsize=(12,8))

 

二、世界地图标注实例

  这一部分主要是解释官方网站的一个例子，为自己绘制地图标注散点做准备。

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
import os
import geopandas as gp
from shapely.geometry import Point

# 官方数据
world = gp.read_file(gp.datasets.get_path('naturalearth_lowres'))
cities = gp.read_file(gp.datasets.get_path('naturalearth_cities'))

#注意world.shape =(177,6)
#  cities.shape = (202,2)

fig, ax = plt.subplots()

ax.set_aspect('equal')

world.plot(ax=ax, color='white', edgecolor='black')
cities.plot(ax=ax, marker='o', color='red', markersize=5)
plt.show()

 

   绘图的机制其实很好理解，就是第一个图层是一个世界地图的多边形，第二个图层是一些城市的经纬度点，把这些经纬度点映射到世界地图上就是这个样了。

三、中国地图标注实例

只需要地图数据的shp文件，还有就是经纬度，就可以实现自己的地图。

china_map=gp.GeoDataFrame.from_file("C:/Users/Administrator/Desktop/map_shp/shipnew
                                      /中国行政区划2004_shp/dijishi_2004.shp",\
                                    encoding = 'gb18030')
data2002com = pd.read_excel(r'C:\Users\Administrator\Desktop\经纬度数据\按公司名称 
                                    \2002.xlsx')

# 几何图形
geo_ploy = china_map['geometry']
# 地图点
geo_point = gp.GeoSeries([Point(x, y) for x, y in zip(lng, lat)])

fig, ax = plt.subplots(figsize=(12,8))

ax.set_aspect('equal')
# 几何图形绘制
geo_ploy.plot(ax=ax, color='white', edgecolor='black')

# 地图点标注
geo_point.plot(ax=ax, marker='o', color='black', markersize=0.1)