2017-04-09 16 views
-1

MODIS由来のNDVI tiff画像のスタックと時系列を作成しました。私は551、674、371374、365(nrow、ncol、ncell、nlayers)のようなラスタ次元を持つラスタ時系列(例えば10年)を持っています。 1.各ピクセルの時系列グラフをプロットする 2.異常な動作を示すピクセルを識別できるように、各ピクセルごとに時系列グラフを1つずつプロットして別のファイルに保存します。ラスター時系列の各ピクセル値をプロットするR

私はRが新しくなって以来、詳細な手順と助けをいただければ幸いです。

ベスト、 クリス

答えて

0

私はコメントでこれを置くつもりだったが、私は少し離れて運ばれ、文字の制限を打ってしまいました。これは、あなたが探している答えではなく、より多くのアドバイスです:

私はdownvoteしませんでしたが、私は最初に私はdownvotted(たくさん)を試みることなくXを行う方法の質問をした最後の時間。私はあなたが最初に試して、あなたがそれらを必要とする特定のものについて質問をすることをお勧めします。例えば

  1. loadまたはread.csvをあなたのデータを取得するために(RへのデータのインポートのためのSO /検索グーグル)
  2. 必要なパッケージをインストールし、:あなたが必要になることがあり

    機能/ステップinstall.packages("ggplot2")。通常、SOの回答がRの場合は、どのパッケージが必要かを示します。

  3. library(dplyr)あなたのデータを読み込んだ後に、データの操作などを行うことがあります。たくさんのデータ計算が必要な場合は、データフレームに非常に似合うようになることをお勧めします。
  4. library(ggplot)画像をプロットする。そこに他のパッケージがあり、これは私の個人的な好みです。
  5. 私は基本的なplot <- ggplot(...)を提案し、連続的に+ geom_point(...)によって(例えばforループ)レイヤーを追加し、それを運のベストggsave(...)

を使用してそれぞれの新しいレイヤーの追加を救います!

+0

ありがとう@Alwin。私は次のコードに別の問題があります。 ファイル< - list.files( ".................."、pattern = "MOD13A1 *") stack1 < - stack (ファイル) stack1 値<-extract(stack1,1:ncell(stack1)) COORD <-xyFromCell(stack1,1:ncell(stack1)) は<-cbind(COORD、値) はwrite.table(コンバイン結合、 "xyvalues.txt") 私は4行目を実行するとNAを取得しますが、ラスタの各セルに値があります。何が間違っているのか、どのようにして値を抽出することができますか? – Kries

関連する問題