void updateCallback(const ros::TimerEvent& e)
{
static bool got_first = false;
latest_dt = (e.current_real - e.last_real).toSec();
latest_jitter = (e.current_real - e.current_expected).toSec();
max_abs_jitter = max(max_abs_jitter, fabs(latest_jitter));
Result::Enum the_result;
bool was_new_frame = true;
if ((not segment_data_enabled) //
and (pose_pub.getNumSubscribers() > 0 or markers_pub.getNumSubscribers() > 0))
{
the_result = MyClient.EnableSegmentData().Result;
if (the_result != Result::Success)
{
ROS_ERROR("Enable segment data call failed");
}
else
{
ROS_ASSERT(MyClient.IsSegmentDataEnabled().Enabled);
ROS_INFO("Segment data enabled");
segment_data_enabled = true;
}
}
if (segment_data_enabled)
{
while (was_new_frame and (the_result = MyClient.GetFrame().Result) == Result::Success)
;
{
now_time = ros::Time::now(); // try to grab as close to getting message as possible
was_new_frame = process_frame();
got_first = true;
}
if (the_result != Result::NoFrame)
{
ROS_ERROR_STREAM("GetFrame() returned " << (int)the_result);
}
if (got_first)
{
max_period_between_updates = max(max_period_between_updates, latest_dt);
last_callback_duration = e.profile.last_duration.toSec();
}
}
diag_updater.update();
}
bool grab_vicon_pose_callback(vicon_mocap::viconGrabPose::Request& req, vicon_mocap::viconGrabPose::Response& resp)
{
ROS_INFO("Got request for a VICON pose");
updateTimer.stop();
tf::StampedTransform transform;
//tf::Quaternion quat;
if (not segment_data_enabled)
{
MyClient.EnableSegmentData();
ROS_ASSERT(MyClient.IsSegmentDataEnabled().Enabled);
segment_data_enabled = true;
}
// Gather data:
int N = req.n_measurements;
double accum_trans[3] = {0, 0, 0};
double accum_quat[4] = {0, 0, 0, 0};
Result::Enum the_result;
int n_success = 0;
for (int k = 0; k < N; k++)
{
ros::Duration(1 / vicon_capture_rate).sleep(); // Wait at least as long as vicon needs to capture the next frame
if ((the_result = MyClient.GetFrame().Result) == Result::Success)
{
try
{
Output_GetSegmentGlobalTranslation trans = MyClient.GetSegmentGlobalTranslation(req.subject_name,
req.segment_name);
Output_GetSegmentGlobalRotationQuaternion quat =
MyClient.GetSegmentGlobalRotationQuaternion(req.subject_name, req.segment_name);
if ((!trans.Occluded) && (!quat.Occluded))
{
accum_trans[0] += trans.Translation[0] / 1000;
accum_trans[1] += trans.Translation[1] / 1000;
accum_trans[2] += trans.Translation[2] / 1000;
accum_quat[0] += quat.Rotation[3];
accum_quat[1] += quat.Rotation[0];
accum_quat[2] += quat.Rotation[1];
accum_quat[3] += quat.Rotation[2];
n_success++;
}
}
catch (tf::TransformException ex)
{
ROS_ERROR("%s", ex.what());
resp.success = false;
return false; // TODO: should we really bail here, or just try again?
}
}
else
{
if (the_result != Result::NoFrame)
{
ROS_ERROR_STREAM("GetFrame() returned " << (int)the_result);
}
}
}
ROS_INFO("Averaging %d measurements", n_success);
// Average the data:
double normalized_quat[4];
double quat_norm = sqrt(
accum_quat[0] * accum_quat[0] + accum_quat[1] * accum_quat[1] + accum_quat[2]
* accum_quat[2] + accum_quat[3] * accum_quat[3]);
for (int i = 0; i < 4; i++)
{
normalized_quat[i] = accum_quat[i] / quat_norm;
}
double normalized_vector[3];
// Copy to inputs:
for (int i = 0; i < 3; i++)
{
normalized_vector[i] = accum_trans[i] / n_success;
}
// copy what we used to service call response:
resp.success = true;
resp.pose.header.stamp = ros::Time::now();
resp.pose.header.frame_id = tf_ref_frame_id;
resp.pose.pose.position.x = normalized_vector[0];
resp.pose.pose.position.y = normalized_vector[1];
resp.pose.pose.position.z = normalized_vector[2];
resp.pose.pose.orientation.w = normalized_quat[0];
resp.pose.pose.orientation.x = normalized_quat[1];
resp.pose.pose.orientation.y = normalized_quat[2];
resp.pose.pose.orientation.z = normalized_quat[3];
updateTimer.start();
return true;
}
// mocap position thread
void DataStreamClient::run(void)
{
while (true) //TODO: exit cleanly?
{
// get frame
while (_vicon_client.GetFrame().Result != Result::Success) {
dbg("Couldn't GetFrame()\n");
usleep(1000);
}
// get timecode and populate subject message
// Output_GetTimecode time_code = _vicon_client.GetTimecode();
vicon::body_t msg;
msg.utime = _timestamp_now();
vicon::vicon_t system;
// get subject count
system.nummodels = _vicon_client.GetSubjectCount().SubjectCount;
// Uh oh, models disappeared - realloc the system
if (system.nummodels < system.models.size())
{
system.models.clear();
system.models.reserve(system.nummodels);
}
if (system.nummodels == 0)
{
std::cout << "no models detected" << std::endl;
continue;
}
// loop through subjects
for (uint subject_index = 0; subject_index < system.nummodels; subject_index++) {
////////////////////////////////////////////////////////////////////////
/// Marker
vicon::model_t& model = findModel(system.models, _vicon_client.GetSubjectName(subject_index).SubjectName);
model.nummarkers = _vicon_client.GetMarkerCount(model.name).MarkerCount;
if (model.nummarkers < model.markers.size())
{
model.markers.clear();
model.markers.reserve(model.nummarkers);
}
for(int32_t j = 0; j < model.nummarkers; j++)
{
vicon::marker_t& marker = findMarker(model.markers, _vicon_client.GetMarkerName(model.name, j).MarkerName);
Output_GetMarkerGlobalTranslation translation = _vicon_client.GetMarkerGlobalTranslation(model.name, marker.name);
marker.o = translation.Occluded;
memcpy(marker.xyz, translation.Translation, 3*sizeof(double));
}
////////////////////////////////////////////////////////////////////////
/// Segments
// get number of segments
model.numsegments = _vicon_client.GetSegmentCount(model.name).SegmentCount;
if (model.numsegments < model.segments.size())
{
model.segments.clear();
model.segments.reserve(model.numsegments);
}
// get subject name
std::string subject_name = _vicon_client.GetSubjectName(subject_index).SubjectName;
// get root segment name
std::string root_segment_name = _vicon_client.GetSubjectRootSegmentName(subject_name).SegmentName;
// loop through segments
for (uint segment_index = 0; segment_index < model.numsegments; segment_index++) {
// get segment name
std::string segment_name = _vicon_client.GetSegmentName(subject_name, segment_index).SegmentName;
vicon::segment_t& segment = findSegment(model.segments, segment_name);
Output_GetSegmentGlobalRotationEulerXYZ A = _vicon_client.GetSegmentGlobalRotationEulerXYZ(model.name, segment.name);
Output_GetSegmentGlobalTranslation T = _vicon_client.GetSegmentGlobalTranslation(model.name, segment.name);
Output_GetSegmentLocalRotationEulerXYZ ba = _vicon_client.GetSegmentLocalRotationEulerXYZ(model.name, segment.name);
Output_GetSegmentLocalTranslation bt = _vicon_client.GetSegmentLocalTranslation(model.name, segment.name);
memcpy(segment.A, A.Rotation, 3*sizeof(double));
memcpy(segment.T, T.Translation, 3*sizeof(double));
memcpy(segment.ba, ba.Rotation, 3*sizeof(double));
memcpy(segment.bt, bt.Translation, 3*sizeof(double));
// check if root segment
if (segment_name == root_segment_name) { //TODO: handle articulated bodies
// get segment translation
// Output_GetSegmentStaticTranslation segment_translation = _vicon_client.GetSegmentStaticTranslation(subject_name, segment_name);
Output_GetSegmentGlobalTranslation segment_translation = _vicon_client.GetSegmentGlobalTranslation(
subject_name, segment_name);
// get segment rotation
//Output_GetSegmentStaticRotationQuaternion segment_rotation = _vicon_client.GetSegmentStaticRotationQuaternion(subject_name, segment_name);
Output_GetSegmentGlobalRotationQuaternion segment_rotation =
_vicon_client.GetSegmentGlobalRotationQuaternion(subject_name, segment_name);
// Output_GetSegmentGlobalRotationEulerXYZ segment_rotation = _vicon_client.GetSegmentGlobalRotationEulerXYZ(subject_name, segment_name);
// populate message with position
for (int i = 0; i < 3; i++)
msg.trans[i] = segment_translation.Translation[i]/1000.0; //vicon data is in mm
msg.quat[0] = segment_rotation.Rotation[3]; //vicon is x,y,z,w
msg.quat[1] = segment_rotation.Rotation[0];
msg.quat[2] = segment_rotation.Rotation[1];
msg.quat[3] = segment_rotation.Rotation[2];
std::string channel = "VICON_" + subject_name;
_lcm.publish(channel.c_str(), &msg);
// break from segment for loop
break;
} // root segment
} // for each segment (model.numsegments)
} // for each subject (system.nummodels)
_lcm.publish("VICON_MARKERS", &system);
}
// disconnect client
dbg("Subject position thread terminated properly, disconnecting client\n");
_vicon_client.Disconnect();
}
