Package {skymodelr}

Convert altitude/azimuth to ENU

Description

Map altitude/azimuth angles (north-referenced) to a unit East-North-Up direction vector.

Usage

altaz_to_enu(alt_rad, azN_rad)

Arguments

Apply Prague RGB correction gain

Description

Apply Prague RGB correction gain

Usage

apply_prague_rgb_gain(x, gain)

Arguments

Attach skymodelr RGB color metadata

Description

Attach skymodelr RGB color metadata

Usage

as_sky_image(image)

Arguments

Convert B-V color index to blackbody temperature

Description

Convert B-V color index to blackbody temperature

Usage

bv_to_temperature_K(bv)

Arguments

Value

Numeric vector of temperatures in Kelvin.

Sample Prague sky radiance at a chosen wavelength.

Description

Evaluate the Prague spectral sky model at arbitrary spherical directions and return radiance at a user-specified wavelength. Use render_mode = "all" for atmosphere + solar disk, "atmosphere" for atmospheric radiance only, or "sun" for the solar disk only.

Usage

calculate_sky_radiance(
  phi,
  theta,
  lambda_nm,
  altitude = 0,
  elevation = 10,
  visibility = 50,
  albedo = 0.5,
  azimuth = 90,
  number_cores = 1,
  wide_spectrum = FALSE,
  render_mode = "all"
)

Arguments

Value

Numeric vector of radiance values at the requested wavelength(s).

Examples

lambda_vals = calculate_sky_radiance(
  phi = c(90, 90),
  theta = c(45, 45),
  lambda_nm = 550,
  altitude = c(0, 10000),
  elevation = 20,
  visibility = 80,
  albedo = 0.1
)
cbind(altitude = c(0, 10000), radiance = lambda_vals)

Sample a direction from the Prague model.

Description

Evaluate the Prague spectral sky model at arbitrary spherical directions without writing an image, returning radiance-only samples.

Usage

calculate_sky_values(
  phi,
  theta,
  altitude = 0,
  elevation = 10,
  visibility = 50,
  albedo = 0.5,
  azimuth = 90,
  number_cores = 1,
  wide_spectrum = FALSE,
  render_mode = "all",
  prague_rgb_correction = TRUE,
  prague_rgb_correction_strength = 1,
  prague_rgb_correction_gain = "auto"
)

Arguments

Value

3-column RGB matrix.

Examples

# Generate a basic atmosphere with the Prague model
value_grid = expand.grid(
  phi = seq(0, 360, by = 30),
  theta = seq(0, 90, by = 10),
  altitude = c(0, 10000)
)
vals = calculate_sky_values(
  phi = value_grid$phi,
  theta = value_grid$theta,
  altitude = value_grid$altitude,
  elevation = 45,
  visibility = 120,
  albedo = 0
)
cbind(value_grid, vals)

Sample sun luminance at the disk center.

Description

Evaluate the Prague or Hosek sky model at the sun center and integrate against the CIE Y curve to return a luminance value. This is useful for relative attenuation (e.g., comparing zenith sun to a low sun).

Usage

calculate_sun_brightness(
  elevation = 10,
  azimuth = 90,
  albedo = 0.5,
  turbidity = 3,
  altitude = 0,
  visibility = 50,
  hosek = TRUE,
  wide_spectrum = FALSE,
  lambda_nm = NULL
)

Arguments

Value

Numeric scalar of sun luminance (CIE Y, relative scale).

Examples

calculate_sun_brightness(elevation = 45, hosek = TRUE)

Clear cached sky data files

Description

Removes files cached by download_sky_data() under tools::R_user_dir("skymodelr", "data").

Usage

clear_sky_data(files = NULL, ask = interactive())

Arguments

Value

Invisibly, the paths successfully removed.

Examples

clear_sky_data()

Compute effective luminous efficacy for a chosen SPD

Description

Compute effective luminous efficacy for a chosen SPD

Usage

compute_K_eff(
  spd_type = c("D65", "BB5778"),
  wavelength_grid = NULL,
  V_lambda = NULL
)

Arguments

Value

Effective luminous efficacy in lm/W.

Compute a unit RGB vector for a chosen SPD

Description

Compute a unit RGB vector for a chosen SPD

Usage

compute_spd_rgb_unit(spd_type = c("D65", "BB5778"))

Arguments

Value

Numeric vector length 3, normalized so sum equals 1.

Download Prague Sky Model Coefficient Data

Description

This model which allows for sun angles below the horizon, wide spectral ranges including infrared, polarization, as well rendering at altitude) need to download a coefficient file. There are three versions of this dataset, listed here in order of increasing size:

Usage

download_sky_data(sea_level = TRUE, wide_spectrum = FALSE)

Arguments

Details

Value

Invisibly, the full path to the data file.

Examples

# Standard (11-channel, sea-level) coefficients
download_sky_data()

# Wide-spectrum sea-level coefficients
download_sky_data(wide_spectrum = TRUE)

# Full altitude-range coefficients
download_sky_data(sea_level = FALSE)

Keep EXR metadata fields supported by libopenexr

Description

Keep EXR metadata fields supported by libopenexr

Usage

filter_supported_exr_metadata(metadata)

Arguments

Generate the moon into a lat-long image array patch

Description

Produce a rasterised moon texture aligned for composition into the sky dome. The Earth phase is inferred from the Sun-Moon geometry, and earthshine is implemented as an emissive term scaled relative to solar irradiance.

Usage

generate_moon_image_latlong(
  datetime,
  lat,
  lon,
  elev_m = 0,
  width = 800,
  height = 800,
  earthshine = TRUE,
  earthshine_albedo = 0.19,
  solar_irradiance_w_m2 = 1300,
  moon_extinction_kV = 0.172,
  verbose = FALSE
)

Arguments

Generate the atmosphere with the moon

Description

Note that this is just a scaled version of generate_sky(), scaled down by the luminance of the moon as compared to the sun. This function takes the phase of the moon into account, along with the increase in luminosity around a full moon (known as opposition surge). Moonlight attenuation uses the Rozenberg/Krisciunas-Schaefer airmass approximation with configurable moon_extinction_kV.

Usage

generate_moon_latlong(
  datetime,
  lat,
  lon,
  filename = NA,
  albedo = 0.5,
  turbidity = 3,
  altitude = 0,
  resolution = 2048,
  number_cores = 1,
  moon_atmosphere = FALSE,
  earthshine = TRUE,
  earthshine_albedo = 0.19,
  solar_irradiance_w_m2 = 1300,
  moon_extinction_kV = 0.172,
  hosek = TRUE,
  wide_spectrum = FALSE,
  visibility = 50,
  moon_texture_width = 801,
  moon_texture_height = 801,
  verbose = FALSE
)

Arguments

Details

The complete lunar disk is calibrated before directions below the horizontal horizon are clipped. The upper limb remains visible after the center sets, with the finite horizon attenuation and tint retained below zero center elevation. The visible segment is not renormalized to the full Moon's brightness. Finite sky-map resolution limits how finely this edge is resolved.

Disk visibility is independent of the atmospheric model's elevation domain: the optional atmospheric component requires center elevation at least zero for Hosek or -4.2 degrees for Prague. Below that domain it is omitted, while any visible lunar limb is still drawn. Horizon depression and refraction are not modeled.

Value

Either the image array, or the array is invisibly returned if a file is written. The array has dimensions ⁠(resolution, 2 * resolution, 4)⁠.

Note

Writing to non-EXR formats will introduce precision loss because HDR data are quantised to the destination format, and low dynamic range outputs like PNG and JPEG files will not represent the true luminosity values encoded in the array.

Examples

# Moonlit sky (Hosek), mid-evening in DC
generate_moon_latlong(
  datetime = as.POSIXct("2025-09-05 19:30:00", tz = "America/New_York"),
  lat = 38.9072,
  lon = -77.0369,
  resolution = 400,
  turbidity = 3,
  verbose = TRUE
) |>
  rayimage::render_exposure(15) |>
  rayimage::plot_image()

Generate a bright-planet image array

Description

Build a planetary luminance map aligned with the sky dome for compositing within generate_sky_latlong().

Usage

generate_planets(
  datetime,
  lon,
  lat,
  filename = NA,
  resolution = 2048,
  turbidity = 3,
  ozone_du = 300,
  altitude = 0,
  color = FALSE,
  planet_width = 1,
  upper_hemisphere_only = TRUE,
  atmosphere_effects = TRUE,
  number_cores = 1,
  verbose = FALSE
)

Arguments

Value

Either the image array, or the array is invisibly returned if a file is written. The array has dimensions ⁠(resolution, 2 * resolution, 4)⁠.

Note

Writing to non-EXR formats will introduce precision loss because HDR data are quantised to the destination format, and low dynamic range outputs like PNG and JPEG files will not represent the true luminosity values encoded in the array.

Examples

# Basic star field over Washington, DC at a fixed time
generate_planets(
  datetime = as.POSIXct("2025-03-21 02:20:00", tz = "America/New_York"),
  lon = -77.0369,
  lat = 38.9072,
  resolution = 400,
  color = TRUE,
  planet_width = 1,
  atmosphere_effects = TRUE,
  upper_hemisphere_only = TRUE,
  number_cores = 2
) |>
  rayimage::plot_image()

Generate a Hosek-Wilkie sky dome array

Description

Evaluate either the Hosek-Wilkie or Prague analytic sky models and return a high-dynamic-range image array for the given solar configuration. An image file is written only when filename is supplied.

Usage

generate_sky(
  filename = NA,
  albedo = 0.1,
  turbidity = 3,
  elevation = 10,
  azimuth = 90,
  altitude = 0,
  resolution = 2048,
  number_cores = 1,
  hosek = TRUE,
  wide_spectrum = FALSE,
  visibility = 50,
  verbose = FALSE,
  render_mode = "all",
  below_horizon = TRUE,
  prague_rgb_correction = TRUE,
  prague_rgb_correction_strength = 1,
  prague_rgb_correction_gain = "auto",
  exr_adopted_white = "D60",
  exr_metadata = TRUE
)

Arguments

Details

For Prague RGB output, prague_rgb_correction = TRUE applies a fixed linear RGB gain to reduce the small magenta / negative-green tint observed in the RGB projection of Prague sky maps. Set prague_rgb_correction = FALSE to recover raw Prague RGB output. The correction is not a white balance and is not applied to wavelength-specific spectral radiance.

EXR metadata: Generated sky maps are linear RGB. When exr_metadata = TRUE, skymodelr tags the image with sRGB/Rec.709 chromaticities and an EXR adoptedNeutral white. The default adopted neutral is D60. This metadata does not convert pixel values; it informs downstream readers how to treat the scene neutral. Use exr_adopted_white = "D65" for D65 adopted neutral, or exr_metadata = FALSE to omit skymodelr EXR metadata. The RGB encoding remains linear sRGB / Rec.709 even when the adopted neutral is D60.

Value

Either the image array, or the array is invisibly returned if a file is written. The array has dimensions ⁠(resolution, 2 * resolution, 4)⁠.

Note

Writing to non-EXR formats will introduce precision loss because HDR data are quantised to the destination format, and low dynamic range outputs like PNG and JPEG files will not represent the true luminosity values encoded in the array.

Examples

sky = generate_sky(
  resolution = 8,
  elevation = 30,
  azimuth = 135,
  render_mode = "atmosphere"
)
dim(sky)


# Hosek model (default): clear morning, Sun SE, with solar disk
generate_sky(
  resolution = 400,
  elevation  = 15,
  azimuth    = 135,
  turbidity  = 3,
  render_mode = "all"
) |>
  rayimage::plot_image()

# Same view but hazier and without the solar disk
generate_sky(
  resolution = 400,
  elevation  = 15,
  azimuth    = 135,
  turbidity  = 6,
  render_mode = "atmosphere"
) |>
  rayimage::plot_image()

# Prague model (requires downloaded coefficients)
generate_sky(
  resolution = 400,
  hosek      = FALSE,
  altitude   = 0,
  visibility = 80,
  albedo     = 0.2,
  elevation  = 5,
  azimuth    = 220,
  number_cores = 2
) |>
  rayimage::plot_image()

generate_sky(
  "raw_prague.exr",
  elevation = 60,
  azimuth = 315,
  hosek = FALSE,
  prague_rgb_correction = FALSE
)

Generate a location and time-specific sky dome (optionally with stars)

Description

Convenience wrapper around generate_sky() that:

1. Computes the Sun's apparent position for datetime, lat, and lon (via Swiss Ephemeris / swephR).

2. Renders the corresponding sky model.

3. Optionally overlays a star field using generate_stars() or moon with generate_moon_latlong().

Usage

generate_sky_latlong(
  datetime,
  lat,
  lon,
  filename = NA,
  albedo = 0.5,
  turbidity = 3,
  altitude = 0,
  resolution = 2048,
  number_cores = 1,
  hosek = TRUE,
  wide_spectrum = FALSE,
  visibility = 50,
  stars = FALSE,
  star_width = 1,
  planets = FALSE,
  moon = FALSE,
  moon_atmosphere = FALSE,
  moon_hosek = TRUE,
  render_mode = "all",
  below_horizon = TRUE,
  prague_rgb_correction = TRUE,
  prague_rgb_correction_strength = 1,
  prague_rgb_correction_gain = "auto",
  exr_adopted_white = "D60",
  exr_metadata = TRUE,
  verbose = FALSE,
  stars_exposure = 0,
  ...
)

Arguments

Details

Solar angles - altitude (degrees above the horizon) and azimuth (degrees clockwise from east, so 90 degrees = south) - are derived internally; you never have to supply them directly.

Black-sky rule - With the Prague model the sky radiance is defined only down to -4.2 degrees, and with the Hosek model it is defined only about 0 degrees. Below that the function skips the sky render and writes only stars when stars = TRUE.

Value

Either the raw data, or the data is invisibly returned if filename is given. The EXR is written to filename.

Examples

sky = generate_sky_latlong(
  datetime = as.POSIXct("2025-03-21 12:00:00", tz = "America/New_York"),
  lat = 38.9072,
  lon = -77.0369,
  resolution = 8
)
dim(sky)


# Morning sunrise on spring solstice over Washington, DC with Prague model
generate_sky_latlong(
  datetime = as.POSIXct("2025-03-21 06:15:00", tz = "America/New_York"),
  lat = 38.9072,
  lon = -77.0369,
  number_cores = 2,
  hosek = FALSE
) |>
  rayimage::plot_image()

generate_sky_latlong(
  datetime = as.POSIXct("2025-03-21 12:00:00", tz = "America/New_York"),
  lat = 38.9072,
  lon = -77.0369,
  number_cores = 2
) |>
  rayimage::plot_image()

generate_sky_latlong(
  datetime = as.POSIXct("2025-03-21 18:00:00", tz = "America/New_York"),
  lat = 38.9072,
  lon = -77.0369,
  number_cores = 2
) |>
  rayimage::plot_image()

generate_sky_latlong(
  datetime = as.POSIXct("2025-03-21 18:30:00", tz = "America/New_York"),
  lat = 38.9072,
  lon = -77.0369,
  number_cores = 2,
  hosek = FALSE,
  verbose = TRUE
) |>
  rayimage::render_exposure(exposure = 2) |>
  rayimage::plot_image()

Generate a star-field array aligned with generate_sky()

Description

Render a star map for a given observer location, time, and atmospheric conditions so it can be composited with generate_sky(). Returns a ⁠(resolution, 2 * resolution, 4)⁠ array with an opaque alpha channel. An image file is written only when filename is supplied.

Usage

generate_stars(
  lon,
  lat,
  datetime,
  filename = NA,
  resolution = 2048,
  turbidity = 3,
  ozone_du = 300,
  altitude = 0,
  color = TRUE,
  star_width = 1,
  upper_hemisphere_only = TRUE,
  atmosphere_effects = TRUE,
  number_cores = 1
)

Arguments

Value

Either the image array, or the array is invisibly returned if a file is written. The array has dimensions ⁠(resolution, 2 * resolution, 4)⁠.

Note

Writing to non-EXR formats will introduce precision loss because HDR data are quantised to the destination format, and low dynamic range outputs like PNG and JPEG files will not represent the true luminosity values encoded in the array.

Examples

star_map = generate_stars(
  resolution = 8,
  lon = -77.0369,
  lat = 38.9072,
  datetime = as.POSIXct("2025-03-21 02:20:00", tz = "UTC"),
  atmosphere_effects = FALSE,
  number_cores = 1
)
dim(star_map)


# Note: exposure has been increased for all examples (via white_point) for
# ease of visibility in documentation

# Basic star field over Washington, DC at a fixed time
generate_stars(
  resolution = 400,
  lon = -77.0369,
  lat = 38.9072,
  datetime = as.POSIXct("2025-03-21 02:20:00", tz = "America/New_York"),
  color = TRUE,
  star_width = 1,
  atmosphere_effects = TRUE,
  upper_hemisphere_only = TRUE,
  number_cores = 2
) |>
  rayimage::plot_image()

# Monochrome stars, no atmospheric extinction/reddening, full sphere
generate_stars(
  resolution = 400,
  lon = -122.4194,
  lat = 37.7749,
  datetime = as.POSIXct("2025-06-01 08:00:00", tz = "UTC"),
  color = FALSE,
  star_width = 1,
  upper_hemisphere_only = FALSE,
  atmosphere_effects = FALSE
) |>
  rayimage::plot_image()

# Sharper stars (smaller PSF) with ozone/turbidity and altitude
generate_stars(
  resolution = 400,
  lon = 10,
  lat = 45,
  datetime = as.POSIXct("2025-12-01 22:00:00", tz = "UTC"),
  star_width = 0.5,
  turbidity = 3.5,
  ozone_du = 320,
  altitude = 1000,
  color = TRUE
) |>
  rayimage::plot_image()

Generate Sun and Moon Disk Radiance Images

Description

Generate a celestial disk directly at texture resolution, independently of an environment map. The Sun uses the Prague solar radiance profile. The Moon uses the existing lunar surface texture, topocentric phase orientation, earthshine, photometry, and atmospheric extinction.

Usage

generate_sun_disk(
  datetime,
  lat,
  lon,
  altitude = 0,
  resolution = 256,
  visibility = 50,
  albedo = 0.5,
  number_cores = 1,
  wide_spectrum = FALSE,
  prague_rgb_correction = TRUE,
  prague_rgb_correction_strength = 1,
  prague_rgb_correction_gain = "auto",
  atmospheric_attenuation = TRUE
)

generate_moon_disk(
  datetime,
  lat,
  lon,
  altitude = 0,
  resolution = 1024,
  visibility = 50,
  albedo = 0.5,
  number_cores = 1,
  wide_spectrum = FALSE,
  prague_rgb_correction = TRUE,
  prague_rgb_correction_strength = 1,
  prague_rgb_correction_gain = "auto",
  earthshine = TRUE,
  earthshine_albedo = 0.19,
  solar_irradiance_w_m2 = 1300,
  moon_extinction_kV = 0.172,
  atmospheric_attenuation = TRUE
)

Arguments

Details

Prague data must be installed with download_sky_data(). The atmosphere is evaluated for one observer altitude and time when atmospheric_attenuation is TRUE. With FALSE, solar radiance comes from the model's intrinsic solar spectrum and lunar radiance uses unattenuated photometry without the atmospheric tint. Solar texture values are then independent of altitude and solar elevation; ephemeris placement and lunar phase still use the requested location and time. Neither mode clips the texture to a geometric horizon. Pair these images with an atmosphere-only sky and disable the corresponding rasterized celestial bodies in that sky to avoid counting their light twice.

The disk is inscribed in the image rectangle. For column j, row i, width w, and height h, pixel centers have coordinates x = 2 * (j - 0.5) / w - 1 and y = 1 - 2 * (i - 0.5) / h. Sample only within x^2 + y^2 <= 1. The corresponding ray direction is normalize(forward + tan(radius) * (x * right + y * up)), where radius = angular_diameter_deg * pi / 360. Image-up is the projection of local vertical into the disk plane; image-right points toward increasing azimuth. At the zenith or nadir, north replaces the vertical reference axis. There is no image reprojection or implicit rotation in the return value.

The solar model has a fixed angular profile, mapped here to the apparent ephemeris diameter without changing radiance. Lunar radiance is normalized over disk solid angle to the existing phase-dependent irradiance, with the solar spectral RGB distribution. With atmospheric_attenuation = TRUE, lunar attenuation and Prague tint are evaluated at the disk center. At or below zero center elevation, their finite horizon values are retained so the upper limb does not disappear prematurely. This continuation preserves a full disk texture; it is not a model of below-horizon visibility or horizon depression at altitude. Geometric horizon clipping is the consumer's responsibility; sample the full disk and discard directions below the horizon when appropriate.

To save an EXR without changing radiance, tag image with rayimage::ray_read_image(image, normalize = FALSE, source_linear = TRUE, assume_colorspace = rayimage::CS_SRGB) and write with clamp = FALSE.

Value

A list containing:

