Description Details Author(s) See Also Examples
Provides easy to use functions to create all-sky grid plots of widely used astronomical coordinate systems (equatorial, ecliptic, galactic) and scatter plots of data on any of these systems including on-the-fly coordinate system conversion. It supports any type of spherical projection defined by the mapproj package.
Package: | SPADAR |
Type: | Package |
Version: | 1.0 |
Date: | 2017-04-29 |
License: | GPL (>= 3) |
The SPADAR package provides simple functions to create spherical projection plots of astronomical data. It enables the creation of coordinate system grids (as backgrounds or overlays) in equatorial, ecliptic and galactic coordinates and also supports the conversion between these coordinate systems (J2000). The plots can adopt multiple spherical projections to the plane (aitoff, mollweide, mercator using the default projection parameters are tested). Coordinate system conversions are performed by a modified version of the euler function from the astrolibR package, and plane projections by the mapproj package.
This package comes with a small dataset containing Ultra High Energy Cosmic Ray events (from the Pierre Auger Observatory 2014 event list), for example and test purposes.
Alberto Krone-Martins
Maintainer: Alberto Krone-Martins <algol@sim.ul.pt>
createAllSkyScatterPlotChart
,
createAllSkyGridChart
,
overplotScatterPlotInAllSkyGridChart
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | # Creates an all sky chart with Equatorial, Ecliptic and Galactic coordinate systems
# superposed. Equatorial as the main system in black, the Ecliptic in blue and the
# Galactic in red
createAllSkyGridChart(mainGrid="equatorial", eqDraw=TRUE, eclDraw=TRUE, galDraq=TRUE,
eqCol="black", eclCol=rgb(0,0,1,0.5), galCol=rgb(1,0,0,0.5), eqLty=1, eclLty=2, galLty=2)
# Creates an all sky scatter plot of the Pierre Auger Observatory 2014 public event list
# (arXiv:1411.6111). Equatorial and Galactic coordinate systems superposed;
# Galactic as the main system in black and the Equatorial in red.
data(uhecrauger2014)
createAllSkyScatterPlotChart(uhecrauger2014$RA, uhecrauger2014$DEC, mainGrid="galactic",
dataCoordSys="equatorial", pointcol="black", pch=19,
pointsize=uhecrauger2014$Eev/max(uhecrauger2014$Eev),
eqDraw=TRUE, eclDraw=FALSE, galDraq=TRUE, galCol="black", eqLty=2, galLty=1,
eqCol=rgb(1,0,0,0.5))
# Creates a grid with the Galactic Coordinate system and overplots data
createAllSkyGridChart(mainGrid="galactic", eqDraw=FALSE, eclDraw=FALSE, galDraq=TRUE,
galCol=rgb(0,0,0,0.5), galLty=1)
overplotScatterPlotInAllSkyGridChart(uhecrauger2014$RA, uhecrauger2014$DEC,
mainGrid="galactic", dataCoordSys="equatorial", pointcol="black", pch=19,
pointsize=uhecrauger2014$Eev/max(uhecrauger2014$Eev))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.