raster

1熱

1答えて

GDALを使ってjp2000をgeotiffに変換するには？

私はgdal_translateを使って、jp2000（.jp2）形式のSentinel-2 imageryをジオトラフ形式に変換しようとしています。ただし、.jp2形式が認識されないようです。 jp2000形式をジオテフに変換するにはどのような方法を使用しますか？ $ gdal_translate B02.jp2 B02.tif ERROR 4: `B02.jp2' not recognise

0熱

1答えて

writeRasterのtiffが開かない

私はいくつかのラスタ計算をしていますが、結果を新しいtiffに書きたいと思います。しかし、RのwriteRasterからの出力.tifは、一般的に使用されているイメージ化ソフトウェアのいずれでも開かない。それを修正するために私のコードに追加すべきことについて何か提案してください。事前におかげで、 setwd() library(sp) library(raster) library(rg

0熱

2答えて

RasterLayerの平均は計算しません（Rラスタパッケージ）

したがって、複数のラスタ（x1、x2、x3、x4、...）を積み重ね、これらのすべて（xmaster）から平均ラスタを計算しました。しかし、私はそのラスタ（xmaster）の平均ピクセル値を求めます。通常、私は要約統計量を表示し、平均値を呼び出します....しかし、平均値は 'xmaster'の要約には現れません！どうしてこんなことが起こっているのか分かりません。私の下のスクリプトを参照してくださ

0熱

1答えて

エクスポート画像ファイル直接使うRのデータフレームに基づいて

I持っRにおけるこの次のデータフレーム： df <- data.frame(a = c(0, 0, 1, 1), b = c(0, 1, 0, 1), col = c("red", "green", "green", "yellow")) 私は1つの画素を含むことになるラスタ画像ファイル（PNG、TIFなど）をエクスポートに興味適切な色で、データフレームの各エン

0熱

1答えて

Landsatデータでgreenbrownパッケージを使用

greenbrownパッケージの 'PhenologyRaster'機能を使用して、自分のスタディエリアの成長期をモデル化しようとしています。しかし、関数を実行するたびに空の出力が得られます（たとえば、SOS.2016のレイヤーはNAと表示されます）。私の質問は次のとおりです。1年のデータで機能を実行しているか、Landsatの時系列が多少不規則なため（つまり1年あたり〜30シーンの頻度）、問題が

0熱

1答えて

Gdal_calcがナノ出力で失敗する

コマンドラインからgdal_calcを使用していくつかのラスタ計算を行いたいが、ナノで満たされた空のラスタを続けておきたい。例えば私はラスタの値を置き換えます。 gdal_calc.py -A input.tif --outfile=output.tif --calc="10*(A==0)+20*(A==1)" この例では0に等しい値が10に置き換える必要がありますし、1に等しい値が20で交

-1熱

1答えて

大きなラスタ周波数テーブル/カウント

freq（）を使用して、Rのラスタのピクセル値の頻度/カウントを計算しようとしています。は比較のための2つの例ラスタを作成します。 library(raster) RastSmall <- raster(nrow=70, ncol=70) RastBig <- raster(nrow=7000, ncol=7000) set.seed(0) RastSmall[] <- round(ru

-2熱

1答えて

点をベクトル線に変換するアルゴリズム。

セグメント（上の画像）からなる線をベクトル線（下の画像）に変換するにはどうしたらいいですか？いくつかのアルゴリズムがありますか？好ましくはC++であるが、擬似コードでさえOKである。

1熱

1答えて

特定の色でRのラスタースタックの最大値をプロットする

私は5ラスターのラスタスタックを持っています（この例では本発明例を作成しました）。このスタックから最大値と最小値の2つのプロットを作成したいと思います。 max/min（a）を使ってかなり簡単です。 r <- raster(nrows=10,ncols=100) r[] <- rnorm(1000) s <- stack(r,r,r,r,r) a<-max(s, na.rm=T) plot