pop.index.plot: Plot of population index

In MetFns: Analysis of Visual Meteor Data

Description

Plots population index of a meteor shower for data containing population index values and given limits on x and y axes.

Usage

pop.index.plot(rdata,xlim1,xlim2,xinc,ylim1,ylim2,yinc,dlim1=NULL,dlim2=NULL,
dinc=NULL,dunit=NULL)

Arguments

rdata	data frame consisting of population index values.
xlim1	numeric vector specifying minimum value on x-axis.
xlim2	numeric vector specifying maximum value on x-axis.
xinc	numeric vector specifying increment between labels on x-axis.
ylim1	numeric vector specifying minimum value on y-axis.
ylim2	numeric vector specifying maximum value on y-axis.
yinc	numeric vector specifying increment between labels on y-axis.
dlim1	character vector specifying start date on secondary x-axis, given in "%Y-%m-%d %H:%M:%S" format (UTC timezone).
dlim2	character vector specifying end date on secondary x-axis, given in "%Y-%m-%d %H:%M:%S" format (UTC timezone).
dinc	numeric vector specifying increment between labels on secondary x-axis.
dunit	character vector specifying unit of increment dinc. Accepted values are "min" for minute, "h" for hour and "day" for day.

Details

rdata represents table of population index values calculated using functions pop.index or pop.index2.

Value

xy plot of population index is made. Solar longitude is on x-axis and population index on y-axis. Population index is represented with black filled circles and 68% confidence intervals. Values of limits (minimum and maximum values) on x and y axis (xlim1,xlim2,ylim1,ylim2), as well as increments between the axis labels (xinc,yinc) should be provided to function call.

Author(s)

Kristina Veljkovic

References

Arlt, R. (2003). Bulletin 19 of the International Leonid Watch: Population index study of the 2002 Leonid meteors. WGN, Journal of the IMO,31:3, 77-87.

See Also

pop.index,pop.index2

Examples

## calculate and make plot of population index for observations of 2015 Orionids,
## time period around maximum, 21th to 25th October, 
## min bin size 0.16 degrees (4 hours), max bin size 1 degree, 
## number of meteors equals 100 
## First select magnitude data for Orionids activity from 21th to 25th October, 
## limiting magnitudes above 5, radiant elevation above 20 degrees, 
## percentage of clouds below 20
ori2015mag<-filter(magn2015,date.start="2015-10-21",date.end="2015-10-25",
shw="ORI",P.up=20,mag.low=5,h.low=20)
oripop<-pop.index2(ori2015mag,date.start="2015-10-21",date.end="2015-10-25",
shw="ORI",kmin=0.16,kmax=1,num=100)
pop.index.plot(oripop,xlim1=207,xlim2=211,xinc=1,ylim1=1.6,ylim2=2.4,yinc=0.2)

## change date tick marks on secondary x-axis, to show date at midnight
pop.index.plot(oripop,xlim1=207,xlim2=211,xinc=1,ylim1=1.5,ylim2=2.5,
yinc=0.2,dlim1="2015-10-21",dlim2="2015-10-25",dinc=1,dunit="day")

MetFns documentation built on May 1, 2019, 6:28 p.m.