問題タブ [sp]

spacetime::SDFDF(完全な時空間グリッドを持つ時空間データ) を時間でサブセット化するにはどうすればよいですか?

ソファ、私は試しました：

しかし、それは私にerror: ncol(i) == 2 is not TRUE.

何か案は？

皆さん、

専門家の意見が必要です。私は一日の大半をこの問題と戦ってきましたが、選択肢がありません。

風景全体で最後の火災からの時間を表すマップを作成しようとしています

明らかに、火災が発生したことのないスポットもあります。これらはここでは 0 として表されますが、やや恣意的です。

私のコードは以下の通りです：

データはこちら

コードを実行したくない場合は、Google Earth で次のように表示されます。

これが主な問題です。

1) 黄褐色の線が正確な境界であるため、ポリゴンの境界情報が保持されているように見えます。

2）「最大」は常にここで報告されているわけではありません。特に、それが表示されるエッジの周りですか？

3）私の目標は、種分布モデリングでこのマップを他のラスターと組み合わせて使用​​ することです（したがって、ラスタースタック、モデル、および予測関数を使用してマップを作成します）。このレイヤーでマップを作成しようとすると、NA があるはずの場所に NA が原因で穴が山積みになってしまいます。

レイヤーを起動してからこの時間なしでモデルを実行すると、すべて問題ありません。それを含めると、結果のマップの NA は通常、最初の写真で黄褐色の「ポリゴン アウトライン」が表示される場所になります。

このラスターレイヤーを鮮明なエッジにし、現在のように混乱させないようにするにはどうすればよいですか???

これは私の最終的なマップの種のマップがどのように見えるかです, 欠落しているデータは火災マップの問題を反映していません.

オブジェクトの正確な地図の縮尺を m または km 単位で調べ、SpatialPolygonsそれをプロットしてspplot、読者が距離を把握できるようにしたいと考えています。???つまり、プロットの に入力する数値を知りたいということです。

多分何かsp::coordinatesまたはImap::gdist？SpatialPolygonsこれを見つけるために必要な事前の再投影はありますか?

別の SpatialPolygon をカバーする SpatialPolygons オブジェクト (Hexagonal グリッド) を作成できるようにしたいと考えています。

すべての六角形の直径を 1 km にし (理想的にはこれを変更できます)、すべての六角形を合わせてオブジェクト全体をカバーしたいと考えています。以下の方法は、そのほんの一部をカバーしているようです...

sp以下は、パッケージを使用した私の試みです：

いつものようにどんな助けも大歓迎です...

注: Edzer Pebesma の提案により、この質問は R-sig-geo (こちら) にクロスポストされ、いくつかの良い反応がありました。

を使用すると、次の予期しない結果が発生しましたcheckPolygonsHoles。

ここで厄介な副作用の 1 つは、元のオブジェクトwrld_simplも変更されることです。この結果は私にはバグのように見えますか、それとも何か見逃していますか?

shape_with_holePS: 以前に編集されたオブジェクトcheckPolygonsHolesは、引き続き奇妙な動作をします: