submit code

Author: beginner-fei

chaoxi/BaseMap/Python_BaseMap1.py

@@ -0,0 +1,96 @@
+import matplotlib.pyplot as plt
+from mpl_toolkits.basemap import Basemap
+import numpy as np
+
+# 画海岸线
+def drawcoast():
+
+    plt.figure(figsize=(12, 8))
+    m = Basemap()   # 创建一个地图
+    m.drawcoastlines()   # 画海岸线
+    plt.show()   # 显示图像
+
+# 地球基本画法
+def draw_basic():
+    map = Basemap(projection='ortho', lat_0=0, lon_0=0)
+    map.drawmapboundary(fill_color='aqua')
+    map.fillcontinents(color='gray',lake_color='aqua')
+    map.drawcoastlines()
+    plt.show()
+
+# 画中国地图
+def draw_china():
+    plt.figure(figsize=(10, 6))
+    m = Basemap(llcrnrlon=77, llcrnrlat=14, urcrnrlon=140, urcrnrlat=51, projection='lcc', lat_1=33, lat_2=45,
+                lon_0=100)
+    m.drawcountries(linewidth=1.5)
+    m.drawcoastlines()
+    plt.show()
+
+# 画地球人口分布图
+def drawearth():
+    names = []
+    pops = []
+    lats = []
+    lons = []
+    countries = []
+    file = open("data/main_city", encoding='utf-8').readlines()
+    for line in file:
+        info = line.split()
+        names.append(info[0])
+        pops.append(float(info[1]))
+        lat = float(info[2][:-1])
+        if info[2][-1] == 'S': lat = -lat
+        lats.append(lat)
+        lon = float(info[3][:-1])
+        if info[3][-1] == 'W': lon = -lon + 360.0
+        lons.append(lon)
+        country = info[4]
+        countries.append(country)
+    # set up map projection with
+    # use low resolution coastlines.
+    map = Basemap(projection='ortho', lat_0=35, lon_0=120, resolution='l')
+    # draw coastlines, country boundaries, fill continents.
+    map.drawcoastlines(linewidth=0.25)
+    map.drawcountries(linewidth=0.25)
+    # draw the edge of the map projection region (the projection limb)
+    map.drawmapboundary(fill_color='#689CD2')
+    # draw lat/lon grid lines every 30 degrees.
+    map.drawmeridians(np.arange(0, 360, 30))
+    map.drawparallels(np.arange(-90, 90, 30))
+    # Fill continent wit a different color
+    map.fillcontinents(color='#BF9E30', lake_color='#689CD2', zorder=0)
+    # compute native map projection coordinates of lat/lon grid.
+    x, y = map(lons, lats)
+    max_pop = max(pops)
+    # Plot each city in a loop.
+    # Set some parameters
+    size_factor = 80.0
+    y_offset = 15.0
+    rotation = 30
+    for i, j, k, name in zip(x, y, pops, names):
+        size = size_factor * k / max_pop
+        cs = map.scatter(i, j, s=size, marker='o', color='#FF5600')
+        plt.text(i, j + y_offset, name, rotation=rotation, fontsize=10)
+
+    plt.title('earth')
+    plt.show()
+
+# 画带投影的地球图片
+def draw_earth1():
+    import matplotlib.pyplot as plt
+    from mpl_toolkits.basemap import Basemap
+    plt.figure(figsize=(8, 8))
+    # 正射投影，投影原点设在了上海周边
+    m = Basemap(projection='ortho', resolution=None, lat_0=30, lon_0=120)
+    # 图像原始分辨率是5400*2700，设置scale = 0.5以后分辨率为2700*1350,如此作图
+    # 迅速不少也不那么占用内存了
+    m.bluemarble(scale=0.5)
+    plt.show()
+
+if __name__ == '__main__':
+    #drawcoast()
+    drawearth()
+    #draw_basic()
+    #draw_china()
+    #draw_earth1()

chaoxi/BaseMap/main_city

@@ -0,0 +1,9 @@
+Shanghai 23019148  31.23N  121.47E  China
+Mumbai   12478447  18.96N  72.82E   India
+Karachi  13050000  24.86N  67.01E   Pakistan
+Delhi    16314838  28.67N  77.21E   India
+Manila   11855975  14.62N  120.97E  Philippines
+Seoul    23616000  37.56N  126.99E  Korea(South)
+Jakarta  28019545   6.18S  106.83E  Indonesia
+Tokyo    35682460  35.67N  139.77E  Japan
+Peking   19612368  39.91N  116.39E  China

chaoxi/README.md

@@ -9,6 +9,8 @@ Python技术 公众号文章代码库
 
 ## 实例代码
 
+[神器-可视化分析之Basemap实战详解（二）](https://github.com/JustDoPython/python-examples/tree/master/chaoxi/BaseMap) 神器-可视化分析之Basemap实战详解（二）
+
 [用 Python 给小表弟画皮卡丘！](https://github.com/JustDoPython/python-examples/tree/master/chaoxi/pkq) 用 Python 给小表弟画皮卡丘！
 
 [惊艳！利用 Python 图像处理绘制专属头像](https://github.com/JustDoPython/python-examples/tree/master/chaoxi/opencv_img) 惊艳！利用 Python 图像处理绘制专属头像