タグ付けされた質問 「shapefile」

ここから地図があります。任意のlon / latペアをマップ上のピクセルに簡単に変換できるようにしたい（逆変換も可能）。マップには.tfwファイルと投影情報が付属しています。 0.02222222222222 0.00000000000000 0.00000000000000 -0.02222222222222 -180.00000000000000 90.00000000000000 と投影情報： Projection: Plate Carree aka Geographic or "LatLong" Earth ellipsoid: Sphere, radius 6370997 m Datum: WGS84 Extent: 180 West to 180 East, 90 North to 90 South Size: 16,200 height samples wide x 8,100 high Resolution: 2.47 km/pixel 私は地図作成関連のことをまったく知らないので、私が理解している限り、まずWGS84（経度/緯度のペア）から地理的投影法に変換する必要があります（同じではありませんか？）。実際には同じように見えますが、上の投影情報の球の半径は6370997で、プレートカリー投影で見つけたspatialreference.comページとは異なります。とにかく、私はDotSpatial.Projectionsライブラリが次のコードでこれを行うことができることを発見しました： var kievCoordinates …

10 shapefile  coordinate-system  c#  latitude-longitude  affine-transformation 

QGISプロジェクトを再投影したかったのですが、そのようなオプションはありません。回答を検索するときogr2ogr、各ベクターレイヤー（この例ではシェープファイル）を再投影することにより、このタスクの提案が表示されます。ただし、期待した結果が得られません。 C:\temp>ogrinfo -al -so misc.shp INFO: Open of `misc.shp' using driver `ESRI Shapefile' successful. Layer name: misc Geometry: Line String Feature Count: 10 Extent: (21.267388, 42.015857) - (21.270225, 42.017470) Layer SRS WKT: PROJCS["MGI_Balkans_zone_7_deprecated", GEOGCS["GCS_MGI", DATUM["Militar_Geographische_Institute", SPHEROID["Bessel_1841",6377397.155,299.1528128]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",21], PARAMETER["scale_factor",0.9999], PARAMETER["false_easting",7500000], PARAMETER["false_northing",0], UNIT["Meter",1]] id: Integer (10.0) C:\temp>ogr2ogr -t_srs …

10 shapefile  coordinate-system  ogr2ogr 

シェープファイルを「反転」する可能性を探しています。私がしたいことは、大陸と島のあるシェープファイルを水のポリゴンを持つシェープファイルに変換することです。Natural Earthの既存の海洋シェープファイルは不正確です。 私はアクティブなOSM寄稿者として、OSMの既存のシェープファイルを使用することでうまくいくと思いました。これらははるかに正確ですが、GISの観点からは非常にバグがあります。しかし、mapserverはそれを処理できるので、これらは問題ないはずだと思いました。これらのシェープファイルは、http：//metro.teczno.com/#coastlineで見つけることができます が、非常に大きいものもあります。 ワールドバウンド（-180〜180、-90〜90）を持つ1つのポリゴンでシェープファイルを作成しました。次に、陸域が他のシェープファイルにある世界境界シェープファイルから領域を切り取ろうとしました。QGIS、ベクターツール->差分関数で試してみました。100％のCPU使用率で12時間後、私はそれを停止しました。常に100％の進行状況を表示しました。RAMは問題ではありませんでした。 「世界のポリゴン」、できればLinuxコマンドラインから陸地を切り取る別の可能性を探しています。

10 shapefile  polygon  clip 

非常に大きなシェープファイル（RIのすべての道路）があります。Pythonで複数のワーカースレッドを使用してファイルを処理できるように、ファイルをチャンクに分割する必要があります。これまでのところ、ベクターデータにオーバーレイされるグリッドをQGISで作成しました。次に、グリッドとシェープファイルを使用して交差アルゴリズムを実行し、ベクトルが分割されて、'id'それらが含まれるボックスにが割り当てられるようにしました。 今、グリッドの特定のボックス内のすべてのベクトルはすべて同じな'id'ので、それらを抽出して独自のシェープファイルに自動的に保存できますか？ですべてのベクトルを選択し、'id'それらを独自のレイヤーにコピーしてからそのレイヤーを保存することは承知していますが、ボックス/ IDが多すぎるため、それぞれをオフに保存するには永久に時間がかかります。QGISを使用して、ID /ボックスごとに新しいシェープファイルを自動的に作成するにはどうすればよいですか？ 基本的に、シェープファイル属性の特定のフィールドの個別の値ごとに新しいシェープファイルを作成します。 元のシェープファイル グリッド交差後のシェープファイル 同じbbox / idを共有するベクトルのセット

