##in the beginning... there were some coordinates
library(grid)
library(raster)
library(rgdal)
library(spatstat)
library(sp) 
##information supplied by user/
x=624115.3
y=4859253.68
rad=1000
##setwd("")
projname<-"Sunriver"
##dir.create(projname)
newproj<-"+proj=utm +zone=10 +ellps=GRS80 +datum=NAD83 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"  
CRS<-CRS(newproj)
radlim=1.5*rad
plot(x,y,xlim=c(x-radlim,x+radlim),ylim=c(y-radlim,y+radlim))
diskbuff<-disc(radius=rad, centre=c(x,y), mask=FALSE, npoly=1200)
arcstart<-data.frame(diskbuff$bdry)
circledf<-data.frame(x=arcstart$x,y=arcstart$y)
center<-data.frame(x=x,y=y) 
circlebind<-(cbind(circledf$x,circledf$y))

##The distance raster (1km)
library(raster)
p1 <- rbind(circlebind,circlebind[1,])
pols <- SpatialPolygons( list( Polygons(list(Polygon(p1)), 1)), proj4string=CRS)
qy <- coordinates(pols)
r <- raster(nrow=100, ncol=100,xmn=x-radlim,xmx=x+radlim,ymn=y-radlim,ymx=y+radlim,crs=CRS)
r <- rasterize(pols, r)
d <- distanceFromPoints(r, qy)
q <- mask(d, r)
dist_1km <- q
writeRaster(dist_1km, filename=file.path(getwd(),projname,"1km_circle.img"), overwrite=TRUE)
plot(dist_1km)
summary(dist_1km)