# ======================================================================== # Red Sun # ======================================================================== # generate pairs of x-y values theta <- seq(0, 2*pi, length = 300) x <- cos(theta) y <- sin(theta) png("red_sun.png", width = 700, height = 400) # set graphical parameters op <- par(bg = "black", mar = rep(0.5, 4)) # open plot plot.new() plot.window(xlim = c(-1, 1), ylim = c(-1, 1), asp = 1) # first solar layer segments(rep(0, 299), rep(0, 299), x[1:299]*runif(299, 0.7), y[1:299]*runif(299,0.7), col = hsv(0.95, 1, 1, runif(299, 0.5)), lwd = 1.5) symbols(0, 0, circles = 0.84, inches = FALSE, fg = "black", bg = "black", add = TRUE) symbols(0, 0, circles = 0.84, inches = FALSE, fg = hsv(0.95, 1, 1, 0.5), bg = hsv(0.95, 1, 1, 0.5), add = TRUE) symbols(0, 0, circles = 0.8, inches = FALSE, fg = hsv(0.95, 1, 1, 0.6), bg = hsv(0.95, 1, 1, 0.6), add = TRUE) symbols(0, 0, circles = 0.76, inches = FALSE, fg = hsv(0.95, 1, 1, 0.7), bg = hsv(0.95, 1, 1, 0.7), add = TRUE) symbols(0, 0, circles = 0.7, inches = FALSE, fg = hsv(0.95, 1, 1, 0.8), bg = hsv(0.95, 1, 1, 0.8), add = TRUE) symbols(0, 0, circles = 0.65, inches = FALSE, fg = hsv(0.95, 1, 1), bg = hsv(0.95, 1, 1), add = TRUE) # second solar layer segments(rep(0, 299), rep(0, 299), x[1:299]*runif(299, 0.7), y[1:299]*runif(299,0.7), col = hsv(0.95, 1, 1, runif(299, 0.5)), lwd = 1.5) symbols(0, 0, circles = 0.74, inches = FALSE, fg = "black", bg = "black", add = TRUE) symbols(0, 0, circles = 0.74, inches = FALSE, fg = hsv(0.95, 1, 1, 0.5), bg = hsv(0.95, 1, 1, 0.5), add = TRUE) symbols(0, 0, circles = 0.7, inches = FALSE, fg = hsv(0.95, 1, 1, 0.6), bg = hsv(0.95, 1, 1, 0.6), add = TRUE) symbols(0, 0, circles = 0.66, inches = FALSE, fg = hsv(0.95, 1, 1, 0.7), bg = hsv(0.95, 1, 1, 0.7), add = TRUE) symbols(0, 0, circles = 0.6, inches = FALSE, fg = hsv(0.95, 1, 1, 0.8), bg = hsv(0.95, 1, 1, 0.8), add = TRUE) symbols(0, 0, circles = 0.55, inches = FALSE, fg = hsv(0.95, 1, 1), bg = hsv(0.95, 1, 1), add = TRUE) # third solar layer segments(rep(0, 299), rep(0, 299), x[1:299]*runif(299, 0.7), y[1:299]*runif(299, 0.7), col = hsv(0.95, 1, 1, runif(299, 0.5)), lwd = 1.5) symbols(0, 0, circles = 0.64, inches = FALSE, fg = "black", bg = "black", add = TRUE) symbols(0, 0, circles = 0.64, inches = FALSE, fg = hsv(.95, 1, 1, 0.5), bg = hsv(0.95, 1, 1, 0.5), add=TRUE) symbols(0, 0, circles = 0.6, inches = FALSE, fg = hsv(0.95, 1, 1, 0.6), bg = hsv(0.95, 1, 1, 0.6), add = TRUE) symbols(0, 0, circles = 0.56, inches = FALSE, fg = hsv(0.95, 1, 1, 0.7), bg = hsv(0.95, 1, 1, 0.7), add = TRUE) symbols(0, 0, circles = 0.5, inches = FALSE, fg = hsv(0.95, 1, 1, 0.8), bg = hsv(0.95, 1, 1, 0.8), add = TRUE) symbols(0, 0, circles = 0.45, inches = FALSE, fg = hsv(0.95, 1, 1), bg = hsv(0.95, 1, 1), add = TRUE) # signature legend("bottomright", legend = "© Gaston Sanchez", bty = "n", text.col = "gray70") # reset par par(op) dev.off()