Question

有关如何将this shapefile剪辑到相应边界框的任何建议吗？特别是，我只想查看这些界限：[（-105.53,39.9）（ - 1044,40.27）]。 site for ogr2ogr有一个使用看起来很有希望的-clipdst选项的示例。我使用Anaconda来管理我的python安装。我安装了GDAL，conda update gdal生成：

Fetching package metadata: ..
# All requested packages already installed.
# packages in environment at /Users/<me>/anaconda:
#
gdal                      1.10.1               np18py27_2

正在运行which ogr2ogr

ogr2ogr is /Users/<me>/anaconda/bin/ogr2ogr

下载并解压缩Natural Earth子单元文件后，我运行以下命令：

ogr2ogr -f GeoJSON -clipdst -105.53 39.9 -104.93 40.27 places.json ne_10m_admin_0_map_subunits.shp

不幸的是，我收到了这个错误：

ERROR 6: GEOS support not enabled.

我已查找此错误，建议here是启用geos（我使用conda install geos安装的软件包）。所以我有geos版本3.3.3。但是，我很遗憾如何启用ogr2ogr来访问geos库。

非常感谢任何帮助。我没有开始使用ogr2ogr，但我没有看到很多替代方案。

更新时间：2015-02-26

重新安装anaconda后，我安装了gdal和geos。这是目前的情况：

$ ogr2ogr -f GeoJSON -clipdst -105.53 39.9 -104.93 40.27 places.json ne_10m_admin_0_map_subunits.shp

FAILURE:
Unable to open datasource `ne_10m_admin_0_map_subunits.shp' with the following drivers.
  -> ESRI Shapefile
  -> MapInfo File
  -> UK .NTF
  -> SDTS
  -> TIGER
  -> S57
  -> DGN
  -> VRT
  -> REC
  -> Memory
  -> BNA
  -> CSV
  -> GML
  -> GPX
  -> KML
  -> GeoJSON
  -> GMT
  -> GPKG
  -> SQLite
  -> WAsP
  -> PostgreSQL
  -> PCIDSK
  -> OpenFileGDB
  -> XPlane
  -> AVCBin
  -> AVCE00
  -> DXF
  -> Geoconcept
  -> GeoRSS
  -> GPSTrackMaker
  -> VFK
  -> PGDump
  -> OSM
  -> GPSBabel
  -> SUA
  -> OpenAir
  -> PDS
  -> WFS
  -> HTF
  -> AeronavFAA
  -> EDIGEO
  -> GFT
  -> GME
  -> SVG
  -> CouchDB
  -> Idrisi
  -> ARCGEN
  -> SEGUKOOA
  -> SEGY
  -> ODS
  -> XLSX
  -> ElasticSearch
  -> PDF
  -> CartoDB
  -> SXF

Answer 1

我遇到了同样的问题并找到了解决方案。解决方法是必须在系统上安装qt库。安装qt rpm软件包后，请安装geos并使用--with-geos选项重新安装gdal。

如何将shapefile剪辑到边界框？

更新时间：2015-02-26

1 个答案: