r-raster

1熱

1答えて

rasterパッケージのrasterize機能を使用することがよくあります。マニュアルに記載されているように、fun引数で使用されるカスタム関数は、na.rm引数を受け入れる必要があります。これは、一般的にすなわち、カスタム関数が「ドット」で書かれていることを意味： funA <- function(x,...)length(x) しかし、第2の提案されたアプローチは、明示的なna.rm引数で

1熱

1答えて

calc（）関数を使用しているエラー：この関数を使用できません。

大きなマルチバンドラスタで計算を実行し、RasterBrickをエクスポートする必要があります。メモリの目的で、ラスタパッケージのcalc()関数を使用します効率。関数は、自分自身で細かい動作しますが、私はcalc()に含めるしようとすると、私はこのエラーを取得しておいてください。 Error in .calcTest(x[1:5], fun, na.rm, forcefun, forceappl

0熱

1答えて

RColorBrewerをプロットするときのエラー

ラスタデータからマップのプロットを作成しようとしています。私はこのコードを使用しています： library(raster) library(rgdal) library(classInt) library(RColorBrewer) NDII = raster("G:\\Sheyenne\\image_differencing\\NDII\\differenced.tif") valu

1熱

1答えて

ラスタからRのポリゴンに領域加重和を抽出するにはどうすればよいですか？

人数Nを保持するラスタ（クラス "RasterLayer"）を持っています。（SpatialPolygonsDataFrameクラスの）管理境界を表すポリゴンがあります。ラスタレイヤーからポリゴンに値を抽出するには、ポリゴンをオーバーレイするラスタータイルの面積の割合で重み付けされます。単一のポリゴンに重なる値（4,8,12,16）を持つ4つのラスタタイルがあるとします。各ラスタの25％がポリゴ

0熱

1答えて

はエラーを返します

ラスタスタックデータから各セルのベクトルにNAを代入または補間する必要があります。私は、NAを補間するために、NAとfun_interpを代入するfun_subという2つの関数を持っています。 fun_subがうまく機能しています。しかし、fun_interpは動作しませんが、私は理由を見つけることができません。私は標準機能を使用して、この道を行くだろう #------EXAMPLE libra

0熱

1答えて

"エラー：ラスタIO時のエラー"を修正するにはどうすればよいですか？

私は種分布モデルの共変量として30m解像度（画像あたり1.4Gb）で12枚の高解像度LANDSAT画像を扱おうとしています。しかし、次のエラーが発生します。“Error: Failure during raster IO”ラスタスタックから値を抽出するとき。 library(raster) preds <- suppressWarnings(raster::stack(files)) env

50熱

1答えて

Rの「ラスタ」パッケージを作成すると、GeoTIFFの正と負の回転行列が区別されますか？

Rのラスタパッケージは、GeoTIFFの正と負の回転を区別していないようです。 Rは回転行列の負の符号を無視しているからです。私は確認するためにrasterソースコードを掘り下げるほど精通していませんが、問題を示す再現可能な例を作成しました： Rのロゴを読み取り、GeoTIFFとして保存してください。 library(raster) b <- brick(system.file("externa

2熱

1答えて

writeRaster出力ファイルのサイズ

私は、ラスターレンガオブジェクトとしてマルチバンド画像を読み込み、さまざまな計算を行うバンドを反復してラスタを新しい.tifとして書き出す関数を持っています。これはうまくいきますが、新しい画像ファイルのファイルサイズは約4倍です（元の画像には4つのバンドがあると仮定しています）。 writeRaster（）関数に私が気づいていないパラメータがあるのか、それとも出力画像が基本的に入力と同じファイ

1熱

1答えて

netCDFから時系列を抽出するR

私はこれを作成しましたfile TRMM_3B42_Daily製品を使用して1998-01-01〜1998-12-31を作成しました。これは私がRで使用されるスクリプトです： lon=seq(-91.875,-86.875,by= 0.25) lat=seq(13.875,16.875,by= 0.25) x_dim <- ncdim_def("lon", "degrees_east", l