10 qgis  shapefile  vector 

私は通常、静的な地形図を作成します。これはほとんどが印刷されます。OSMとラスターファイル（主に異なるDEM）をよく使用しています。私がマッピングしているエリアが20.000km²を超えることはめったにありません。QGIS、SAGA、GRASSなどの排他的なオープンソースソフトウェアを使用しています。 シェープファイルの問題は知っていますが、処理は簡単です。そして、私はpostgisやSpatiaLiteにはあまり詳しくありません。私はpostgisを試しましたが、それを処理するのに多くの問題がありましたが、SpatiaLiteの方が簡単に思えます。 どう思いますか：シェープファイルを使い続けるか、SpatiaLiteやpostgisなどのデータベースに切り替えますか？私の場合の長所と短所は何ですか？

10 postgis  shapefile  spatialite 

特定のポリゴンにボロノイ図を作成しようとしていますが、エッジを処理するときにいくつかの問題があります。ここに私がそれをする方法があります： ボロノイのポイントとダイアグラムのポリゴンの準備 ボルノイを作成 ポリゴンを使用してエッジをカットします ユニオンを使用して、カットされたボロノイとポリゴンを結合します（これが最善の方法かどうかはわかりません）そして、ここでエッジの問題が発生します。1）元のボロノイマップからフレームを削除（青でマーク）2）エッジをエリアポリゴンのエッジ（赤でマーク）まで延長し、エッジ側のサービスエリアが現実的になるようにします。 これどうやってするの？

9 qgis  shapefile  polygon  editing  voronoi-thiessen 

シェープファイルをQGISにインポートし、それに対して空間SQLクエリを実行する方法はありますか？ 私は明らかにogr2ogr -sqlこれを行うために使用できますが、それを見るには毎回QGISに出力をインポートする必要があり、これは苦痛です。 シェープファイルをQGISにインポートし、それに対してSQLを実行して、結果を即座に確認できればすばらしいと思います。

9 qgis  shapefile 

非常に大きなESRIシェープファイル（> 8000ポリゴン、.shpファイルは> 32MB）があります。これをWindows 7のRに読み込みます。 私はrgdalとreadOGR（）に精通しています。ただし、このサイズのシェープファイルでは、シェープファイル全体をメモリに読み込むのは非常に時間がかかります。さらに、シェープファイルに含まれる8,000のいくつかのクラスター化されたポリゴン（<100）のみに関心があります。これを実行する必要があるプロセスは非常に反復的です。つまり、シェープファイルを何度もロードする必要があるため、現在は非常に低速です。 したがって、プロセスを加速するために、シェープファイルの一部のみを読み取ることができるかどうか疑問に思っています。このような大量のデータを読み込んで、そのほとんどを破棄するのは実際的ではないようです。

9 shapefile  r  rgdal 

ggmapを使用して、いくつかのロケーションポイントを含むマップに、シェープファイルの自治体の境界線（ポリゴン）を含めたいと思います。このスクリプトは、ポリゴンのプロット以外のすべてを実行します。 library(rgdal) library(ggmap) # Get shapefile with Drammen municipality borders tmpzip<-tempfile() tmpdir<-tempfile() dir.create(tmpdir) download.file("http://www.kartverket.no/Documents/Kart/N50-N5000%20Kartdata/33_N5000_shape.zip",tmpzip) unzip(tmpzip, exdir=tmpdir) kommune <- readOGR(dsn=tmpdir, layer="NO_AdminOmrader_pol") kommune<-kommune[kommune$NAVN=="Drammen",] kommune<-spTransform(kommune, CRS("+init=epsg:4326")) # Get location point data subscr<-data.frame(lon=c(10.1237,10.2161,10.2993),lat=c(59.7567,59.7527,59.6863), pop=c(58,12,150)) coordinates(subscr)<-~lon+lat proj4string(subscr)<-CRS("+init=epsg:4326") lon <- c(10.0937,10.3293) lat <- c(59.7916,59.6563) map <- get_map(location = c(lon[1], lat[2], lon[2], lat[1]), maptype = "roadmap", source …

9 shapefile  r  ggmap 

シェープファイルをディゾルブすると同時に属性フィールドを保持する方法はありますか？ これは私がフィールドを保持したい方法です： シェープファイルに保存したい2つのフィールドがあります。最初の列はlongタイプのNet_IDで、2番目の列は文字列タイプのGeologyです。 NET_IDフィールドを使用してシェープファイルを解体したいと思います。また、下図のように地質学分野を守りたいと思います。

9 arcgis-desktop  shapefile  arcgis-10.2  dissolve 

