python画地形地貌图_opencv检测瑕疵python

我们可以使用basemap这个工具包来实现中国地图的绘制

首先需要加载一些包：

import numpy as np

import matplotlib.pyplot as plt

from mpl_toolkits.basemap import basemap

basemap包就是气象画图的利器，现在我们就可以愉快的画图了！

plt.figure(1)

map=basemap()

map.drawcoastlines()

plt.title(r’$world\ map$’,fontsize=24)

plt.show()

第2行创建一个地图，第3行添加海岸线，这样一个世界地图就出来了，怎么样，很简单吧。(plt.show()这行代码是用来显示图片的)

我们发现这只是海岸线图，那么怎么将国界线添加上去呢？很简单，只要添加一行代码就可以了。

map.drawcountries()

那么怎么添加河流呢？可能有些同学已经猜到了，就是drawrivers()

map.drawrivers(color=’blue’,linewidth=0.3)

好了，现在我们可以开始画中国地图了！

其实只要在创建地图时指定一下范围就可以了，查阅资料发现，中国的经纬度范围是东经135度2分30秒-东经73度40分，北纬3度52分-北纬53度33分。

map=basemap(llcrnrlon=70,llcrnrlat=3,urcrnrlon=139,urcrnrlat=54)

好了，一个中国地图就出来了！但是我们发现，好像少了点什么，没错就是省界。我们可以在下载中国大陆和台湾省的行政区域的shape文件，下载后解压，然后加入下面的代码。

chn=’g:\python_material\mapofchina’

chn的值就是解压后的地图文件所在的地址。

下面我们就可以加入省界了！

map.readshapefile(chn+’\gadm36_chn_shp\gadm36_chn_1′,

‘states’,drawbounds=true)

别忘了把台湾省加上去

map.readshapefile(chn+’\gadm36_twn_shp\gadm36_twn_1′,

‘taiwan’,drawbounds=true)

还可以在地图上加上经纬度，比如我们要画5条经纬线，可以这么做：

parallels = np.linspace(3,55,5)

map.drawparallels(parallels,labels=[true,false,false,false])

meridians = np.linspace(70,140,5)

map.drawmeridians(meridians,labels=[false,false,false,true])

大功告成！但是，emmm，我们发现好像有点歪？我们可以在创建地图时选择投影参数。

map=basemap(llcrnrlon=70,llcrnrlat=3,urcrnrlon=137,urcrnrlat=54,

projection = ‘lcc’, lat_1 = 33, lat_2 = 45, lon_0 = 100)

这回正式的完成了！

附上所有代码：

import numpy as np

import matplotlib.pyplot as plt

from mpl_toolkits.basemap import basemap

plt.figure(1)

map=basemap(llcrnrlon=70,llcrnrlat=3,urcrnrlon=139,urcrnrlat=54,

projection = ‘lcc’, lat_1 = 33, lat_2 = 45, lon_0 = 100)

map.drawcoastlines()

map.drawcountries()

map.drawrivers(color=’blue’,linewidth=0.3)

chn=’g:\python_material\mapofchina’

map.readshapefile(chn+’\gadm36_chn_shp\gadm36_chn_1′,

‘states’,drawbounds=true)

map.readshapefile(chn+’\gadm36_twn_shp\gadm36_twn_1′,

‘taiwan’,drawbounds=true)

parallels = np.linspace(3,55,5)

map.drawparallels(parallels,labels=[true,false,false,false])

meridians = np.linspace(70,140,5)

map.drawmeridians(meridians,labels=[false,false,false,true])

plt.title(r’$china\ map$’,fontsize=24)

plt.show()

以上就是利用python绘制中国地图的详细内容，更多关于python 绘制地图的资料请关注萬仟网其它相关文章！

希望与广大网友互动？？

点此进行留言吧！