| Title: | Plot Objects Moving in Orbits |
|---|---|
| Description: | Visualize the objects in orbits in 2D and 3D. The packages is under developing to plot the orbits of objects in polar coordinate system. See the examples in demo. |
| Authors: | Lele Shu [aut, cre] |
| Maintainer: | Lele Shu <[email protected]> |
| License: | GPL (>= 3) |
| Version: | 0.0.1 |
| Built: | 2025-02-14 05:12:14 UTC |
| Source: | https://github.com/shulele/roundandround |
Calculate c in Focus (c, 0)
ab2c(a = 1, ab)ab2c(a = 1, ab)
a |
Semi-major (Ellipse) or Radium (Ring). |
ab |
Semi-major over semi-minor. ab=1 for a Ring. |
c in Focus (c, 0)
x1=PCS2CCS(a=10, ab=1.5) x2=PCS2CCS(a=9, ab=1.2) c1 = ab2c(a=10, ab=1.5) c2 = ab2c(a=9, ab=1.2) plot(x1, type='n', xlim=c(-10,10), ylim=c(-10,10), asp=1) abline(h=0, v=0, asp=1, lty=2) lines(x1, col=2); points(c1, 0, col=2) lines(x2, col=3); points(c2, 0, col=3)x1=PCS2CCS(a=10, ab=1.5) x2=PCS2CCS(a=9, ab=1.2) c1 = ab2c(a=10, ab=1.5) c2 = ab2c(a=9, ab=1.2) plot(x1, type='n', xlim=c(-10,10), ylim=c(-10,10), asp=1) abline(h=0, v=0, asp=1, lty=2) lines(x1, col=2); points(c1, 0, col=2) lines(x2, col=3); points(c2, 0, col=3)
Add arrows in Polar Coordinate System
Arrow.pcs(theta, r1 = 0, r2 = 1e+06, o1 = c(0, 0), o2 = o1, ab1 = 1, ab2 = ab1, ...)Arrow.pcs(theta, r1 = 0, r2 = 1e+06, o1 = c(0, 0), o2 = o1, ab1 = 1, ab2 = ab1, ...)
theta |
Angle in polar coordinate system |
r1, r2
|
Radius of start and end points of the arrow. |
o1, o2
|
Origin |
ab1, ab2
|
Semi-major over semi-minor. ab=1 for a Ring. |
... |
More options for graphics::arrows function. |
x1=PCS2CCS(a=10, ab=1.5) c1 = ab2c(a=10, ab=1.5) plot(x1, type='n', xlim=c(-10,10), ylim=c(-10,10), asp=1) abline(h=0, v=0, asp=1, lty=2) graphics::lines(x1, col=2); points(c1, 0, col=2) # focus Arrow.pcs(theta = 1:12 * 30, r1=0, r2=10, ab1=1.5, length=.1, col=2, o1 = c(c1,0), o2=c(0,0))x1=PCS2CCS(a=10, ab=1.5) c1 = ab2c(a=10, ab=1.5) plot(x1, type='n', xlim=c(-10,10), ylim=c(-10,10), asp=1) abline(h=0, v=0, asp=1, lty=2) graphics::lines(x1, col=2); points(c1, 0, col=2) # focus Arrow.pcs(theta = 1:12 * 30, r1=0, r2=10, ab1=1.5, length=.1, col=2, o1 = c(c1,0), o2=c(0,0))
Arrow3D
Plot 3D Arrow axis.
Arrow3D
Arrow3D(len = 10, orig = c(0, 0, 0), cols = c(2:4), ...)Arrow3D(len = 10, orig = c(0, 0, 0), cols = c(2:4), ...)
len |
Length of the arrow. |
orig |
Origin of the axis. |
cols |
Colors of axis. |
... |
More options of arrow3d(). |
Degree to Radian
d2r(x)d2r(x)
x |
Degree |
Radian
r = (1:100)/100 * 4 * pi d = r2d(r) rr = d2r(d) plot(d, sin(rr)); abline(h=0 ) abline(v = 360)r = (1:100)/100 * 4 * pi d = r2d(r) rr = d2r(d) plot(d, sin(rr)); abline(h=0 ) abline(v = 360)
Orbit.location
Calculate location of a planet
Orbit.location
Orbit.location(t, p.orb, a = 1, theta = 0, orig = c(0, 0), ab = 1)Orbit.location(t, p.orb, a = 1, theta = 0, orig = c(0, 0), ab = 1)
t |
Time (day). |
p.orb |
Period of the orbit. |
a |
Radius or Semi-major of the orbit. |
theta |
angle in PCS. |
orig |
Reference orgin. |
ab |
Semi-major over semi-minor. ab=1 for a Ring. |
(x,y) in Cartesian Coordinate System.
tday = seq(0, 365, 30) x=Orbit.location(t=tday, p.orb = 365, a=10) plot(PCS2CCS(0:360, a=10), type='l') plotplanet(orig=x, rad = .51) grid()tday = seq(0, 365, 30) x=Orbit.location(t=tday, p.orb = 365, a=10) plot(PCS2CCS(0:360, a=10), type='l') plotplanet(orig=x, rad = .51) grid()
Give the orbit the parameter
orbit.parameter(a, b = NULL, ab = NULL)orbit.parameter(a, b = NULL, ab = NULL)
a |
Semi-major axis |
b |
Semi-minor axis |
ab |
Semi-major over semi-minor. ab=1 for a Ring. |
orbit.parameter(a=1, b=1.5)orbit.parameter(a=1, b=1.5)
Convert Polar Coordinate System to Cartesian Coordinate System.
PCS2CCS(theta = 0:360, a = 1, ab = 1, orig = c(0, 0), rotation = 0, clockwise = FALSE)PCS2CCS(theta = 0:360, a = 1, ab = 1, orig = c(0, 0), rotation = 0, clockwise = FALSE)
theta |
angle in PCS. |
a |
Semi-major (Ellipse) or Radium (Ring). |
ab |
Semi-major over semi-minor. ab=1 for a Ring. |
orig |
Reference orgin. Default = c(0, 0) |
rotation |
Rotation of the theta=0 |
clockwise |
Whether clockwise, Default = FALSE |
(x,y) in Cartesian Coordinate System.
x1=PCS2CCS(a=10, ab=1.5) x2=PCS2CCS(a=9, ab=1.2) c1 = ab2c(a=10, ab=1.5) c2 = ab2c(a=9, ab=1.2) plot(x1, type='n', xlim=c(-10,10), ylim=c(-10,10), asp=1) abline(h=0, v=0, asp=1, lty=2) lines(x1, col=2); points(c1, 0, col=2) lines(x2, col=3); points(c2, 0, col=3) # Test 2 x1=PCS2CCS(a=10, ab=1.5, clockwise = FALSE, rotation=0); x2=PCS2CCS(a=8, ab=1.5, clockwise = FALSE, rotation=45); plot(x1, asp=1, col=terrain.colors(nrow(x1)), pch=19) points(x2, asp=1, col=terrain.colors(nrow(x1)))x1=PCS2CCS(a=10, ab=1.5) x2=PCS2CCS(a=9, ab=1.2) c1 = ab2c(a=10, ab=1.5) c2 = ab2c(a=9, ab=1.2) plot(x1, type='n', xlim=c(-10,10), ylim=c(-10,10), asp=1) abline(h=0, v=0, asp=1, lty=2) lines(x1, col=2); points(c1, 0, col=2) lines(x2, col=3); points(c2, 0, col=3) # Test 2 x1=PCS2CCS(a=10, ab=1.5, clockwise = FALSE, rotation=0); x2=PCS2CCS(a=8, ab=1.5, clockwise = FALSE, rotation=45); plot(x1, asp=1, col=terrain.colors(nrow(x1)), pch=19) points(x2, asp=1, col=terrain.colors(nrow(x1)))
Plot a clock
plotclock(time = c(as.numeric(format(Sys.time(), format = "%H")), as.numeric(format(Sys.time(), format = "%M")), as.numeric(format(Sys.time(), format = "%S"))), rad = 1, ab = 1, orig = c(0, 0), val = 1:12, angle = c((time[1] + ifelse(is.na(time[2]), 0, time[2]/60) + ifelse(is.na(time[3]), 0, time[3]/3600)) * 30, (time[2] + ifelse(is.na(time[3]), 0, time[3]/60)) * 6, time[3] * 6) * (-1) + 90, val.arg = list(col = "blue", cex = 1), arr.arg = list(col = c(1, 3, 2), lwd = ring.arg$lwd/(1:3) * 1.5, lty = rep(1, 3), arrlen = rep(0.1, 3), length = c(0.5, 0.8, 0.9) * rad), ring.arg = list(col = "gold", type = "l", lwd = 4, len.tick = rad * 0.05), fun.plot = lines, add = F, ...)plotclock(time = c(as.numeric(format(Sys.time(), format = "%H")), as.numeric(format(Sys.time(), format = "%M")), as.numeric(format(Sys.time(), format = "%S"))), rad = 1, ab = 1, orig = c(0, 0), val = 1:12, angle = c((time[1] + ifelse(is.na(time[2]), 0, time[2]/60) + ifelse(is.na(time[3]), 0, time[3]/3600)) * 30, (time[2] + ifelse(is.na(time[3]), 0, time[3]/60)) * 6, time[3] * 6) * (-1) + 90, val.arg = list(col = "blue", cex = 1), arr.arg = list(col = c(1, 3, 2), lwd = ring.arg$lwd/(1:3) * 1.5, lty = rep(1, 3), arrlen = rep(0.1, 3), length = c(0.5, 0.8, 0.9) * rad), ring.arg = list(col = "gold", type = "l", lwd = 4, len.tick = rad * 0.05), fun.plot = lines, add = F, ...)
time |
vector, c(h, m, s) |
rad |
Radius of the clock |
ab |
Semi-major over semi-minor. ab=1 for the planet |
orig |
Origin of the clock. |
val |
The labels on clock edge |
angle |
Angle of the three Arrows, default angles are calculated from the time |
val.arg |
Arguments for plot values on clock. |
arr.arg |
Arguments for plot three arrows on clock. |
ring.arg |
Arguments for plot ring on clock. |
fun.plot |
Plot function |
add |
Whether add plot to existing plot. |
... |
More options in plot function. |
plot(0, type='n', xlim=c(-1,1)*1.5, ylim=c(-1,1)*1.5, asp=1) plotclock(add=TRUE,fun.plot = lines, orig=c(0.5,0), rad=.25, val=rep('', 60), time=c(NA, NA, as.numeric(format(Sys.time(), '%S')))) plotclock(add=TRUE,fun.plot = lines, orig=c(0, .5), rad=.25, val=rep('', 60), time=c(NA, as.numeric(format(Sys.time(), '%M')), NA) ) plotclock(add=TRUE,fun.plot = lines, orig=c(-.5,0), rad=.25, val=rep('', 60), time=c(as.numeric(format(Sys.time(), '%H')), NA, NA)) plotclock(add=TRUE, fun.plot=lines)plot(0, type='n', xlim=c(-1,1)*1.5, ylim=c(-1,1)*1.5, asp=1) plotclock(add=TRUE,fun.plot = lines, orig=c(0.5,0), rad=.25, val=rep('', 60), time=c(NA, NA, as.numeric(format(Sys.time(), '%S')))) plotclock(add=TRUE,fun.plot = lines, orig=c(0, .5), rad=.25, val=rep('', 60), time=c(NA, as.numeric(format(Sys.time(), '%M')), NA) ) plotclock(add=TRUE,fun.plot = lines, orig=c(-.5,0), rad=.25, val=rep('', 60), time=c(as.numeric(format(Sys.time(), '%H')), NA, NA)) plotclock(add=TRUE, fun.plot=lines)
Plot a fan
plotfan(theta, a = 1, ab = 1, orig = c(0, 0), fun = graphics::plot, ...)plotfan(theta, a = 1, ab = 1, orig = c(0, 0), fun = graphics::plot, ...)
theta |
Angle in polar coordinate system |
a |
Radius of start and end points of the arrow. |
ab |
Semi-major over semi-minor. ab=1 for a Ring. |
orig |
Origin |
fun |
Plot function. default = plot |
... |
More options in plot function |
n=50 theta =0:n plotpcs(theta = theta, a=1, ab=1, asp=1) plotfan(theta = theta * 10, a=1, ab=1, fun=polygon)n=50 theta =0:n plotpcs(theta = theta, a=1, ab=1, asp=1) plotfan(theta = theta * 10, a=1, ab=1, fun=polygon)
Plot in polar coordinate system
plotpcs(theta, a = 1, ab = 1, orig = c(0, 0), fun = graphics::plot, ...)plotpcs(theta, a = 1, ab = 1, orig = c(0, 0), fun = graphics::plot, ...)
theta |
Angle in polar coordinate system |
a |
Radius of start and end points of the arrow. |
ab |
Semi-major over semi-minor. ab=1 for a Ring. |
orig |
Origin |
fun |
Plot function. default = plot |
... |
More options in plot function |
n=50 par(mfrow=c(2,1)) plotpcs(theta = 1:n * 15, a=1:n/10, ab=1, type='l', asp=1) plotpcs(theta = 1:n * 10, a=1:n/10, ab=1, type='l', asp=1) xy = PCS2CCS(theta = 1:n * 10, a=1:n/10, ab=1) xy[,1]=xy[,1]+1 points(xy, pch=19, col=terrain.colors(nrow(xy)))n=50 par(mfrow=c(2,1)) plotpcs(theta = 1:n * 15, a=1:n/10, ab=1, type='l', asp=1) plotpcs(theta = 1:n * 10, a=1:n/10, ab=1, type='l', asp=1) xy = PCS2CCS(theta = 1:n * 10, a=1:n/10, ab=1) xy[,1]=xy[,1]+1 points(xy, pch=19, col=terrain.colors(nrow(xy)))
Plot a planet
plotplanet(orig = c(0, 0), rad = 1, theta = 0, fun = graphics::lines, cols = "gray", ab = 1, arrow = TRUE, arrow.len = 0.1, ...)plotplanet(orig = c(0, 0), rad = 1, theta = 0, fun = graphics::lines, cols = "gray", ab = 1, arrow = TRUE, arrow.len = 0.1, ...)
orig |
Origin |
rad |
Radius of the planet |
theta |
Angle of the Arrow inside of the planet |
fun |
Function to plot the planet |
cols |
Color of planet and arrow. |
ab |
Semi-major over semi-minor. ab=1 for the planet |
arrow |
Whether plot the arrow. |
arrow.len |
Length in arrow function. |
... |
More options in plot function. |
a = 10; ab=1.5 x1=PCS2CCS(a=a, ab=ab) c1 = ab2c(a=a, ab=ab) plot(x1, type='l', xlim=c(-10,10), ylim=c(-10,10), asp=1, col='gray') Arrow.pcs(theta = 1:12 * 30, r1=0, r2=a, ab1=ab, length=.1, col=2, o1 = c(c1,0), o2=c(0,0)) pos = PCS2CCS(theta = 1:12 * 30, a=a, ab=ab) plotplanet(orig = pos, arrow.len=0.1)a = 10; ab=1.5 x1=PCS2CCS(a=a, ab=ab) c1 = ab2c(a=a, ab=ab) plot(x1, type='l', xlim=c(-10,10), ylim=c(-10,10), asp=1, col='gray') Arrow.pcs(theta = 1:12 * 30, r1=0, r2=a, ab1=ab, length=.1, col=2, o1 = c(c1,0), o2=c(0,0)) pos = PCS2CCS(theta = 1:12 * 30, a=a, ab=ab) plotplanet(orig = pos, arrow.len=0.1)
Radian to degree
r2d(x)r2d(x)
x |
Radian |
Degree
r = (1:100)/100 * 4 * pi d = r2d(r) rr = d2r(d) plot(d, sin(rr)); abline(h=0 ) abline(v = 360)r = (1:100)/100 * 4 * pi d = r2d(r) rr = d2r(d) plot(d, sin(rr)); abline(h=0 ) abline(v = 360)
Roate a (x,y) with an Angle.
rotate(x, theta, orig = c(0, 0))rotate(x, theta, orig = c(0, 0))
x |
Coordinate (x,y) |
theta |
An angle (degree) |
orig |
Pivot of rotation |
Rotated (x,y)
plot(PCS2CCS(a=2), asp=1) grid() x=rbind(c(1,1)) for(i in 1:12){ points(rotate(x, 30*i, orig=c(1,0))) }plot(PCS2CCS(a=2), asp=1) grid() x=rbind(c(1,1)) for(i in 1:12){ points(rotate(x, 30*i, orig=c(1,0))) }
SpaceObject
Class of planet
SpaceObject
Class of SpaceObject
shapePloting function of the shape
radiusRadius for sphere
Period.Rotatedata.frame 1*3 c(Period.Rotate, Period.Orbit, Period.Synodic)
Orbit
Class of Orbit
Orbit
Class of SpaceOrbit
abShape of the object, ab=1 Sphere, ab!=1 Ellipsoid
eeccentric of the orbit
radiusRadius for sphere (ab=1), or Semi-major axis for Ellipsoid (ab!=1)
perioddata.frame 1*3 c(Period.Rotate, Period.Orbit, Period.Synodic)
InclinationInclination.
CenterObjectCentral Object.
Status.planet
Calculate the status of planet
Status.planet
Status.planet(t, p.orb, ab = 1, r.orb = 1, orig = c(0, 0))Status.planet(t, p.orb, ab = 1, r.orb = 1, orig = c(0, 0))
t |
Time (day). |
p.orb |
Orbital Period. |
ab |
Semi-major over semi-minor. ab=1 for a Ring. |
r.orb |
Radius of the orbit. |
orig |
Reference orgin. |
(x,y) in Cartesian Coordinate System.
tday = seq(0, 365, 30) x=Status.planet(t=tday, p.orb = 365, r.orb=10) plot(PCS2CCS(0:360, a=10), type='l') plotplanet(orig=x[,-1], rad = .51) grid()tday = seq(0, 365, 30) x=Status.planet(t=tday, p.orb = 365, r.orb=10) plot(PCS2CCS(0:360, a=10), type='l') plotplanet(orig=x[,-1], rad = .51) grid()