問題タブ [sp]

マップの拡大部分を元のマップに追加し、元のマップと拡大/ズーム部分の両方を示す 1 つのマップを最終製品として作成したいと考えています。meuse例としてデータセットを使用します。

rasterまたはrasterVISパッケージに、拡大されたラスターの一部を元のマップに追加できるコマンドがあるかどうかはわかりません。機能を試してみましたparが、うまくいきません。どんな提案も歓迎します。

を使用してデジタル標高モデルからルートの標高データを抽出しています

• new.points は、mygrid の CRS で変換された約 7000 ポイントの座標 (経度/緯度) を持つ SpatialPoints オブジェクトです。
• mygrid は、800 万を超える要素を持つ SpatialGridDataFrame です
  (詳細については、前の質問を参照してください) 。

my.elev にいくつかの NA 値があるので、コードをデバッグしたところ、new.points のほぼすべてのポイントが複数回繰り返されていることがわかりました (私のルートでは、いくつかのセグメントが 2 回交差しています)
。 .elev 値が正しい
- 2 番目の値は NA 値 (場合によってはまったく異なる値) です。

new.points で重複した値を削除する問題は簡単に解決できますが、over 関数が同じポイントに対して同じ値を返さないのはなぜでしょうか。

この質問は他の場所で回答されていると確信していますが、検索しても見つけられませんでした。

各都市の人口とともに、国内の都市を表すポイントがあります。郡のポリゴン ファイルもあります。各郡内で最大の都市の場所を見つけたいです。

これはどのように行うことができますか？

ここにいくつかのデータがあります

