StateTimeSeries
and StateVector
¶>>> from gwpy.timeseries import (StateTimeSeries, StateVector)
A large quantity of important data from gravitational-wave detectors can be distilled into simple boolean (True
or False
) statements informing something about the state of the instrument at a given time.
These statements can be used to identify times during which a particular control system was active, or when the signal in a seismometer was above an alarming threshold, for example.
In GWpy, these data are represented by special cases (sub-classes
) of the TimeSeries
object:
StateTimeSeries |
Boolean array representing a good/bad state determination of some data. |
StateVector |
Binary array representing good/bad state determinations of some data. |
The example of a threshold on signal time-series is the core of a large amount of low-level data quality information, used in searches for gravitational waves, and detector characterisation, and is described by the StateTimeSeries
object, a specific type of TimeSeries
containing only boolean values.
These arrays can be generated from simple arrays of booleans, as follows:
>>> from gwpy.timeseries import StateTimeSeries
>>> state = StateTimeSeries([True, True, False, False, False, True, False],
sample_rate=1, epoch=1064534416)
>>> state
<StateTimeSeries([ True, True, False, False, False, True, False], dtype=bool
name=None
unit=Unit(dimensionless)
epoch=<Time object: scale='tai' format='gps' value=1064534416.0>
channel=None
sample_rate=<Quantity 1 Hz>)>
Alternatively, applying a standard mathematical comparison to a regular TimeSeries
will return a StateTimeSeries
:
>>> from gwpy.timeseries import TimeSeries
>>> seisdata = TimeSeries.fetch('L1:HPI-BS_BLRMS_Z_3_10', 1064534416, 1064538016)
>>> seisdata.unit = 'nm/s'
>>> highseismic = seisdata > 400
>>> highseismic
<StateTimeSeries([False, False, False, ..., False, False, False], dtype=bool
name='L1:HPI-BS_BLRMS_Z_3_10 > 400 nm / s'
unit=Unit("nm / s")
epoch=<Time object: scale='tai' format='gps' value=1064534416.0>
channel=Channel("L1:HPI-BS_BLRMS_Z_3_10")
sample_rate=<Quantity 16.0 Hz>)>
The StateTimeSeries
includes a handy StateTimeSeries.to_dqflag()
method to convert the boolean array into a DataQualityFlag
:
>>> segments = highseismic.to_dqflag(round=True)
>>> segments
<DataQualityFlag(valid=[[1064534416 ... 1064538016)],
active=[[1064535295 ... 1064535296)
[1064535896 ... 1064535897)
[1064536969 ... 1064536970)
[1064537086 ... 1064537088)
[1064537528 ... 1064537529)],
ifo=None,
name=None,
version=None,
comment='L1:HPI-BS_BLRMS_Z_3_10 > 400 nm / s')>
Which the StateTimeSeries
represents a single True
/False
statement about the state of a system, the StateVector
gives a grouping of these, with a binary bitmask mapping bits in a binary word to descriptions of multiple states in a given compound system.
For example, the state of the full laser interferometer was described in Initial LIGO by a combination of separate states, including:
- operator set to go to ‘science mode’
- EPICS control system record (conlog) OK
- instrument locked in all resonant cavities
- no signal injections
- no unauthorised excitations
Additionall, the higher bits 5-15 were set ‘ON’ at all times, such that the word 0xffff
indicated ‘science mode’ operation of the instrument.
This StateVector
can be read from the S6 frames as:
>>> from gwpy.timeseries import StateVector
>>> state = StateVector.fetch('H1:IFO-SV_STATE_VECTOR', 968631675, 968632275,
['science', 'conlog', 'up', 'no injection', 'no excitation'])
>>> print(state)
StateVector([65528 65528 65528 ..., 65535 65535 65535],
name=H1:IFO-SV_STATE_VECTOR,
unit=,
epoch=968631675.0,
channel=H1:IFO-SV_STATE_VECTOR,
sample_rate=16.0 Hz,
bitmask=BitMask(0: science
1: conlog
2: up
3: no injection
4: no excitation,
channel=H1:IFO-SV_STATE_VECTOR,
epoch=968631675.0))
As can be seen, the input bitmask (['science', 'conlog', 'up', 'no injection', 'no excitation']
) is represented through the BitMask
class, recording the bits as a list with some metdata about their purpose.
The StateVector
fetched in the above example can then be parsed into a series of DataQualityFlag
objects, recording the active segments for that bit in the vector:
>>> flags = state.to_dqflags(round=True)
>>> print(flags[0])
<DataQualityFlag(valid=[[968631675 ... 968632275)],
active=[[968632248 ... 968632275)],
ifo=None,
name='science',
version=None,
comment='H1:IFO-SV_STATE_VECTOR bit 0')>
This reference contains the following Class
entries:
StateVector |
Binary array representing good/bad state determinations of some data. |
StateTimeSeries |
Boolean array representing a good/bad state determination of some data. |
StateVectorDict |
gwpy.timeseries.
StateVector
[source]¶Bases: gwpy.timeseries.core.TimeSeriesBase
Binary array representing good/bad state determinations of some data.
Each binary bit represents a single boolean condition, with the
definitions of all the bits stored in the StateVector.bits
attribute.
Parameters: | value : array-like
unit :
epoch :
bits :
sample_rate :
times :
name :
dtype :
copy :
subok :
|
---|
Notes
Key methods:
fetch (channel, start, end[, bits, host, ...]) |
Fetch data from NDS into a StateVector . |
read (*args, **kwargs) |
Read data into a StateVector |
write (data, *args, **kwargs) |
Write this StateVector to a file |
to_dqflags ([bits, minlen, dtype, round]) |
Convert this StateVector into a DataQualityDict |
plot ([format, bits]) |
Plot the data for this StateVector |
Attributes Summary
T |
Same as self.transpose(), except that self is returned if self.ndim < 2. |
base |
Base object if memory is from some other object. |
bits |
The list of bit names for this StateVector . |
boolean |
A mapping of this StateVector to a 2-D array containing all binary bits as booleans, for each time point. |
cgs |
Returns a copy of the current Quantity instance with CGS units. |
channel |
Instrumental channel associated with these data |
ctypes |
An object to simplify the interaction of the array with the ctypes module. |
data |
|
dt |
X-axis sample separation |
dtype |
Data-type of the array’s elements. |
duration |
Duration of this series in seconds |
dx |
X-axis sample separation |
epoch |
GPS epoch for these data. |
equivalencies |
A list of equivalencies that will be applied by default during unit conversions. |
flags |
Information about the memory layout of the array. |
flat |
A 1-D iterator over the Quantity array. |
imag |
The imaginary part of the array. |
info |
Container for meta information like name, description, format. |
isscalar |
True if the value of this quantity is a scalar, or False if it is an array-like object. |
itemsize |
Length of one array element in bytes. |
name |
Name for this data set |
nbytes |
Total bytes consumed by the elements of the array. |
ndim |
Number of array dimensions. |
real |
The real part of the array. |
sample_rate |
Data rate for this TimeSeries in samples per second (Hertz). |
shape |
Tuple of array dimensions. |
si |
Returns a copy of the current Quantity instance with SI units. |
size |
Number of elements in the array. |
span |
X-axis [low, high) segment encompassed by these data |
strides |
Tuple of bytes to step in each dimension when traversing an array. |
times |
Series of GPS times for each sample |
unit |
The physical unit of these data |
value |
The numerical value of this quantity. |
x0 |
X-axis value of the first data point |
xindex |
Positions of the data on the x-axis |
xspan |
X-axis [low, high) segment encompassed by these data |
xunit |
Unit of x-axis index |
Methods Summary
all ([axis, out, keepdims]) |
Returns True if all elements evaluate to True. |
any ([axis, out, keepdims]) |
Returns True if any of the elements of a evaluate to True. |
append (other[, gap, inplace, pad, resize]) |
Connect another series onto the end of the current one. |
argmax ([axis, out]) |
Return indices of the maximum values along the given axis. |
argmin ([axis, out]) |
Return indices of the minimum values along the given axis of a . |
argpartition (kth[, axis, kind, order]) |
Returns the indices that would partition this array. |
argsort ([axis, kind, order]) |
Returns the indices that would sort this array. |
astype (dtype[, order, casting, subok, copy]) |
Copy of the array, cast to a specified type. |
byteswap (inplace) |
Swap the bytes of the array elements |
choose (choices[, out, mode]) |
Use an index array to construct a new array from a set of choices. |
clip ([min, max, out]) |
Return an array whose values are limited to [min, max] . |
compress (condition[, axis, out]) |
Return selected slices of this array along given axis. |
conj () |
Complex-conjugate all elements. |
conjugate () |
Return the complex conjugate, element-wise. |
copy ([order]) |
Return a copy of the array. |
copy_metadata () |
Return a deepcopy of the metadata for this array |
crop ([start, end, copy]) |
Crop this series to the given x-axis extent. |
cumprod ([axis, dtype, out]) |
Return the cumulative product of the elements along the given axis. |
cumsum ([axis, dtype, out]) |
Return the cumulative sum of the elements along the given axis. |
decompose ([bases]) |
Generates a new Quantity with the units decomposed. |
diagonal ([offset, axis1, axis2]) |
Return specified diagonals. |
diff ([n, axis]) |
Calculate the n-th discrete difference along given axis. |
dot (b[, out]) |
Dot product of two arrays. |
dump (file) |
Dump a pickle of the array to the specified file. |
dumps () |
Returns the pickle of the array as a string. |
ediff1d ([to_end, to_begin]) |
|
fetch (channel, start, end[, bits, host, ...]) |
Fetch data from NDS into a StateVector . |
fetch_open_data (ifo, start, end[, name, host]) |
Fetch open-access data from the LIGO Open Science Center |
fill (value) |
Fill the array with a scalar value. |
find (channel, start, end[, frametype, pad, ...]) |
Find and read data from frames for a channel |
flatten ([order]) |
Return a copy of the array collapsed into one dimension. |
from_hdf5 (*args, **kwargs) |
Read an array from the given HDF file. |
from_lal (*args, **kwargs) |
Generate a new TimeSeries from a LAL TimeSeries of any type. |
from_nds2_buffer (*args, **kwargs) |
Construct a new TimeSeries from an nds2.buffer object |
from_pycbc (ts) |
Convert a pycbc.types.timeseries.TimeSeries into a TimeSeries |
get (channel, start, end[, bits]) |
Get data for this channel from frames or NDS |
get_bit_series ([bits]) |
Get the StateTimeSeries for each bit of this StateVector . |
getfield (dtype[, offset]) |
Returns a field of the given array as a certain type. |
insert (obj, values[, axis]) |
Insert values along the given axis before the given indices and return a new Quantity object. |
is_compatible (other) |
Check whether this series and other have compatible metadata |
is_contiguous (other[, tol]) |
Check whether other is contiguous with self. |
item (*args) |
Copy an element of an array to a standard Python scalar and return it. |
itemset (*args) |
Insert scalar into an array (scalar is cast to array’s dtype, if possible) |
max ([axis, out]) |
Return the maximum along a given axis. |
mean ([axis, dtype, out, keepdims]) |
Returns the average of the array elements along given axis. |
median ([axis]) |
Compute the median along the specified axis. |
min ([axis, out, keepdims]) |
Return the minimum along a given axis. |
nansum ([axis, out, keepdims]) |
|
newbyteorder ([new_order]) |
Return the array with the same data viewed with a different byte order. |
nonzero () |
Return the indices of the elements that are non-zero. |
override_unit (unit[, parse_strict]) |
Forcefully reset the unit of these data |
pad (pad_width, **kwargs) |
Pad this series to a new size |
partition (kth[, axis, kind, order]) |
Rearranges the elements in the array in such a way that value of the element in kth position is in the position it would be in a sorted array. |
plot ([format, bits]) |
Plot the data for this StateVector |
prepend (other[, gap, inplace, pad, resize]) |
Connect another series onto the start of the current one. |
prod ([axis, dtype, out, keepdims]) |
Return the product of the array elements over the given axis |
ptp ([axis, out]) |
Peak to peak (maximum - minimum) value along a given axis. |
put (indices, values[, mode]) |
Set a.flat[n] = values[n] for all n in indices. |
ravel ([order]) |
Return a flattened array. |
read (*args, **kwargs) |
Read data into a StateVector |
repeat (repeats[, axis]) |
Repeat elements of an array. |
resample (rate) |
Resample this StateVector to a new rate |
reshape (shape[, order]) |
Returns an array containing the same data with a new shape. |
resize (new_shape[, refcheck]) |
Change shape and size of array in-place. |
round ([decimals, out]) |
Return a with each element rounded to the given number of decimals. |
searchsorted (v[, side, sorter]) |
Find indices where elements of v should be inserted in a to maintain order. |
setfield (val, dtype[, offset]) |
Put a value into a specified place in a field defined by a data-type. |
setflags ([write, align, uic]) |
Set array flags WRITEABLE, ALIGNED, and UPDATEIFCOPY, respectively. |
sort ([axis, kind, order]) |
Sort an array, in-place. |
squeeze ([axis]) |
Remove single-dimensional entries from the shape of a . |
std ([axis, dtype, out, ddof, keepdims]) |
Returns the standard deviation of the array elements along given axis. |
sum ([axis, dtype, out, keepdims]) |
Return the sum of the array elements over the given axis. |
swapaxes (axis1, axis2) |
Return a view of the array with axis1 and axis2 interchanged. |
take (indices[, axis, out, mode]) |
Return an array formed from the elements of a at the given indices. |
to (unit[, equivalencies]) |
Returns a new Quantity object with the specified units. |
to_dqflags ([bits, minlen, dtype, round]) |
Convert this StateVector into a DataQualityDict |
to_hdf5 (*args, **kwargs) |
Convert this array to a h5py.Dataset . |
to_lal (*args, **kwargs) |
Bogus function inherited from superclass, do not use. |
to_pycbc (*args, **kwargs) |
Convert this TimeSeries into a PyCBC |
tobytes ([order]) |
Construct Python bytes containing the raw data bytes in the array. |
tofile (fid[, sep, format]) |
Write array to a file as text or binary (default). |
tolist () |
Return the array as a (possibly nested) list. |
tostring ([order]) |
Construct Python bytes containing the raw data bytes in the array. |
trace ([offset, axis1, axis2, dtype, out]) |
Return the sum along diagonals of the array. |
transpose (*axes) |
Returns a view of the array with axes transposed. |
update (other[, inplace]) |
Update this series by appending new data from an other and dropping the same amount of data off the start. |
value_at (x) |
Return the value of this Series at the given xindex value |
var ([axis, dtype, out, ddof, keepdims]) |
Returns the variance of the array elements, along given axis. |
view ([dtype, type]) |
New view of array with the same data. |
write (data, *args, **kwargs) |
Write this StateVector to a file |
zip () |
Zip the xindex and value arrays of this Series |
Attributes Documentation
T
¶Same as self.transpose(), except that self is returned if self.ndim < 2.
Examples
>>> x = np.array([[1.,2.],[3.,4.]])
>>> x
array([[ 1., 2.],
[ 3., 4.]])
>>> x.T
array([[ 1., 3.],
[ 2., 4.]])
>>> x = np.array([1.,2.,3.,4.])
>>> x
array([ 1., 2., 3., 4.])
>>> x.T
array([ 1., 2., 3., 4.])
base
¶Base object if memory is from some other object.
Examples
The base of an array that owns its memory is None:
>>> x = np.array([1,2,3,4])
>>> x.base is None
True
Slicing creates a view, whose memory is shared with x:
>>> y = x[2:]
>>> y.base is x
True
bits
¶The list of bit names for this StateVector
.
Type: | Bits |
---|
boolean
¶A mapping of this StateVector
to a 2-D array containing all
binary bits as booleans, for each time point.
cgs
¶Returns a copy of the current Quantity
instance with CGS units. The
value of the resulting object will be scaled.
ctypes
¶An object to simplify the interaction of the array with the ctypes module.
This attribute creates an object that makes it easier to use arrays when calling shared libraries with the ctypes module. The returned object has, among others, data, shape, and strides attributes (see Notes below) which themselves return ctypes objects that can be used as arguments to a shared library.
Parameters: | None : |
---|---|
Returns: | c : Python object
|
See also
numpy.ctypeslib
Notes
Below are the public attributes of this object which were documented in “Guide to NumPy” (we have omitted undocumented public attributes, as well as documented private attributes):
Be careful using the ctypes attribute - especially on temporary
arrays or arrays constructed on the fly. For example, calling
(a+b).ctypes.data_as(ctypes.c_void_p)
returns a pointer to memory
that is invalid because the array created as (a+b) is deallocated
before the next Python statement. You can avoid this problem using
either c=a+b
or ct=(a+b).ctypes
. In the latter case, ct will
hold a reference to the array until ct is deleted or re-assigned.
If the ctypes module is not available, then the ctypes attribute of array objects still returns something useful, but ctypes objects are not returned and errors may be raised instead. In particular, the object will still have the as parameter attribute which will return an integer equal to the data attribute.
Examples
>>> import ctypes
>>> x
array([[0, 1],
[2, 3]])
>>> x.ctypes.data
30439712
>>> x.ctypes.data_as(ctypes.POINTER(ctypes.c_long))
<ctypes.LP_c_long object at 0x01F01300>
>>> x.ctypes.data_as(ctypes.POINTER(ctypes.c_long)).contents
c_long(0)
>>> x.ctypes.data_as(ctypes.POINTER(ctypes.c_longlong)).contents
c_longlong(4294967296L)
>>> x.ctypes.shape
<numpy.core._internal.c_long_Array_2 object at 0x01FFD580>
>>> x.ctypes.shape_as(ctypes.c_long)
<numpy.core._internal.c_long_Array_2 object at 0x01FCE620>
>>> x.ctypes.strides
<numpy.core._internal.c_long_Array_2 object at 0x01FCE620>
>>> x.ctypes.strides_as(ctypes.c_longlong)
<numpy.core._internal.c_longlong_Array_2 object at 0x01F01300>
data
¶dtype
¶Data-type of the array’s elements.
Parameters: | None : |
---|---|
Returns: | d : numpy dtype object |
See also
Examples
>>> x
array([[0, 1],
[2, 3]])
>>> x.dtype
dtype('int32')
>>> type(x.dtype)
<type 'numpy.dtype'>
duration
¶Duration of this series in seconds
equivalencies
¶A list of equivalencies that will be applied by default during unit conversions.
flags
¶Information about the memory layout of the array.
Notes
The flags
object can be accessed dictionary-like (as in a.flags['WRITEABLE']
),
or by using lowercased attribute names (as in a.flags.writeable
). Short flag
names are only supported in dictionary access.
Only the UPDATEIFCOPY, WRITEABLE, and ALIGNED flags can be changed by
the user, via direct assignment to the attribute or dictionary entry,
or by calling ndarray.setflags
.
The array flags cannot be set arbitrarily:
False
.True
if the data is truly aligned.True
if the array owns its own memory
or the ultimate owner of the memory exposes a writeable buffer
interface or is a string.Arrays can be both C-style and Fortran-style contiguous simultaneously. This is clear for 1-dimensional arrays, but can also be true for higher dimensional arrays.
Even for contiguous arrays a stride for a given dimension
arr.strides[dim]
may be arbitrary if arr.shape[dim] == 1
or the array has no elements.
It does not generally hold that self.strides[-1] == self.itemsize
for C-style contiguous arrays or self.strides[0] == self.itemsize
for
Fortran-style contiguous arrays is true.
Attributes
flat
¶A 1-D iterator over the Quantity array.
This returns a QuantityIterator
instance, which behaves the same
as the flatiter
instance returned by flat
,
and is similar to, but not a subclass of, Python’s built-in iterator
object.
imag
¶The imaginary part of the array.
Examples
>>> x = np.sqrt([1+0j, 0+1j])
>>> x.imag
array([ 0. , 0.70710678])
>>> x.imag.dtype
dtype('float64')
info
¶Container for meta information like name, description, format. This is required when the object is used as a mixin column within a table, but can be used as a general way to store meta information.
isscalar
¶True if the value
of this quantity is a scalar, or False if it
is an array-like object.
Note
This is subtly different from numpy.isscalar
in that
numpy.isscalar
returns False for a zero-dimensional array
(e.g. np.array(1)
), while this is True for quantities,
since quantities cannot represent true numpy scalars.
itemsize
¶Length of one array element in bytes.
Examples
>>> x = np.array([1,2,3], dtype=np.float64)
>>> x.itemsize
8
>>> x = np.array([1,2,3], dtype=np.complex128)
>>> x.itemsize
16
nbytes
¶Total bytes consumed by the elements of the array.
Notes
Does not include memory consumed by non-element attributes of the array object.
Examples
>>> x = np.zeros((3,5,2), dtype=np.complex128)
>>> x.nbytes
480
>>> np.prod(x.shape) * x.itemsize
480
ndim
¶Number of array dimensions.
Examples
>>> x = np.array([1, 2, 3])
>>> x.ndim
1
>>> y = np.zeros((2, 3, 4))
>>> y.ndim
3
real
¶The real part of the array.
See also
numpy.real
Examples
>>> x = np.sqrt([1+0j, 0+1j])
>>> x.real
array([ 1. , 0.70710678])
>>> x.real.dtype
dtype('float64')
sample_rate
¶Data rate for this TimeSeries
in samples per second (Hertz).
This attribute is stored internally by the dx
attribute
Type: | Quantity scalar |
---|
shape
¶Tuple of array dimensions.
Notes
May be used to “reshape” the array, as long as this would not require a change in the total number of elements
Examples
>>> x = np.array([1, 2, 3, 4])
>>> x.shape
(4,)
>>> y = np.zeros((2, 3, 4))
>>> y.shape
(2, 3, 4)
>>> y.shape = (3, 8)
>>> y
array([[ 0., 0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0.]])
>>> y.shape = (3, 6)
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
ValueError: total size of new array must be unchanged
si
¶Returns a copy of the current Quantity
instance with SI units. The
value of the resulting object will be scaled.
size
¶Number of elements in the array.
Equivalent to np.prod(a.shape)
, i.e., the product of the array’s
dimensions.
Examples
>>> x = np.zeros((3, 5, 2), dtype=np.complex128)
>>> x.size
30
>>> np.prod(x.shape)
30
strides
¶Tuple of bytes to step in each dimension when traversing an array.
The byte offset of element (i[0], i[1], ..., i[n])
in an array a
is:
offset = sum(np.array(i) * a.strides)
A more detailed explanation of strides can be found in the “ndarray.rst” file in the NumPy reference guide.
See also
numpy.lib.stride_tricks.as_strided
Notes
Imagine an array of 32-bit integers (each 4 bytes):
x = np.array([[0, 1, 2, 3, 4],
[5, 6, 7, 8, 9]], dtype=np.int32)
This array is stored in memory as 40 bytes, one after the other
(known as a contiguous block of memory). The strides of an array tell
us how many bytes we have to skip in memory to move to the next position
along a certain axis. For example, we have to skip 4 bytes (1 value) to
move to the next column, but 20 bytes (5 values) to get to the same
position in the next row. As such, the strides for the array x
will be
(20, 4)
.
Examples
>>> y = np.reshape(np.arange(2*3*4), (2,3,4))
>>> y
array([[[ 0, 1, 2, 3],
[ 4, 5, 6, 7],
[ 8, 9, 10, 11]],
[[12, 13, 14, 15],
[16, 17, 18, 19],
[20, 21, 22, 23]]])
>>> y.strides
(48, 16, 4)
>>> y[1,1,1]
17
>>> offset=sum(y.strides * np.array((1,1,1)))
>>> offset/y.itemsize
17
>>> x = np.reshape(np.arange(5*6*7*8), (5,6,7,8)).transpose(2,3,1,0)
>>> x.strides
(32, 4, 224, 1344)
>>> i = np.array([3,5,2,2])
>>> offset = sum(i * x.strides)
>>> x[3,5,2,2]
813
>>> offset / x.itemsize
813
times
¶Series of GPS times for each sample
value
¶The numerical value of this quantity.
Methods Documentation
all
(axis=None, out=None, keepdims=False)¶Returns True if all elements evaluate to True.
Refer to numpy.all
for full documentation.
See also
numpy.all
any
(axis=None, out=None, keepdims=False)¶Returns True if any of the elements of a
evaluate to True.
Refer to numpy.any
for full documentation.
See also
numpy.any
append
(other, gap='raise', inplace=True, pad=0.0, resize=True)[source]¶Connect another series onto the end of the current one.
Parameters: | other :
gap :
inplace :
pad :
resize :
|
---|---|
Returns: | series :
|
argmax
(axis=None, out=None)¶Return indices of the maximum values along the given axis.
Refer to numpy.argmax
for full documentation.
See also
numpy.argmax
argmin
(axis=None, out=None)¶Return indices of the minimum values along the given axis of a
.
Refer to numpy.argmin
for detailed documentation.
See also
numpy.argmin
argpartition
(kth, axis=-1, kind='introselect', order=None)¶Returns the indices that would partition this array.
Refer to numpy.argpartition
for full documentation.
New in version 1.8.0.
See also
numpy.argpartition
argsort
(axis=-1, kind='quicksort', order=None)¶Returns the indices that would sort this array.
Refer to numpy.argsort
for full documentation.
See also
numpy.argsort
astype
(dtype, order='K', casting='unsafe', subok=True, copy=True)¶Copy of the array, cast to a specified type.
Parameters: | dtype : str or dtype
order : {‘C’, ‘F’, ‘A’, ‘K’}, optional
casting : {‘no’, ‘equiv’, ‘safe’, ‘same_kind’, ‘unsafe’}, optional
subok : bool, optional
copy : bool, optional
|
---|---|
Returns: | arr_t : ndarray |
Raises: | ComplexWarning :
|
Notes
Starting in NumPy 1.9, astype method now returns an error if the string dtype to cast to is not long enough in ‘safe’ casting mode to hold the max value of integer/float array that is being casted. Previously the casting was allowed even if the result was truncated.
Examples
>>> x = np.array([1, 2, 2.5])
>>> x
array([ 1. , 2. , 2.5])
>>> x.astype(int)
array([1, 2, 2])
byteswap
(inplace)¶Swap the bytes of the array elements
Toggle between low-endian and big-endian data representation by returning a byteswapped array, optionally swapped in-place.
Parameters: | inplace : bool, optional
|
---|---|
Returns: | out : ndarray
|
Examples
>>> A = np.array([1, 256, 8755], dtype=np.int16)
>>> map(hex, A)
['0x1', '0x100', '0x2233']
>>> A.byteswap(True)
array([ 256, 1, 13090], dtype=int16)
>>> map(hex, A)
['0x100', '0x1', '0x3322']
Arrays of strings are not swapped
>>> A = np.array(['ceg', 'fac'])
>>> A.byteswap()
array(['ceg', 'fac'],
dtype='|S3')
choose
(choices, out=None, mode='raise')¶Use an index array to construct a new array from a set of choices.
Refer to numpy.choose
for full documentation.
See also
numpy.choose
clip
(min=None, max=None, out=None)¶Return an array whose values are limited to [min, max]
.
One of max or min must be given.
Refer to numpy.clip
for full documentation.
See also
numpy.clip
compress
(condition, axis=None, out=None)¶Return selected slices of this array along given axis.
Refer to numpy.compress
for full documentation.
See also
numpy.compress
conj
()¶Complex-conjugate all elements.
Refer to numpy.conjugate
for full documentation.
See also
numpy.conjugate
conjugate
()¶Return the complex conjugate, element-wise.
Refer to numpy.conjugate
for full documentation.
See also
numpy.conjugate
copy
(order='C')[source]¶Return a copy of the array.
Parameters: | order : {‘C’, ‘F’, ‘A’, ‘K’}, optional
|
---|
See also
Examples
>>> x = np.array([[1,2,3],[4,5,6]], order='F')
>>> y = x.copy()
>>> x.fill(0)
>>> x
array([[0, 0, 0],
[0, 0, 0]])
>>> y
array([[1, 2, 3],
[4, 5, 6]])
>>> y.flags['C_CONTIGUOUS']
True
crop
(start=None, end=None, copy=False)[source]¶Crop this series to the given x-axis extent.
Parameters: | start :
end :
copy :
|
---|---|
Returns: | series :
|
Notes
If either start
or end
are outside of the original
Series
span, warnings will be printed and the limits will
be restricted to the xspan
cumprod
(axis=None, dtype=None, out=None)¶Return the cumulative product of the elements along the given axis.
Refer to numpy.cumprod
for full documentation.
See also
numpy.cumprod
cumsum
(axis=None, dtype=None, out=None)¶Return the cumulative sum of the elements along the given axis.
Refer to numpy.cumsum
for full documentation.
See also
numpy.cumsum
decompose
(bases=[])¶Generates a new Quantity
with the units
decomposed. Decomposed units have only irreducible units in
them (see astropy.units.UnitBase.decompose
).
Parameters: | bases : sequence of UnitBase, optional
|
---|---|
Returns: | newq :
|
diagonal
(offset=0, axis1=0, axis2=1)¶Return specified diagonals. In NumPy 1.9 the returned array is a read-only view instead of a copy as in previous NumPy versions. In a future version the read-only restriction will be removed.
Refer to numpy.diagonal()
for full documentation.
See also
numpy.diagonal
diff
(n=1, axis=-1)[source]¶Calculate the n-th discrete difference along given axis.
The first difference is given by
out[n] = a[n+1] - a[n]
along the given axis, higher differences are calculated by usingdiff
recursively.
Parameters: | a : array_like
|
---|---|
Returns: | diff : ndarray
. : |
Examples
>>> x = np.array([1, 2, 4, 7, 0])
>>> np.diff(x)
array([ 1, 2, 3, -7])
>>> np.diff(x, n=2)
array([ 1, 1, -10])
>>> x = np.array([[1, 3, 6, 10], [0, 5, 6, 8]])
>>> np.diff(x)
array([[2, 3, 4],
[5, 1, 2]])
>>> np.diff(x, axis=0)
array([[-1, 2, 0, -2]])
dot
(b, out=None)¶Dot product of two arrays.
Refer to numpy.dot
for full documentation.
See also
numpy.dot
Examples
>>> a = np.eye(2)
>>> b = np.ones((2, 2)) * 2
>>> a.dot(b)
array([[ 2., 2.],
[ 2., 2.]])
This array method can be conveniently chained:
>>> a.dot(b).dot(b)
array([[ 8., 8.],
[ 8., 8.]])
dump
(file)¶Dump a pickle of the array to the specified file. The array can be read back with pickle.load or numpy.load.
Parameters: | file : str
|
---|
dumps
()[source]¶Returns the pickle of the array as a string. pickle.loads or numpy.loads will convert the string back to an array.
Parameters: | None : |
---|
ediff1d
(to_end=None, to_begin=None)¶fetch
(channel, start, end, bits=None, host=None, port=None, verbose=False, connection=None, type=102)[source]¶Fetch data from NDS into a StateVector
.
Parameters: |
start :
end :
bits :
host :
port :
verify :
connection :
verbose :
type :
dtype :
|
---|
fetch_open_data
(ifo, start, end, name='quality/simple', host='https://losc.ligo.org')[source]¶Fetch open-access data from the LIGO Open Science Center
Parameters: | ifo :
start :
end :
name :
host :
|
---|
fill
(value)¶Fill the array with a scalar value.
Parameters: | value : scalar
|
---|
Examples
>>> a = np.array([1, 2])
>>> a.fill(0)
>>> a
array([0, 0])
>>> a = np.empty(2)
>>> a.fill(1)
>>> a
array([ 1., 1.])
find
(channel, start, end, frametype=None, pad=None, dtype=None, nproc=1, verbose=False, **readargs)[source]¶Find and read data from frames for a channel
Parameters: |
start :
end :
frametype :
pad :
nproc :
dtype : verbose :
**readargs :
|
---|
flatten
(order='C')¶Return a copy of the array collapsed into one dimension.
Parameters: | order : {‘C’, ‘F’, ‘A’, ‘K’}, optional
|
---|---|
Returns: | y : ndarray
|
Examples
>>> a = np.array([[1,2], [3,4]])
>>> a.flatten()
array([1, 2, 3, 4])
>>> a.flatten('F')
array([1, 3, 2, 4])
from_hdf5
(*args, **kwargs)[source]¶Read an array from the given HDF file.
This method has been deprecated in favour of the unified I/O method:
Class.write(..., format=’hdf5’)
from_nds2_buffer
(*args, **kwargs)[source]¶Construct a new TimeSeries
from an nds2.buffer
object
Parameters: | buffer_ :
**metadata :
|
---|---|
Returns: | timeseries :
|
Notes
This classmethod requires the nds2-client package
from_pycbc
(ts)[source]¶Convert a pycbc.types.timeseries.TimeSeries
into a TimeSeries
Parameters: | ts :
|
---|---|
Returns: | timeseries :
|
get
(channel, start, end, bits=None, **kwargs)[source]¶Get data for this channel from frames or NDS
Parameters: |
start :
end :
bits :
pad :
dtype : nproc :
verbose :
**kwargs other keyword arguments to pass to either : |
---|
See also
StateVector.fetch
StateVector.find
get_bit_series
(bits=None)[source]¶Get the StateTimeSeries
for each bit of this StateVector
.
Parameters: | bits :
|
---|---|
Returns: | bitseries :
|
getfield
(dtype, offset=0)¶Returns a field of the given array as a certain type.
A field is a view of the array data with a given data-type. The values in the view are determined by the given type and the offset into the current array in bytes. The offset needs to be such that the view dtype fits in the array dtype; for example an array of dtype complex128 has 16-byte elements. If taking a view with a 32-bit integer (4 bytes), the offset needs to be between 0 and 12 bytes.
Parameters: | dtype : str or dtype
offset : int
|
---|
Examples
>>> x = np.diag([1.+1.j]*2)
>>> x[1, 1] = 2 + 4.j
>>> x
array([[ 1.+1.j, 0.+0.j],
[ 0.+0.j, 2.+4.j]])
>>> x.getfield(np.float64)
array([[ 1., 0.],
[ 0., 2.]])
By choosing an offset of 8 bytes we can select the complex part of the array for our view:
>>> x.getfield(np.float64, offset=8)
array([[ 1., 0.],
[ 0., 4.]])
insert
(obj, values, axis=None)¶Insert values along the given axis before the given indices and return
a new Quantity
object.
This is a thin wrapper around the numpy.insert
function.
Parameters: | obj : int, slice or sequence of ints
values : array-like
axis : int, optional
|
---|---|
Returns: | out :
|
Examples
>>> import astropy.units as u
>>> q = [1, 2] * u.m
>>> q.insert(0, 50 * u.cm)
<Quantity [ 0.5, 1., 2.] m>
>>> q = [[1, 2], [3, 4]] * u.m
>>> q.insert(1, [10, 20] * u.m, axis=0)
<Quantity [[ 1., 2.],
[ 10., 20.],
[ 3., 4.]] m>
>>> q.insert(1, 10 * u.m, axis=1)
<Quantity [[ 1., 10., 2.],
[ 3., 10., 4.]] m>
is_compatible
(other)[source]¶Check whether this series and other have compatible metadata
This method tests that the sample size
, and the
unit
match.
is_contiguous
(other, tol=3.814697265625e-06)[source]¶Check whether other is contiguous with self.
Parameters: | other :
tol :
|
---|---|
Returns: | 1 :
-1 :
0 :
|
Notes
if a raw numpy.ndarray
is passed as other, with no metadata, then
the contiguity check will always pass
item
(*args)¶Copy an element of an array to a standard Python scalar and return it.
Parameters: | *args : Arguments (variable number and type)
|
---|---|
Returns: | z : Standard Python scalar object
|
Notes
When the data type of a
is longdouble or clongdouble, item() returns
a scalar array object because there is no available Python scalar that
would not lose information. Void arrays return a buffer object for item(),
unless fields are defined, in which case a tuple is returned.
item
is very similar to a[args], except, instead of an array scalar,
a standard Python scalar is returned. This can be useful for speeding up
access to elements of the array and doing arithmetic on elements of the
array using Python’s optimized math.
Examples
>>> x = np.random.randint(9, size=(3, 3))
>>> x
array([[3, 1, 7],
[2, 8, 3],
[8, 5, 3]])
>>> x.item(3)
2
>>> x.item(7)
5
>>> x.item((0, 1))
1
>>> x.item((2, 2))
3
itemset
(*args)¶Insert scalar into an array (scalar is cast to array’s dtype, if possible)
There must be at least 1 argument, and define the last argument
as item. Then, a.itemset(*args)
is equivalent to but faster
than a[args] = item
. The item should be a scalar value and args
must select a single item in the array a
.
Parameters: | *args : Arguments
|
---|
Notes
Compared to indexing syntax, itemset
provides some speed increase
for placing a scalar into a particular location in an ndarray
,
if you must do this. However, generally this is discouraged:
among other problems, it complicates the appearance of the code.
Also, when using itemset
(and item
) inside a loop, be sure
to assign the methods to a local variable to avoid the attribute
look-up at each loop iteration.
Examples
>>> x = np.random.randint(9, size=(3, 3))
>>> x
array([[3, 1, 7],
[2, 8, 3],
[8, 5, 3]])
>>> x.itemset(4, 0)
>>> x.itemset((2, 2), 9)
>>> x
array([[3, 1, 7],
[2, 0, 3],
[8, 5, 9]])
max
(axis=None, out=None)¶Return the maximum along a given axis.
Refer to numpy.amax
for full documentation.
See also
numpy.amax
mean
(axis=None, dtype=None, out=None, keepdims=False)¶Returns the average of the array elements along given axis.
Refer to numpy.mean
for full documentation.
See also
numpy.mean
median
(axis=None, **kwargs)[source]¶Compute the median along the specified axis.
Returns the median of the array elements.
Parameters: | a : array_like
axis : {int, sequence of int, None}, optional
out : ndarray, optional
overwrite_input : bool, optional
keepdims : bool, optional
|
---|---|
Returns: | median : ndarray
|
See also
mean
, percentile
Notes
Given a vector V
of length N
, the median of V
is the
middle value of a sorted copy of V
, V_sorted
- i
e., V_sorted[(N-1)/2]
, when N
is odd, and the average of the
two middle values of V_sorted
when N
is even.
Examples
>>> a = np.array([[10, 7, 4], [3, 2, 1]])
>>> a
array([[10, 7, 4],
[ 3, 2, 1]])
>>> np.median(a)
3.5
>>> np.median(a, axis=0)
array([ 6.5, 4.5, 2.5])
>>> np.median(a, axis=1)
array([ 7., 2.])
>>> m = np.median(a, axis=0)
>>> out = np.zeros_like(m)
>>> np.median(a, axis=0, out=m)
array([ 6.5, 4.5, 2.5])
>>> m
array([ 6.5, 4.5, 2.5])
>>> b = a.copy()
>>> np.median(b, axis=1, overwrite_input=True)
array([ 7., 2.])
>>> assert not np.all(a==b)
>>> b = a.copy()
>>> np.median(b, axis=None, overwrite_input=True)
3.5
>>> assert not np.all(a==b)
min
(axis=None, out=None, keepdims=False)¶Return the minimum along a given axis.
Refer to numpy.amin
for full documentation.
See also
numpy.amin
nansum
(axis=None, out=None, keepdims=False)¶newbyteorder
(new_order='S')¶Return the array with the same data viewed with a different byte order.
Equivalent to:
arr.view(arr.dtype.newbytorder(new_order))
Changes are also made in all fields and sub-arrays of the array data type.
Parameters: | new_order : string, optional
|
---|---|
Returns: | new_arr : array
|
nonzero
()¶Return the indices of the elements that are non-zero.
Refer to numpy.nonzero
for full documentation.
See also
numpy.nonzero
override_unit
(unit, parse_strict='raise')[source]¶Forcefully reset the unit of these data
Use of this method is discouraged in favour of to()
,
which performs accurate conversions from one unit to another.
The method should really only be used when the original unit of the
array is plain wrong.
Parameters: |
|
---|---|
Raises: | ValueError :
|
pad
(pad_width, **kwargs)[source]¶Pad this series to a new size
Parameters: | pad_width :
**kwargs :
|
---|---|
Returns: | series :
|
See also
numpy.pad
partition
(kth, axis=-1, kind='introselect', order=None)¶Rearranges the elements in the array in such a way that value of the element in kth position is in the position it would be in a sorted array. All elements smaller than the kth element are moved before this element and all equal or greater are moved behind it. The ordering of the elements in the two partitions is undefined.
New in version 1.8.0.
Parameters: | kth : int or sequence of ints
axis : int, optional
kind : {‘introselect’}, optional
order : str or list of str, optional
|
---|
See also
numpy.partition
argpartition
sort
Notes
See np.partition
for notes on the different algorithms.
Examples
>>> a = np.array([3, 4, 2, 1])
>>> a.partition(a, 3)
>>> a
array([2, 1, 3, 4])
>>> a.partition((1, 3))
array([1, 2, 3, 4])
plot
(format='segments', bits=None, **kwargs)[source]¶Plot the data for this StateVector
Parameters: | format :
bits :
**kwargs :
|
---|---|
Returns: | plot :
|
prepend
(other, gap='raise', inplace=True, pad=0.0, resize=True)[source]¶Connect another series onto the start of the current one.
Parameters: | other :
gap :
inplace :
pad :
|
---|---|
Returns: | series :
|
prod
(axis=None, dtype=None, out=None, keepdims=False)¶Return the product of the array elements over the given axis
Refer to numpy.prod
for full documentation.
See also
numpy.prod
ptp
(axis=None, out=None)¶Peak to peak (maximum - minimum) value along a given axis.
Refer to numpy.ptp
for full documentation.
See also
numpy.ptp
put
(indices, values, mode='raise')¶Set a.flat[n] = values[n]
for all n
in indices.
Refer to numpy.put
for full documentation.
See also
numpy.put
ravel
([order])¶Return a flattened array.
Refer to numpy.ravel
for full documentation.
See also
numpy.ravel
ndarray.flat
read
(*args, **kwargs)¶Read data into a StateVector
Parameters: |
start :
end :
bits :
format :
nproc :
gap :
pad :
|
---|
Notes
The available built-in formats are:
Format | Read | Write | Auto-identify |
---|---|---|---|
csv | Yes | Yes | Yes |
framecpp | Yes | Yes | No |
gwf | Yes | No | Yes |
hdf | Yes | Yes | No |
hdf5 | Yes | Yes | Yes |
lalframe | Yes | No | No |
losc | Yes | No | No |
txt | Yes | Yes | Yes |
repeat
(repeats, axis=None)¶Repeat elements of an array.
Refer to numpy.repeat
for full documentation.
See also
numpy.repeat
resample
(rate)[source]¶Resample this StateVector
to a new rate
Because of the nature of a state-vector, downsampling is done by taking the logical ‘and’ of all original samples in each new sampling interval, while upsampling is achieved by repeating samples.
Parameters: | rate :
|
---|---|
Returns: | vector :
|
reshape
(shape, order='C')¶Returns an array containing the same data with a new shape.
Refer to numpy.reshape
for full documentation.
See also
numpy.reshape
resize
(new_shape, refcheck=True)¶Change shape and size of array in-place.
Parameters: | new_shape : tuple of ints, or
refcheck : bool, optional
|
---|---|
Returns: | None : |
Raises: | ValueError :
SystemError :
|
See also
resize
Notes
This reallocates space for the data area if necessary.
Only contiguous arrays (data elements consecutive in memory) can be resized.
The purpose of the reference count check is to make sure you
do not use this array as a buffer for another Python object and then
reallocate the memory. However, reference counts can increase in
other ways so if you are sure that you have not shared the memory
for this array with another Python object, then you may safely set
refcheck
to False.
Examples
Shrinking an array: array is flattened (in the order that the data are stored in memory), resized, and reshaped:
>>> a = np.array([[0, 1], [2, 3]], order='C')
>>> a.resize((2, 1))
>>> a
array([[0],
[1]])
>>> a = np.array([[0, 1], [2, 3]], order='F')
>>> a.resize((2, 1))
>>> a
array([[0],
[2]])
Enlarging an array: as above, but missing entries are filled with zeros:
>>> b = np.array([[0, 1], [2, 3]])
>>> b.resize(2, 3) # new_shape parameter doesn't have to be a tuple
>>> b
array([[0, 1, 2],
[3, 0, 0]])
Referencing an array prevents resizing...
>>> c = a
>>> a.resize((1, 1))
Traceback (most recent call last):
...
ValueError: cannot resize an array that has been referenced ...
Unless refcheck
is False:
>>> a.resize((1, 1), refcheck=False)
>>> a
array([[0]])
>>> c
array([[0]])
round
(decimals=0, out=None)¶Return a
with each element rounded to the given number of decimals.
Refer to numpy.around
for full documentation.
See also
numpy.around
searchsorted
(v, side='left', sorter=None)¶Find indices where elements of v should be inserted in a to maintain order.
For full documentation, see numpy.searchsorted
See also
numpy.searchsorted
setfield
(val, dtype, offset=0)¶Put a value into a specified place in a field defined by a data-type.
Place val
into a
‘s field defined by dtype
and beginning offset
bytes into the field.
Parameters: | val : object
dtype : dtype object
offset : int, optional
|
---|---|
Returns: | None : |
See also
Examples
>>> x = np.eye(3)
>>> x.getfield(np.float64)
array([[ 1., 0., 0.],
[ 0., 1., 0.],
[ 0., 0., 1.]])
>>> x.setfield(3, np.int32)
>>> x.getfield(np.int32)
array([[3, 3, 3],
[3, 3, 3],
[3, 3, 3]])
>>> x
array([[ 1.00000000e+000, 1.48219694e-323, 1.48219694e-323],
[ 1.48219694e-323, 1.00000000e+000, 1.48219694e-323],
[ 1.48219694e-323, 1.48219694e-323, 1.00000000e+000]])
>>> x.setfield(np.eye(3), np.int32)
>>> x
array([[ 1., 0., 0.],
[ 0., 1., 0.],
[ 0., 0., 1.]])
setflags
(write=None, align=None, uic=None)¶Set array flags WRITEABLE, ALIGNED, and UPDATEIFCOPY, respectively.
These Boolean-valued flags affect how numpy interprets the memory
area used by a
(see Notes below). The ALIGNED flag can only
be set to True if the data is actually aligned according to the type.
The UPDATEIFCOPY flag can never be set to True. The flag WRITEABLE
can only be set to True if the array owns its own memory, or the
ultimate owner of the memory exposes a writeable buffer interface,
or is a string. (The exception for string is made so that unpickling
can be done without copying memory.)
Parameters: | write : bool, optional
align : bool, optional
uic : bool, optional
|
---|
Notes
Array flags provide information about how the memory area used for the array is to be interpreted. There are 6 Boolean flags in use, only three of which can be changed by the user: UPDATEIFCOPY, WRITEABLE, and ALIGNED.
WRITEABLE (W) the data area can be written to;
ALIGNED (A) the data and strides are aligned appropriately for the hardware (as determined by the compiler);
UPDATEIFCOPY (U) this array is a copy of some other array (referenced by .base). When this array is deallocated, the base array will be updated with the contents of this array.
All flags can be accessed using their first (upper case) letter as well as the full name.
Examples
>>> y
array([[3, 1, 7],
[2, 0, 0],
[8, 5, 9]])
>>> y.flags
C_CONTIGUOUS : True
F_CONTIGUOUS : False
OWNDATA : True
WRITEABLE : True
ALIGNED : True
UPDATEIFCOPY : False
>>> y.setflags(write=0, align=0)
>>> y.flags
C_CONTIGUOUS : True
F_CONTIGUOUS : False
OWNDATA : True
WRITEABLE : False
ALIGNED : False
UPDATEIFCOPY : False
>>> y.setflags(uic=1)
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
ValueError: cannot set UPDATEIFCOPY flag to True
sort
(axis=-1, kind='quicksort', order=None)¶Sort an array, in-place.
Parameters: | axis : int, optional
kind : {‘quicksort’, ‘mergesort’, ‘heapsort’}, optional
order : str or list of str, optional
|
---|
See also
numpy.sort
argsort
lexsort
searchsorted
partition
Notes
See sort
for notes on the different sorting algorithms.
Examples
>>> a = np.array([[1,4], [3,1]])
>>> a.sort(axis=1)
>>> a
array([[1, 4],
[1, 3]])
>>> a.sort(axis=0)
>>> a
array([[1, 3],
[1, 4]])
Use the order
keyword to specify a field to use when sorting a
structured array:
>>> a = np.array([('a', 2), ('c', 1)], dtype=[('x', 'S1'), ('y', int)])
>>> a.sort(order='y')
>>> a
array([('c', 1), ('a', 2)],
dtype=[('x', '|S1'), ('y', '<i4')])
squeeze
(axis=None)¶Remove single-dimensional entries from the shape of a
.
Refer to numpy.squeeze
for full documentation.
See also
numpy.squeeze
std
(axis=None, dtype=None, out=None, ddof=0, keepdims=False)¶Returns the standard deviation of the array elements along given axis.
Refer to numpy.std
for full documentation.
See also
numpy.std
sum
(axis=None, dtype=None, out=None, keepdims=False)¶Return the sum of the array elements over the given axis.
Refer to numpy.sum
for full documentation.
See also
numpy.sum
swapaxes
(axis1, axis2)¶Return a view of the array with axis1
and axis2
interchanged.
Refer to numpy.swapaxes
for full documentation.
See also
numpy.swapaxes
take
(indices, axis=None, out=None, mode='raise')¶Return an array formed from the elements of a
at the given indices.
Refer to numpy.take
for full documentation.
See also
numpy.take
to
(unit, equivalencies=[])¶Returns a new Quantity
object with the specified
units.
Parameters: | unit : equivalencies : list of equivalence pairs, optional
|
---|
to_dqflags
(bits=None, minlen=1, dtype=<type 'float'>, round=False)[source]¶Convert this StateVector
into a DataQualityDict
The StateTimeSeries
for each bit is converted into a
DataQualityFlag
with the bits combined into a dict.
Parameters: | minlen :
bits :
|
---|---|
Returns: | DataQualityFlag list :
|
See also
StateTimeSeries.to_dqflag()
StateVector
bitsto_hdf5
(*args, **kwargs)[source]¶Convert this array to a h5py.Dataset
.
This method has been deprecated in favour of the unified I/O method:
Class.write(..., format=’hdf5’)
to_pycbc
(*args, **kwargs)[source]¶Convert this TimeSeries
into a PyCBC
TimeSeries
Parameters: | copy :
|
---|---|
Returns: | timeseries :
|
tobytes
(order='C')¶Construct Python bytes containing the raw data bytes in the array.
Constructs Python bytes showing a copy of the raw contents of data memory. The bytes object can be produced in either ‘C’ or ‘Fortran’, or ‘Any’ order (the default is ‘C’-order). ‘Any’ order means C-order unless the F_CONTIGUOUS flag in the array is set, in which case it means ‘Fortran’ order.
New in version 1.9.0.
Parameters: | order : {‘C’, ‘F’, None}, optional
|
---|---|
Returns: | s : bytes
|
Examples
>>> x = np.array([[0, 1], [2, 3]])
>>> x.tobytes()
b'\x00\x00\x00\x00\x01\x00\x00\x00\x02\x00\x00\x00\x03\x00\x00\x00'
>>> x.tobytes('C') == x.tobytes()
True
>>> x.tobytes('F')
b'\x00\x00\x00\x00\x02\x00\x00\x00\x01\x00\x00\x00\x03\x00\x00\x00'
tofile
(fid, sep="", format="%s")¶Write array to a file as text or binary (default).
Data is always written in ‘C’ order, independent of the order of a
.
The data produced by this method can be recovered using the function
fromfile().
Parameters: | fid : file or str
sep : str
format : str
|
---|
Notes
This is a convenience function for quick storage of array data. Information on endianness and precision is lost, so this method is not a good choice for files intended to archive data or transport data between machines with different endianness. Some of these problems can be overcome by outputting the data as text files, at the expense of speed and file size.
tolist
()¶Return the array as a (possibly nested) list.
Return a copy of the array data as a (nested) Python list. Data items are converted to the nearest compatible Python type.
Parameters: | none : |
---|---|
Returns: | y : list
|
Notes
The array may be recreated, a = np.array(a.tolist())
.
Examples
>>> a = np.array([1, 2])
>>> a.tolist()
[1, 2]
>>> a = np.array([[1, 2], [3, 4]])
>>> list(a)
[array([1, 2]), array([3, 4])]
>>> a.tolist()
[[1, 2], [3, 4]]
tostring
(order='C')[source]¶Construct Python bytes containing the raw data bytes in the array.
Constructs Python bytes showing a copy of the raw contents of data memory. The bytes object can be produced in either ‘C’ or ‘Fortran’, or ‘Any’ order (the default is ‘C’-order). ‘Any’ order means C-order unless the F_CONTIGUOUS flag in the array is set, in which case it means ‘Fortran’ order.
This function is a compatibility alias for tobytes. Despite its name it returns bytes not strings.
Parameters: | order : {‘C’, ‘F’, None}, optional
|
---|---|
Returns: | s : bytes
|
Examples
>>> x = np.array([[0, 1], [2, 3]])
>>> x.tobytes()
b'\x00\x00\x00\x00\x01\x00\x00\x00\x02\x00\x00\x00\x03\x00\x00\x00'
>>> x.tobytes('C') == x.tobytes()
True
>>> x.tobytes('F')
b'\x00\x00\x00\x00\x02\x00\x00\x00\x01\x00\x00\x00\x03\x00\x00\x00'
trace
(offset=0, axis1=0, axis2=1, dtype=None, out=None)¶Return the sum along diagonals of the array.
Refer to numpy.trace
for full documentation.
See also
numpy.trace
transpose
(*axes)¶Returns a view of the array with axes transposed.
For a 1-D array, this has no effect. (To change between column and
row vectors, first cast the 1-D array into a matrix object.)
For a 2-D array, this is the usual matrix transpose.
For an n-D array, if axes are given, their order indicates how the
axes are permuted (see Examples). If axes are not provided and
a.shape = (i[0], i[1], ... i[n-2], i[n-1])
, then
a.transpose().shape = (i[n-1], i[n-2], ... i[1], i[0])
.
Parameters: | axes : None, tuple of ints, or
|
---|---|
Returns: | out : ndarray
|
See also
ndarray.T
Examples
>>> a = np.array([[1, 2], [3, 4]])
>>> a
array([[1, 2],
[3, 4]])
>>> a.transpose()
array([[1, 3],
[2, 4]])
>>> a.transpose((1, 0))
array([[1, 3],
[2, 4]])
>>> a.transpose(1, 0)
array([[1, 3],
[2, 4]])
update
(other, inplace=True)[source]¶Update this series by appending new data from an other and dropping the same amount of data off the start.
This is a convenience method that just calls append
with
resize=False
.
value_at
(x)[source]¶Return the value of this Series
at the given xindex
value
Parameters: |
|
---|---|
Returns: | y :
|
var
(axis=None, dtype=None, out=None, ddof=0, keepdims=False)¶Returns the variance of the array elements, along given axis.
Refer to numpy.var
for full documentation.
See also
numpy.var
view
(dtype=None, type=None)¶New view of array with the same data.
Parameters: | dtype : data-type or ndarray sub-class, optional
type : Python type, optional
|
---|
Notes
a.view()
is used two different ways:
a.view(some_dtype)
or a.view(dtype=some_dtype)
constructs a view
of the array’s memory with a different data-type. This can cause a
reinterpretation of the bytes of memory.
a.view(ndarray_subclass)
or a.view(type=ndarray_subclass)
just
returns an instance of ndarray_subclass
that looks at the same array
(same shape, dtype, etc.) This does not cause a reinterpretation of the
memory.
For a.view(some_dtype)
, if some_dtype
has a different number of
bytes per entry than the previous dtype (for example, converting a
regular array to a structured array), then the behavior of the view
cannot be predicted just from the superficial appearance of a
(shown
by print(a)
). It also depends on exactly how a
is stored in
memory. Therefore if a
is C-ordered versus fortran-ordered, versus
defined as a slice or transpose, etc., the view may give different
results.
Examples
>>> x = np.array([(1, 2)], dtype=[('a', np.int8), ('b', np.int8)])
Viewing array data using a different type and dtype:
>>> y = x.view(dtype=np.int16, type=np.matrix)
>>> y
matrix([[513]], dtype=int16)
>>> print(type(y))
<class 'numpy.matrixlib.defmatrix.matrix'>
Creating a view on a structured array so it can be used in calculations
>>> x = np.array([(1, 2),(3,4)], dtype=[('a', np.int8), ('b', np.int8)])
>>> xv = x.view(dtype=np.int8).reshape(-1,2)
>>> xv
array([[1, 2],
[3, 4]], dtype=int8)
>>> xv.mean(0)
array([ 2., 3.])
Making changes to the view changes the underlying array
>>> xv[0,1] = 20
>>> print(x)
[(1, 20) (3, 4)]
Using a view to convert an array to a recarray:
>>> z = x.view(np.recarray)
>>> z.a
array([1], dtype=int8)
Views share data:
>>> x[0] = (9, 10)
>>> z[0]
(9, 10)
Views that change the dtype size (bytes per entry) should normally be avoided on arrays defined by slices, transposes, fortran-ordering, etc.:
>>> x = np.array([[1,2,3],[4,5,6]], dtype=np.int16)
>>> y = x[:, 0:2]
>>> y
array([[1, 2],
[4, 5]], dtype=int16)
>>> y.view(dtype=[('width', np.int16), ('length', np.int16)])
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
ValueError: new type not compatible with array.
>>> z = y.copy()
>>> z.view(dtype=[('width', np.int16), ('length', np.int16)])
array([[(1, 2)],
[(4, 5)]], dtype=[('width', '<i2'), ('length', '<i2')])
write
(data, *args, **kwargs)¶Write this StateVector
to a file
Parameters: | outfile :
|
---|
Notes
The available built-in formats are:
Format | Read | Write | Auto-identify |
---|---|---|---|
csv | Yes | Yes | Yes |
framecpp | Yes | Yes | No |
hdf | Yes | Yes | No |
hdf5 | Yes | Yes | Yes |
txt | Yes | Yes | Yes |
zip
()[source]¶Zip the xindex
and value
arrays of this Series
Returns: | stacked : 2-d |
---|
Examples
>>> a = Series([0, 2, 4, 6, 8], xindex=[-5, -4, -3, -2, -1])
>>> a.zip()
array([[-5., 0.],
[-4., 2.],
[-3., 4.],
[-2., 6.],
[-1., 8.]])
gwpy.timeseries.
StateTimeSeries
[source]¶Bases: gwpy.timeseries.core.TimeSeriesBase
Boolean array representing a good/bad state determination of some data.
Parameters: | value : array-like
unit :
epoch :
sample_rate :
times :
name :
dtype :
copy :
subok :
|
---|
Notes
The input data array is cast to the bool
data type upon creation of
this series.
Attributes Summary
T |
Same as self.transpose(), except that self is returned if self.ndim < 2. |
base |
Base object if memory is from some other object. |
cgs |
Returns a copy of the current Quantity instance with CGS units. |
channel |
Instrumental channel associated with these data |
ctypes |
An object to simplify the interaction of the array with the ctypes module. |
data |
|
dt |
X-axis sample separation |
dtype |
Data-type of the array’s elements. |
duration |
Duration of this series in seconds |
dx |
X-axis sample separation |
epoch |
GPS epoch for these data. |
equivalencies |
A list of equivalencies that will be applied by default during unit conversions. |
flags |
Information about the memory layout of the array. |
flat |
A 1-D iterator over the Quantity array. |
imag |
The imaginary part of the array. |
info |
Container for meta information like name, description, format. |
isscalar |
True if the value of this quantity is a scalar, or False if it is an array-like object. |
itemsize |
Length of one array element in bytes. |
name |
Name for this data set |
nbytes |
Total bytes consumed by the elements of the array. |
ndim |
Number of array dimensions. |
real |
The real part of the array. |
sample_rate |
Data rate for this TimeSeries in samples per second (Hertz). |
shape |
Tuple of array dimensions. |
si |
Returns a copy of the current Quantity instance with SI units. |
size |
Number of elements in the array. |
span |
X-axis [low, high) segment encompassed by these data |
strides |
Tuple of bytes to step in each dimension when traversing an array. |
times |
Series of GPS times for each sample |
unit |
The physical unit of these data |
value |
The numerical value of this quantity. |
x0 |
X-axis value of the first data point |
xindex |
Positions of the data on the x-axis |
xspan |
X-axis [low, high) segment encompassed by these data |
xunit |
Unit of x-axis index |
Methods Summary
all ([axis, out, keepdims]) |
Returns True if all elements evaluate to True. |
any ([axis, out, keepdims]) |
Returns True if any of the elements of a evaluate to True. |
append (other[, gap, inplace, pad, resize]) |
Connect another series onto the end of the current one. |
argmax ([axis, out]) |
Return indices of the maximum values along the given axis. |
argmin ([axis, out]) |
Return indices of the minimum values along the given axis of a . |
argpartition (kth[, axis, kind, order]) |
Returns the indices that would partition this array. |
argsort ([axis, kind, order]) |
Returns the indices that would sort this array. |
astype (dtype[, order, casting, subok, copy]) |
Copy of the array, cast to a specified type. |
byteswap (inplace) |
Swap the bytes of the array elements |
choose (choices[, out, mode]) |
Use an index array to construct a new array from a set of choices. |
clip ([min, max, out]) |
Return an array whose values are limited to [min, max] . |
compress (condition[, axis, out]) |
Return selected slices of this array along given axis. |
conj () |
Complex-conjugate all elements. |
conjugate () |
Return the complex conjugate, element-wise. |
copy ([order]) |
Return a copy of the array. |
copy_metadata () |
Return a deepcopy of the metadata for this array |
crop ([start, end, copy]) |
Crop this series to the given x-axis extent. |
cumprod ([axis, dtype, out]) |
Return the cumulative product of the elements along the given axis. |
cumsum ([axis, dtype, out]) |
Return the cumulative sum of the elements along the given axis. |
decompose ([bases]) |
Generates a new Quantity with the units decomposed. |
diagonal ([offset, axis1, axis2]) |
Return specified diagonals. |
diff ([n, axis]) |
Calculate the n-th discrete difference along given axis. |
dot (b[, out]) |
Dot product of two arrays. |
dump (file) |
Dump a pickle of the array to the specified file. |
dumps () |
Returns the pickle of the array as a string. |
ediff1d ([to_end, to_begin]) |
|
fetch (*args, **kwargs) |
Fetch data from NDS into a TimeSeries . |
fetch_open_data (ifo, start, end[, name, ...]) |
Fetch open-access data from the LIGO Open Science Center |
fill (value) |
Fill the array with a scalar value. |
find (channel, start, end[, frametype, pad, ...]) |
Find and read data from frames for a channel |
flatten ([order]) |
Return a copy of the array collapsed into one dimension. |
from_hdf5 (*args, **kwargs) |
Read an array from the given HDF file. |
from_lal (*args, **kwargs) |
Generate a new TimeSeries from a LAL TimeSeries of any type. |
from_nds2_buffer (*args, **kwargs) |
Construct a new TimeSeries from an nds2.buffer object |
from_pycbc (ts) |
Convert a pycbc.types.timeseries.TimeSeries into a TimeSeries |
get (channel, start, end[, pad, dtype, verbose]) |
Get data for this channel from frames or NDS |
getfield (dtype[, offset]) |
Returns a field of the given array as a certain type. |
insert (obj, values[, axis]) |
Insert values along the given axis before the given indices and return a new Quantity object. |
is_compatible (other) |
Check whether this series and other have compatible metadata |
is_contiguous (other[, tol]) |
Check whether other is contiguous with self. |
item (*args) |
Copy an element of an array to a standard Python scalar and return it. |
itemset (*args) |
Insert scalar into an array (scalar is cast to array’s dtype, if possible) |
max ([axis, out]) |
Return the maximum along a given axis. |
mean ([axis, dtype, out, keepdims]) |
Returns the average of the array elements along given axis. |
median ([axis]) |
Compute the median along the specified axis. |
min ([axis, out, keepdims]) |
Return the minimum along a given axis. |
nansum ([axis, out, keepdims]) |
|
newbyteorder ([new_order]) |
Return the array with the same data viewed with a different byte order. |
nonzero () |
Return the indices of the elements that are non-zero. |
override_unit (unit[, parse_strict]) |
Forcefully reset the unit of these data |
pad (pad_width, **kwargs) |
Pad this series to a new size |
partition (kth[, axis, kind, order]) |
Rearranges the elements in the array in such a way that value of the element in kth position is in the position it would be in a sorted array. |
plot (**kwargs) |
Plot the data for this TimeSeries |
prepend (other[, gap, inplace, pad, resize]) |
Connect another series onto the start of the current one. |
prod ([axis, dtype, out, keepdims]) |
Return the product of the array elements over the given axis |
ptp ([axis, out]) |
Peak to peak (maximum - minimum) value along a given axis. |
put (indices, values[, mode]) |
Set a.flat[n] = values[n] for all n in indices. |
ravel ([order]) |
Return a flattened array. |
read (*args, **kwargs) |
Read in data |
repeat (repeats[, axis]) |
Repeat elements of an array. |
reshape (shape[, order]) |
Returns an array containing the same data with a new shape. |
resize (new_shape[, refcheck]) |
Change shape and size of array in-place. |
round ([decimals, out]) |
Return a with each element rounded to the given number of decimals. |
searchsorted (v[, side, sorter]) |
Find indices where elements of v should be inserted in a to maintain order. |
setfield (val, dtype[, offset]) |
Put a value into a specified place in a field defined by a data-type. |
setflags ([write, align, uic]) |
Set array flags WRITEABLE, ALIGNED, and UPDATEIFCOPY, respectively. |
sort ([axis, kind, order]) |
Sort an array, in-place. |
squeeze ([axis]) |
Remove single-dimensional entries from the shape of a . |
std ([axis, dtype, out, ddof, keepdims]) |
Returns the standard deviation of the array elements along given axis. |
sum ([axis, dtype, out, keepdims]) |
Return the sum of the array elements over the given axis. |
swapaxes (axis1, axis2) |
Return a view of the array with axis1 and axis2 interchanged. |
take (indices[, axis, out, mode]) |
Return an array formed from the elements of a at the given indices. |
to (unit[, equivalencies]) |
Returns a new Quantity object with the specified units. |
to_dqflag ([name, minlen, dtype, round, ...]) |
Convert this StateTimeSeries into a |
to_hdf5 (*args, **kwargs) |
Convert this array to a h5py.Dataset . |
to_lal (*args, **kwargs) |
Bogus function inherited from superclass, do not use. |
to_pycbc (*args, **kwargs) |
Convert this TimeSeries into a PyCBC |
tobytes ([order]) |
Construct Python bytes containing the raw data bytes in the array. |
tofile (fid[, sep, format]) |
Write array to a file as text or binary (default). |
tolist () |
Return the array as a (possibly nested) list. |
tostring ([order]) |
Construct Python bytes containing the raw data bytes in the array. |
trace ([offset, axis1, axis2, dtype, out]) |
Return the sum along diagonals of the array. |
transpose (*axes) |
Returns a view of the array with axes transposed. |
update (other[, inplace]) |
Update this series by appending new data from an other and dropping the same amount of data off the start. |
value_at (x) |
Return the value of this Series at the given xindex value |
var ([axis, dtype, out, ddof, keepdims]) |
Returns the variance of the array elements, along given axis. |
view ([dtype, type]) |
New view of array with the same data. |
write (data, *args, **kwargs) |
Write out data |
zip () |
Zip the xindex and value arrays of this Series |
Attributes Documentation
T
¶Same as self.transpose(), except that self is returned if self.ndim < 2.
Examples
>>> x = np.array([[1.,2.],[3.,4.]])
>>> x
array([[ 1., 2.],
[ 3., 4.]])
>>> x.T
array([[ 1., 3.],
[ 2., 4.]])
>>> x = np.array([1.,2.,3.,4.])
>>> x
array([ 1., 2., 3., 4.])
>>> x.T
array([ 1., 2., 3., 4.])
base
¶Base object if memory is from some other object.
Examples
The base of an array that owns its memory is None:
>>> x = np.array([1,2,3,4])
>>> x.base is None
True
Slicing creates a view, whose memory is shared with x:
>>> y = x[2:]
>>> y.base is x
True
cgs
¶Returns a copy of the current Quantity
instance with CGS units. The
value of the resulting object will be scaled.
ctypes
¶An object to simplify the interaction of the array with the ctypes module.
This attribute creates an object that makes it easier to use arrays when calling shared libraries with the ctypes module. The returned object has, among others, data, shape, and strides attributes (see Notes below) which themselves return ctypes objects that can be used as arguments to a shared library.
Parameters: | None : |
---|---|
Returns: | c : Python object
|
See also
numpy.ctypeslib
Notes
Below are the public attributes of this object which were documented in “Guide to NumPy” (we have omitted undocumented public attributes, as well as documented private attributes):
Be careful using the ctypes attribute - especially on temporary
arrays or arrays constructed on the fly. For example, calling
(a+b).ctypes.data_as(ctypes.c_void_p)
returns a pointer to memory
that is invalid because the array created as (a+b) is deallocated
before the next Python statement. You can avoid this problem using
either c=a+b
or ct=(a+b).ctypes
. In the latter case, ct will
hold a reference to the array until ct is deleted or re-assigned.
If the ctypes module is not available, then the ctypes attribute of array objects still returns something useful, but ctypes objects are not returned and errors may be raised instead. In particular, the object will still have the as parameter attribute which will return an integer equal to the data attribute.
Examples
>>> import ctypes
>>> x
array([[0, 1],
[2, 3]])
>>> x.ctypes.data
30439712
>>> x.ctypes.data_as(ctypes.POINTER(ctypes.c_long))
<ctypes.LP_c_long object at 0x01F01300>
>>> x.ctypes.data_as(ctypes.POINTER(ctypes.c_long)).contents
c_long(0)
>>> x.ctypes.data_as(ctypes.POINTER(ctypes.c_longlong)).contents
c_longlong(4294967296L)
>>> x.ctypes.shape
<numpy.core._internal.c_long_Array_2 object at 0x01FFD580>
>>> x.ctypes.shape_as(ctypes.c_long)
<numpy.core._internal.c_long_Array_2 object at 0x01FCE620>
>>> x.ctypes.strides
<numpy.core._internal.c_long_Array_2 object at 0x01FCE620>
>>> x.ctypes.strides_as(ctypes.c_longlong)
<numpy.core._internal.c_longlong_Array_2 object at 0x01F01300>
data
¶dtype
¶Data-type of the array’s elements.
Parameters: | None : |
---|---|
Returns: | d : numpy dtype object |
See also
Examples
>>> x
array([[0, 1],
[2, 3]])
>>> x.dtype
dtype('int32')
>>> type(x.dtype)
<type 'numpy.dtype'>
duration
¶Duration of this series in seconds
equivalencies
¶A list of equivalencies that will be applied by default during unit conversions.
flags
¶Information about the memory layout of the array.
Notes
The flags
object can be accessed dictionary-like (as in a.flags['WRITEABLE']
),
or by using lowercased attribute names (as in a.flags.writeable
). Short flag
names are only supported in dictionary access.
Only the UPDATEIFCOPY, WRITEABLE, and ALIGNED flags can be changed by
the user, via direct assignment to the attribute or dictionary entry,
or by calling ndarray.setflags
.
The array flags cannot be set arbitrarily:
False
.True
if the data is truly aligned.True
if the array owns its own memory
or the ultimate owner of the memory exposes a writeable buffer
interface or is a string.Arrays can be both C-style and Fortran-style contiguous simultaneously. This is clear for 1-dimensional arrays, but can also be true for higher dimensional arrays.
Even for contiguous arrays a stride for a given dimension
arr.strides[dim]
may be arbitrary if arr.shape[dim] == 1
or the array has no elements.
It does not generally hold that self.strides[-1] == self.itemsize
for C-style contiguous arrays or self.strides[0] == self.itemsize
for
Fortran-style contiguous arrays is true.
Attributes
flat
¶A 1-D iterator over the Quantity array.
This returns a QuantityIterator
instance, which behaves the same
as the flatiter
instance returned by flat
,
and is similar to, but not a subclass of, Python’s built-in iterator
object.
imag
¶The imaginary part of the array.
Examples
>>> x = np.sqrt([1+0j, 0+1j])
>>> x.imag
array([ 0. , 0.70710678])
>>> x.imag.dtype
dtype('float64')
info
¶Container for meta information like name, description, format. This is required when the object is used as a mixin column within a table, but can be used as a general way to store meta information.
isscalar
¶True if the value
of this quantity is a scalar, or False if it
is an array-like object.
Note
This is subtly different from numpy.isscalar
in that
numpy.isscalar
returns False for a zero-dimensional array
(e.g. np.array(1)
), while this is True for quantities,
since quantities cannot represent true numpy scalars.
itemsize
¶Length of one array element in bytes.
Examples
>>> x = np.array([1,2,3], dtype=np.float64)
>>> x.itemsize
8
>>> x = np.array([1,2,3], dtype=np.complex128)
>>> x.itemsize
16
nbytes
¶Total bytes consumed by the elements of the array.
Notes
Does not include memory consumed by non-element attributes of the array object.
Examples
>>> x = np.zeros((3,5,2), dtype=np.complex128)
>>> x.nbytes
480
>>> np.prod(x.shape) * x.itemsize
480
ndim
¶Number of array dimensions.
Examples
>>> x = np.array([1, 2, 3])
>>> x.ndim
1
>>> y = np.zeros((2, 3, 4))
>>> y.ndim
3
real
¶The real part of the array.
See also
numpy.real
Examples
>>> x = np.sqrt([1+0j, 0+1j])
>>> x.real
array([ 1. , 0.70710678])
>>> x.real.dtype
dtype('float64')
sample_rate
¶Data rate for this TimeSeries
in samples per second (Hertz).
This attribute is stored internally by the dx
attribute
Type: | Quantity scalar |
---|
shape
¶Tuple of array dimensions.
Notes
May be used to “reshape” the array, as long as this would not require a change in the total number of elements
Examples
>>> x = np.array([1, 2, 3, 4])
>>> x.shape
(4,)
>>> y = np.zeros((2, 3, 4))
>>> y.shape
(2, 3, 4)
>>> y.shape = (3, 8)
>>> y
array([[ 0., 0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0.]])
>>> y.shape = (3, 6)
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
ValueError: total size of new array must be unchanged
si
¶Returns a copy of the current Quantity
instance with SI units. The
value of the resulting object will be scaled.
size
¶Number of elements in the array.
Equivalent to np.prod(a.shape)
, i.e., the product of the array’s
dimensions.
Examples
>>> x = np.zeros((3, 5, 2), dtype=np.complex128)
>>> x.size
30
>>> np.prod(x.shape)
30
strides
¶Tuple of bytes to step in each dimension when traversing an array.
The byte offset of element (i[0], i[1], ..., i[n])
in an array a
is:
offset = sum(np.array(i) * a.strides)
A more detailed explanation of strides can be found in the “ndarray.rst” file in the NumPy reference guide.
See also
numpy.lib.stride_tricks.as_strided
Notes
Imagine an array of 32-bit integers (each 4 bytes):
x = np.array([[0, 1, 2, 3, 4],
[5, 6, 7, 8, 9]], dtype=np.int32)
This array is stored in memory as 40 bytes, one after the other
(known as a contiguous block of memory). The strides of an array tell
us how many bytes we have to skip in memory to move to the next position
along a certain axis. For example, we have to skip 4 bytes (1 value) to
move to the next column, but 20 bytes (5 values) to get to the same
position in the next row. As such, the strides for the array x
will be
(20, 4)
.
Examples
>>> y = np.reshape(np.arange(2*3*4), (2,3,4))
>>> y
array([[[ 0, 1, 2, 3],
[ 4, 5, 6, 7],
[ 8, 9, 10, 11]],
[[12, 13, 14, 15],
[16, 17, 18, 19],
[20, 21, 22, 23]]])
>>> y.strides
(48, 16, 4)
>>> y[1,1,1]
17
>>> offset=sum(y.strides * np.array((1,1,1)))
>>> offset/y.itemsize
17
>>> x = np.reshape(np.arange(5*6*7*8), (5,6,7,8)).transpose(2,3,1,0)
>>> x.strides
(32, 4, 224, 1344)
>>> i = np.array([3,5,2,2])
>>> offset = sum(i * x.strides)
>>> x[3,5,2,2]
813
>>> offset / x.itemsize
813
times
¶Series of GPS times for each sample
value
¶The numerical value of this quantity.
Methods Documentation
all
(axis=None, out=None, keepdims=False)¶Returns True if all elements evaluate to True.
Refer to numpy.all
for full documentation.
See also
numpy.all
any
(axis=None, out=None, keepdims=False)¶Returns True if any of the elements of a
evaluate to True.
Refer to numpy.any
for full documentation.
See also
numpy.any
append
(other, gap='raise', inplace=True, pad=0.0, resize=True)[source]¶Connect another series onto the end of the current one.
Parameters: | other :
gap :
inplace :
pad :
resize :
|
---|---|
Returns: | series :
|
argmax
(axis=None, out=None)¶Return indices of the maximum values along the given axis.
Refer to numpy.argmax
for full documentation.
See also
numpy.argmax
argmin
(axis=None, out=None)¶Return indices of the minimum values along the given axis of a
.
Refer to numpy.argmin
for detailed documentation.
See also
numpy.argmin
argpartition
(kth, axis=-1, kind='introselect', order=None)¶Returns the indices that would partition this array.
Refer to numpy.argpartition
for full documentation.
New in version 1.8.0.
See also
numpy.argpartition
argsort
(axis=-1, kind='quicksort', order=None)¶Returns the indices that would sort this array.
Refer to numpy.argsort
for full documentation.
See also
numpy.argsort
astype
(dtype, order='K', casting='unsafe', subok=True, copy=True)¶Copy of the array, cast to a specified type.
Parameters: | dtype : str or dtype
order : {‘C’, ‘F’, ‘A’, ‘K’}, optional
casting : {‘no’, ‘equiv’, ‘safe’, ‘same_kind’, ‘unsafe’}, optional
subok : bool, optional
copy : bool, optional
|
---|---|
Returns: | arr_t : ndarray |
Raises: | ComplexWarning :
|
Notes
Starting in NumPy 1.9, astype method now returns an error if the string dtype to cast to is not long enough in ‘safe’ casting mode to hold the max value of integer/float array that is being casted. Previously the casting was allowed even if the result was truncated.
Examples
>>> x = np.array([1, 2, 2.5])
>>> x
array([ 1. , 2. , 2.5])
>>> x.astype(int)
array([1, 2, 2])
byteswap
(inplace)¶Swap the bytes of the array elements
Toggle between low-endian and big-endian data representation by returning a byteswapped array, optionally swapped in-place.
Parameters: | inplace : bool, optional
|
---|---|
Returns: | out : ndarray
|
Examples
>>> A = np.array([1, 256, 8755], dtype=np.int16)
>>> map(hex, A)
['0x1', '0x100', '0x2233']
>>> A.byteswap(True)
array([ 256, 1, 13090], dtype=int16)
>>> map(hex, A)
['0x100', '0x1', '0x3322']
Arrays of strings are not swapped
>>> A = np.array(['ceg', 'fac'])
>>> A.byteswap()
array(['ceg', 'fac'],
dtype='|S3')
choose
(choices, out=None, mode='raise')¶Use an index array to construct a new array from a set of choices.
Refer to numpy.choose
for full documentation.
See also
numpy.choose
clip
(min=None, max=None, out=None)¶Return an array whose values are limited to [min, max]
.
One of max or min must be given.
Refer to numpy.clip
for full documentation.
See also
numpy.clip
compress
(condition, axis=None, out=None)¶Return selected slices of this array along given axis.
Refer to numpy.compress
for full documentation.
See also
numpy.compress
conj
()¶Complex-conjugate all elements.
Refer to numpy.conjugate
for full documentation.
See also
numpy.conjugate
conjugate
()¶Return the complex conjugate, element-wise.
Refer to numpy.conjugate
for full documentation.
See also
numpy.conjugate
copy
(order='C')[source]¶Return a copy of the array.
Parameters: | order : {‘C’, ‘F’, ‘A’, ‘K’}, optional
|
---|
See also
Examples
>>> x = np.array([[1,2,3],[4,5,6]], order='F')
>>> y = x.copy()
>>> x.fill(0)
>>> x
array([[0, 0, 0],
[0, 0, 0]])
>>> y
array([[1, 2, 3],
[4, 5, 6]])
>>> y.flags['C_CONTIGUOUS']
True
crop
(start=None, end=None, copy=False)[source]¶Crop this series to the given x-axis extent.
Parameters: | start :
end :
copy :
|
---|---|
Returns: | series :
|
Notes
If either start
or end
are outside of the original
Series
span, warnings will be printed and the limits will
be restricted to the xspan
cumprod
(axis=None, dtype=None, out=None)¶Return the cumulative product of the elements along the given axis.
Refer to numpy.cumprod
for full documentation.
See also
numpy.cumprod
cumsum
(axis=None, dtype=None, out=None)¶Return the cumulative sum of the elements along the given axis.
Refer to numpy.cumsum
for full documentation.
See also
numpy.cumsum
decompose
(bases=[])¶Generates a new Quantity
with the units
decomposed. Decomposed units have only irreducible units in
them (see astropy.units.UnitBase.decompose
).
Parameters: | bases : sequence of UnitBase, optional
|
---|---|
Returns: | newq :
|
diagonal
(offset=0, axis1=0, axis2=1)¶Return specified diagonals. In NumPy 1.9 the returned array is a read-only view instead of a copy as in previous NumPy versions. In a future version the read-only restriction will be removed.
Refer to numpy.diagonal()
for full documentation.
See also
numpy.diagonal
diff
(n=1, axis=-1)[source]¶Calculate the n-th discrete difference along given axis.
The first difference is given by
out[n] = a[n+1] - a[n]
along the given axis, higher differences are calculated by usingdiff
recursively.
Parameters: | a : array_like
|
---|---|
Returns: | diff : ndarray
. : |
Examples
>>> x = np.array([1, 2, 4, 7, 0])
>>> np.diff(x)
array([ 1, 2, 3, -7])
>>> np.diff(x, n=2)
array([ 1, 1, -10])
>>> x = np.array([[1, 3, 6, 10], [0, 5, 6, 8]])
>>> np.diff(x)
array([[2, 3, 4],
[5, 1, 2]])
>>> np.diff(x, axis=0)
array([[-1, 2, 0, -2]])
dot
(b, out=None)¶Dot product of two arrays.
Refer to numpy.dot
for full documentation.
See also
numpy.dot
Examples
>>> a = np.eye(2)
>>> b = np.ones((2, 2)) * 2
>>> a.dot(b)
array([[ 2., 2.],
[ 2., 2.]])
This array method can be conveniently chained:
>>> a.dot(b).dot(b)
array([[ 8., 8.],
[ 8., 8.]])
dump
(file)¶Dump a pickle of the array to the specified file. The array can be read back with pickle.load or numpy.load.
Parameters: | file : str
|
---|
dumps
()[source]¶Returns the pickle of the array as a string. pickle.loads or numpy.loads will convert the string back to an array.
Parameters: | None : |
---|
ediff1d
(to_end=None, to_begin=None)¶fetch
(*args, **kwargs)[source]¶Fetch data from NDS into a TimeSeries
.
Parameters: |
start :
end :
host :
port :
verify :
connection :
verbose :
type :
dtype :
|
---|
fetch_open_data
(ifo, start, end, name='strain/Strain', sample_rate=4096, host='https://losc.ligo.org')[source]¶Fetch open-access data from the LIGO Open Science Center
Parameters: | ifo :
start :
end :
name :
sample_rate :
host :
|
---|
fill
(value)¶Fill the array with a scalar value.
Parameters: | value : scalar
|
---|
Examples
>>> a = np.array([1, 2])
>>> a.fill(0)
>>> a
array([0, 0])
>>> a = np.empty(2)
>>> a.fill(1)
>>> a
array([ 1., 1.])
find
(channel, start, end, frametype=None, pad=None, dtype=None, nproc=1, verbose=False, **readargs)[source]¶Find and read data from frames for a channel
Parameters: |
start :
end :
frametype :
pad :
nproc :
dtype : verbose :
**readargs :
|
---|
flatten
(order='C')¶Return a copy of the array collapsed into one dimension.
Parameters: | order : {‘C’, ‘F’, ‘A’, ‘K’}, optional
|
---|---|
Returns: | y : ndarray
|
Examples
>>> a = np.array([[1,2], [3,4]])
>>> a.flatten()
array([1, 2, 3, 4])
>>> a.flatten('F')
array([1, 3, 2, 4])
from_hdf5
(*args, **kwargs)[source]¶Read an array from the given HDF file.
This method has been deprecated in favour of the unified I/O method:
Class.write(..., format=’hdf5’)
from_nds2_buffer
(*args, **kwargs)[source]¶Construct a new TimeSeries
from an nds2.buffer
object
Parameters: | buffer_ :
**metadata :
|
---|---|
Returns: | timeseries :
|
Notes
This classmethod requires the nds2-client package
from_pycbc
(ts)[source]¶Convert a pycbc.types.timeseries.TimeSeries
into a TimeSeries
Parameters: | ts :
|
---|---|
Returns: | timeseries :
|
get
(channel, start, end, pad=None, dtype=None, verbose=False, **kwargs)[source]¶Get data for this channel from frames or NDS
This method dynamically accesses either frames on disk, or a remote NDS2 server to find and return data for the given interval
Parameters: |
start :
end :
pad :
dtype : nproc :
verbose :
**kwargs other keyword arguments to pass to either : |
---|
See also
TimeSeries.fetch
TimeSeries.find
getfield
(dtype, offset=0)¶Returns a field of the given array as a certain type.
A field is a view of the array data with a given data-type. The values in the view are determined by the given type and the offset into the current array in bytes. The offset needs to be such that the view dtype fits in the array dtype; for example an array of dtype complex128 has 16-byte elements. If taking a view with a 32-bit integer (4 bytes), the offset needs to be between 0 and 12 bytes.
Parameters: | dtype : str or dtype
offset : int
|
---|
Examples
>>> x = np.diag([1.+1.j]*2)
>>> x[1, 1] = 2 + 4.j
>>> x
array([[ 1.+1.j, 0.+0.j],
[ 0.+0.j, 2.+4.j]])
>>> x.getfield(np.float64)
array([[ 1., 0.],
[ 0., 2.]])
By choosing an offset of 8 bytes we can select the complex part of the array for our view:
>>> x.getfield(np.float64, offset=8)
array([[ 1., 0.],
[ 0., 4.]])
insert
(obj, values, axis=None)¶Insert values along the given axis before the given indices and return
a new Quantity
object.
This is a thin wrapper around the numpy.insert
function.
Parameters: | obj : int, slice or sequence of ints
values : array-like
axis : int, optional
|
---|---|
Returns: | out :
|
Examples
>>> import astropy.units as u
>>> q = [1, 2] * u.m
>>> q.insert(0, 50 * u.cm)
<Quantity [ 0.5, 1., 2.] m>
>>> q = [[1, 2], [3, 4]] * u.m
>>> q.insert(1, [10, 20] * u.m, axis=0)
<Quantity [[ 1., 2.],
[ 10., 20.],
[ 3., 4.]] m>
>>> q.insert(1, 10 * u.m, axis=1)
<Quantity [[ 1., 10., 2.],
[ 3., 10., 4.]] m>
is_compatible
(other)[source]¶Check whether this series and other have compatible metadata
This method tests that the sample size
, and the
unit
match.
is_contiguous
(other, tol=3.814697265625e-06)[source]¶Check whether other is contiguous with self.
Parameters: | other :
tol :
|
---|---|
Returns: | 1 :
-1 :
0 :
|
Notes
if a raw numpy.ndarray
is passed as other, with no metadata, then
the contiguity check will always pass
item
(*args)¶Copy an element of an array to a standard Python scalar and return it.
Parameters: | *args : Arguments (variable number and type)
|
---|---|
Returns: | z : Standard Python scalar object
|
Notes
When the data type of a
is longdouble or clongdouble, item() returns
a scalar array object because there is no available Python scalar that
would not lose information. Void arrays return a buffer object for item(),
unless fields are defined, in which case a tuple is returned.
item
is very similar to a[args], except, instead of an array scalar,
a standard Python scalar is returned. This can be useful for speeding up
access to elements of the array and doing arithmetic on elements of the
array using Python’s optimized math.
Examples
>>> x = np.random.randint(9, size=(3, 3))
>>> x
array([[3, 1, 7],
[2, 8, 3],
[8, 5, 3]])
>>> x.item(3)
2
>>> x.item(7)
5
>>> x.item((0, 1))
1
>>> x.item((2, 2))
3
itemset
(*args)¶Insert scalar into an array (scalar is cast to array’s dtype, if possible)
There must be at least 1 argument, and define the last argument
as item. Then, a.itemset(*args)
is equivalent to but faster
than a[args] = item
. The item should be a scalar value and args
must select a single item in the array a
.
Parameters: | *args : Arguments
|
---|
Notes
Compared to indexing syntax, itemset
provides some speed increase
for placing a scalar into a particular location in an ndarray
,
if you must do this. However, generally this is discouraged:
among other problems, it complicates the appearance of the code.
Also, when using itemset
(and item
) inside a loop, be sure
to assign the methods to a local variable to avoid the attribute
look-up at each loop iteration.
Examples
>>> x = np.random.randint(9, size=(3, 3))
>>> x
array([[3, 1, 7],
[2, 8, 3],
[8, 5, 3]])
>>> x.itemset(4, 0)
>>> x.itemset((2, 2), 9)
>>> x
array([[3, 1, 7],
[2, 0, 3],
[8, 5, 9]])
max
(axis=None, out=None)¶Return the maximum along a given axis.
Refer to numpy.amax
for full documentation.
See also
numpy.amax
mean
(axis=None, dtype=None, out=None, keepdims=False)¶Returns the average of the array elements along given axis.
Refer to numpy.mean
for full documentation.
See also
numpy.mean
median
(axis=None, **kwargs)[source]¶Compute the median along the specified axis.
Returns the median of the array elements.
Parameters: | a : array_like
axis : {int, sequence of int, None}, optional
out : ndarray, optional
overwrite_input : bool, optional
keepdims : bool, optional
|
---|---|
Returns: | median : ndarray
|
See also
mean
, percentile
Notes
Given a vector V
of length N
, the median of V
is the
middle value of a sorted copy of V
, V_sorted
- i
e., V_sorted[(N-1)/2]
, when N
is odd, and the average of the
two middle values of V_sorted
when N
is even.
Examples
>>> a = np.array([[10, 7, 4], [3, 2, 1]])
>>> a
array([[10, 7, 4],
[ 3, 2, 1]])
>>> np.median(a)
3.5
>>> np.median(a, axis=0)
array([ 6.5, 4.5, 2.5])
>>> np.median(a, axis=1)
array([ 7., 2.])
>>> m = np.median(a, axis=0)
>>> out = np.zeros_like(m)
>>> np.median(a, axis=0, out=m)
array([ 6.5, 4.5, 2.5])
>>> m
array([ 6.5, 4.5, 2.5])
>>> b = a.copy()
>>> np.median(b, axis=1, overwrite_input=True)
array([ 7., 2.])
>>> assert not np.all(a==b)
>>> b = a.copy()
>>> np.median(b, axis=None, overwrite_input=True)
3.5
>>> assert not np.all(a==b)
min
(axis=None, out=None, keepdims=False)¶Return the minimum along a given axis.
Refer to numpy.amin
for full documentation.
See also
numpy.amin
nansum
(axis=None, out=None, keepdims=False)¶newbyteorder
(new_order='S')¶Return the array with the same data viewed with a different byte order.
Equivalent to:
arr.view(arr.dtype.newbytorder(new_order))
Changes are also made in all fields and sub-arrays of the array data type.
Parameters: | new_order : string, optional
|
---|---|
Returns: | new_arr : array
|
nonzero
()¶Return the indices of the elements that are non-zero.
Refer to numpy.nonzero
for full documentation.
See also
numpy.nonzero
override_unit
(unit, parse_strict='raise')[source]¶Forcefully reset the unit of these data
Use of this method is discouraged in favour of to()
,
which performs accurate conversions from one unit to another.
The method should really only be used when the original unit of the
array is plain wrong.
Parameters: |
|
---|---|
Raises: | ValueError :
|
pad
(pad_width, **kwargs)[source]¶Pad this series to a new size
Parameters: | pad_width :
**kwargs :
|
---|---|
Returns: | series :
|
See also
numpy.pad
partition
(kth, axis=-1, kind='introselect', order=None)¶Rearranges the elements in the array in such a way that value of the element in kth position is in the position it would be in a sorted array. All elements smaller than the kth element are moved before this element and all equal or greater are moved behind it. The ordering of the elements in the two partitions is undefined.
New in version 1.8.0.
Parameters: | kth : int or sequence of ints
axis : int, optional
kind : {‘introselect’}, optional
order : str or list of str, optional
|
---|
See also
numpy.partition
argpartition
sort
Notes
See np.partition
for notes on the different algorithms.
Examples
>>> a = np.array([3, 4, 2, 1])
>>> a.partition(a, 3)
>>> a
array([2, 1, 3, 4])
>>> a.partition((1, 3))
array([1, 2, 3, 4])
plot
(**kwargs)[source]¶Plot the data for this TimeSeries
prepend
(other, gap='raise', inplace=True, pad=0.0, resize=True)[source]¶Connect another series onto the start of the current one.
Parameters: | other :
gap :
inplace :
pad :
|
---|---|
Returns: | series :
|
prod
(axis=None, dtype=None, out=None, keepdims=False)¶Return the product of the array elements over the given axis
Refer to numpy.prod
for full documentation.
See also
numpy.prod
ptp
(axis=None, out=None)¶Peak to peak (maximum - minimum) value along a given axis.
Refer to numpy.ptp
for full documentation.
See also
numpy.ptp
put
(indices, values, mode='raise')¶Set a.flat[n] = values[n]
for all n
in indices.
Refer to numpy.put
for full documentation.
See also
numpy.put
ravel
([order])¶Return a flattened array.
Refer to numpy.ravel
for full documentation.
See also
numpy.ravel
ndarray.flat
read
(*args, **kwargs)¶Read in data
The arguments passed to this method depend on the format
The available built-in formats are:
Format | Read | Write | Auto-identify |
---|---|---|---|
hdf | Yes | No | No |
hdf5 | Yes | Yes | Yes |
repeat
(repeats, axis=None)¶Repeat elements of an array.
Refer to numpy.repeat
for full documentation.
See also
numpy.repeat
reshape
(shape, order='C')¶Returns an array containing the same data with a new shape.
Refer to numpy.reshape
for full documentation.
See also
numpy.reshape
resize
(new_shape, refcheck=True)¶Change shape and size of array in-place.
Parameters: | new_shape : tuple of ints, or
refcheck : bool, optional
|
---|---|
Returns: | None : |
Raises: | ValueError :
SystemError :
|
See also
resize
Notes
This reallocates space for the data area if necessary.
Only contiguous arrays (data elements consecutive in memory) can be resized.
The purpose of the reference count check is to make sure you
do not use this array as a buffer for another Python object and then
reallocate the memory. However, reference counts can increase in
other ways so if you are sure that you have not shared the memory
for this array with another Python object, then you may safely set
refcheck
to False.
Examples
Shrinking an array: array is flattened (in the order that the data are stored in memory), resized, and reshaped:
>>> a = np.array([[0, 1], [2, 3]], order='C')
>>> a.resize((2, 1))
>>> a
array([[0],
[1]])
>>> a = np.array([[0, 1], [2, 3]], order='F')
>>> a.resize((2, 1))
>>> a
array([[0],
[2]])
Enlarging an array: as above, but missing entries are filled with zeros:
>>> b = np.array([[0, 1], [2, 3]])
>>> b.resize(2, 3) # new_shape parameter doesn't have to be a tuple
>>> b
array([[0, 1, 2],
[3, 0, 0]])
Referencing an array prevents resizing...
>>> c = a
>>> a.resize((1, 1))
Traceback (most recent call last):
...
ValueError: cannot resize an array that has been referenced ...
Unless refcheck
is False:
>>> a.resize((1, 1), refcheck=False)
>>> a
array([[0]])
>>> c
array([[0]])
round
(decimals=0, out=None)¶Return a
with each element rounded to the given number of decimals.
Refer to numpy.around
for full documentation.
See also
numpy.around
searchsorted
(v, side='left', sorter=None)¶Find indices where elements of v should be inserted in a to maintain order.
For full documentation, see numpy.searchsorted
See also
numpy.searchsorted
setfield
(val, dtype, offset=0)¶Put a value into a specified place in a field defined by a data-type.
Place val
into a
‘s field defined by dtype
and beginning offset
bytes into the field.
Parameters: | val : object
dtype : dtype object
offset : int, optional
|
---|---|
Returns: | None : |
See also
Examples
>>> x = np.eye(3)
>>> x.getfield(np.float64)
array([[ 1., 0., 0.],
[ 0., 1., 0.],
[ 0., 0., 1.]])
>>> x.setfield(3, np.int32)
>>> x.getfield(np.int32)
array([[3, 3, 3],
[3, 3, 3],
[3, 3, 3]])
>>> x
array([[ 1.00000000e+000, 1.48219694e-323, 1.48219694e-323],
[ 1.48219694e-323, 1.00000000e+000, 1.48219694e-323],
[ 1.48219694e-323, 1.48219694e-323, 1.00000000e+000]])
>>> x.setfield(np.eye(3), np.int32)
>>> x
array([[ 1., 0., 0.],
[ 0., 1., 0.],
[ 0., 0., 1.]])
setflags
(write=None, align=None, uic=None)¶Set array flags WRITEABLE, ALIGNED, and UPDATEIFCOPY, respectively.
These Boolean-valued flags affect how numpy interprets the memory
area used by a
(see Notes below). The ALIGNED flag can only
be set to True if the data is actually aligned according to the type.
The UPDATEIFCOPY flag can never be set to True. The flag WRITEABLE
can only be set to True if the array owns its own memory, or the
ultimate owner of the memory exposes a writeable buffer interface,
or is a string. (The exception for string is made so that unpickling
can be done without copying memory.)
Parameters: | write : bool, optional
align : bool, optional
uic : bool, optional
|
---|
Notes
Array flags provide information about how the memory area used for the array is to be interpreted. There are 6 Boolean flags in use, only three of which can be changed by the user: UPDATEIFCOPY, WRITEABLE, and ALIGNED.
WRITEABLE (W) the data area can be written to;
ALIGNED (A) the data and strides are aligned appropriately for the hardware (as determined by the compiler);
UPDATEIFCOPY (U) this array is a copy of some other array (referenced by .base). When this array is deallocated, the base array will be updated with the contents of this array.
All flags can be accessed using their first (upper case) letter as well as the full name.
Examples
>>> y
array([[3, 1, 7],
[2, 0, 0],
[8, 5, 9]])
>>> y.flags
C_CONTIGUOUS : True
F_CONTIGUOUS : False
OWNDATA : True
WRITEABLE : True
ALIGNED : True
UPDATEIFCOPY : False
>>> y.setflags(write=0, align=0)
>>> y.flags
C_CONTIGUOUS : True
F_CONTIGUOUS : False
OWNDATA : True
WRITEABLE : False
ALIGNED : False
UPDATEIFCOPY : False
>>> y.setflags(uic=1)
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
ValueError: cannot set UPDATEIFCOPY flag to True
sort
(axis=-1, kind='quicksort', order=None)¶Sort an array, in-place.
Parameters: | axis : int, optional
kind : {‘quicksort’, ‘mergesort’, ‘heapsort’}, optional
order : str or list of str, optional
|
---|
See also
numpy.sort
argsort
lexsort
searchsorted
partition
Notes
See sort
for notes on the different sorting algorithms.
Examples
>>> a = np.array([[1,4], [3,1]])
>>> a.sort(axis=1)
>>> a
array([[1, 4],
[1, 3]])
>>> a.sort(axis=0)
>>> a
array([[1, 3],
[1, 4]])
Use the order
keyword to specify a field to use when sorting a
structured array:
>>> a = np.array([('a', 2), ('c', 1)], dtype=[('x', 'S1'), ('y', int)])
>>> a.sort(order='y')
>>> a
array([('c', 1), ('a', 2)],
dtype=[('x', '|S1'), ('y', '<i4')])
squeeze
(axis=None)¶Remove single-dimensional entries from the shape of a
.
Refer to numpy.squeeze
for full documentation.
See also
numpy.squeeze
std
(axis=None, dtype=None, out=None, ddof=0, keepdims=False)¶Returns the standard deviation of the array elements along given axis.
Refer to numpy.std
for full documentation.
See also
numpy.std
sum
(axis=None, dtype=None, out=None, keepdims=False)¶Return the sum of the array elements over the given axis.
Refer to numpy.sum
for full documentation.
See also
numpy.sum
swapaxes
(axis1, axis2)¶Return a view of the array with axis1
and axis2
interchanged.
Refer to numpy.swapaxes
for full documentation.
See also
numpy.swapaxes
take
(indices, axis=None, out=None, mode='raise')¶Return an array formed from the elements of a
at the given indices.
Refer to numpy.take
for full documentation.
See also
numpy.take
to
(unit, equivalencies=[])¶Returns a new Quantity
object with the specified
units.
Parameters: | unit : equivalencies : list of equivalence pairs, optional
|
---|
to_dqflag
(name=None, minlen=1, dtype=<type 'float'>, round=False, label=None, description=None)[source]¶Convert this StateTimeSeries
into a
DataQualityFlag
Each contiguous set of True
values are grouped as a
Segment
running from the GPS time the first
found True
, to the GPS time of the next False
(or the end
of the series)
Parameters: | minlen : dtype :
round :
|
---|---|
Returns: | dqflag :
|
to_hdf5
(*args, **kwargs)[source]¶Convert this array to a h5py.Dataset
.
This method has been deprecated in favour of the unified I/O method:
Class.write(..., format=’hdf5’)
to_pycbc
(*args, **kwargs)[source]¶Convert this TimeSeries
into a PyCBC
TimeSeries
Parameters: | copy :
|
---|---|
Returns: | timeseries :
|
tobytes
(order='C')¶Construct Python bytes containing the raw data bytes in the array.
Constructs Python bytes showing a copy of the raw contents of data memory. The bytes object can be produced in either ‘C’ or ‘Fortran’, or ‘Any’ order (the default is ‘C’-order). ‘Any’ order means C-order unless the F_CONTIGUOUS flag in the array is set, in which case it means ‘Fortran’ order.
New in version 1.9.0.
Parameters: | order : {‘C’, ‘F’, None}, optional
|
---|---|
Returns: | s : bytes
|
Examples
>>> x = np.array([[0, 1], [2, 3]])
>>> x.tobytes()
b'\x00\x00\x00\x00\x01\x00\x00\x00\x02\x00\x00\x00\x03\x00\x00\x00'
>>> x.tobytes('C') == x.tobytes()
True
>>> x.tobytes('F')
b'\x00\x00\x00\x00\x02\x00\x00\x00\x01\x00\x00\x00\x03\x00\x00\x00'
tofile
(fid, sep="", format="%s")¶Write array to a file as text or binary (default).
Data is always written in ‘C’ order, independent of the order of a
.
The data produced by this method can be recovered using the function
fromfile().
Parameters: | fid : file or str
sep : str
format : str
|
---|
Notes
This is a convenience function for quick storage of array data. Information on endianness and precision is lost, so this method is not a good choice for files intended to archive data or transport data between machines with different endianness. Some of these problems can be overcome by outputting the data as text files, at the expense of speed and file size.
tolist
()¶Return the array as a (possibly nested) list.
Return a copy of the array data as a (nested) Python list. Data items are converted to the nearest compatible Python type.
Parameters: | none : |
---|---|
Returns: | y : list
|
Notes
The array may be recreated, a = np.array(a.tolist())
.
Examples
>>> a = np.array([1, 2])
>>> a.tolist()
[1, 2]
>>> a = np.array([[1, 2], [3, 4]])
>>> list(a)
[array([1, 2]), array([3, 4])]
>>> a.tolist()
[[1, 2], [3, 4]]
tostring
(order='C')[source]¶Construct Python bytes containing the raw data bytes in the array.
Constructs Python bytes showing a copy of the raw contents of data memory. The bytes object can be produced in either ‘C’ or ‘Fortran’, or ‘Any’ order (the default is ‘C’-order). ‘Any’ order means C-order unless the F_CONTIGUOUS flag in the array is set, in which case it means ‘Fortran’ order.
This function is a compatibility alias for tobytes. Despite its name it returns bytes not strings.
Parameters: | order : {‘C’, ‘F’, None}, optional
|
---|---|
Returns: | s : bytes
|
Examples
>>> x = np.array([[0, 1], [2, 3]])
>>> x.tobytes()
b'\x00\x00\x00\x00\x01\x00\x00\x00\x02\x00\x00\x00\x03\x00\x00\x00'
>>> x.tobytes('C') == x.tobytes()
True
>>> x.tobytes('F')
b'\x00\x00\x00\x00\x02\x00\x00\x00\x01\x00\x00\x00\x03\x00\x00\x00'
trace
(offset=0, axis1=0, axis2=1, dtype=None, out=None)¶Return the sum along diagonals of the array.
Refer to numpy.trace
for full documentation.
See also
numpy.trace
transpose
(*axes)¶Returns a view of the array with axes transposed.
For a 1-D array, this has no effect. (To change between column and
row vectors, first cast the 1-D array into a matrix object.)
For a 2-D array, this is the usual matrix transpose.
For an n-D array, if axes are given, their order indicates how the
axes are permuted (see Examples). If axes are not provided and
a.shape = (i[0], i[1], ... i[n-2], i[n-1])
, then
a.transpose().shape = (i[n-1], i[n-2], ... i[1], i[0])
.
Parameters: | axes : None, tuple of ints, or
|
---|---|
Returns: | out : ndarray
|
See also
ndarray.T
Examples
>>> a = np.array([[1, 2], [3, 4]])
>>> a
array([[1, 2],
[3, 4]])
>>> a.transpose()
array([[1, 3],
[2, 4]])
>>> a.transpose((1, 0))
array([[1, 3],
[2, 4]])
>>> a.transpose(1, 0)
array([[1, 3],
[2, 4]])
update
(other, inplace=True)[source]¶Update this series by appending new data from an other and dropping the same amount of data off the start.
This is a convenience method that just calls append
with
resize=False
.
value_at
(x)[source]¶Return the value of this Series
at the given xindex
value
Parameters: |
|
---|---|
Returns: | y :
|
var
(axis=None, dtype=None, out=None, ddof=0, keepdims=False)¶Returns the variance of the array elements, along given axis.
Refer to numpy.var
for full documentation.
See also
numpy.var
view
(dtype=None, type=None)¶New view of array with the same data.
Parameters: | dtype : data-type or ndarray sub-class, optional
type : Python type, optional
|
---|
Notes
a.view()
is used two different ways:
a.view(some_dtype)
or a.view(dtype=some_dtype)
constructs a view
of the array’s memory with a different data-type. This can cause a
reinterpretation of the bytes of memory.
a.view(ndarray_subclass)
or a.view(type=ndarray_subclass)
just
returns an instance of ndarray_subclass
that looks at the same array
(same shape, dtype, etc.) This does not cause a reinterpretation of the
memory.
For a.view(some_dtype)
, if some_dtype
has a different number of
bytes per entry than the previous dtype (for example, converting a
regular array to a structured array), then the behavior of the view
cannot be predicted just from the superficial appearance of a
(shown
by print(a)
). It also depends on exactly how a
is stored in
memory. Therefore if a
is C-ordered versus fortran-ordered, versus
defined as a slice or transpose, etc., the view may give different
results.
Examples
>>> x = np.array([(1, 2)], dtype=[('a', np.int8), ('b', np.int8)])
Viewing array data using a different type and dtype:
>>> y = x.view(dtype=np.int16, type=np.matrix)
>>> y
matrix([[513]], dtype=int16)
>>> print(type(y))
<class 'numpy.matrixlib.defmatrix.matrix'>
Creating a view on a structured array so it can be used in calculations
>>> x = np.array([(1, 2),(3,4)], dtype=[('a', np.int8), ('b', np.int8)])
>>> xv = x.view(dtype=np.int8).reshape(-1,2)
>>> xv
array([[1, 2],
[3, 4]], dtype=int8)
>>> xv.mean(0)
array([ 2., 3.])
Making changes to the view changes the underlying array
>>> xv[0,1] = 20
>>> print(x)
[(1, 20) (3, 4)]
Using a view to convert an array to a recarray:
>>> z = x.view(np.recarray)
>>> z.a
array([1], dtype=int8)
Views share data:
>>> x[0] = (9, 10)
>>> z[0]
(9, 10)
Views that change the dtype size (bytes per entry) should normally be avoided on arrays defined by slices, transposes, fortran-ordering, etc.:
>>> x = np.array([[1,2,3],[4,5,6]], dtype=np.int16)
>>> y = x[:, 0:2]
>>> y
array([[1, 2],
[4, 5]], dtype=int16)
>>> y.view(dtype=[('width', np.int16), ('length', np.int16)])
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
ValueError: new type not compatible with array.
>>> z = y.copy()
>>> z.view(dtype=[('width', np.int16), ('length', np.int16)])
array([[(1, 2)],
[(4, 5)]], dtype=[('width', '<i2'), ('length', '<i2')])
write
(data, *args, **kwargs)¶Write out data
The arguments passed to this method depend on the format
The available built-in formats are:
Format | Read | Write | Auto-identify |
---|---|---|---|
hdf | Yes | Yes | No |
hdf5 | Yes | Yes | Yes |
zip
()[source]¶Zip the xindex
and value
arrays of this Series
Returns: | stacked : 2-d |
---|
Examples
>>> a = Series([0, 2, 4, 6, 8], xindex=[-5, -4, -3, -2, -1])
>>> a.zip()
array([[-5., 0.],
[-4., 2.],
[-3., 4.],
[-2., 6.],
[-1., 8.]])
gwpy.timeseries.
StateVectorDict
(*args, **kwds)[source]¶Bases: gwpy.timeseries.core.TimeSeriesBaseDict
Methods Summary
append (other[, copy]) |
|
clear (() -> None. Remove all items from od.) |
|
copy () |
|
crop ([start, end, copy]) |
Crop each entry of this TimeSeriesBaseDict . |
fetch (*args, **kwargs) |
Fetch data from NDS for a number of channels. |
find (channels, start, end[, frametype, pad, ...]) |
Find and read data from frames for a number of channels. |
fromkeys ((S[, ...) |
If not specified, the value defaults to None. |
get (channels, start, end[, pad, dtype, ...]) |
Retrieve data for multiple channels from frames or NDS |
has_key ((k) -> True if D has a key k, else False) |
|
items (() -> list of (key, value) pairs in od) |
|
iteritems () |
od.iteritems -> an iterator over the (key, value) pairs in od |
iterkeys (() -> an iterator over the keys in od) |
|
itervalues () |
od.itervalues -> an iterator over the values in od |
keys (() -> list of keys in od) |
|
plot ([label]) |
Plot the data for this TimeSeriesBaseDict . |
pop ((k[,d]) -> v, ...) |
value. If key is not found, d is returned if given, otherwise KeyError |
popitem (() -> (k, v), ...) |
Pairs are returned in LIFO order if last is true or FIFO order if false. |
prepend (other, **kwargs) |
|
read (*args, **kwargs) |
Read data for multiple bit vector channels into a StateVectorDict |
resample (rate, **kwargs) |
Resample items in this dict. |
setdefault ((k[,d]) -> od.get(k,d), ...) |
|
update (([E, ...) |
If E present and has a .keys() method, does: for k in E: D[k] = E[k] |
values (() -> list of values in od) |
|
viewitems (...) |
|
viewkeys (...) |
|
viewvalues (...) |
|
write (data, *args, **kwargs) |
Write this TimeSeriesDict to a file |
Methods Documentation
clear
() → None. Remove all items from od.¶crop
(start=None, end=None, copy=False)[source]¶Crop each entry of this TimeSeriesBaseDict
.
This method calls the crop()
method of all entries and
modifies this dict in place.
Parameters: | start :
end :
|
---|
See also
TimeSeries.crop
fetch
(*args, **kwargs)[source]¶Fetch data from NDS for a number of channels.
Parameters: | channels :
start :
end :
host :
port :
verify :
verbose :
connection :
dtype :
|
---|---|
Returns: | data :
|
find
(channels, start, end, frametype=None, pad=None, dtype=None, nproc=1, verbose=False, allow_tape=True, observatory=None, **readargs)[source]¶Find and read data from frames for a number of channels.
Parameters: | channels :
start :
end :
frametype :
pad :
dtype :
nproc :
verbose :
allow_tape :
**readargs :
|
---|
fromkeys
(S[, v]) → New ordered dictionary with keys from S.¶If not specified, the value defaults to None.
get
(channels, start, end, pad=None, dtype=None, verbose=False, allow_tape=False, **kwargs)[source]¶Retrieve data for multiple channels from frames or NDS
This method dynamically accesses either frames on disk, or a remote NDS2 server to find and return data for the given interval
Parameters: | channels :
start :
end :
frametype :
pad :
dtype :
nproc :
verbose :
allow_tape :
**kwargs :
|
---|
has_key
(k) → True if D has a key k, else False¶items
() → list of (key, value) pairs in od¶iteritems
()¶od.iteritems -> an iterator over the (key, value) pairs in od
iterkeys
() → an iterator over the keys in od¶itervalues
()¶od.itervalues -> an iterator over the values in od
keys
() → list of keys in od¶plot
(label='key', **kwargs)[source]¶Plot the data for this TimeSeriesBaseDict
.
Parameters: | label :
**kwargs :
|
---|
pop
(k[, d]) → v, remove specified key and return the corresponding¶value. If key is not found, d is returned if given, otherwise KeyError is raised.
popitem
() → (k, v), return and remove a (key, value) pair.¶Pairs are returned in LIFO order if last is true or FIFO order if false.
read
(*args, **kwargs)¶Read data for multiple bit vector channels into a StateVectorDict
Parameters: |
channels :
start :
end :
bits : format :
nproc :
gap :
pad :
|
---|---|
Returns: | statevectordict :
|
Notes
The available built-in formats are:
Format | Read | Write | Auto-identify |
---|---|---|---|
framecpp | Yes | No | No |
gwf | Yes | No | Yes |
lalframe | Yes | No | No |
resample
(rate, **kwargs)[source]¶Resample items in this dict.
This operation over-writes items inplace.
Parameters: |
kwargs :
|
---|
setdefault
(k[, d]) → od.get(k,d), also set od[k]=d if k not in od¶update
([E, ]**F) → None. Update D from mapping/iterable E and F.¶If E present and has a .keys() method, does: for k in E: D[k] = E[k] If E present and lacks .keys() method, does: for (k, v) in E: D[k] = v In either case, this is followed by: for k, v in F.items(): D[k] = v
values
() → list of values in od¶viewitems
() → a set-like object providing a view on od's items¶viewkeys
() → a set-like object providing a view on od's keys¶viewvalues
() → an object providing a view on od's values¶write
(data, *args, **kwargs)¶Write this TimeSeriesDict
to a file
Notes
The available built-in formats are:
Format | Read | Write | Auto-identify |
---|---|---|---|
framecpp | Yes | Yes | No |