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

強力なデータ分析パッケージ内で美しいグラフを作成するためのggplot2の利点を享受して、Rでコロプレスマッピング手法を試してきました。 マッピングに関しては、QGISから日常的に得られる結果に匹敵する結果を出すことができませんでした。シェフィールドのフローのデータ（ここでは完全に複製可能なデータと.qgsファイル）を使用して、QGISは次のものを簡単に作成しました。 Rを使用して（ここで説明するコードを使用して）作成できる最高のものは、次の2つのマップでした。1つ目はGISToolsを使用し、2つ目はggplot2を使用しました。 障害のある行の問題（おそらくfortify（）コマンドの問題、またはここで説明したreadOGR（）を使用してシェープファイルを読み取れない問題）を解決できれば、ggplot2オプションははるかに魅力的に思えます。 質問は2倍です。ggplot2オプションはRの最適なコロプレスマッピングソリューションですか。もしそうなら、どのようにして白い線の問題を解決できますか。 私がやったことを見つけるための複製可能なコードはここにあります。 編集-choropleth（）出力オプションをより魅力的にしたので： このオプションのさは、高解像度でエクスポートし、凡例を削除することで軽減できます（他のR凡例が利用可能です）。それでも線の太さを変更することはできませんでしたが、改善されています！赤い線は、5000人以上を雇用するゾーンへのフローを表します。

18 r  thematic-map 

SpatialPolygonsDataFrame プロットをより効率的にするために、Rでaを一般化したいと思います。私が持っているポリゴンは、読み込みが遅く、ドキュメントに統合するのが難しい巨大なpdfプロットになります。 パッケージのgSimplify関数を使用しようとしましたrgeosが、残念ながらポリゴンのトポロジは保存されません。 使用できる代替手段はありますか？

18 r  topology  generalization 

RI では、パッケージとカスタムマップを使用して、アニメーション化されたコロプレスマップをrMaps作成します。ここで説明した Diego Valleの例に従って、データマップにカスタムマップ情報を使用しようとしました。 ただし、マップは表示されません。ただし、適切な塗りつぶしキーを持つ凡例は、年のアニメーションバーと共に表示されます。topoJsonファイルは問題ないようです。 library(rCharts);library(rMaps);library(rgdal) map = readOGR("https://raw.githubusercontent.com/gupierre/GHA_vol/master/Data/GHA_adm/GHA_adm1.topo.json", layer="admin_1",stringsAsFactors=F) #Data: SDLOG<-read.table("https://raw.githubusercontent.com/gupierre/GHA_vol/master/Data/sdlog.csv", header=T,stringsAsFactors=F, sep =",") head(SDLOG) year name sdlog type >1 2004 Greater Accra 0.20030687 urban >2 2005 Greater Accra 0.15322712 urban >3 2006 Greater Accra 0.10782677 urban #Map: library(rCharts);library(rMaps) d1 <- ichoropleth(sdlog ~ name, data = SDLOG[SDLOG$type=="rural",], ncuts = …

17 r  topojson  choropleth 