• image: a numeric array with dimensions c(height, width, 3) containing linear sRGB radiance, with no exposure, tone mapping, or alpha channel. Antialiased lunar coverage is already applied: do not multiply by alpha. Out-of-gamut solar RGB values may be negative and should be preserved.

• azimuth_deg: disk-center azimuth, clockwise from true north, in [0, 360).

• elevation_deg: disk-center elevation above the local horizon, in degrees.

• angular_diameter_deg: topocentric apparent angular diameter in degrees.

• atmospheric_attenuation: whether atmospheric filtering is included.

• projection: "rectilinear", using the disk mapping described below.

Examples

time = as.POSIXct("2026-01-28 21:00:00", tz = "Pacific/Auckland")
moon = generate_moon_disk(time, lat = -36.87593, lon = 174.7647)
moon[c("azimuth_deg", "elevation_deg", "angular_diameter_deg")]
# Let a renderer evaluate Earth's atmosphere at each interaction instead.
intrinsic_moon = generate_moon_disk(time, lat = -36.87593, lon = 174.7647,
  atmospheric_attenuation = FALSE)

Prepare Location and Time Metadata for Native Prague Sky Queries

Description

Resolve the full-altitude Prague dataset and the Sun's topocentric position without loading coefficients or generating an image. This lets rendering libraries evaluate the model natively while using skymodelr's ephemeris and RGB calibration through a public interface.

Usage

get_prague_sky_metadata(
  datetime,
  lat,
  lon,
  altitude = 0,
  visibility = 50,
  albedo = 0.5,
  prague_rgb_correction = TRUE,
  prague_rgb_correction_strength = 1,
  prague_rgb_correction_gain = "auto"
)

Arguments

Details

Install the data explicitly with download_sky_data(sea_level = FALSE). This function never downloads files or prompts. It supports the visible-spectrum, full-altitude dataset and returns the actual Sun position even below the Prague model's minimum elevation of -4.2 degrees. Native consumers should use zero solar radiance below that limit while retaining transmission for other light sources. This function does not evaluate radiance, attenuation, or atmospheric refraction.

Value

A list with filename (the installed full-altitude coefficient file), elevation_deg, azimuth_deg (clockwise from north), angular_diameter_deg, rgb_gain (the applied R, G, B gains), and the reference altitude, visibility, and albedo. No native pointers are included; the result can be serialized. The coefficient path is local to this computer and should be resolved again on another machine.

Examples

info = get_prague_sky_metadata(
  as.POSIXct("2026-06-21 20:00:00", tz = "America/New_York"),
  lat = 40.7, lon = -74
)
info[c("elevation_deg", "azimuth_deg", "rgb_gain")]

Get skymodelr EXR adopted white

Description

Get skymodelr EXR adopted white

Usage

get_skymodelr_adopted_white(exr_adopted_white)

Arguments

Convert POSIXct timestamps to Julian Date

Description

Transform a POSIXct time value into an astronomical Julian Date (days since noon UTC 1 Jan 4713 BCE).

Usage

jd_utc(time_utc)

Arguments

List cached sky data files

Description

Lists files cached by download_sky_data() under tools::R_user_dir("skymodelr", "data").

Usage

list_sky_data()

Value

A data frame with columns file, path, size, and modified.

Examples

list_sky_data()

Compute luminous efficacy for a blackbody SPD

Description

Compute luminous efficacy for a blackbody SPD

Usage

luminous_efficacy_blackbody(
  temperature_K,
  wavelengths_nm = NULL,
  V_lambda = NULL
)

Arguments

Value

Numeric vector of K_eff values in lm/W.

Convert photometric illuminance to radiometric irradiance

Description

Convert photometric illuminance to radiometric irradiance

Usage

lux_to_radiometric_irradiance(E_v_lux, K_eff)

Arguments

Value

Irradiance in W/m^2.

Moon V-band magnitude from phase angle

Description

Moon V-band magnitude from phase angle

Usage

moon_mag(phase_deg, dist_km = 384400, mean_km = 384400)

Arguments

Value

Apparent V magnitude m.

Moon disk solid angle in steradians

Description

Moon disk solid angle in steradians

Usage

moon_solid_angle_sr(diameter_deg = 0.533)

Arguments

Value

Solid angle in steradians.

Normalize a vector to unit length

Description

Scale a numeric vector so it has unit Euclidean length.

Usage

normalize(v)

