我正在尝试使用较小的值转换多边形shapefile。 colunm 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是一个64位双精度栅格，所有像素值都等于0.
此多边形的要素之间没有重叠，只是边界触摸，所以我没有使用fun作为rasterize的参数。
在rasterize之后，当我从r_pol检查minValue和maxValue时，我获得了0.000002和0.119419，而不是0.08687903和0.1140689 }}
我不知道问题出在哪里。你能帮帮我吗？