- 积分
- 55955
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-6-21
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
发布MeteoInfo 1.3.2。引入了JTS Topolygo Suite代码,JTS提供了一套空间数据操作的核心算法,是很多开源GIS软件的核心部件。引入JTS可以极大的增强MeteoInfo的空间数据操作、分析功能。新版本增强了图层编辑功能,增加了剪切、复制、粘贴图元要素,给多边形增加洞,分隔和合并要素的功能。增加了空间分析菜单,目前只提供了图层缓冲区分析(buffer)和裁剪分析,其它空间分析功能会在后续版本添加。
在脚本中增加了topology模块,引入该模块可以进行要素的空间计算:intersect, union, difference, symdifference, buffer, convexhull, split。
- import mipylib.topology as tp
- axesm()
- lworld = shaperead('D:/Temp/map/country1.shp')
- geoshow(lworld, edgecolor='k', facecolor='g')
- #Add line
- lat = [15, 0, -45, -25]
- lon = [-100, 0, 70, 110]
- line1 = geoshow(lat, lon, size=2, color='r')
- buf1 = tp.buffer(line1,5)
- geoshow(buf1, color='y')
- geoshow(lat, lon, size=2, color='r')
- #Add polygon
- lat = array([30, 0, 18, 48, 30])
- lon = array([60, 70, 130, 120, 60])
- g1 = geoshow(lat, lon, displaytype='polygon', color=[150,230,230,230], edgecolor='b', size=2)
- lat = lat + 10
- lon = lon + 10
- g2 = geoshow(lat, lon, displaytype='polygon', color=[150,230,230,230], edgecolor='b', size=2)
- g3 = tp.intersect(g1, g2)
- geoshow(g3, color='r')
- #Set extent
- xlim(-180, 180)
- ylim(-90, 90)
- xticks(arange(-180, 181, 30))
- yticks(arange(-90, 91, 30))
- title('Buffer and intersection')
|
|