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

次のようなcsv形式のデータセットがあります。 latitude, longitude, value -45, 45, 10 .... ArcGISまたはR + grassを使用して、これらのデータから郡レベルの平均を計算したいと思います。 ArcGIS USA Countiesレイヤーがありますが、空間結合を正常に完了することができませんでした。 これどうやってするの？

9 arcgis-desktop  data  geoprocessing  r 

ポリゴンの値に基づいてラスターをマスクしたい：つまり、ポリゴンでカバーされていないラスターのすべての値をNAに設定します。ただし、raster::mask()部分的に覆われている場合でも、セルをマスクする（つまり、NAに設定する）ようです。 これは質問に似ています：ポリゴンでラスターをトリミングする-部分的にポリゴンの外側にあるセルが欠落していますが、とは異なりcrop()、オプションmask()がありませんsnap="out"。 チェック：三角形の正しい点がセルに接触していますが、そのセルはマスクされています。 library(raster) library(sp) r <- raster(xmn=1, xmx=5, ymn=1, ymx=5, nrows=4, ncols=4) r[] <- 1:length(r) Sr1 = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2))) Sr2 = Polygon(cbind(c(4.1,4,2),c(2,3,2))) SpP = SpatialPolygons(list(Polygons(list(Sr1), "s1"), Polygons(list(Sr2), "s2")), 1:2) plot(mask(r, SpP)) plot(SpP, add=TRUE)

9 raster  r  masking 

で測地線バッファリングの理解、Esriのジオプロセシング開発チームは、ユークリッドと測地線バッファリングを区別します。「投影されたフィーチャクラスで実行されるユークリッドバッファリングは、誤解を招く技術的に不正確なバッファを生成する可能性があります。ただし、測地線バッファは、投影座標系によって導入される歪みの影響を受けないため、常に地理的に正確な結果を生成します。」 ポイントグローバルデータセットを使用する必要があり、座標は投影されません（+proj=longlat +ellps=WGS84 +datum=WGS84）。メートル単位で幅が指定されている場合、Rに測地線バッファーを作成する関数はありますか？パッケージgBufferから承知しておりrgeosます。この関数は、使用される空間オブジェクト（例）の単位でバッファーを作成するため、座標を投影して、目的のX kmのバッファーを作成できるようにする必要があります。投影してから、gBuffer実際にユークリッドバッファーを作成する手段を適用します。以下は、私の懸念を示すためのコードです。 require(rgeos) require(sp) require(plotKML) # Generate a random grid-points for a (almost) global bounding box b.box <- as(raster::extent(120, -120, -60, 60), "SpatialPolygons") proj4string(b.box) <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" set.seed(2017) pts <- sp::spsample(b.box, n=100, type="regular") plot(pts@coords) # Project to Mollweide to be able to apply buffer with …

9 coordinate-system  r  buffer  rgeos  geodesic 

多くのポリゴンを持つシェープファイルを使用しています。 「area_sqkm」という名前のフィールドをもう1つ追加して、シェープファイルの各ポリゴンの面積を計算するにはどうすればよいですか？

9 shapefile  r  area 

R maptoolsライブラリを使用して、シェープファイルを解析してポリゴンのリストを作成しています。関数readShapeSpatial shp <- readShapeSpatial("<path to my shapefile>") 私に与えSpatailPolygonDataFrameのオブジェクトを。私の例では、SpatialPolygonsDataFrameに次の列があります。 > names(shp) [1] "AREA" "PERIMETER" "COMAREA_" "COMAREA_ID" "AREA_NUMBE" [6] "COMMUNITY" "AREA_NUM_1" "SHAPE_AREA" "SHAPE_LEN 特定のポリゴンをで削除できることを知っていますrow.id。たとえば、 shp.dropI <- shp[-i, ] shp.subset <- shp[i %in% c(1,2,3),] 次に、特定のポリゴンを条件付きでドロップしたいとしAREA > 10ます。これをエレガントに実装するにはどうすればよいですか？私が今持っている唯一の方法は、すべての行を反復処理して対応するを見つけることrow.idです。

9 shapefile  r  sp 

CSVファイルをシェープファイルに結合したいと思いました。 結合を容易にするフィールドを特定したとき、たとえば、 map <- spChFIDs(map, as.character(map$ID)) 戻った spChFIDs（SP、x）のエラー：長さが異なります 誰かアドバイスできますか？

9 shapefile  r  csv  spatial-join  sp 

この質問での回答のおかげで、イギリスの一部、この場合はペンブロークシャーの選挙区のサブセットを作成して地図を描くことができました。結果のデータフレームは大きく、Ordnance Surveyデータが含まれているため、ここに投稿することは困難ですが、データフレームは次のようになります。 > str(bar) 'data.frame': 134609 obs. of 7 variables: $ long : num 214206 214203 214202 214198 214187 ... $ lat : num 207320 207333 207339 207347 207357 ... $ order: int 1 2 3 4 5 6 7 8 9 10 ... $ hole : logi FALSE FALSE …

9 shapefile  r 

推定年齢のポイントのセットが与えられた場合、予測される移動の方向を視覚化するための最良の方法について質問があります。 コミュニティを表す一連の地理座標があるとします。それぞれの点について、そのコミュニティの年齢の推定値があります。時間の経過に伴う一般的な移動の最も可能性の高い方向を示す矢印をマップに描画したいと思います。また、矢印の大きさによって、グラデーションの一貫性や、時間と空間の間に意味のある整合があることをどの程度確信しているかを示します。移動のプロセスの理論については質問していません。矢印の角度、位置、大きさをどのように決定するかについてです。 たとえば、以下の赤い点は緑の点より古いです。最初のボックスでは、移行が南東であったという賢明な結論のように見えます。2番目のボックスには明確なパターンがないため、矢印は小さくなっています。 これを行うための標準的で原則的な方法はありますか？私が考えた方法の1つは、最も古い点と2番目に古い点の間、次に2番目に古い点と3番目に古い点の間にベクトルを描画することでした。次に、ベクトルを合計します（ただし、開始位置を決定する方法は？）。それとも、すべてのポイントのペアの間のベクトルを取り、それらを相対年齢で重み付けする必要がありますか？または、コンターを計算してから、最高点から始まり下り坂に向かう経路を描く必要がありますか？ Rベースのソリューションが理想的です。

