地理情報システム

私は巨大な.kmlファイル（最大10 Gb）で作業しており、それらをRに読み込む効率的な方法が必要です。これまでは、QGISを介してシェープファイルに変換し、その後readShapePolyとreadOGR（後者ちなみに、前者よりも1000倍高速です）。面倒で時間がかかるため、QGIS中間段階を切り取りたいと考えています。 .kmlファイルを直接読み込む方法は？ 私は、これはまた、readOGRで行うことができます参照します。残念ながら、実際の例を実装する方法はわかりません（.kmlファイルを長時間準備した後xx <- readOGR(paste(td, "cities.kml", sep="/"), "cities")）。ここの「都市」は空間オブジェクトの名前のようです。 Roger Bivandは、「OGRのKMLドライバーはファイルにアクセスするためにこの名前を必要とするため、この名前の発見方法は明らかではありません。1つの可能性は次のとおりです。 system(paste("ogrinfo", paste(td, "cities.kml", sep="/")), intern=TRUE) 」 しかし、これは私にとってもうまくいきません。試用するテスト.kmlファイルを次に示します。作業ディレクトリにあると、readOGR("x.kml", "id")このエラーメッセージが生成されます。 Error in ogrInfo(dsn = dsn, layer = layer, encoding = encoding, use_iconv = use_iconv) : Cannot open layer . そしてsystem(paste("ogrinfo", "x.kml"), intern=TRUE)生成： [1] "Had to open data source read-only." "INFO: Open …

42 kml  r 

QGIS 2.16の[ベクター]メニューにある通常の処理ツール（ベクターオーバーレイ、解析など）がすべてありません。それらは処理ツールボックスに表示されますが、ベクトルメニューには表示されません。誰もこの問題に遭遇しましたか？

41 qgis  qgis-processing 

Rの空間データを分析するには、いくつかの素晴らしい可能性があります。現在のプロジェクトのコンテキストでは、Rをより頻繁に使用して分析したいと思います。 これまで、データをプロットおよび探索するための実用的なツールを多数提供するggplot2パッケージを使用してマップをプロットしました。それでも、私のマップは、ArcGISでプロットするために使用したマップほど美しくありません。 それで、出版の目的のために、少し実験して地図をもっと見栄え良くするために使用できる良いチュートリアル、本、実用的なヒント、パッケージなどがあるのだろうか？

41 r  cartography 

この質問の答えを考え出すのに少し時間を費やしました。Google検索からすぐに明らかではないので、ここに回答を投稿すると便利だと思います。連続していないポリゴンに関する追加の質問もあります。 簡単な答え：次のコマンドを使用します。 centroids <- getSpPPolygonsLabptSlots(polys) （これは、R、spの包括的な空間パッケージのSpatialPolygonsDataFrame Rデータクラスのクラスの説明で見つかりました） これはまったく同じことをするようです cents <- SpatialPointsDataFrame(coords=cents, data=sids@data, proj4string=CRS("+proj=longlat +ellps=clrk66")) 次のコードでは、どのRインストールでも複製可能であるはずです（試してみてください！） #Rcentroids install.packages("GISTools") library(GISTools) sids <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1], proj4string=CRS("+proj=longlat +ellps=clrk66")) class(sids) plot(sids) writeSpatialShape(sids, "sids") cents <- coordinates(sids) cents <- SpatialPointsDataFrame(coords=cents, data=sids@data, proj4string=CRS("+proj=longlat +ellps=clrk66")) points(cents, col = "Blue") writeSpatialShape(cents, "cents") centroids <- getSpPPolygonsLabptSlots(sids) points(centroids, pch = 3, …

41 r  centroids 

これは私にはそのような単純な質問のように思えます（おそらくそうです）が、答えを与える例を見つけることができないようです。PostGISを使用して、ポリゴンの外側にあるポイントを選択したいだけです。最終的に、これはST_Intersectsの逆であり、私が見る限りです。 例：タックスロットレイヤーと住所ポイントレイヤーがあります。ST_Intersectsを使用する必要があると思いますが、逆選択を行うようにするにはどうすればよいですか？次のコードの前にNOTステートメントを追加することも考えられましたが、うまくいきませんでした。 CREATE table t_intersect AS SELECT hp.gid, hp.st_address, hp.city, hp.st_num, hp.the_geom FROM public.parcel as par, public.housepoints as hp WHERE ST_Intersects(hp.the_geom,par.the_geom);

41 postgis  geoprocessing  sql  intersection 

同じ地理的領域内に2つのポイントパターン分布がある場合、それら2つの分布を視覚的および定量的に比較するにはどうすればよいですか？ また、小さな領域内に多くのポイントがあるため、ピンマップを表示するだけでは情報が得られないと想定します。

41 point  spatial-statistics  visualisation  gis-principle 

私が取り組んでいるプロジェクトでは、方向性を持つ1,000ゾーンを超えるネットワークフローを表す必要があります。 情報をすばらしく非圧倒的な方法で表示するための提案はありますか？ 重い都市部でネットワークフローを表示する例はありますか？ 色の提案と厚さのスキーム？ ネットワークフローを表すときに「やるべきこと」はありませんか？

