空間ポリゴン オブジェクトと空間ポイント オブジェクトがあります。後者は、データフレーム内の xy 緯度データ ベクトル (それぞれ緯度と経度と呼ばれます) から作成されましたが、前者は rgdal を使用して直接 R に読み込まれました。私のコードは次のとおりです。
boros <- readOGR(dsn = ".", "nybb")
rats <- read.csv("nycrats_missing_latlong_removed_4.2.14.csv", header = TRUE)
coordinates(rats) <- ~longitude + latitude
この時点では、どちらの空間オブジェクトも投影されません。これらのオブジェクトを次のように投影すると:
proj4string(boros) <- CRS("+proj=lcc")
proj4string(rats) <- CRS("+proj=lcc")
両方のオブジェクトが投影され、両方とも次のように plot() 関数で正常にマップされます。
plot(boros)
plot(rats)
ただし、それらを一緒にプロットしようとすると:
plot(boros)
plot(rats, add = TRUE)
私は最初のプロットのみを取得し、rats オブジェクトを boros に重ねません。ただし、これは大きな問題です。エラー メッセージが表示されないため、これら 2 つの空間オブジェクトが互いに通信できる間の切断がどこにあるのかを判断できませんでした。どちらのコマンドもエラーや警告なしでスムーズに実行されますが、プロットは 1 つしか残っていません。そして、proj4string() で各オブジェクトの射影をチェックすると、各オブジェクトに対して同じ射影が返されます。
plot() を使用して一緒にマッピングできるように、CRS と投影が一致する 2 つの空間オブジェクトを作成するさまざまな方法を試すために、数日間にわたって何時間も費やしてきました。ちなみに、私がとったアプローチの 1 つは、ArcGIS でラット オブジェクトのシェープファイルを作成することでした。これは、ファイルを作成するのにうまく機能しました。しかし、R で 2 つの空間オブジェクトを連携させることができないという問題がまだ残っています。どんな助けでも大歓迎です。上記の 2 つのデータ ファイルを含むドロップボックス リンクも含めました: https://www.dropbox.com/sh/x0urdo6guprnw8y/tQdfzSZ384