8 coordinate-system  r  vector 

として与えられた点から複数の線を作成したいと思いますsf objects。 私がいくつかのポイントを持っている場合 library(sf) pts <- st_multipoint(matrix(c(10, 10, 15, 20, 30, 30), nrow = 3, byrow = TRUE), dim = "XY") 私はst_castそれらのラインを作成するために使用しています lines <- st_cast(pts, "MULTILINESTRING") 私は常にsf object複数のセグメントを持つ1つを取得しますが、取得したいのは複数の行（この例では2つ）です。

8 r  sf 

直線と一連の点を考えると、sf各点が直線のどちら側にあるかを特定する方法を理解できません。 別の質問を基にした、再現可能な小さな例を以下に示します。 # Load Libraries ---------------------------------------------------------- library('sf') # Test data --------------------------------------------------------------- points.df <- data.frame( 'x' = c(-53.50000, -54.15489, -54.48560, -52.00000, -52.57810, -49.22097, -48.00000), 'y' = c(-38.54859, -41.00000, -38.80000, -38.49485, -38.00000, -40.50000, -37.74859) ) line.df <- data.frame( 'x' = c(-54.53557, -52.00000, -50.00000, -48.00000, -46.40190), 'y' = c(-39.00000, -38.60742, -38.08149, -38.82503, -37.00000) …

8 r  sf 

解像度と範囲が異なる2つのラスターがあります。 > res(Elevation) [1] 0.002083333 0.002083333 > res(Ann_precip) [1] 0.008333333 0.008333333 > extent(Elevation) class : Extent xmin : -15.07722 xmax : -7.641806 ymin : 7.193611 ymax : 12.67694 > extent(Ann_precip) class : Extent xmin : -15.075 xmax : -7.641667 ymin : 7.191667 ymax : 12.675 私の質問は、これらの2つのラスターが一致する解像度と範囲を持つためには、 A）raster::aggregate関数を使用する > 0.008333333/0.002083333 [1] …

8 raster  r  extents  resolution 

SpatialPolygonsDataFrameとSpatialPointsDataFrameを空間結合したいのですが。 コードの下で、複数のポリゴン内の1つのポイントデータのみを取得します。 sp::over(meuse, srdf) 一部のデータは複数のポリゴンにポイントされます。ポイントデータを取得したいのですが、すべてのポリゴンデータを保持します。 例： ポリゴンaとb内のポイントA。 今得ることができます： | point | polygon | ------------------- | A | a | 私はこのようになります： | point | polygon | ------------------- | A | a | | A | b | データを使用： coordinates(meuse) = ~x+y sr1=Polygons(list(Polygon(cbind(c(180114, 180553, 181127, 181477, 181294, 181007, 180409, 180162, 180114), c(332349, …

8 r  sp 

レーダーから得たビッグデータの扱い方がわかりません。 1日あたり144の画像.tifがあり、各ピクセルの1日の平均値を抽出し、1つの画像.tifにエクスポートする必要があります。 144枚すべての画像をスタックすることを考えましたが、続行する方法と結果を1つの画像にエクスポートする方法がわかりません。 私はこれを始めました： set("D:/R/Radar/Rasday1") f <- list.files(getwd()) ras <- lapply(f,raster) STACK1 <- stack(ras) ## with dimensions :1006, 804, 808824 (nrow, ncol, ncell) 次に、ラスターを作成して次のようなことを実行してもよいかどうかわかりません： r1<- raster(ncol=804, nrow=1006) media<-mean(STACK1, r1)

8 r  extract 

サンプリングサイトのマップを作成する必要があります。ジャーナルは、このマップが世界のどこにあるかを示す差し込みマップを作成するように要求されました。これが私のマップコードです。 library(maps) library(GISTools) map('state', fill = FALSE, xlim = c(-125, -114), ylim = c(32.2, 42.5), xlab = "lon", ylab = "lat") map.axes(cex.axis=0.8) points(-121.6945, 39.36708, bg = "black", pch = 21) maps::map.scale(x=-124, y=34, ratio=FALSE, relwidth=0.3) north.arrow(xb=-116, yb=41, len=0.22, lab="N") 差し込みマップを作成する方法を見つけようとしましたが、今のところうまくいきません。

8 r  maps 

道路を使用してGPS位置を比較し、データにエラーがないことを確認する必要があります。 ローマ市の情報があるサイトからmap.osmをダウンロードし、Rの空間線に変換しました。 src <- osmsource_file(file = "rome.osm") bb <- center_bbox(41.9, 12.3, 1000, 1000) rome <- get_osm(bb, src) ways <- find(rome, way(tags(k == "highway"))) ways <- find_down(rome, way(ways)) ways <- subset(rome, ids = ways) hw_lines <- as_sp(ways, "lines") したがって、空間線オブジェクトでdistm関数を使用しようとしましたが、空間オブジェクトを使用できない（適切な形式ではない）と非難されています。 OSMを使用して、RのGPSデータと比較できる道路のみの情報を取得する別の方法はありますか？

8 openstreetmap  r  gps  distance  road