00001 //================================================================================================= 00002 // Copyright (c) 2011, Johannes Meyer, TU Darmstadt 00003 // All rights reserved. 00004 00005 // Redistribution and use in source and binary forms, with or without 00006 // modification, are permitted provided that the following conditions are met: 00007 // * Redistributions of source code must retain the above copyright 00008 // notice, this list of conditions and the following disclaimer. 00009 // * Redistributions in binary form must reproduce the above copyright 00010 // notice, this list of conditions and the following disclaimer in the 00011 // documentation and/or other materials provided with the distribution. 00012 // * Neither the name of the Flight Systems and Automatic Control group, 00013 // TU Darmstadt, nor the names of its contributors may be used to 00014 // endorse or promote products derived from this software without 00015 // specific prior written permission. 00016 00017 // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND 00018 // ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED 00019 // WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE 00020 // DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER BE LIABLE FOR ANY 00021 // DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES 00022 // (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; 00023 // LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND 00024 // ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT 00025 // (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS 00026 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 00027 //================================================================================================= 00028 00029 #include <hector_pose_estimation/measurements/baro.h> 00030 #include <hector_pose_estimation/pose_estimation.h> 00031 00032 #include <boost/bind.hpp> 00033 00034 namespace hector_pose_estimation { 00035 00036 BaroModel::BaroModel() 00037 { 00038 stddev_ = 1.0; 00039 qnh_ = 1013.25; 00040 parameters().add("qnh", qnh_); 00041 } 00042 00043 BaroModel::~BaroModel() {} 00044 00045 ColumnVector BaroModel::ExpectedValueGet() const 00046 { 00047 y_(1) = qnh_ * pow(1.0 - (0.0065 * (x_(POSITION_Z) + getElevation())) / 288.15, 5.255); 00048 return y_; 00049 } 00050 00051 Matrix BaroModel::dfGet(unsigned int i) const 00052 { 00053 C_(1,POSITION_Z) = qnh_ * 5.255 * pow(1.0 - (0.0065 * (x_(POSITION_Z) + getElevation())) / 288.15, 4.255) * (-0.0065 / 288.15); 00054 return C_; 00055 } 00056 00057 double BaroModel::getAltitude(const BaroUpdate& update) 00058 { 00059 return 288.15 / 0.0065 * (1.0 - pow(update.getVector()(1) / qnh_, 1.0/5.255)); 00060 } 00061 00062 BaroUpdate::BaroUpdate() : qnh_(0) {} 00063 BaroUpdate::BaroUpdate(double pressure) : qnh_(0) { setValue(pressure); } 00064 BaroUpdate::BaroUpdate(double pressure, double qnh) : qnh_(qnh) { setValue(pressure); } 00065 00066 void Baro::onReset() 00067 { 00068 HeightBaroCommon::onReset(); 00069 } 00070 00071 bool Baro::beforeUpdate(PoseEstimation &estimator, const Baro::Update &update) { 00072 if (update.qnh() != 0) setQnh(update.qnh()); 00073 setElevation(resetElevation(estimator, boost::bind(&BaroModel::getAltitude, getModel(), update))); 00074 return true; 00075 } 00076 00077 } // namespace hector_pose_estimation