問題タブ [maptools]

1900 年以降にフロリダと交差したハリケーン/熱帯暴風雨を示すマップの小さな倍数グリッドをプロットしようとしています。いくつかの空間クエリを使用して、このプロジェクトのすべての大西洋の暴風雨のデータベースをサブセット化しました。

私は現在、限られた数のハリケーン トラックのラインシェープファイルを、フロリダ州、隣接するいくつかの州、フロリダのいくつかの主要都市、そしてもちろんオキチョビー湖のポリゴンの上にプロットしています。簡単なコードは次のとおりです。

私が遭遇している3つの大きな問題：

1) ループは、各嵐のマップを私に与えています。できればラベル付きで、各年 (嵐のない年でも) とその年のすべての嵐のグリッドにフロリダ/サウス マップをコードで作成することをお勧めします。

2) ループの特定の行だけでなく、すべての嵐の中で風速の色を設定したいと思います。そうすれば、強い嵐 (1992 年のアンドリューのような) は、その年の唯一の嵐であっても、より暗く表示されます。おそらく、それに応じてスタイルを設定できるカテゴリ (H1、H2 など) 変数を再コーディングすることで、これを解決できます。

3) No. 1 を把握できると仮定すると、各ストーム パスでラベルをレンダリングするのに問題があります。maptools のドキュメントは良くありません。

とにかく、ここまでの出力は次のとおりです (タイトルはシェープファイルの 2 つのフィールドを連結したものです)。

本当の問題はその 1 です。よろしくお願いします。

R のシェープファイルの EU NUTS コードはどこにありますか?

NUTS = 統計のための地域単位の命名法

各 NUTS リージョンには NUTS コードがあります。

.

詳細な NUTS コードは、この Excel ファイルでも確認できます。

NUTS シェープファイルをダウンロードすると、この NUTS コードが見つかりません。

シェープファイルNUTS_2013_20M_SH.zipを次からダウンロードしました。

ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units/nuts#nuts13

ダウンロード後、次のコードを使用してシェープファイルを R にアップロードします。

ノート：

• dsn =シェープファイルを保存したフォルダーへのパスです。
• layer =は、ファイル拡張子を除いたファイルの名前です (例: .shp)

を使用するsummary(NUTS_shape)と、「NUTS_BN_ID」のみが見つかりますが、NUTS コードは見つかりません。またstr(NUTS_shape[4461,])、example-item 4461 のリスト構造を見てみると、NUTS コードはないようです。

NUTS シェープファイルに NUTS コードが含まれているかどうか、または NUTS_BN_ID が何らかの形でこのコードにリンクされているかどうかを知っている人はいますか?

queryResult という都市のいくつかの場所のすべての緯度と経度を含むテーブルがあり、次のことを行います。

1 - 都市のラスター マップを取得します [たとえば、ブラックプール]

これにより、マップ上の次のポイントが返されます

ここで、これらすべての緯度と経度の外側の境界 [都市の境界線] を、次の期待される結果のように描画したいと考えています。

更新 1 : 入力データの提供と提案されたahullソリューションの適用:

@spacedman と @cuttlefish44 によって提案された ahull ソリューションを適用したところ、予想されるポリゴンとは大きく異なる次の結果が得られました。

すべての緯度と経度を含む .csv ファイルを次のリンクからダウンロードできます: Blackpool Lat,Lon

Google が推奨するエリア境界は次のようになります。

生成された Web マップにaddPolylines特定の値の等高線をオーバーレイするために使用していますが、提供されているすべてのデータをプロットしているわけではありません。最初の画像のように、値が125の少なくとも 2 つの領域が明らかですが、リーフレットで生成されたマップには 1 つの領域しか表示されません。データはここからダウンロードでき、SpatialLinesDataFrame オブジェクトを生成する再現可能なスクリプトと、次のようにマップを生成できます。leafletaddPolylines
addPolylinesleaflet

c.linesSP3 本の線の座標がありますが、addPolylines1 つの領域のみの等高線をプロットします

addPolylines最初に提供されたすべての座標をプロットしない理由は明らかではありません。提案を大歓迎します。