humotion

/*

* This file is part of humotion

*

* Copyright(c) sschulz <AT> techfak.uni-bielefeld.de

* http://opensource.cit-ec.de/projects/humotion

*

* This file may be licensed under the terms of the

* GNU Lesser General Public License Version 3 (the ``LGPL''),

* or (at your option) any later version.

*

* Software distributed under the License is distributed

* on an ``AS IS'' basis, WITHOUT WARRANTY OF ANY KIND, either

* express or implied. See the LGPL for the specific language

* governing rights and limitations.

*

* You should have received a copy of the LGPL along with this

* program. If not, go to http://www.gnu.org/licenses/lgpl.html

* or write to the Free Software Foundation, Inc.,

* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.

*

* The development of this software was supported by the

* Excellence Cluster EXC 277 Cognitive Interaction Technology.

* The Excellence Cluster EXC 277 is a grant of the Deutsche

* Forschungsgemeinschaft (DFG) in the context of the German

* Excellence Initiative.

*/

#include "server/eye_motion_generator.h"

#include "server/server.h"

using namespace std;

using namespace humotion;

using namespace humotion::server;

//saccade detection threshold in deg/s

const float EyeMotionGenerator::SACCADE_SPEED_THRESHOLD = 15.0;

//! constructor

EyeMotionGenerator::EyeMotionGenerator(JointInterface *j) : GazeMotionGenerator(j, 3, 1.0/Server::MOTION_UPDATERATE){

}

//! destructor

EyeMotionGenerator::~EyeMotionGenerator(){

}

//! set up an eyemotion profile

//! this will use speed and acceleration calc formulas from literature:

//! \param dof id of joint

//! \param target angle

//! \param current angle

void EyeMotionGenerator::setup_eyemotion(int dof, float target, float now){

    //get distance to target:

    float distance_abs = fabs(target - now);

    //get max speed: factor can be found in encyc britannica: "linear.... being 300° per second for 10° and 500° per second for 30°" (max=700)

    float max_speed = fmin(700.0, 10.0*distance_abs + 200.0);

    //max accel: use data from ""Speed and Accuracy of Saccadic Eye Movements: Characteristics of Impulse Variability in the Oculomotor System"

    //                 http://www-personal.umich.edu/~kornblum/files/journal_exp_psych_HPP_15-3.pdf

    //           gathered from table2

    float max_accel = fmin(80000.0, 1526.53*distance_abs + 10245.4);

    //feed reflexxes api with data

    reflexxes_set_input(dof, target, now, max_speed, max_accel);

}

//! calculate joint targets

void EyeMotionGenerator::calculate_targets(){

    //use the target values for a faster response

    float neck_pan_now  = joint_interface->get_target_position(JointInterface::ID_NECK_PAN);

    float neck_tilt_now = joint_interface->get_target_position(JointInterface::ID_NECK_TILT);

    //calculate target angles for the eyes:

    //right eye is dominant -> direct output

    float eye_pan_r_target = (requested_gaze_state.pan + requested_gaze_state.vergence/2.0) - (neck_pan_now);

    //left eye is non dominant -> filtered output: TODO: activate low pass filtered output

    //FIXME: USE LOWPASS FILTER HERE --> output_angle[angle_names::EYE_PAN_L] + 0.1 * (eye_pan_l_target - output_angle[angle_names::EYE_PAN_L]) etc;

    float eye_pan_l_target = (requested_gaze_state.pan - requested_gaze_state.vergence/2.0) - (neck_pan_now);

    //tilt:

    float eye_tilt_target = requested_gaze_state.tilt - (neck_tilt_now);

    if (0){

        eye_pan_r_target = 0.0;

        eye_pan_l_target = 0.0;

        eye_tilt_target  = 0.0;

    }

    //check and take care of limits

    eye_pan_l_target = limit_target(JointInterface::ID_EYES_LEFT_LR, eye_pan_l_target);

    eye_pan_r_target = limit_target(JointInterface::ID_EYES_RIGHT_LR, eye_pan_r_target);

    eye_tilt_target = limit_target(JointInterface::ID_EYES_BOTH_UD, eye_tilt_target);

    //fetch current angles:

    float eye_pan_l_now = get_current_position(JointInterface::ID_EYES_LEFT_LR);

    float eye_pan_r_now = get_current_position(JointInterface::ID_EYES_RIGHT_LR);

    float eye_tilt_now  = get_current_position(JointInterface::ID_EYES_BOTH_UD);

    //pass paramaters to reflexxes api:

    setup_eyemotion(0, eye_pan_l_target, eye_pan_l_now);

    setup_eyemotion(1, eye_pan_r_target, eye_pan_r_now);

    setup_eyemotion(2, eye_tilt_target, eye_tilt_now);

    cout << "EYE MOTION 2 " << eye_tilt_target << " now=" << eye_tilt_now << "\n";

    //call reflexxes to handle profile calculation:

    reflexxes_calculate_profile();

    //tell the joint about the new values:

    joint_interface->set_target_position(JointInterface::ID_EYES_LEFT_LR,

                                         reflexxes_position_output->NewPositionVector->VecData[0],

                                         reflexxes_position_output->NewVelocityVector->VecData[0]);

    joint_interface->set_target_position(JointInterface::ID_EYES_RIGHT_LR,

                                         reflexxes_position_output->NewPositionVector->VecData[1],

                                         reflexxes_position_output->NewVelocityVector->VecData[1]);

    joint_interface->set_target_position(JointInterface::ID_EYES_BOTH_UD,

                                         reflexxes_position_output->NewPositionVector->VecData[2],

                                         reflexxes_position_output->NewVelocityVector->VecData[2]);

    cout << reflexxes_position_output->NewPositionVector->VecData[2] << " " << reflexxes_position_output->NewVelocityVector->VecData[2] ;

}

//! pubish the calculated targets to the joint subsystem

void EyeMotionGenerator::publish_targets(){

    //publish values if there is an active gaze input within the last timerange

    if (gaze_target_input_active()){

        joint_interface->publish_target_position(JointInterface::ID_EYES_LEFT_LR);

        joint_interface->publish_target_position(JointInterface::ID_EYES_RIGHT_LR);

        joint_interface->publish_target_position(JointInterface::ID_EYES_BOTH_UD);

    }

}