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

python画地形地貌图_opencv检测瑕疵python我们可以使用basemap这个工具包来实现中国地图的绘制首先需要加载一些包:importnumpyasnpimportmatplotlib.pyplotaspltfrommpl_toolkits.basemapimportbasemapbasemap包就是气象画图的利器,现在我们就可以愉快的画图了!plt.figure(1)map=basemap()map.drawcoastli…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

我们可以使用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()这行代码是用来显示图片的)

6d92e9eda134fc6693384ebc994c7fa6.png

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

map.drawcountries()

13cdc6718ce865f6f21eb19e90c2c093.png

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

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

b37d83d80de182404191ea055a454446.png

好了,现在我们可以开始画中国地图了!

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

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

6ecacd3f7fef4a6de3228029fbfe7d57.png

好了,一个中国地图就出来了!但是我们发现,好像少了点什么,没错就是省界。我们可以在下载中国大陆和台湾省的行政区域的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)

5e69b0422c952cdc2bb632e24f0e86c6.png

还可以在地图上加上经纬度,比如我们要画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])

38c937f0a24d2632ea829d27528517fc.png

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

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

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

98dc34f1065a68493b510f00096b7863.png

这回正式的完成了!

附上所有代码:

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 绘制地图的资料请关注萬仟网其它相关文章!

希望与广大网友互动??

点此进行留言吧!

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/227279.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 如何将ER图转化为关系模型(超详细、含例题)「建议收藏」

    如何将ER图转化为关系模型(超详细、含例题)「建议收藏」本文转载自:http://blog.csdn.net/HaoDaWang/article/details/78098937?locationNum=4&fps=1超级感谢博主分享本篇博文中的例题只涉及到N:M和N:M:P的转换,大家看完后可以再去看下百度文库里面的这个例题,这个例题几个转换准则都用到了。https://wenku.baidu.com/view/33849a73a…

    2022年9月22日
    1
  • 不用go-(xxoo)agent的也有可能是程序员

    不用go-(xxoo)agent的也有可能是程序员记得是去年某个时候,在一个群里问大家是怎么翻(什么)墙的,发现很多人是用go-(xxoo)agent的,而我当时完全没听说过go-agent,因为我已经N多年不翻(什么)墙了,于是,有人就说:“连go-(xxoo)agent都没用过你还是程序员吗?”说实话,我也不希望我是程序员啊!可我真的没用过go-(xxoo)agent啊!想起我一位前同事,工作很多年了,却从来没用过linux!(也不知道他

    2022年7月17日
    21
  • Java: Hook技术[通俗易懂]

    Java: Hook技术[通俗易懂]一、什么是Hook  Hook翻译成中文就是勾子的意思,在java中它表示在事件到达终点前进行拦截或监控的一种行为。 二、Hook的实现实现hook我们必须要知道java的反射和动态代理。 1、反射  反射是java的查看、检测、修改自身的一种行为。  在编译阶段,编译器将我们编写的java文件编译成.class文件。而在运行期,jvm又将.class文件通过…

    2022年5月9日
    145
  • Ext.apply(src,apply) 和 Ext.applyIf(src,apply)比较(转)

    Ext.apply(src,apply) 和 Ext.applyIf(src,apply)比较(转)Ext.onReady(function(){/**Ext.apply(src,apply)和Ext.applyIf(src,apply)两个方法的使用和区别比较*///Ext.app

    2022年7月2日
    22
  • 【用python编写一个简单的单线程wifi暴力破解工具】

    源代码a.txt:密码文件crack.py:wifi破解模块main.py:主模块scan.py:wifi扫描模块scan.pyimportpywifiimporttime#WiFi扫描模块defwifi_scan():#初始化wifiwifi=pywifi.PyWiFi()#使用第一个无线网卡interface=wifi.interfaces()[0]#开始扫描interface.scan()

    2022年4月9日
    49
  • 深入浅出理解卷积运算

    深入浅出理解卷积运算提起卷积运算相信大家都不陌生 这是一种很常见的运算 我们在学习 信号与系统 时就一直在和卷积打交道 在后来的一些课程中也有卷积运算的身影 比如 自动控制原理现代部分 中的卷积定理等 在学习 信号与系统 时我们知道了卷积的定义 对于两个函数 f x f x f x 和 g x g x g x 他们的卷积 f g n f g n f g n 的公式如下 连续形式 f nbsp g n f g n d mathrm f mathrm g n int infty

    2025年8月27日
    4

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号