2017-11-22 9 views
-1

SDMの環境ラスタレイヤをRにアップロードしようとしています。それらはすべて同じ範囲、ピクセルサイズなどであり、ほとんどのレイヤーはアップロードされていますが、私のBIOCLIMレイヤーはアップロードされていません。環境ラスタレイヤをアップロードして積み重ねる際に問題が発生するR

私は各層をプロットすることができたし、それがうまく見えたが、私はそれらを一緒に積み重ねたときに、私はこのエラーを取得:.makeRasterList(RLIST)で

無視データのない層(複数可)

私のスタックを見ると、私の19の層のうちの1つだけがそこにあります。このエラーは、ピクセルの一部にデータがないことを意味していると仮定していますが、これを修正する方法はわかりません。

編集: 私はそれは(私はそれが私の他の非気候の層が働いていたので、かつてのかもしれ信じる)のArcMapまたはRに問題があるかはわかりません。ここでは19層の一つのための私のコードは次のとおりです。

ann_mean_temp <- 'ann_mean_temp3.tif' 
ann_mean_temp=brick(ann_mean_temp) 
plot(ann_mean_temp) 

ann_mean_temp 
class  : RasterBrick 
dimensions : 785, 887, 696295, 1 (nrow, ncol, ncell, nlayers) 
resolution : 1000, 1000 (x, y) 
extent  : 936315.5, 1823316, -341835.1, 443164.9 (xmin, xmax, ymin, 
ymax) 
coord. ref. : +proj=aea +lat_1=20 +lat_2=60 +lat_0=40 +lon_0=-96 +x_0=0 
+y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0 
data source : C:\Users\Anna\Documents\ArcGIS\R\ann_mean_temp3.tif 
names  : ann_mean_temp3 
min values :    57 
max values :   155 

climate <- stack(ann_mean_temp, diurnal_range, Isothermality, 
      temp_seasonality, max_temp_warmest, min_temp_coldest, 
      temp_range, temp_wettest, temp_driest, 
      mean_temp_warmest, mean_temp_coldest, ann_precip, 
      precip_wettest_m, precip_driest_m, precip_seasonality, 
      precip_wettest_q, precip_driest_q, precip_warmest, 
      precip_coldest) 

私は今、このエラーを取得しています: compareRasterでエラー(X):異なる数または列

私も中にいることを言及する必要がありますarcmapはBIOCLIMから各レイヤーをダウンロードし、それを投影し、マスクで抽出して各レイヤーを私の研究領域にクリップし、tifファイルとしてエクスポートすることでした。

+0

メッセージには、関連付けられた値を持たないRasterLayersを使用しているRasterStackを作成していることが示唆されています。おそらく2つのレイヤーを持つ例を挙げることができますか? – RobertH

+0

こんにちは。応答していただきありがとうございます。 私は積み重ねるためにRにインポートしようとしている19個のBIOCLIMレイヤを持っています。例えば、2つの層は年降水量と年平均気温です。 私は各レイヤーを個別にインポートしてプロットすることができましたが、一緒に積み重ねるとこのメッセージが表示されます。 – Raspberria

+0

はい、あなたは何か間違っているようです。あなたがコードを表示しないと、それが何であるかを指摘することはできません。できるだけ問題を単純化してみてください。空の(値が関連付けられていない)RasterLayerを追加しているようです。 – RobertH

答えて

0

all I did was download each layer from BIOCLIM, project it, extract by mask to clip each layer to my study area, and export it as a tif file.

そして、それは問題ないはずですが、あなたはArcMapには、(環境設定を確認してください)行うのは難しいことができ、ラスタの原点と解像度を変更しないことを確認する必要があります。また

は、現在、いくつかのコーヒーを持って行く

f <- list.files(pattern="bio") 
s <- stack(f) 

e <- extent(c(936315.5, 1823316, -341835.1, 443164.9)) 
r <- raster(nrow=758, ncol=887, ext=e, res=1000, crs="+proj=aea +lat_1=20 +lat_2=60 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83") 

y <- projectRaster(s, r) 

R.

内のすべての手順を実行します。このプロセスは、あなたが帰ってきたときに行うべきです。

関連する問題