Arguments

Normalize Prague RGB correction option

Description

Normalize Prague RGB correction option

Usage

normalize_prague_rgb_correction(prague_rgb_correction)

Arguments

Prepare Prague RGB correction gain

Description

Prepare Prague RGB correction gain

Usage

prepare_prague_rgb_gain(gain = "auto", strength = 1)

Arguments

Build EXR metadata for skymodelr RGB output

Description

Build EXR metadata for skymodelr RGB output

Usage

sky_exr_metadata(image = NULL)

Arguments

Compute radiometric star amplitude and RGB unit vector

Description

Compute radiometric star amplitude and RGB unit vector

Usage

star_radiometric_amplitude(mV, bv, wavelengths_nm = NULL, V_lambda = NULL)

Arguments

Value

List with E_v (lux), E_e (W/m^2), K_eff (lm/W), and rgb_unit.

Star catalog used for sky rendering

Description

A data frame of stellar positions and photometric/color information used by generate_stars() to render star fields.

Usage

data(stars)

Format

A data.frame with 9,110 rows and 8 variables:

bsc_number

Numeric identifier from the source catalog.

ra_rad

Right ascension in radians.

dec_rad

Declination in radians.

v_mag

Apparent visual magnitude (V band).

spec

Spectral type string.

r

Relative red channel weight derived from spectral type.

g

Relative green channel weight derived from spectral type.

b

Relative blue channel weight derived from spectral type.

Source

Derived from the Bright Star Catalogue, 5th Revised Edition (BSC5), also known as the Yale Bright Star Catalogue. The source catalogue is commonly identified as Hoffleit, D. and Warren, W. H. Jr. (1991), "The Bright Star Catalogue, 5th Revised Ed.", Yale University Observatory, and is distributed as machine-readable catalogue V/50 by the Centre de Donnees astronomiques de Strasbourg (CDS): https://cdsarc.cds.unistra.fr/viz-bin/cat/V/50.

Compute topocentric moon/sun directions

Description

Query Swiss Ephemeris for topocentric sun and moon direction vectors and magnitudes for a given observer location.

Usage

swe_dirs_topo_moon_sun(
  datetime,
  lat,
  lon,
  elev_m = 0,
  moon_extinction_kV = 0.172
)

Arguments

Compute bright-planet positions via Swiss Ephemeris

Description

Produce a data frame of planetary apparent positions and magnitudes for the observer.

Usage

swe_dirs_topo_planets_df(datetime, lat, lon, elev_m = 0)

Arguments

Compute sun/moon positions via Swiss Ephemeris

Description

Produce a data frame of apparent positions and magnitudes for the observer.

Usage

swe_dirs_topo_sunmoon_df(datetime, lat, lon, elev_m = 0)

Arguments

Convert Swiss ephemeris data to a starfield frame

Description

Build a data frame compatible with make_starfield_rcpp() from Swiss Ephemeris results.

Usage

sweph_info_to_stars_df(sweph_info)

Arguments

Tag generated sky EXR metadata

Description

Tag generated sky EXR metadata

Usage

tag_generated_sky_exr_metadata(sky, hosek, exr_metadata, exr_adopted_white)

Arguments

Tag skymodelr EXR metadata

Description

Tag skymodelr EXR metadata

Usage

tag_skymodelr_exr_metadata(
  sky,
  adopted_white_xy = .skymodelr_d60_xy,
  adopted_white_name = .skymodelr_exr_adopted_neutral_label,
  rgb_colorspace = rayimage::CS_SRGB,
  model_name = NULL,
  prague_rgb_correction = NULL
)

Arguments

Validate Prague RGB correction gain

Description

Validate Prague RGB correction gain

Usage

validate_prague_rgb_gain(gain)

Arguments

Convert V magnitude to illuminance in lux

Description

Convert V magnitude to illuminance in lux

Usage

vmag_to_ev_lux(mV, zero_point = 1)

Arguments

Value

Numeric vector of photometric illuminance in lux.

Write a skymodelr image with EXR color metadata

Description

Write a skymodelr image with EXR color metadata

Usage

write_sky_image(image, filename, ...)

Arguments

Convert xy chromaticity to XYZ with Y equal to 1

Description

Convert xy chromaticity to XYZ with Y equal to 1

Usage

xy_to_xyz_y1(xy)

Arguments