私が持っているSpatialPointsDataFrameいくつかの追加データで。ポリゴン内のこれらのポイントを抽出すると同時に、SPDFオブジェクトとそれに対応するデータを保持したいと思います。 これまでのところ、私はほとんど運がなく、共通のIDを照合してマージすることに頼りましたが、これは、個々のIDSでデータをグリッド化したためにのみ機能します。 ここに簡単な例があります、私は赤い四角の中のポイントを探しています。 library(sp) set.seed(357) pts <- data.frame(x = rnorm(100), y = rnorm(100), var1 = runif(100), var2 = sample(letters, 100, replace = TRUE)) coordinates(pts) <- ~ x + y class(pts) plot(pts) axis(1); axis(2) ply <- matrix(c(-1,-1, 1,-1, 1,1, -1,1, -1,-1), ncol = 2, byrow = TRUE) ply <- SpatialPolygons(list(Polygons(list(Polygon(ply)), ID …

17 r  overlay 

私は空間分析に慣れていないので、以下に概説するように、私が試みているプロジェクトの一般的な方向性に感謝します（最初から始めています）。 目標：ほとんどの世帯を中央ビジネス地区（CBD）に接続するために、故郷に2000フィートの歩道を設置するのに最適な場所を見つけること。既存の構造（世帯）、道路、歩道（既に設置済み）を示すシェープファイルがあります。 ここに私の提案された解決策/思考プロセスがあります： インプレースの歩道ネットワークを、重み（距離）で接続されたノードのデータベースに変換します。すべての交差点をクリックして、QGIS（または他のプログラム）でこれを直接行う方法はありますか？ pgRoutingなどのルーティング機能を使用して、中央ビジネス地区から1.2マイル以内にある世帯数（緯度経度のポイントまたはポリゴンなど）を計算します。これが基本ケースの「世帯アクセス」値になります。 道路レイヤーをガイドとして使用して、さらに2000フィート（たとえば10フィートのセグメント）の歩道をランダムに歩道レイヤーに配置します。これは、任意の数の新しい歩道を構築することと同等です。 （1）のように新しい歩行者ネットワークを使用してノードと重みを再計算し、次に（2）のようにCBDから1.2マイル以内にある世帯の数を再計算します。追加の歩道で増加するはずです。追加の歩道の場所と関連する「世帯アクセス」値をファイル（スプレッドシートなど）に保存します。 モンテカルロシミュレーションと同様に、ステップ（3）と（4）を10000回繰り返します。データポイントの10000セットを使用して、CBDから1.2マイル以内の世帯数を最大化する歩道配置場所を選択します。 この思考プロセスは現実的に聞こえますか？誰か提案はありますか？ -QGISとRのいくつかの組み合わせを使用してこれを達成したいのですが、目標を達成するためにPostGISやPython（または他の何か）を学ぶことは自由です。

17 qgis  distance  pgrouting  r  simulation 

この投稿はコミュニティWikiです。この投稿を改善するには、既存の回答を編集してください。現在、新しい回答は受け付けていません。 私は地理統計（および3 tのその他の多くの単語:)に非常に興味があります。 私はあまりプログラム指向ではありませんが、このテーマでいくつかのリストに従っています。 1つはstackexchangeのRフィルターです（質問のダイジェストメールを時々受け取ります）。 統計全般に関連する非常に複雑な言語があることがわかりました。 私は、この言語の一部が何を意味するのか、それを日常のGISにどのように適用するのかを、素人の視点から学ぶことに興味を持っています。 地理統計の知識/理解のために現在のソースをリストしてください。 ここに役立つ情報の例があります... オンラインのライブチャンス をトレーニングするための旅行 編集：リクエストごと。私はオンライントレーニングを楽しんでいますが、教室やステップバイステップブック（またはpdf）でより多くを学んでいます。

17 r  training 

のレイヤーがありcoord. ref. : NAます。 私はどのように割り当てることができますcoord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0するにはcoord. ref. : NA？ class : RasterLayer dimensions : 127, 146, 18542 (nrow, ncol, ncell) resolution : 0.01, 0.01 (x, y) extent : -122.39, -120.93, 36.72, 37.99 (xmin, xmax, ymin, ymax) coord. ref. : NA data source : C:\2012Tif\2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc.tif …

16 raster  coordinate-system  r 

関数で構築した空間バッファー（30000バッファー）のリストがありますlapply。 buff.pts <- lapply(1:nrow(pts.prj), FUN=function(l){ buff <- gBuffer(pts.prj[l,], width=1000) ## 1km return(buff) })) > head(buff.pts) [[1]] class : SpatialPolygons features : 1 extent : 307941.8, 311941.8, 4994518, 4998518 (xmin, xmax, ymin, ymax) coord. ref. : +proj=tmerc +lat_0=0 +lon_0=-73.5 +k=0.9999 +x_0=304800 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0 [[2]] class : SpatialPolygons …

16 polygon  r  sp  rgeos 

一連の単一境界ポリゴンを使用して土壌ポリゴンのRで真のGISクリップを実行したいのですが、適切に実行するR関数が見つかりません。clipESRIのArcMapの関数と同じように機能するはずです。パッケージover内のメソッドを試しましたが、sppolys over polysでは機能しないようです。 1つの提案は、次のコードを使用して、gIntersectionin rgeosパッケージをクリップとして使用することです。 #------------------------------------ library(rgeos) library(maptools) #Read layers as SpatialPolygonsDataFrame (both the same Albers projection) Soils_poly = readShapePoly("Soils_polygons") #Note - Has 400 polygons clipper_poly = readShapePoly("clipper_polygon") #Note - Has 1 polygon #Try gintersection as clip Clipped_polys = gIntersection(Clipper_Tile_poly, Soils_poly) #----------------------------------- これには、実行に5分かかり（速度が遅すぎる）、次のエラーが発生します。 RGEOSBinTopoFunc（spgeom1、spgeom2、byid、id、drop_not_poly、 "rgeos_intersection"）のエラー：TopologyException：-721459.77681285271 2009506.5980877089で発信dirEdgeが見つかりません また、重複をチェックするためにこのコードを試しました。 gIntersects(Clipper_Tile_poly, Soils_poly) そして結果は真でした。clipESRI …

16 r  polygon  layers  vector  clip 

私はQGISを初めて使用し、QGISとRのテーブルを使用できるデータベースを作成する簡単な方法を探しています。統計分析に使用されるExcelスプレッドシートに多くのデータが保存されています。 Rで タスクの詳細：最初のタスクは、テーブルをシェープファイルに結合してデータを可視化することです。次に、QGIS内でシェープファイルを使用して空間変換を行い、新しく作成したデータを統計分析で使用します。したがって、このデータをデータベースに結合する必要があります。 最後になりましたが、データベースは他のコンテキストで使用されるため、QGISおよびRの将来の分析のためにデータが適切な方法で保存されることを保証したくありません。 MS Accessがそうする良い方法であると私に言われました。誰かアドバイスをいただけますか？おそらく文学やそのようなものについても。

16 qgis  r  spatial-database  database  ms-access 

現在、次のスクリプトを使用して、テーブルからいくつかの属性データを多数の個々のシェープファイルに追加しています。 library(rgdal) specieslist <- read.csv("SpeciesList1.txt", header=F) attdata <- read.table("TestAtt.csv", sep = ",", header=T) for (n in 1:dim(specieslist)[1]) { speciesname <- specieslist[n,1] shp <- readOGR("Mesoamerica_modified_polygons", speciesname) shp$ENGL_NAME<-attdata[n,2] writeOGR(shp, "PolygonsV2", speciesname, driver="ESRI Shapefile") } 最後に次の警告が表示されます。 1: In writeOGR(shp, "PolygonsV2", speciesname, driver = "ESRI Shapefile") : Field names abbreviated for ESRI Shapefile driver …

16 shapefile  gdal  r  ogr  attribute-table 

2つのレイヤーがあります。多くのタイルを含むポリゴンシェイプレイヤーと、カラーマップの多くのカテゴリを含むCORINE 2006土地被覆を含むラスターレイヤー。シェープレイヤーの各ポリゴンについて、ラスターレイヤーの各土地被覆カテゴリの合計を取得します。 たとえば、IDが '2'のポリゴンがあり、このポリゴンに対して次のような属性が必要です（パーセントまたは平方メートル）。 耕地：15％ フォレスト：11％ 通り：2％（...など） 私は草、qgis（機能なし）、saga（すべてを合計値に合計する）r（合計）でそれをしようとしましたが、まだ解決策が見つかりませんでした。ほとんどのプラグイン（qgisのゾーン統計）は、0-1ラスターレイヤーのみをサポートします。v.rast.statsも役に立たなかった。Iamは、優れたスマートソリューションにオープンです。たぶん私は間違ったアプローチを使用したか、間違いを犯した。 Arcgisでは、このタスクは非常に簡単です。覚えていれば正しいですが、私はあなたの日常のLinuxユーザーにとって良い解決策がまだありません。 私はDebian Linuxシステムを実行しているため、このOS用のプログラムしか使用できません。 編集：この質問にはまだ多くのビューと訪問者がいるので、QGISプラグインを作成しました。これは、ラスターレイヤーの土地被覆を計算することもできます。ポリゴンオーバーレイはまだコーディングしていませんが、間違いなく計画されています。ここでプラグインを見つけて、最初にScipyライブラリをインストールします。

16 qgis  raster  polygon  r  statistics 

コードをspパッケージから新しいsfパッケージに移行しています。以前のコードでは、SpatialDataFrame（censimentoMap）とSpatialPointDataFrame（indirizzi.sp）のポリゴンがあり、以下の手順で各ポイントのポリゴンセルID（Cell110）を取得しました。 points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"]) 実際に、2つのsfオブジェクトを作成しました。 shape_sf <- st_read(dsn = shape_dsn) shape_sf <- st_transform(x=shape_sf, crs=crs_string) そして indirizzi_sf = st_as_sf(df, coords = c("lng", "lat"), crs = crs_string) そして、上記の命令に相当するsfを探しています... ids<-sapply(st_intersects(x=indirizzi_sf,y=shshape_sfpeCrif), function(z) if (length(z)==0) NA_integer_ else z[1]) cell_ids <- shape_sf[ids,"Cell110"]

16 r  overlay  sp  sf 

Bing、Yahoo、OpenStreetmapを使用してバッチジオコーディングを実行するRのパッケージまたはRのソースコードを見つける必要があります。何か提案はありますか？

15 r  geocoding  batch  bing-maps 

今日、オーストラリアの道路形状ファイルを探していました。最終的にGeofabrikに行きました。Open Street Mapには.osm形式のデータがあります。しかし、彼らは.shpオーストラリアのためにそれを持っていません。 可能であれば、Rを使用するようにosmファイルを変換しますshapefiles。検索しましたが、まだ解決策が見つかりませんでした。これに何か方法はありますか？あるいは、osmファイルをRに読み込み、データをデータフレームに変換して、図を使用して図形を描画することは可能ggplot2ですか？ との2つのシェープファイルを使用して、NZの次の図を作成maptoolsしましたggplot2。理想的には、オーストラリアの同様の地図を作成したいと思います。

15 shapefile  openstreetmap  r  ggplot2