このシェープファイルを適切な境界ボックスにクリップする方法について何か提案はありますか? 特に、[(-105.53 , 39.9) (-104.93,40.27)] という境界のみを表示したいと思います。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.
このエラーを調べましたが、ここでの提案は有効にすることです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