使用基于功能输出的值创建栅格

Question

我有一个函数，可根据输入斜率和距离计算价格。我想将价格写入栅格作为栅格值。我怎么做？ OpenSource和ArcMap解决方案可以工作。

slopeRaster = "slope.tif" 
emptyRaster = "emptyraster.tif" # How do I create an empty raster? 
road = "road.shp" 

for cell in emptyraster: 
    # get slope from sloperaster at cell location 
    ... 
    slope = ... 

    # get distance to nearest road from center of cell 
    ... 
    distance = ... 

    # calculate price for cell 
    price = pricefunc(slope, distance) 

    # write price to cell as value # How do I write a value to a raster 

Answer 1

你可以很容易地在R中做到这一点。我建议你download and install it（它是免费的和开源的）。你唯一需要做的就是研究如何在R中编写你的价格函数，这就是为什么我建议你发布代码。一旦你定义了pricefunc，你就可以从R命令行运行这些命令。

# Install required packages 
install.packages(c("raster" , "spatstat" , "sp" , "rgdal") , dep = TRUE) 

# Load required packages 
require(raster) 
require(spatstat) 
require(sp) 
require(rgdal) 

# Read in your data files (you might have to alter the directory paths here, the R default is to look in your $USERHOME$ directory R uses/not \ to delimit directories 
slp <- raster("slope.tif") 
roads <- readShapeLines("road.shp") 


# Create point segment pattern from Spatial Lines 
distPSP <- as.psp(roads) 


# Create point pattern from slope raster values 
slpPPP <- as.ppp(values(slp)) 


# Calculate distances from lines for each cell 
distances <- nncross(slpPPP , distPSP) 


# Create raster with calcualted distances 
rDist <- raster(slp) 
values(rDist) <- distances 


# Define your princefunc() here. It should take two input values, slope and distance and return one value, which I have called price 
pricefunc <- function(slp , dist){ 
    ...my code 
     ... more code 
    ...more code 
    return(price) 
} 


# Calculate price raster using your price function and save as output.tif 
rPrice <- overlay(slp , rDist , fun = function(x , y){ pricefunc(x , y) } , filename = "output.tif")