structure(list(Country = c("us", "us", "us", "us", "us", "us", "us", "us", "us", "us", "us")) 「、
「私たち」、「私たち」、「私たち」、「私たち」、「私たち」、「私たち」、「私たち」、「私たち」、「私たち」、「私たち」、「私たち」、「私たち」、 "us", "us"), City = c("cabarrus", "cox store", "cal-vel", "briarwood townhouses", "barker heights", "davie"
crossroads」、「crab point village」、「azalea」、「chesterfield」、「charlesmont」、「connor」、「clover garden」、「corriher heights」、「callisons」、「crestview acres」、「clegg」、「canaan」 park", "chantilly", "belgrade", "brices crossroads", "bluff", "butner", "bottom", "bandy", "bostian heights"), AccentCity = c("Cabarrus", "Cox Store" 、「Cal-Vel」、「Briarwood Townhouses」、「Barker Heights」、「Davie Crossroads」、「Crab Point Village」、「Azalea」、「Chesterfield」、「Charlesmont」、「Connor」、「Clover Garden」、「コリハー ハイツ」、「カリソンズ」、"Crestview Acres"、"Clegg"、"Canaan Park"、"Chantilly"、"Belgrade"、"Brices Crossroads"、"Bluff"、"Butner"、"Bottom"、"Bandy"、"Bostian Heights")、地域= c("NC", "NC", "NC", "NC", "NC", "NC", "NC", "NC", "NC", "NC", "NC", "NC" 、"NC"、"NC"、"NC"、"NC"、"NC"、"NC"、"NC"、"NC"、"NC"、"NC"、"NC"、"NC"、" NC")、人口 = c(NA_integer_、NA_integer_、NA_integer_、NA_integer_、NA_integer_、NA_integer_、NA_integer_、NA_integer_、NA_integer_、NA_integer_、NA_integer_、NA_integer_、A_integer_、NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_), Latitude = (35.2369444, 35.275, 36.4291667, 35.295, 35.3111111, 35.8319444, 34.7602778, 35.58, 35.81, 5.9341667, 35.7419444, 36.1883333, 35.5605556, 35.0841667, 35.0213889, 35.8655556, 36.2761111, 36.3016667, 34.88, 34.8186111, 35.8377778, 36.1319444, 36.4747222, 35.6419444, 35.7544444), Longitude = c(-80.5419444, -82.0352778, -78.9694444, -81.5238889, -82.4441667, -80.535 , -76.7305556, -82.4713889, -81.6611111, -81.5127778, -78.1486111, -79.4630556, -80.635, -76.7255556, -80.5427778, -78.8497222, -79.7852778, -76.1711111, -77.2352778, -78.1016667, -82.8580556, -78.7569444, - 80.7741667, -81.09, -80.9294444)), .Names = c("国", "都市", "アクセント都市","地域", "人口", "緯度", "経度"), row.names = c(544L, 889L, 551L, 434L, 190L, 975L, 894L, 147L, 717L, 700L, 831L, 773L, 862L, 559L, 915L、753L、584L、695L、262L、437L、372L、537L、406L、178L、02L)、クラス = "data.frame")

ノースカロライナで読むべきいくつかのコード

各郡内で人口が最大の都市を見つけたいと考えています。再現可能な例がなくてすみません。私だったら答えが出ます！

大規模なデータセット (>1500 obs) に空間コンポーネントを追加するのに苦労している初心者の R プログラマー

与えられたリスト: 観測、対応する x/y ポイント (緯度/経度ではない)、および属性。問題の観測の特定の半径 (10) 内にあるすべての観測の変数 "Z" の合計を返すアルゴリズムを作成するにはどうすればよいですか。観測ごとにこれを行う必要があります。

単純化されたデータ ヘッダーは次のように読み取られます。

任意の支援をいただければ幸いです。ありがとう。

約 1000 人ごとにたくさんの場所があります。以前のデータセットの合計は約 250 万で、私の処理スクリプトの実行には約 20 時間かかりました。しかし、今では 2,400 万件の観測があり、コードをクリーンアップし、可能であれば並列処理を使用する必要があると考えています。

各ポイントについて、最も近いポリゴン (ほとんどのポイントはポリゴン内にありません) とそのポリゴンの距離を見つけたいと考えています。ポイントは主に海での観測であり、ポリゴンはポイントに最も近い沿岸 (米国) の郡です。

私は rgeos パッケージの gDistance 関数を使用してこれを行っており、一連のループを実行して (わかっています、わかっています)、個人ごとの処理を分割しています。これをどうにかして plyr/dplyr 構文に移動する方法を見つけようと多くの時間を費やしましたが、うまくいきません。私の問題の一部として、オブジェクト クラスが SpatialPoint と SpatialPoylgonDataFrames であることに関係していると思います。

人々が持っているかもしれない提案に感謝します。

R データにいくつかの形状ファイルが保存されています。私は最初に以下を実行します

library(sp)

次に、次のコードを実行します。

このコードで必要なのは、複数の形状ファイルが複数のパネルにプロットされる PDF ファイルを用意することだけです。ただし、このコードは出力を提供していません。コードのどこが間違っていますか?

多くの同様のトピックがありますが (たとえば、 hereとhere )、SO で見た他の質問とは少し異なる最終目標があります。参考までに、R v 3.1.0 を使用しています。

私は2つのマトリックスを持っています。それぞれにポイントの座標が含まれています。1 つ目 ( A ) には 2,107,377 ポイントが含まれ、2 つ目 ( B ) には 26,577 ポイントが含まれます。

Aの各点が最も近いBの点を見つけたい。つまり、Aのポイント 1 とBの各ポイントの間の距離(26,577 距離) を計算し、最小値を保存します。A のすべての点(最小値 2,107,377)に対してこれを実行したいと思います。目標は、最も近いBの点に基づいてAの点をグループ化することです。したがって、Bの一部の点は割り当てられません。他の (多くの) はA の複数のポイントに割り当てられます。

私が試してみました：

しかし、メモリ割り当ての問題が発生しました (16 Mb を超えるベクトルを割り当てることができませんでした)。

私は今 for ループを実行しています:

しかし、これは同じ結果につながると私は予想していますが、もっとゆっくりです。

まったく別のアプローチ（おそらくパッケージを学習する）を試す前にraster、誰かが何かアイデアを持っているかどうかを確認したいと思いました。

RI には、重複したポイント (座標と属性) を持つ SpatialPointsDataFrame があり、同じデータを持つすべてのポイントを削除したいと思います ...

spパッケージ内に関数が見つかりましたが、そのremove.duplicates()場所でのみ削除されるようです...別の方法はありますか?

ありがとうございました

E.