小さな値を持つポリゴン シェープファイルを変換しようとしています。列の値のpropEmp範囲は0.000002～0.119419です。これは私の試みです：
# Load shapefile
emp_planejado <- shapefile("./planejado/7_wgs/emp_planejado.shp")
# Load raster model
r_bioma<- raster("./_GRID/grid_caatinga_disol_64bit.tif")
# List values from tipologia field
tipologia<-unique(emp_planejado$tipologia)
for (tp in tipologia){
# Select features for each value in tipologia
tipo<- emp_planejado[emp_planejado$tipologia==tp,]
# Rasterize
r_pol <- rasterize(tipo,r_bioma,field="propEmp",background=NA,mask=F)
# Merge
raster_merge <- merge(r_pol,r_bioma)
# Save raster
writeRaster(raster_merge,filename= paste0("./planejado/8_raster/",tp,"_planejado"),format="GTiff",NAflag=-9999,overwrite=TRUE)
}
r_biomaすべてのピクセル値が 0 に等しい倍精度ラスターの 64 ビットです。
このポリゴンのフィーチャ間にオーバーラップはなく、境界に接しているだけなのでfun、 の引数として使用しませんでしrasterizeた。
の後rasterize、 から minValue と maxValue をチェックすると、 and を取得r_polする代わりに、0.000002and0.119419を取得0.08687903し0.1140689ます。
どこに問題があるのかわからない。私たちを手伝ってくれますか？