您可以在R
. 我建议您下载并安装它(它是免费且开源的)。您唯一需要做的就是弄清楚如何在 R 中编写价格函数,这就是我建议您发布该代码的原因。一旦你定义了价格函数,你就可以从 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" )