Skip to main content

Sparse sensor placement

Project description

Build Documentation Status PyPI Codecov Binder

PySensors is a Scikit-learn style Python package for the sparse placement of sensors, either for reconstruction or classification tasks.

Sparse sensor placement

Sparse sensor placement concerns the problem of selecting a small subset of sensor or measurement locations in a way that allows one to perform some task nearly as well as if one had access to measurements at every location.

PySensors provides objects designed for the tasks of reconstruction and classification.

Reconstruction

Reconstruction deals with predicting the values of a quantity of interest at different locations other than those where sensors are located. For example, one might predict the temperature at a point in the middle of a lake based on temperature readings taken at various other positions in the lake.

PySensors provides the SensorSelector class to aid in the solution of reconstruction problems.

Take representative examples of the types of data to be reconstructed (in this case polynomials)

x = numpy.linspace(0, 1, 1001)
data = numpy.vander(x, 11).T  # Select

and feed them into a SensorSelector

model = pysensors.Sensorselector(n_sensors=10)
model.fit(x)

Use the predict method to reconstruct a new function sampled at the chosen sensor locations:

f = numpy.abs(x[method.selected_sensors]**2 - 0.5)
f_pred = model.predict(f)
A plot showing the function to be reconstructed, the learned sensor locations, and the reconstruction.

Classification

Classification is the problem of predicting which category an example belongs to, given a set of training data (e.g. determining whether digital photos are of dogs or cats). The SSPOC class is used to solve classification problems. Users familiar with Scikit-learn will find it intuitive:

model = pysensors.classification.SSPOC()
model.fit(x, y)  # Learn sensor locations and fit a linear classifier
y_pred = model.predict(x_test[:, model.selected_sensors])  #  Get predictions

See our set of classification examples for more information.

Bases

The basis in which measurement data are represented can have a dramatic effect on performance. PySensors implements the three bases most commonly used for sparse sensor placement: raw measurements, SVD/POD/PCA modes, and random projections. Bases can be easily incorporated into SensorSelector and SSPOC classes:

basis = pysensors.basis.SVD(n_basis_modes=20)
recon_model = pysensors.SensorSelector(basis=basis)
class_model = pysensors.classification.SSPOC(basis=basis)

See this example for further discussion of these options.

Installation

Installing with pip

If you are using Linux or macOS you can install PySensors with pip:

pip install python-sensors

Note: the name you type in here is python-sensors and is not pysensors.

Installing from source

First clone this repository:

git clone https://github.com/dynamicslab/pysensors.git

Then, to install the package, run

cd pysensors
pip install .

If you do not have pip you can instead use

python setup.py install

If you do not have root access, you should add the --user option to the install commands above.

Features

The primary PySensors object is the SensorSelector class, which is used to choose sensor locations optimized for reconstruction tasks. Other implemented objects include

  • SSPOC - a class for sparse sensor selection for classification problems

  • basis - submodule implementing different bases in which to represent data - Identity - use raw measurement data - SVD - efficiently compute first k left singular vectors - RandomProjection - Gaussian random projections of measurements

  • Convenience functions to aid in the analysis of error as number of sensors or basis modes are varied

Documentation

PySensors has a documentation site hosted by readthedocs. Examples are available online, as static Jupyter notebooks and as interactive notebooks.

Community guidelines

Contributing examples

If you have used PySensors to solve an interesting problem, please consider submitting an example Jupyter notebook showcasing your work!

Contributing code

We welcome contributions to PySensors. To contribute a new feature please submit a pull request. To get started we recommend installing the packages in requirements-dev.txt via

pip install -r requirements-dev.txt

This will allow you to run unit tests and automatically format your code. To be accepted your code should conform to PEP8 and pass all unit tests. Code can be tested by invoking

pytest

We recommend using pre-commit to format your code. Once you have staged changes to commit

git add path/to/changed/file.py

you can run the following to automatically reformat your staged code

pre-commit

Note that you will then need to re-stage any changes pre-commit made to your code.

Reporting issues or bugs

If you find a bug in the code or want to request a new feature, please open an issue.

References

  • Manohar, Krithika, Bingni W. Brunton, J. Nathan Kutz, and Steven L. Brunton. “Data-driven sparse sensor placement for reconstruction: Demonstrating the benefits of exploiting known patterns.” IEEE Control Systems Magazine 38, no. 3 (2018): 63-86. [DOI]

  • Clark, Emily, Travis Askham, Steven L. Brunton, and J. Nathan Kutz. “Greedy sensor placement with cost constraints.” IEEE Sensors Journal 19, no. 7 (2018): 2642-2656. [DOI]

  • Brunton, Bingni W., Steven L. Brunton, Joshua L. Proctor, and J Nathan Kutz.

    “Sparse sensor placement optimization for classification.” SIAM Journal on Applied Mathematics 76.5 (2016): 2099-2122. [DOI]

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

python-sensors-0.3.0.tar.gz (1.3 MB view hashes)

Uploaded Source

Built Distribution

python_sensors-0.3.0-py3-none-any.whl (25.0 kB view hashes)

Uploaded Python 3

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page