41 network  transportation  flow-map 

QGIS googleレイヤープラグインを使用して、区画をシェープファイルとしてデジタル化します。後でPostGISにインポートしたいと思います。グーグルマップは特別な「グーグルメルカトル」投影法を使用していることを知っています。この投影法の正しいコードは何ですか：900913または3857？または、これら2つのコードは同一の投影方程式を記述していますか？

41 google-maps  coordinate-system  epsg 

実際、SinnottがHaversine式を公開したとき、計算精度は制限されていました。現在、JavaScript（および最新のコンピューターと言語のほとんど）はIEEE 754 64ビット浮動小数点数を使用しており、15桁の精度を実現しています。この精度により、余弦公式の単純な 球面法則（cos c = cos a cos b + sin a sin b cos C）により、約1メートルの距離まで良好な条件の結果が得られます。 これを考慮すると、ほとんどの場合、単純なコサインの法則、またはHaversineよりもより正確な楕円形のVincenty式を使用する価値があります。（球体モデルの精度の制限に関する以下の注意事項に留意してください）。 ソース：http : //www.movable-type.co.uk/scripts/latlong.html 余弦の法則がより望ましい理由は何ですか？ 注：引用されたテキストは、によって更新されたその作者に述べたように下回ります。

41 distance  algorithm  spherical-geometry  haversine 

路面を表すポリゴンの平均幅を調べることに興味があります。また、道路の中心線をベクトルとして持っています（正確には中心にない場合があります）。この例では、道路の中心線は赤で、ポリゴンは青です。 私が考えていたブルートフォースアプローチの1つは、ラインを少しずつバッファリングし、バッファをフィッシュネットグリッドと交差させ、道路ポリゴンとフィッシュネットグリッドを交差させ、両方の交差測定の交差面積を計算し、これをエラーはわずかです。しかし、これは大雑把なアプローチであり、よりエレガントなソリューションがあるかどうか疑問に思っています。さらに、これにより、大きな道路と小さな道路の幅が隠されます。 ArcGIS 10、PostGIS 2.0、またはQGISソフトウェアを使用するソリューションに興味があります。この質問を見て、ArcGIS10用のDan Pattersonのツールをダウンロードしましたが、必要なものを計算できませんでした。 ArcGIS 10で最小境界ジオメトリツールを発見したところ、次の緑のポリゴンを作成できます。 これは、グリッドに沿った道路に適したソリューションのように見えますが、それ以外では機能しません。そのため、他の提案にも興味があります。

41 qgis  arcgis-desktop  arcgis-10.0  postgis  polygon 

ArcMapライセンスがなくてもPythonを使用してシェープファイルの内容を見ることができるかどうか疑問に思いました。状況は、ESRIソフトウェアだけでなく、多くの異なるアプリケーションからシェープファイルを作成できることです。空間参照、フィーチャタイプ、属性名と定義、およびシェープファイルのフィールドの内容をチェックし、それらを許容値のセットと比較するPythonスクリプトを作成したいと思います。組織にESRIライセンスがない場合でも、このスクリプトが機能するようにします。このようなことを行うには、ArcPyを使用する必要がありますか、それともArcPyを使用せずにシェープファイルを掘り下げることができますか？

40 python  shapefile 

いくつかのgdal APIを使用gdal_translateしてPythonコードから呼び出すことは可能ですか？ファイルシステムからgdal_translate.exeを実行するだけではなく、コード内でなんらかの方法で呼び出すので、gdal_translate実行可能ファイルが存在する正確なディレクトリを知る必要はありません。

40 python  gdal 

複数のポリゴンがあり、すべて英国の沖合にあります。UTMゾーンで分割せずに、これらのポリゴンの面積を平方キロメートルで計算しようとしています。 これを正確に行うことは可能ですか？私はワールドメルカトルを使用していますが、このスケールでは大きな歪みがあります。たとえば、ポリゴンが東海岸に限定されている場合、UTM31Nを使用します。しかし、それらは英国海域の至る所にあるので、私は何を使うべきか分かりません。 この領域をカバーする座標系はありますか？ブリティッシュナショナルグリッドは、使用できるほど海まで十分に拡張されていません。

40 coordinate-system  utm  united-kingdom 

誰かからtiffファイルを受け取りましたが、ジオリファレンスされているかどうかわかりません。これを確認する簡単な方法はありますか？

40 geotiff-tiff  georeferencing 

フィールドに郵便番号データを含むCSVファイルがありますが、それは非空間的です。たとえば、EastingとNorthingはありません。 フィールドに郵便番号といくつかの追加データを含むシェープファイルがあります。 属性（郵便番号フィールド）に基づいてCSVをシェープファイルに結合したい。ArcMapでこれを実行できることは知っていますが、QGISで同じことをどのように実現できますか？

40 qgis  csv  attribute-joins