休業。この質問には、より焦点を当てる必要があります。現在、回答を受け付けていません。 この質問を改善してみませんか？質問を更新して、この投稿を編集するだけで1つの問題に焦点を当てます。 2年前休業。 OGR / GDALを使用してPythonでスクリプトを作成しています。 シェープファイルのセットとGeoTiffラスターファイルのセットがあります。 シェープファイルがラスター領域と交差しない場合、スクリプトに無視させたいのですが。 シェープファイルは長方形ではないため、layer.GetExtent（）から返されたxmin / xmax、ymin / ymaxの値を単純に比較することはできません。全体的な形状を表す実際のポリゴンが必要です。次に、そのポリゴンがラスタースクエアと交差するかどうかを判断する方法が必要です。 どういうわけか、シェープファイル内のすべてのポリゴンを1つのフィーチャにマージし、そのフィーチャのジオメトリを読み取って、その情報をラスター範囲と比較できると考えていました。ただし、具体的にどのように実行するかはわかりません。 シェープファイルからボーダーポリゴン情報を抽出する方法は？ そのポリゴンが特定の正方形の領域と交差するかどうかを判断するにはどうすればよいですか？

9 python  raster  shapefile  gdal  ogr 

私はデンマークの自治体のシェープファイルの重心をいくつか見つけて、その後それらの間の運転時間を見つけようとしています。私が使っRさんreadShapePolyからの機能maptoolsと組み合わせるgCentroidから機能rgeos、およびすべての作品を。しかし、私は次のような空間点を取得します SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) arguments: NA これは明らかにGoogleで移動時間を取得するために使用できるものではありません。これらの数値を経度と緯度に変換する方法を探していますが、方法がわかりません。 ライブラリreadOGRからを使用してデータを読み取ると、rgdal同じ座標が取得されますが、これは投影法であると想定していることについて次のことを教えてくれます（ただし、座標は同じです） Slot "proj4string": CRS arguments: +proj=utm +zone=32 +ellps=intl +units=m +no_defs 再現可能な例：https : //github.com/sebastianbarfort/shapefilesの例のデータをここに配置しました： これは問題を再現するはずです： library(maptools) library(rgdal) library(rgeos) map = readShapePoly("~/Downloads/shapefiles-master/kommuner1983.shp") centroid = gCentroid(map) centroid

9 shapefile  r  coordinate-system 

私が働いている町のすべての建物と家を含む大きなシェープファイルがあります（約90,000個のフィーチャ）。建物/家屋のデータは町の測量技師によって保存され、悪い習慣とそのデータへのさまざまな測量士のアクセスにより、多くの建物/家屋が2回保存され、重複してマップに表示されています。 それらの一部は正確に複製され（一方が他方の上に表示されます）、もう一方は2つのオブジェクト間にスペースを空けて複製されます（一方のオブジェクトが他方の中にあるように-添付のスクリーンショットを参照）。 そのデータを整理して、町に正しい建物や家だけができるようにしたいので、私の質問は次のとおりです。 すべての複製されたフィーチャ（正確なものと他の内部にあるものの両方）を見つけるために実行できるGIS分析またはSQL式はありますか？私はArcGISとQGISの両方を持っているので、あなたの提案をすべて受け入れます。

9 qgis  arcgis-desktop  shapefile  sql  analysis 

私はマイクボストックのチュートリアル「マップを作ろう」に従っています。彼はogr2ogr、機能の属性に基づいてシェープファイルをフィルタリングする方法を説明しています。QGISでシェープファイルを開き、属性テーブルを見て、それらの機能を確認できます。コマンドラインからその情報を見つける方法はありますか？

9 shapefile  gdal  ogr2ogr 

私はGISが初めてです。Pythonツールを使用してコロプレスを作成する方法のチュートリアルを見つけました。フォローしたいのですが、自分のデータと地図を使います。ただし、チュートリアルではUSマップのSVGファイルを使用しており、CAの郵便番号マップを使用する必要があります。米国国勢調査のシェープファイルで私が見つけた唯一のマップ。kartograph.pyを使用してSVGに変換しようとしましたが、そうするとファイル内の郵便番号が失われます。ファイルをSVGに変換し、各パスを適切な郵便番号としてidしたままにするにはどうすればよいですか。 私は無料のツールのチュートリアルを使用して米国郡の主題図を作成する方法に従っています。これで、私が何をしようとしているのかがわかります。私はSVG形式のマップを持っていません。 私はjavascriptユーティリティへの参照を見てきました。私はJavaScriptを使用したことがありませんが、簡単なものを使用できると仮定します。

9 python  shapefile  svg  choropleth