![]() ![]() #Write out final km/km2 road density raster #Rescale step1 output with respect to cell size(30m) and radius of a circleįinal_rdDensity <- (step1*0.03)/3.14159265 #multiply km2/km2 density by number of cells in the moving window ![]() #Read-in an already created raster mask (cells are all set to 0) I then convert the km2/km2 output to km/km2. I currently execute a backwards approach using raster::focal function, calculating a density of burned in road features. Here is a link on how line density works: Īnd this is how to use it in a python (arcpy) script: ArcMap has a quick and easy tool that handles this, but I need a pure R solution. Essentially I need to calculate a km/km2 road density within a 500m pixel search radius. I need to calculate the magnitude-per-unit area of polylines that fall within a radius around each cell.
0 Comments
Leave a Reply. |