ogr

I'm trying to extract data from an OSM.PBF file using the python GDAL/OGR module. Currently my code looks like this: import gdal, ogr osm = ogr.Open('file.osm.pbf') ## Select multipolygon from the layer layer = osm.GetLayer(3) # Create list to store pubs pubs = [] for feat in layer: if feat.GetField('amenity') == 'pub': pubs.append(feat) While this little bit of code works fine with small.pbf files (15mb). However, when parsing files larger than 50mb I get the following error: ERROR 1: Too many features have accumulated in points layer. Use OGR_INTERLEAVED_READING=YES MODE When I turn this

How can I get the installed GDAL/OGR version from python? I aware of the gdal-config program and are currently using the following: In [3]: import commands In [4]: commands.getoutput('gdal-config --version') Out[4]: '1.7.2' However, I suspect there is a way to do this using the python API itself. Any dice? gdal.VersionInfo() does what I want: >>> osgeo.gdal.VersionInfo() '1604' This works on both my Windows box and Ubuntu install. gdal.__version__ gives an error on my Windows installation, although it works on my Ubuntu installation: >>> import osgeo.gdal >>> print osgeo.gdal.__version__

I have a raster file and a WGS84 lat/lon point. I would like to know what value in the raster corresponds with the point. My feeling is that I should use GetSpatialRef() on the raster object or one of its bands and then apply a ogr.osr.CoordinateTransformation() to the point to map it to the raster's space. My hope would then be that I could simply ask the rasters' bands what is at that point. However, the raster object doesn't seem to have a GetSpatialRef() or a way to access a geo-located point, so I'm somewhat at a loss for how to do this. Any thoughts? Say i have a geotiff file test.tif.

I am seeking a way to take an existing ESRI Shapefile that has a Feature set of 200 countries. Each country Feature has an attribute of "NAME." My objective is to create a Python script that adds an arbitrary (for now) additional attribute, say, "POPULATION". Of course I have the OSGeo and GeoDjango modules installed. I'm as far as: from osgeo import ogr infile = ogr.Open('sample.shp', 1) #'sample.shp' is a pre-existing ESRI shapefile described above inlyr = infile.GetLayerByIndex(0) Am I missing an OGR function that will allow me to insert Feature attribute fields into an existing Shapefile?

OK, After a bit of fiddling, I've tweaked a script from the site hyperlink in the second comment line. The purpose of the script is to clip/mask a LARGE raster (i.e. that cannot fit into a 32-bit Python 2.7.5 application) in GTiff format with a shapefile with multiple polygons (each with a "Name" record) and save the clipped rasters into a "clip" sub-directory, where each masked grid is named after each polygon's "Name". Like the original script, it assumes that the GTiff and shapefile are in the same projection and overlap correctly, and it processes ~100 masks/sec. However, I have modified

说明： 这是一份来自网上的学习笔记，教程来自于犹他州立大学，课程名称为：Geoprocessing with Python using Open Source GIS 。 这份教程在网上有几个网站进行过发布，也无从判断哪份才是原创。作为一份学习笔记，里面很多技术细节都语焉不详，但是不管怎样，这份资料都是很难得的。 开放地理空间实验室进行了整理，修改了一些细节，并进行重新发布。感谢原作者的工作。 地理数据处理软件包GDAL简介 GDAL OGR GDAL python教程（1）——用OGR读写矢量数据 为什么用open source？ Open source RS/GIS模块 相关模块 导入库： 读取数据层 释放内存 读完了再说怎么写 GDAL python教程（2）——几何形状geometry与投影projection 建立新的几何形状 创建复合几何形状multi geometry GDAL python教程（3）——过滤器，简单的空间分析，函数和模块 属性过滤器Attribute filters 空间过滤器Spatial filters 统计每种cover类型各有多少个Feature 空间操作 下面看看简单的地理数据处理geoprocessing GDAL python教程（4）——用GDAL读取栅格数据 导入GDAL支持库 GDAL python教程（5）—