データピクセルがないためラスターを合計することはできません

毎日の土壌水分データを含む一連のラスター（.vrt）があります。月ごとの測定値を得るために、すべてのラスターのピクセルを合計します。ただし、毎日全世界がカバーされていないため、他の日にメジャーが存在する場所ではデータ値がありません。

私がしたいことは、各ラスターの値を合計することです。ただし、nodataピクセルが合計に含まれるたびに、結果のピクセルはnodataとして直接分類されるようです。私は逆に、すべてのnodata値を無視し、残りを合計します。

私は問題を解決する2つの方法を考えました：

1. nodata値を無視してラスターを合計する
2. nodataピクセルを値0に変換してから、すべてのラスターを合計する

残念ながら、これを行うためのツールが見つかりません。

誰も私を助けることができますか？

ArcGISにアクセスできる場合、セル統計ツールには、ラスタをsum \ mean \ min \ maxしたnodataを無視するオプションがあります。

まず、gdal_calc.pyを使用して、すべての-9999を0に変更し、NoData値を0に設定できます。

例えば：

gdal_calc.py -A input.tif --outfile=input_with_NoData.tif --calc="A+9999*(A==-9999)" --NoDataValue=0

次に、gdal_translateで-a_nodataオプションの後にnoneを使用して、NoData値を無視できます。

-a_nodata 値：

指定されたnodata値を出力バンドに割り当てます。GDAL 1.8.0以降、noneに設定して、ソースファイルにnodata値が存在する場合に出力ファイルにnodata値を設定しないようにすることができます。

例：

gdal_translate -a_nodata none input_with_NoData.tif output_without_NoData.tif

Rで：

library(raster)
s <- stack('file1', 'file2', ...)
ss <- sum(s, na.rm=TRUE)

Spatial Analystの条件付きツールは、null値をゼロに変換する簡単な方法です。その後、合計を取ることは簡単です。

私は少し前に同じ問題を抱えていましたが、それをなんとか解決しました。

QGISでのnodata値の表示についてのメモと同様に、nodata値は、実際に値が-9999、0、NULLなどであるかどうかに関係なく、マップウィンドウおよび「オブジェクト情報」ツールを使用するときにnodata値として常に表示されます。したがって、nodata値を0に再分類した後、レイヤー設定に移動すると、メタデータでnodataの値が0であることがわかります。

私はついに、クエストに最適なツール、GRASSを見つけましたr.series。[集計操作]設定[最大]および[NULLの伝播]をオフにすると、ツールは要件に合ったレイヤーを作成します。

私はしばらく前にこの問題を抱えていたことを覚えています。解決策として、ラスター内のすべてのセルに+1し、ラスターを一緒に追加してから、最後の減算を行いました（ただし、追加した+1の数が多い）。それはかなり粗野ですが、仕事が終われば...！

または、QGISの[処理]ツールボックスで[グリッド値の再分類]を使用して、null値をゼロに変換します。