#GPS in Decimal Degrees in the form longitude latitude (raw data)
library(maptools)
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)
#making the points into a map projection (latlong in WGS84)
RM215.sp <- SpatialPoints(RM215, proj4string=CRS("+proj=longlat
+datum=WGS84"))
#start time for the sequence
d060101 <- as.POSIXct("2006-01-01", tz="EST")
#sequence in days
study_seq <- seq(from=d060101, length.out=761, by="days")
#sunrise
up.215 <- sunriset(RM215.sp, study_seq, direction="sunrise",
POSIXct.out=TRUE)
#sunset
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)
# put all of the observations side by side sunrise next to sunset then put
all sites in columns
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))
#name columns the site and then sunrise (st) and sunset (ss)
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")
On Tue, Jul 15, 2008 at 6:10 PM, jim holtman <[EMAIL PROTECTED]> wrote:
> I am reasonably sure that it could be more efficient, but it would be
> very helpfulto have some comments in the code and a description of
> "what is the problem you are trying to solve". It looks like you can
> setup some lists and process them more compactly.
>
> On Tue, Jul 15, 2008 at 3:21 PM, stephen sefick <[EMAIL PROTECTED]> wrote:
> > # 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]]
> >
> > ______________________________________________
> > R-help@r-project.org 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.
> >
>
>
>
> --
> Jim Holtman
> Cincinnati, OH
> +1 513 646 9390
>
> What is the problem you are trying to solve?
>
--
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]]
______________________________________________
R-help@r-project.org 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.
