0

皆さん、

専門家の意見が必要です。私は一日の大半をこの問題と戦ってきましたが、選択肢がありません。

風景全体で最後の火災からの時間を表すマップを作成しようとしています

明らかに、火災が発生したことのないスポットもあります。これらはここでは 0 として表されますが、やや恣意的です。

私のコードは以下の通りです：

データはこちら

#example for SO
library(rgdal)
library(raster)
library(plotKML)

##load Data
Fire <- readOGR(dsn="***YOURPATHHERE******/H_Fire_Poly_Jan3_2105", "Fire_poly_KootClip")


##build raster1
r<-raster(res=300,extent(Fire))

pop<-rasterize(Fire,r,"FIRE_YEAR",fun="max",na.rm=TRUE,background=2015)
pop<-projectRaster(pop,res=300, crs="+proj=utm +zone=11 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs")

plot(pop)


##build raster 2
background<-raster(extent(pop),res=300,vals=2015)
projection(background)<-CRS(proj4string(pop))

##make sure all is good
plot(background)
plot(pop,add=T)

##algebra
final<-background-pop

##plot
plot(final)

##export
plotKML(final)

コードを実行したくない場合は、Google Earth で次のように表示されます。 バツ

これが主な問題です。

1) 黄褐色の線が正確な境界であるため、ポリゴンの境界情報が保持されているように見えます。

2）「最大」は常にここで報告されているわけではありません。特に、それが表示されるエッジの周りですか？

3）私の目標は、種分布モデリングでこのマップを他のラスターと組み合わせて使用​​ することです（したがって、ラスタースタック、モデル、および予測関数を使用してマップを作成します）。このレイヤーでマップを作成しようとすると、NA があるはずの場所に NA が原因で穴が山積みになってしまいます。

レイヤーを起動してからこの時間なしでモデルを実行すると、すべて問題ありません。それを含めると、結果のマップの NA は通常、最初の写真で黄褐色の「ポリゴン アウトライン」が表示される場所になります。

このラスターレイヤーを鮮明なエッジにし、現在のように混乱させないようにするにはどうすればよいですか???

これは私の最終的なマップの種のマップがどのように見えるかです, 欠落しているデータは火災マップの問題を反映していません.

c

4

1 に答える 1

0

あなたの混乱は、計算によるものではなく、Google Earth での視覚化によるものだと思います。これを R でプロットすると問題ないように見えますが、projectRaster は使用しないでください。

library(rgdal)
library(raster)
Fire <- shapefile("Fire_poly_KootClip")
r <- raster(res=300,extent(Fire))
pop <- rasterize(Fire, r, "FIRE_YEAR", fun="max", na.rm=TRUE, background=2015)
x <- pop-2015
plot(x)

より簡単に行うことです

Fire$since <- Fire$FIRE_YEAR - 2015
x <- rasterize(Fire, r, "FIRE_YEAR", fun="max", na.rm=TRUE, background=-100)
plot(x)

-100 ではありません。火災がないことは、「今」ではなく「昔」に似ているからです。

于 2015-03-21T21:51:53.327 に答える