問題タブ [sp]

私のデータは、緯度と経度を度、分、秒の形式で投影したものです。これは、テキスト形式で保存されています ( character-class、 は無視してrow.namesください)。

メーリングリストで提供された回答により、データを10進数形式に変換しました...

その後、 をオブジェクトに変換しdata.frameました。SpatialPoints

次に、次の方法で投影を変換しましたspTransform

この Web サイトを使用して、例の最初の座標を変換し(座標は南と西を指すため、負の値を指定するようにしてください)、次の座標を提供します: 633726.46 mE 8468757.51 mN Zone 18L、これは正しいです。対照的spに、変換後のオブジェクトは次のとおりです。

マニュアルではspTransform、メタデータを正しく提供する必要があると具体的に述べています。私のデータは WGS84 楕円体を使用しており、18L ゾーンにあります。より正確にはepsg:32718です。どこが間違っているのですか？

このrasterパッケージを使用して、ASCII ラスターと ESRI シェープファイルの 2 つのデータセットを読み込みました。ラスターの (水温) データを、湖の海岸線であるシェープファイルの全範囲に抽出したいと考えています。

SpatialPolygonsDataFrame関数を使用して読み込むと、ESRI シェープファイルは として扱われshapefile()ます。

shapefile <- shapefile("shore.shp",verbose=TRUE)

関数を使用raster()して ASCII ラスターを読み込みました。

raster <- raster("1995_001.asc",native=TRUE,crs="+proj=merc +datum=WGS84 +ellps=WGS84 +units=m")

シェープファイルの座標参照情報は次のとおりです。

+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0

ラスターのそれ (つまりcrs、関数内の引数を使用して次のように強制されraster()ます):

+proj=merc +datum=WGS84 +ellps=WGS84 +units=m +towgs84=0,0,0

次にspTransform()、パッケージ内の関数を使用しrgdalて、シェープファイルの空間参照をラスターの空間参照に強制しました。

spTransform(shapefile, CRS(projection(raster)))

最後に、以下を提出しました。

extract(raster,shapefile,method="simple",fun=mean,small=TRUE,na.rm=TRUE,df=FALSE)

ただし、型のオブジェクトをextract()返すだけです。この問題は、座標参照の明示的な強制から生じていると思います。NULLlist

さらに、show()各データセットで関数を使用した結果は次のとおりです。

> show(raster) class : RasterLayer dimensions : 1024, 1024, 1048576 (nrow, ncol, ncell) resolution : 1800, 1800 (x, y) extent : -10288022, -8444822, 4675974, 6519174 (xmin, xmax, ymin, ymax) coord. ref. : NA data source : in memory names : layer values : -9999, 8.97 (min, max)

> show(shapefile) class : SpatialPolygonsDataFrame features : 1 extent : 597568.5, 998261.6, 278635.3, 668182.2 (xmin, xmax, ymin, ymax) coord. ref. : +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0 variables : 3 names : AREA, PERIMETER, HECTARES min values : 59682523455.695, 5543510.075, 5968252.346 max values : 59682523455.695, 5543510.075, 5968252.346

これらのフォーラムで解決策のない同様の質問を多数検索しました。誰か（仮想の）手を貸してくれませんか？

事前にどうもありがとうございました。

インドの 357 グリッドのデータがあります。すべてのグリッドに何らかの値があります。Rでプロットしたい.次の行を使用する

しかし、私は非常に悪い地図を手に入れました。 このイメージを改善するにはどうすればよいですか？R Plot Filled Longitude-Latitude Grid Cells on Mapでいくつかの良い方法を見てき ましたが 、残念ながらこれらの方法は私の場合はうまくいきませんでした。いずれも高く評価されます。

SpatialPolygons(DataFrame) オブジェクトがあります。SpP

ポリゴンごとに、その重心 ( coordinates(SpP)) と同じポリゴン内の他のポイントとの間の最大距離を計算したいと思います (この最も遠いポイントは端にあると思いますか?)。

誰かがそれを行う方法を教えてもらえますか?

もともと SpatialPolygons オブジェクトに含まれていた多数のポリゴンから座標を抽出しようとしています:

SpP オブジェクトから Sr1 と Sr2 の座標を抽出しようとしています。スタック交換の他の場所でこのコードを見ました:

角括弧内のインデックスの不一致を理解できませんが、実行されます。しかし、出力は、Sr1 または Sr2 で指定した座標と一致しません。インデックスのすべての組み合わせを試しましたが、探している答えが得られません!

bubbleパッケージの関数を使用してデータをプロットし、レイヤーspを追加したいと思います。SpatialPolygonsを使用してデータを簡単にプロットできますspplotが、何らかの理由で では機能しませんbubble。例えば：

を に変換することでこれを行ったSpatialPolygons人もSpatialLinesいますが、ggplot2 を介して実行できるように見えますが (投稿を参照)、なぜバブル プロットが で機能しないのか混乱していSpatialPolygonsます。それを機能させるためのトリックはありますか？

名前付きのポリゴン形状ファイルCROPS.shpと次のキーがあります-

rgdalパッケージを使用して形状ファイルを という名前の変数に読み込みますmy_crops。現在、、、 、、のmy_crops@data6 つの異なるフィールドがあります。これらの各フィールドには、1 から 4 までの値があります。Y1Y2Y3Y4Y5Y6

たとえば、 の最初の数行は次のcrops@dataようになります。

私は何をしたいですか？

最初のプロットは field に基づいてY1おり、最後のプロットは field に基づいていますY6。プロット (ポリゴン形状) の色は、各シェープファイル プロットの値 1 ～ 4 に基づいている必要があります。

プロットには、ページの下部に単一の凡例が付随する必要があります。

次を使用して、drawpoly からポリゴンを取得しています。

しかし、ポリゴンから座標を取得しようとしていますが、できません:

poligon@polygons [[1]] クラス「Polygons」スロット「Polygons」のオブジェクト: [[1]] クラス「Polygon」スロット「labpt」のオブジェクト: [1] -46.37327 -23.91955

前もって感謝します

spplotコマンドを使用して、spの単純なデフォルトのコロプレススタイルの色塗りつぶしを使用して、格子状に2つのマップを並べてプロットしています。私が使用している構文は、次の例と同じです。

マップは、意図したとおりに並べて完全かつ正確にプロットされます。ただし、次のように各マップの上に警告が表示されます。

「パケット 1 の長さを使用するエラーは「ユニット」オブジェクトでなければなりません」(チャート 1)

「パケット 2 の長さを使用するエラーは「ユニット」オブジェクトでなければなりません」(チャート 2)

トラブルシューティングのためにここでデータを提供できないことをお詫び申し上げます。しかし、誰かが何が間違っているのかを理解する正しい方向に私を向けることができることを望んでいました.

注: 両方のマップは、次のものと同等のものを使用して、個々のマップとして完全に良好にチャート化されます。

明らかに問題は格子ビューにありますが、私は (Sarkar ch8 を調べようとしたにもかかわらず) 何が問題なのかを知ることができるほどパケットに精通していません。

別の方法として、チャートが必要に応じてマッピングされている場合、これらの警告を簡単にオフにする方法はありますか?

助けてくれてありがとう