FilesExpand file tree

 master

/

Python_BaseMap1.py

Copy path

Blame

More file actions

Blame

More file actions

Latest commit

 

History

History

History

96 lines (88 loc) · 3.04 KB

 master

/

Python_BaseMap1.py

Top

File metadata and controls

• Code
• Blame

96 lines (88 loc) · 3.04 KB

Raw

Copy raw file

Download raw file

Open symbols panel

Edit and raw actions

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

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()