ImageVerifierCode 换一换
格式:DOC , 页数:3 ,大小:100.01KB ,
资源ID:11670016      下载积分:10 金币
快捷下载
登录下载
邮箱/手机:
温馨提示:
快捷下载时,用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)。 如填写123,账号就是123,密码也是123。
特别说明:
请自助下载,系统不会自动发送文件的哦; 如果您已付费,想二次下载,请登录后访问:我的下载记录
支付方式: 支付宝    微信支付   
验证码:   换一换

加入VIP,免费下载
 

温馨提示:由于个人手机设置不同,如果发现不能下载,请复制以下地址【https://www.docduoduo.com/d-11670016.html】到电脑端继续下载(重复下载不扣费)。

已注册用户请登录:
账号:
密码:
验证码:   换一换
  忘记密码?
三方登录: 微信登录   QQ登录   微博登录 

下载须知

1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。
2: 试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。
3: 文件的所有权益归上传用户所有。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 本站仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

版权提示 | 免责声明

本文(Python与开源GIS:根据空间过滤器选择要素.doc)为本站会员(HR专家)主动上传,道客多多仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知道客多多(发送邮件至docduoduo@163.com或直接QQ联系客服),我们立即给予删除!

Python与开源GIS:根据空间过滤器选择要素.doc

1、Python与开源GIS:根据空间过滤器选择要素2012-12-03 21:12作者:bukunu来源:本站浏览: 98 views我要评论字号:大中小摘要: 空间过滤器Spatial filters 如果说按照属性进行要素的选择,还是带有数据库的特征的话,那么,根据空间位置进行选择,则是地地道道的GIS功能。在OGR中,使用了Spatial filters(空间过滤)这一术语表征这一功能。 OGR提供的空间过滤功能.空间过滤器Spatial filters如果说按照属性进行要素的选择,还是带有数据库的特征的话,那么,根据空间位置进行选择,则是地地道道的GIS功能。在OGR中,使用了Spati

2、al filters(空间过滤)这一术语表征这一功能。OGR提供的空间过滤功能有两种,一种是SetSpatialFilter(),过滤某一类型的Feature,例如参数中使用Polygon,就是选出Layer中的所有Polygon所覆盖的要素(注意,只要相交即可,不必完全包含)。下面这段代码用了两套数据。world_borders是全球国界数据,cover.shp则是覆盖了非洲南部地区的一个多边形。下面先定义一个根据图层直接生成shape文件的函数,方便后面调用。12345678from osgeo import ogrdef create_shp_by_layer(shp, layer):o

3、utputfile = shpif os.access(outputfile, os.F_OK):driver.DeleteDataSource(outputfile)newds = driver. CreateDataSource ( outputfile )pt_layer = newds.CopyLayer ( layer, )newds.Destroy ()下面代码是使用cover.shp中的多边形来选择全球国界数据:1234567891011121314driver = ogr.GetDriverByName(ESRI Shapefile)world_shp = /gdata/wor

4、ld_borders.shpcover_shp = /gdata/cover.shpworld_ds = ogr.Open(world_shp)cover_ds = ogr.Open(cover_shp)world_layer = world_ds.GetLayer(0)cover_layer = cover_ds.GetLayer(0)print(world_layer.GetFeatureCount()cover_feats = cover_layer.GetNextFeature()poly = cover_feats.GetGeometryRef()world_layer.SetSpa

5、tialFilter(poly)out_shp = /gdata/world_cover.shpcreate_shp_by_layer(out_shp, world_layer)结果可以通过下面的图来看。另外还有SetSpatialFilterRect(minx, miny, maxx, maxy),参数输入四个坐标,可以选中矩形内的要素。123world_layer.SetSpatialFilterRect(50, 60, 25, 35)out_shp = /gdata/world_spatial_filter.shpcreate_shp_by_layer(out_shp, world_layer)同样可以打开GIS软件来查看结果。SetSpatialFilter(None)则是清空空间属性过滤器。这个可以通过打印图层中要素的数目来查看。

本站链接:文库   一言   我酷   合作


客服QQ:2549714901微博号:道客多多官方知乎号:道客多多

经营许可证编号: 粤ICP备2021046453号世界地图

道客多多©版权所有2020-2025营业执照举报