# I am sure that I could be more efficient than this but how? Thanks in
advance.
#GPS in Decimal Degrees in the form longitude latitude
RM215 <- matrix(c(-82.1461363, 33.5959109), nrow=1)
SC <- matrix(c(-82.025888, 33.606454), nrow=1)
RM202 <- matrix(c(-81.9906723, 33.5027653), nrow=1)
RM198 <- matrix(c(-81.926823, 33.4634678), nrow=1)
HC <- matrix(c(-81.920505, 33.46192), nrow=1)
RM190 <- matrix(c(-81.9317347, 33.3839097), nrow=1)
BC <- matrix(c(-81.948189, 33.373043), nrow=1)
RM185 <- matrix(c(-81.941, 33.3453), nrow=1)
RM179 <- matrix(c(-81.890929, 33.317914), nrow=1)
RM148 <- matrix(c(-81.7547337, 33.1514072), nrow=1)
RM119 <- matrix(c(-81.501919, 32.94038), nrow=1)
RM61 <- matrix(c(-81.262388, 32.524739), nrow=1)
RM215.sp <- SpatialPoints(RM215, proj4string=CRS("+proj=longlat
+datum=WGS84"))
d060101 <- as.POSIXct("2006-01-01", tz="EST")
study_seq <- seq(from=d060101, length.out=761, by="days")
up.215 <- sunriset(RM215.sp, study_seq, direction="sunrise",
POSIXct.out=TRUE)
down.215 <- sunriset(RM215.sp, study_seq, direction="sunset",
POSIXct.out=TRUE)
SC.sp <- SpatialPoints(SC, proj4string=CRS("+proj=longlat +datum=WGS84"))
up.SC <- sunriset(SC.sp, study_seq, direction="sunrise", POSIXct.out=TRUE)
down.SC <- sunriset(SC.sp, study_seq, direction="sunset", POSIXct.out=TRUE)
RM202.sp <- SpatialPoints(RM202, proj4string=CRS("+proj=longlat
+datum=WGS84"))
up.202 <- sunriset(RM202.sp, study_seq, direction="sunrise",
POSIXct.out=TRUE)
down.202 <- sunriset(RM202.sp, study_seq, direction="sunset",
POSIXct.out=TRUE)
RM198.sp <- SpatialPoints(RM198, proj4string=CRS("+proj=longlat
+datum=WGS84"))
up.198 <- sunriset(RM198.sp, study_seq, direction="sunrise",
POSIXct.out=TRUE)
down.198 <- sunriset(RM198.sp, study_seq, direction="sunset",
POSIXct.out=TRUE)
HC.sp <- SpatialPoints(HC, proj4string=CRS("+proj=longlat +datum=WGS84"))
up.HC <- sunriset(HC.sp, study_seq, direction="sunrise", POSIXct.out=TRUE)
down.HC <- sunriset(HC.sp, study_seq, direction="sunset", POSIXct.out=TRUE)
RM190.sp <- SpatialPoints(RM190, proj4string=CRS("+proj=longlat
+datum=WGS84"))
up.190 <- sunriset(RM190.sp, study_seq, direction="sunrise",
POSIXct.out=TRUE)
down.190 <- sunriset(RM190.sp, study_seq, direction="sunset",
POSIXct.out=TRUE)
BC.sp <- SpatialPoints(BC, proj4string=CRS("+proj=longlat +datum=WGS84"))
up.BC <- sunriset(BC.sp, study_seq, direction="sunrise", POSIXct.out=TRUE)
down.BC <- sunriset(BC.sp, study_seq, direction="sunset", POSIXct.out=TRUE)
RM185.sp <- SpatialPoints(RM185, proj4string=CRS("+proj=longlat
+datum=WGS84"))
up.185 <- sunriset(RM185.sp, study_seq, direction="sunrise",
POSIXct.out=TRUE)
down.185 <- sunriset(RM185.sp, study_seq, direction="sunset",
POSIXct.out=TRUE)
RM179.sp <- SpatialPoints(RM179, proj4string=CRS("+proj=longlat
+datum=WGS84"))
up.179 <- sunriset(RM179.sp, study_seq, direction="sunrise",
POSIXct.out=TRUE)
down.179 <- sunriset(RM179.sp, study_seq, direction="sunset",
POSIXct.out=TRUE)
RM148.sp <- SpatialPoints(RM148, proj4string=CRS("+proj=longlat
+datum=WGS84"))
up.148 <- sunriset(RM148.sp, study_seq, direction="sunrise",
POSIXct.out=TRUE)
down.148 <- sunriset(RM148.sp, study_seq, direction="sunset",
POSIXct.out=TRUE)
RM119.sp <- SpatialPoints(RM119, proj4string=CRS("+proj=longlat
+datum=WGS84"))
up.119 <- sunriset(RM119.sp, study_seq, direction="sunrise",
POSIXct.out=TRUE)
down.119 <- sunriset(RM119.sp, study_seq, direction="sunset",
POSIXct.out=TRUE)
RM61.sp <- SpatialPoints(RM61, proj4string=CRS("+proj=longlat
+datum=WGS84"))
up.61 <- sunriset(RM61.sp, study_seq, direction="sunrise", POSIXct.out=TRUE)
down.61 <- sunriset(RM61.sp, study_seq, direction="sunset",
POSIXct.out=TRUE)
sunriseset <- cbind(as.data.frame(up.215$time),
as.data.frame(down.215$time), as.data.frame(up.SC$time),
as.data.frame(down.SC$time),
as.data.frame(up.202$time), as.data.frame(down.202$time),
as.data.frame(up.198$time), as.data.frame(down.198$time),
as.data.frame(up.HC$time), as.data.frame(down.HC$time),
as.data.frame(up.190$time), as.data.frame(down.190$time),
as.data.frame(up.BC$time),
as.data.frame(down.BC$time), as.data.frame(up.185$time),
as.data.frame(down.185$time),
as.data.frame(up.179$time), as.data.frame(down.179$time),
as.data.frame(up.148$time), as.data.frame(down.148$time),
as.data.frame(up.119$time), as.data.frame(down.119$time),
as.data.frame(up.61$time), as.data.frame(down.61$time))
colnames(sunriseset) <- c("RM215sr", "RM215ss","SCsr", "SCss","RM202sr",
"RM202ss" ,"RM198sr", "RM198sr",
"HCsr", "HCsr", "RM190sr", "RM190ss","BCsr","BCss","RM185sr",
"RM185ss","RM179sr", "RM179sss",
"RM148sr", "RM148ss","RM119sr", "RM119ss","RM61sr", "RM61sr")
--
Let's not spend our time and resources thinking about things that are so
little or so large that all they really do for us is puff us up and make us
feel like gods. We are mammals, and have not exhausted the annoying little
problems of being mammals.
-K. Mullis
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
