制作一个中国蔬菜产量图

微博上看到的数据

import geopandas as gpd

china = gpd.read_file("MLgis/feifuli3_geometry/china_provinces.json")

albers_proj = '+proj=aea +lat_1=25 +lat_2=47 +lon_0=105'

veg_product = pd.read_json('{"2020\\u5e7431\\u4e2a\\u7701\\u4efd\\u852c\\u83dc\\u4ea7\\u91cf":{"0":"\\u5c71\\u4e1c","1":"\\u6cb3\\u5357","2":"\\u6c5f\\u82cf","3":"\\u6cb3\\u5317","4":"\\u56db\\u5ddd","5":"\\u6e56\\u5317","6":"\\u6e56\\u5357","7":"\\u5e7f\\u897f","8":"\\u5e7f\\u4e1c","9":"\\u8d35\\u5dde","10":"\\u4e91\\u5357","11":"\\u5b89\\u5fbd","12":"\\u91cd\\u5e86","13":"\\u8fbd\\u5b81","14":"\\u9655\\u897f","15":"\\u6d59\\u6c5f","16":"\\u65b0\\u7586","17":"\\u6c5f\\u897f","18":"\\u798f\\u5efa","19":"\\u7518\\u8083","20":"\\u5185\\u8499\\u53e4","21":"\\u5c71\\u897f","22":"\\u9ed1\\u9f99\\u6c5f","23":"\\u6d77\\u5357","24":"\\u5b81\\u590f","25":"\\u5409\\u6797","26":"\\u5929\\u6d25","27":"\\u4e0a\\u6d77","28":"\\u9752\\u6d77","29":"\\u5317\\u4eac","30":"\\u897f\\u85cf"},"\\u5355\\u4f4d\\uff1a\\u4e07\\u5428":{"0":8434.7,"1":7612.4,"2":5728.1,"3":5198.2,"4":4813.4,"5":4119.4,"6":4110.1,"7":3830.8,"8":3706.8,"9":2990.9,"10":2507.9,"11":2330.9,"12":2092.6,"13":1960.0,"14":1957.7,"15":1945.5,"16":1714.9,"17":1642.7,"18":1630.2,"19":1478.5,"20":1075.1,"21":861.2,"22":674.3,"23":572.8,"24":566.4,"25":464.9,"26":266.5,"27":252.9,"28":151.4,"29":137.9,"30":84.3}}')


veg_product.columns = ['prov', 'product']

china['prov_name'] = china['NAME'].str.replace('省', '')
china['prov_name'] = china['prov_name'].str.replace('市', '')
china['prov_name'] = china['prov_name'].str.replace('自治区', '')
china['prov_name'] = china['prov_name'].str.replace(r'壮族|回族|维吾尔', '')
china['prov_name']

data_with_geometry = pd.merge(left = veg_product,
                              right = china,
                              left_on='prov',
                              right_on = 'prov_name',
                              how = 'left')

data_with_geometry = gpd.GeoDataFrame(data_with_geometry, crs = 'EPSG:4326')
fig, ax = plt.subplots(figsize=(12, 8))
ax = data_with_geometry.to_crs(albers_proj).plot(ax=ax,
                                                 column = 'product',
                                                 cmap = 'Greens',
                                                 legend = True,
                                                 scheme = 'NaturalBreaks',
                                                 k =10,
                                                 legend_kwds = {
                                                     "loc":"lower left",
                                                     "title":'中国蔬菜产量',
                                                     "shadow":True
                                                 })

成图效果如下：