- 积分
- 55951
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-6-21
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
利用geoshow()函数可以在地图上添加自定义点、线和多边形。
脚本程序:
- axesm(tickfontsize=12)
- lworld = shaperead('D:/Temp/map/country1.shp')
- geoshow(lworld, edgecolor='k', facecolor='g')
- #Add point
- geoshow(15, -100, size=14, color='r', marker='S')
- #Add line
- lat = [15, 0, -45, -25, nan, 15, -30, -60, -65]
- lon = [-100, 0, 70, 110, nan, -100, -20, 100, 150]
- geoshow(lat, lon, size=2, color='b')
- #Add line and polygon
- lat = array([30, 0, 18, 48, 30])
- lon = array([60, 70, 130, 120, 60])
- geoshow(lat, lon, color='r', size=2, linestyle=':')
- lon = lon - 100
- geoshow(lat, lon, displaytype='polygon', color=[150,230,230,230], edgecolor='b', size=2)
- #Set extent
- xlim(-180, 180)
- ylim(-90, 90)
- xticks(arange(-180, 181, 30))
- yticks(arange(-90, 91, 30))
- title('Add shapes')
|
|