Shulou(Shulou.com)06/01 Report--

Editor to share with you how to load shp files in matplotlib Basemap. I hope you will get something after reading this article. Let's discuss it together.

Matplotlib is a drawing plug-in in python

Matplotlib supports the effect of two-dimensional map, but also supports the effect of three-dimensional map. In the production of big data, you can take the map as the bottom and add the effect of a bar chart next to it to achieve a three-dimensional effect. There are many ready-made libraries in python, which is more convenient in development and can be directly used by import.

Matplotlib's Basemap contains parsing commonly used shp format vector data, by reading the shp file, it is very convenient to draw on Basemap.

Here is an example of how to read shp and draw the data in shp on a map based on an example of Basemap:

# reference the corresponding class library

Import numpy as np

Import matplotlib.pyplot as plt

From mpl_toolkits.basemap import Basemap as Basemap

# set the coordinate system and coordinate display range of the map

M = Basemap (llcrnrlon=-100.,llcrnrlat=0.,urcrnrlon=-20.,urcrnrlat=57.

Projection='lcc',lat_1=20.,lat_2=40.,lon_0=-60.

Resolution = 'lumbago mills area1 000.

Fig=plt.figure ()

# read the shp file, note that the second parameter is the name of the shp file in the map, which will be used next.

Shp_info = m.readshapefile ('huralll020','hurall',drawbounds=False)

Print (shp_info)

# find the shp file in the map and read the attributes in the shp file

Names = []

For shapedict in m.hurall_info:

Cat = shapedict ['CATEGORY']

Name = shapedict ['NAME']

If cat in ['H4pm] H5' and name not in names:

If name! = 'NOT NAMED': names.append (name)

Print (names)

Print (len (names))

# draw on the map according to the attributes

For shapedict,shape in zip (m. Hurall):

Name = shapedict ['NAME']

Cat = shapedict ['CATEGORY']

If name in names:

Xx,yy = zip (* shape)

# plot for drawing on a map

If cat in ['H4pm and H5']:

M.plot (xx,yy,linewidth=1.5,color='r')

Elif cat in ['H _ 1, H _ 2, and H _ 3']:

M.plot (xx,yy,color='k')

# draw other auxiliary elements on the map, national boundaries, latitude and longitude lines, etc.

M.drawcoastlines ()

M.drawcountries ()

M.drawmapboundary (fill_color='#99ffff')

M.fillcontinents (color='#cc9966',lake_color='#99ffff')

M.drawparallels (np.arange (10, 10, 70, 20), labels= [1, 1, 1, 1, 1, 0, 0, 0])

M.drawmeridians (np.arange (- 100d0jue 20), labels= [0je 0rem 0je 1])

# set the title of the diagram

Plt.title ('Atlantic Hurricane Tracks (Storms Reaching Category 4, 1851-2004)')

Plt.show ()

After reading this article, I believe you have a certain understanding of "how to load shp files in matplotlib Basemap". If you want to know more about it, you are welcome to follow the industry information channel. Thank you for reading!

Welcome to subscribe "Shulou Technology Information " to get latest news, interesting things and hot topics in the IT industry, and controls the hottest and latest Internet news, technology news and IT industry trends.