python动态演示数据gdp_利用Python制作中国GDP分布图和动态演示

利用Python制作中国GDP分布图和动态演示

数据读取

## 导入相关模块

import pandas as pd

import geopandas as gpd

import numpy as np

import matplotlib.pyplot as plt

from shapely.geometry import Point

import matplotlib.patches as mpatches

from mpl_toolkits.basemap import Basemap

from matplotlib_scalebar.scalebar import ScaleBar

plt.rcParams['font.family'] = ['sans-serif']

plt.rcParams['font.sans-serif'] = ['SimHei']# 替换sans-serif字体为黑体

plt.rcParams['axes.unicode_minus'] = False # 解决坐标轴负数的负号显示问题

regibns = gpd.GeoDataFrame.from_file("F:\ArcGIS\ArcGIS文件\全国\中国地图\省级行政区.shp")

regibns.plot()

## 查看坐标系

regibns.crs

## 绘制国界线

nine_lines = gpd.read_file('F:\ArcGIS\ArcGIS文件\全国\中国地图\国界线.shp')

nine_lines.plot()

## 经纬网

jingwei = gpd.GeoDataFrame.from_file('F:\ArcGIS\ArcGIS文件\全国\中国地图\经纬网.shp')

fig = plt.figure(figsize=(8,8)) #设置画布大小

ax = plt.gca()

regibns.plot(ax=ax)

jingwei.plot(ax=ax)

nine_lines.plot(ax=ax)

taiwan = gpd.GeoDataFrame.from_file("F:\ArcGIS\ArcGIS文件\全国\中国地图\省级行政区.shp")

taiwan = taiwan[regibns['NAME'].isin(['台湾'])]

taiwan.plot()

数据清洗

regibns = regibns[['GDP_1994(','GDP_1997(','GDP_1998(','GDP_1999(&