3atv精品不卡视频,97人人超碰国产精品最新,中文字幕av一区二区三区人妻少妇,久久久精品波多野结衣,日韩一区二区三区精品

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁(yè) > 编程资源 > 编程问答 >内容正文

编程问答

ROS学习笔记之——amcl源码的解读

發(fā)布時(shí)間:2023/12/20 编程问答 29 豆豆
生活随笔 收集整理的這篇文章主要介紹了 ROS学习笔记之——amcl源码的解读 小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

之前博文《ROS學(xué)習(xí)筆記之——gmapping與amcl 》已經(jīng)介紹過gmapping與amcl了。本博文詳細(xì)的看一下amcl的代碼。如需要修改amcl,用新的包,可以選擇到/opt/ros/melodic/include/amcl目錄下,把舊的amcl刪掉,然后再編譯新的,當(dāng)然也可以通過更改名字(更改包的名字)。個(gè)人認(rèn)為,更改名字是比較好的方法了哈~

?

目錄

AMCL源碼

TF

TF的數(shù)據(jù)類型

發(fā)布TF

訂閱TF

關(guān)于TF的一些測(cè)試

參考資料


?

從之前博文以及實(shí)驗(yàn)中發(fā)現(xiàn),amcl定位有兩個(gè)關(guān)鍵點(diǎn)。

1,需要一個(gè)初始的姿態(tài)估計(jì)。初始的位姿估計(jì)并不需要很準(zhǔn)確,大致的位置和orientation對(duì)了就可以了。

2,在移動(dòng)的過程中,如果把a(bǔ)mcl的點(diǎn)可視化出來,會(huì)發(fā)現(xiàn)隨著機(jī)器人的移動(dòng),會(huì)逐漸收斂。(在之前博客中已經(jīng)做過類似的仿真了,后來實(shí)驗(yàn)也可以得到類似的效果,但收斂的速度和最終收斂的大小,相對(duì)沒有仿真的這么好)

amcl使用粒子濾波器在已知地圖的情況下定位機(jī)器人位姿。它根據(jù)粒子集合的統(tǒng)計(jì)數(shù)據(jù)插入一定數(shù)量的新粒子,來解決全局定位失效或者說是機(jī)器人綁架的問題。 并根據(jù)KLD采樣自適應(yīng)的調(diào)節(jié)粒子數(shù)量,可以在定位精度和計(jì)算代價(jià)之間進(jìn)行權(quán)衡。

?

AMCL源碼

打開AMCL的源碼,會(huì)發(fā)現(xiàn)有大量的文檔,但是關(guān)鍵的節(jié)點(diǎn)應(yīng)該就是amcl_node.cpp文件

?

關(guān)于源碼的解讀見下面注釋

https://gaoyichao.com/Xiaotu/?book=turtlebot&title=index#part6

https://gaoyichao.com/Xiaotu/?book=turtlebot&title=amcl%E7%9A%84%E6%80%BB%E4%BD%93%E4%B8%9A%E5%8A%A1%E9%80%BB%E8%BE%91

https://gaoyichao.com/Xiaotu/?book=turtlebot&title=amcl_demo%E4%B9%8B%E5%AE%9A%E4%BD%8D%E5%99%A8

通過turtlbot navigation中的amcl.launch文件定義了運(yùn)行amcl的配置參數(shù),

?

下面來看看amcl_node.cpp文件

/** Copyright (c) 2008, Willow Garage, Inc.* All rights reserved.** This library is free software; you can redistribute it and/or* modify it under the terms of the GNU Lesser General Public* License as published by the Free Software Foundation; either* version 2.1 of the License, or (at your option) any later version.** This library is distributed in the hope that it will be useful,* but WITHOUT ANY WARRANTY; without even the implied warranty of* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU* Lesser General Public License for more details.** You should have received a copy of the GNU Lesser General Public* License along with this library; if not, write to the Free Software* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA**//* Author: Brian Gerkey */#include <algorithm> #include <vector> #include <map> #include <cmath> #include <memory>#include <boost/bind.hpp> #include <boost/thread/mutex.hpp>// Signal handling #include <signal.h>#include "amcl/map/map.h" #include "amcl/pf/pf.h" #include "amcl/sensors/amcl_odom.h" #include "amcl/sensors/amcl_laser.h" #include "portable_utils.hpp"#include "ros/assert.h"// roscpp #include "ros/ros.h"// Messages that I need #include "sensor_msgs/LaserScan.h" #include "geometry_msgs/PoseWithCovarianceStamped.h" #include "geometry_msgs/PoseArray.h" #include "geometry_msgs/Pose.h" #include "geometry_msgs/PoseStamped.h" #include "nav_msgs/GetMap.h" #include "nav_msgs/SetMap.h" #include "std_srvs/Empty.h"// For transform support #include "tf2/LinearMath/Transform.h" #include "tf2/convert.h" #include "tf2/utils.h" #include "tf2_geometry_msgs/tf2_geometry_msgs.h" #include "tf2_ros/buffer.h" #include "tf2_ros/message_filter.h" #include "tf2_ros/transform_broadcaster.h" #include "tf2_ros/transform_listener.h" #include "message_filters/subscriber.h"// Dynamic_reconfigure #include "dynamic_reconfigure/server.h" #include "amcl/AMCLConfig.h"// Allows AMCL to run from bag file #include <rosbag/bag.h> #include <rosbag/view.h> #include <boost/foreach.hpp>// For monitoring the estimator #include <diagnostic_updater/diagnostic_updater.h>#define NEW_UNIFORM_SAMPLING 1using namespace amcl;// Pose hypothesis 姿態(tài)的假設(shè) typedef struct {// Total weight (weights sum to 1)double weight;//姿態(tài)的權(quán)重。每個(gè)姿態(tài)有對(duì)應(yīng)的權(quán)重來衡量其確信度// Mean of pose esimatepf_vector_t pf_pose_mean;//粒子估計(jì)姿態(tài)的均值// Covariance of pose estimatepf_matrix_t pf_pose_cov;//估計(jì)的協(xié)方差} amcl_hyp_t;static double normalize(double z) {return atan2(sin(z),cos(z)); } static double angle_diff(double a, double b) {double d1, d2;a = normalize(a);b = normalize(b);d1 = a-b;d2 = 2*M_PI - fabs(d1);if(d1 > 0)d2 *= -1.0;if(fabs(d1) < fabs(d2))return(d1);elsereturn(d2); }static const std::string scan_topic_ = "scan";//激光雷達(dá)的數(shù)據(jù)/* This function is only useful to have the whole code work* with old rosbags that have trailing slashes for their frames*/ inline std::string stripSlash(const std::string& in) {std::string out = in;if ( ( !in.empty() ) && (in[0] == '/') )out.erase(0,1);return out; }class AmclNode {public:AmclNode();~AmclNode();/*** @brief Uses TF and LaserScan messages from bag file to drive AMCL instead* @param in_bag_fn input bagfile* @param trigger_global_localization whether to trigger (觸發(fā)) global localization* before starting to process the bagfile*/void runFromBag(const std::string &in_bag_fn, bool trigger_global_localization = false);int process();void savePoseToServer();private:std::shared_ptr<tf2_ros::TransformBroadcaster> tfb_;std::shared_ptr<tf2_ros::TransformListener> tfl_;std::shared_ptr<tf2_ros::Buffer> tf_;bool sent_first_transform_;tf2::Transform latest_tf_;bool latest_tf_valid_;// Pose-generating function used to uniformly distribute particles over// the mapstatic pf_vector_t uniformPoseGenerator(void* arg); #if NEW_UNIFORM_SAMPLINGstatic std::vector<std::pair<int,int> > free_space_indices; #endif// Callbacksbool globalLocalizationCallback(std_srvs::Empty::Request& req,std_srvs::Empty::Response& res);bool nomotionUpdateCallback(std_srvs::Empty::Request& req,std_srvs::Empty::Response& res);bool setMapCallback(nav_msgs::SetMap::Request& req,nav_msgs::SetMap::Response& res);void laserReceived(const sensor_msgs::LaserScanConstPtr& laser_scan);void initialPoseReceived(const geometry_msgs::PoseWithCovarianceStampedConstPtr& msg);void handleInitialPoseMessage(const geometry_msgs::PoseWithCovarianceStamped& msg);void mapReceived(const nav_msgs::OccupancyGridConstPtr& msg);void handleMapMessage(const nav_msgs::OccupancyGrid& msg);void freeMapDependentMemory();map_t* convertMap( const nav_msgs::OccupancyGrid& map_msg );void updatePoseFromServer();void applyInitialPose();//parameter for what odom to usestd::string odom_frame_id_;//paramater to store latest odom posegeometry_msgs::PoseStamped latest_odom_pose_;//parameter for what base to usestd::string base_frame_id_;std::string global_frame_id_;bool use_map_topic_;bool first_map_only_;ros::Duration gui_publish_period;ros::Time save_pose_last_time;ros::Duration save_pose_period;geometry_msgs::PoseWithCovarianceStamped last_published_pose;map_t* map_;char* mapdata;int sx, sy;double resolution;message_filters::Subscriber<sensor_msgs::LaserScan>* laser_scan_sub_;tf2_ros::MessageFilter<sensor_msgs::LaserScan>* laser_scan_filter_;ros::Subscriber initial_pose_sub_;std::vector< AMCLLaser* > lasers_;std::vector< bool > lasers_update_;std::map< std::string, int > frame_to_laser_;// Particle filterpf_t *pf_;double pf_err_, pf_z_;bool pf_init_;pf_vector_t pf_odom_pose_;double d_thresh_, a_thresh_;int resample_interval_;int resample_count_;double laser_min_range_;double laser_max_range_;//Nomotion update controlbool m_force_update; // used to temporarily let amcl update samples even when no motion occurs...AMCLOdom* odom_;AMCLLaser* laser_;ros::Duration cloud_pub_interval;ros::Time last_cloud_pub_time;// For slowing play-back when reading directly from a bag fileros::WallDuration bag_scan_period_;void requestMap();// Helper to get odometric pose from transform systembool getOdomPose(geometry_msgs::PoseStamped& pose,double& x, double& y, double& yaw,const ros::Time& t, const std::string& f);//time for tolerance on the published transform,//basically defines how long a map->odom transform is good forros::Duration transform_tolerance_;ros::NodeHandle nh_;ros::NodeHandle private_nh_;ros::Publisher pose_pub_;ros::Publisher particlecloud_pub_;ros::ServiceServer global_loc_srv_;ros::ServiceServer nomotion_update_srv_; //to let amcl update samples without requiring motionros::ServiceServer set_map_srv_;ros::Subscriber initial_pose_sub_old_;ros::Subscriber map_sub_;diagnostic_updater::Updater diagnosic_updater_;void standardDeviationDiagnostics(diagnostic_updater::DiagnosticStatusWrapper& diagnostic_status);double std_warn_level_x_;double std_warn_level_y_;double std_warn_level_yaw_;amcl_hyp_t* initial_pose_hyp_;bool first_map_received_;bool first_reconfigure_call_;boost::recursive_mutex configuration_mutex_;dynamic_reconfigure::Server<amcl::AMCLConfig> *dsrv_;amcl::AMCLConfig default_config_;ros::Timer check_laser_timer_;int max_beams_, min_particles_, max_particles_;double alpha1_, alpha2_, alpha3_, alpha4_, alpha5_;double alpha_slow_, alpha_fast_;double z_hit_, z_short_, z_max_, z_rand_, sigma_hit_, lambda_short_;//beam skip related paramsbool do_beamskip_;double beam_skip_distance_, beam_skip_threshold_, beam_skip_error_threshold_;double laser_likelihood_max_dist_;odom_model_t odom_model_type_;double init_pose_[3];double init_cov_[3];laser_model_t laser_model_type_;bool tf_broadcast_;bool selective_resampling_;void reconfigureCB(amcl::AMCLConfig &config, uint32_t level);ros::Time last_laser_received_ts_;ros::Duration laser_check_interval_;void checkLaserReceived(const ros::TimerEvent& event); };#if NEW_UNIFORM_SAMPLING std::vector<std::pair<int,int> > AmclNode::free_space_indices; #endif#define USAGE "USAGE: amcl"boost::shared_ptr<AmclNode> amcl_node_ptr;//定義一個(gè)amcl節(jié)點(diǎn)的指針void sigintHandler(int sig) {// Save latest pose as we're shutting down.amcl_node_ptr->savePoseToServer();ros::shutdown(); }//主函數(shù) int main(int argc, char** argv) {ros::init(argc, argv, "amcl");ros::NodeHandle nh;// Override default sigint handlersignal(SIGINT, sigintHandler);// Make our node available to sigintHandleramcl_node_ptr.reset(new AmclNode());//通過指針來創(chuàng)建amcl類(AmclNode)if (argc == 1){// run using ROS inputros::spin();}else if ((argc >= 3) && (std::string(argv[1]) == "--run-from-bag")){if (argc == 3){amcl_node_ptr->runFromBag(argv[2]);}else if ((argc == 4) && (std::string(argv[3]) == "--global-localization")){amcl_node_ptr->runFromBag(argv[2], true);}}// Without this, our boost locks are not shut down nicelyamcl_node_ptr.reset();// To quote Morgan, Hooray!return(0); }AmclNode::AmclNode() :sent_first_transform_(false),latest_tf_valid_(false),map_(NULL),pf_(NULL),resample_count_(0),odom_(NULL),laser_(NULL),private_nh_("~"),initial_pose_hyp_(NULL),first_map_received_(false),first_reconfigure_call_(true) {boost::recursive_mutex::scoped_lock l(configuration_mutex_);// 讀入一系列的參數(shù)// Grab params off the param serverprivate_nh_.param("use_map_topic", use_map_topic_, false);private_nh_.param("first_map_only", first_map_only_, false);double tmp;private_nh_.param("gui_publish_rate", tmp, -1.0);gui_publish_period = ros::Duration(1.0/tmp);private_nh_.param("save_pose_rate", tmp, 0.5);save_pose_period = ros::Duration(1.0/tmp);private_nh_.param("laser_min_range", laser_min_range_, -1.0);private_nh_.param("laser_max_range", laser_max_range_, -1.0);private_nh_.param("laser_max_beams", max_beams_, 30);private_nh_.param("min_particles", min_particles_, 100);private_nh_.param("max_particles", max_particles_, 5000);private_nh_.param("kld_err", pf_err_, 0.01);private_nh_.param("kld_z", pf_z_, 0.99);private_nh_.param("odom_alpha1", alpha1_, 0.2);private_nh_.param("odom_alpha2", alpha2_, 0.2);private_nh_.param("odom_alpha3", alpha3_, 0.2);private_nh_.param("odom_alpha4", alpha4_, 0.2);private_nh_.param("odom_alpha5", alpha5_, 0.2);private_nh_.param("do_beamskip", do_beamskip_, false);private_nh_.param("beam_skip_distance", beam_skip_distance_, 0.5);private_nh_.param("beam_skip_threshold", beam_skip_threshold_, 0.3);if (private_nh_.hasParam("beam_skip_error_threshold_")){private_nh_.param("beam_skip_error_threshold_", beam_skip_error_threshold_);}else{private_nh_.param("beam_skip_error_threshold", beam_skip_error_threshold_, 0.9);}private_nh_.param("laser_z_hit", z_hit_, 0.95);private_nh_.param("laser_z_short", z_short_, 0.1);private_nh_.param("laser_z_max", z_max_, 0.05);private_nh_.param("laser_z_rand", z_rand_, 0.05);private_nh_.param("laser_sigma_hit", sigma_hit_, 0.2);private_nh_.param("laser_lambda_short", lambda_short_, 0.1);private_nh_.param("laser_likelihood_max_dist", laser_likelihood_max_dist_, 2.0);std::string tmp_model_type;private_nh_.param("laser_model_type", tmp_model_type, std::string("likelihood_field"));if(tmp_model_type == "beam")laser_model_type_ = LASER_MODEL_BEAM;else if(tmp_model_type == "likelihood_field")laser_model_type_ = LASER_MODEL_LIKELIHOOD_FIELD;else if(tmp_model_type == "likelihood_field_prob"){laser_model_type_ = LASER_MODEL_LIKELIHOOD_FIELD_PROB;}else{ROS_WARN("Unknown laser model type \"%s\"; defaulting to likelihood_field model",tmp_model_type.c_str());laser_model_type_ = LASER_MODEL_LIKELIHOOD_FIELD;}private_nh_.param("odom_model_type", tmp_model_type, std::string("diff"));if(tmp_model_type == "diff")odom_model_type_ = ODOM_MODEL_DIFF;else if(tmp_model_type == "omni")odom_model_type_ = ODOM_MODEL_OMNI;else if(tmp_model_type == "diff-corrected")odom_model_type_ = ODOM_MODEL_DIFF_CORRECTED;else if(tmp_model_type == "omni-corrected")odom_model_type_ = ODOM_MODEL_OMNI_CORRECTED;else{ROS_WARN("Unknown odom model type \"%s\"; defaulting to diff model",tmp_model_type.c_str());odom_model_type_ = ODOM_MODEL_DIFF;}private_nh_.param("update_min_d", d_thresh_, 0.2);private_nh_.param("update_min_a", a_thresh_, M_PI/6.0);private_nh_.param("odom_frame_id", odom_frame_id_, std::string("odom"));private_nh_.param("base_frame_id", base_frame_id_, std::string("base_link"));private_nh_.param("global_frame_id", global_frame_id_, std::string("map"));private_nh_.param("resample_interval", resample_interval_, 2);private_nh_.param("selective_resampling", selective_resampling_, false);double tmp_tol;private_nh_.param("transform_tolerance", tmp_tol, 0.1);private_nh_.param("recovery_alpha_slow", alpha_slow_, 0.001);private_nh_.param("recovery_alpha_fast", alpha_fast_, 0.1);private_nh_.param("tf_broadcast", tf_broadcast_, true);// For diagnostics(診斷)private_nh_.param("std_warn_level_x", std_warn_level_x_, 0.2);private_nh_.param("std_warn_level_y", std_warn_level_y_, 0.2);private_nh_.param("std_warn_level_yaw", std_warn_level_yaw_, 0.1);transform_tolerance_.fromSec(tmp_tol);{double bag_scan_period;private_nh_.param("bag_scan_period", bag_scan_period, -1.0);bag_scan_period_.fromSec(bag_scan_period);}odom_frame_id_ = stripSlash(odom_frame_id_);base_frame_id_ = stripSlash(base_frame_id_);global_frame_id_ = stripSlash(global_frame_id_);updatePoseFromServer();cloud_pub_interval.fromSec(1.0);tfb_.reset(new tf2_ros::TransformBroadcaster());tf_.reset(new tf2_ros::Buffer());tfl_.reset(new tf2_ros::TransformListener(*tf_));pose_pub_ = nh_.advertise<geometry_msgs::PoseWithCovarianceStamped>("amcl_pose", 2, true);//發(fā)布的位姿信息particlecloud_pub_ = nh_.advertise<geometry_msgs::PoseArray>("particlecloud", 2, true);//發(fā)布的粒子濾波的信息//注冊(cè)三個(gè)服務(wù)global_loc_srv_ = nh_.advertiseService("global_localization", &AmclNode::globalLocalizationCallback,this);//用于獲取機(jī)器人的全局定位nomotion_update_srv_= nh_.advertiseService("request_nomotion_update", &AmclNode::nomotionUpdateCallback, this);//用于手動(dòng)的觸發(fā)粒子更新并發(fā)布新的位姿估計(jì)set_map_srv_= nh_.advertiseService("set_map", &AmclNode::setMapCallback, this);//設(shè)定機(jī)器人位姿和地圖信息,調(diào)用setMapCallback//構(gòu)建激光傳感器的消息過濾器對(duì)象和tf2的過濾器,并注冊(cè)回調(diào)函數(shù)laserReceived // 這里的message_filter為ROS系統(tǒng)提供了一些通用的消息過濾方法 //它對(duì)接收到的消息進(jìn)行緩存,只有當(dāng)滿足過濾條件后才輸出,在需要消息同步的時(shí)候應(yīng)用比較多。 //這里主要是同時(shí)監(jiān)聽激光掃描消息和里程計(jì)坐標(biāo)變換,同步兩者的輸出。 laser_scan_sub_ = new message_filters::Subscriber<sensor_msgs::LaserScan>(nh_, scan_topic_, 100);//訂閱激光雷達(dá)的信息laser_scan_filter_ = new tf2_ros::MessageFilter<sensor_msgs::LaserScan>(*laser_scan_sub_,*tf_,odom_frame_id_,100,nh_);//調(diào)用laserReceived函數(shù) ,在里面實(shí)現(xiàn)了將激光雷達(dá)信息進(jìn)行tf變換 laser_scan_filter_->registerCallback(boost::bind(&AmclNode::laserReceived,this, _1));訂閱初始位姿,并調(diào)用initialPoseReceived函數(shù) initial_pose_sub_ = nh_.subscribe("initialpose", 2, &AmclNode::initialPoseReceived, this);//進(jìn)行位姿的初始化//通過topic"initialpose"來對(duì)機(jī)器人的位姿做初始化估計(jì)if(use_map_topic_) {map_sub_ = nh_.subscribe("map", 1, &AmclNode::mapReceived, this);ROS_INFO("Subscribed to map topic.");} else {requestMap();}m_force_update = false;dsrv_ = new dynamic_reconfigure::Server<amcl::AMCLConfig>(ros::NodeHandle("~"));dynamic_reconfigure::Server<amcl::AMCLConfig>::CallbackType cb = boost::bind(&AmclNode::reconfigureCB, this, _1, _2);dsrv_->setCallback(cb);// 15s timer to warn on lack of receipt of laser scans, #5209laser_check_interval_ = ros::Duration(15.0);check_laser_timer_ = nh_.createTimer(laser_check_interval_, boost::bind(&AmclNode::checkLaserReceived, this, _1));diagnosic_updater_.setHardwareID("None");diagnosic_updater_.add("Standard deviation", this, &AmclNode::standardDeviationDiagnostics); }void AmclNode::reconfigureCB(AMCLConfig &config, uint32_t level) {boost::recursive_mutex::scoped_lock cfl(configuration_mutex_);//we don't want to do anything on the first call//which corresponds to startupif(first_reconfigure_call_){first_reconfigure_call_ = false;default_config_ = config;return;}if(config.restore_defaults) {config = default_config_;//avoid loopingconfig.restore_defaults = false;}d_thresh_ = config.update_min_d;a_thresh_ = config.update_min_a;resample_interval_ = config.resample_interval;laser_min_range_ = config.laser_min_range;laser_max_range_ = config.laser_max_range;gui_publish_period = ros::Duration(1.0/config.gui_publish_rate);save_pose_period = ros::Duration(1.0/config.save_pose_rate);transform_tolerance_.fromSec(config.transform_tolerance);max_beams_ = config.laser_max_beams;alpha1_ = config.odom_alpha1;alpha2_ = config.odom_alpha2;alpha3_ = config.odom_alpha3;alpha4_ = config.odom_alpha4;alpha5_ = config.odom_alpha5;z_hit_ = config.laser_z_hit;z_short_ = config.laser_z_short;z_max_ = config.laser_z_max;z_rand_ = config.laser_z_rand;sigma_hit_ = config.laser_sigma_hit;lambda_short_ = config.laser_lambda_short;laser_likelihood_max_dist_ = config.laser_likelihood_max_dist;if(config.laser_model_type == "beam")laser_model_type_ = LASER_MODEL_BEAM;else if(config.laser_model_type == "likelihood_field")laser_model_type_ = LASER_MODEL_LIKELIHOOD_FIELD;else if(config.laser_model_type == "likelihood_field_prob")laser_model_type_ = LASER_MODEL_LIKELIHOOD_FIELD_PROB;if(config.odom_model_type == "diff")odom_model_type_ = ODOM_MODEL_DIFF;else if(config.odom_model_type == "omni")odom_model_type_ = ODOM_MODEL_OMNI;else if(config.odom_model_type == "diff-corrected")odom_model_type_ = ODOM_MODEL_DIFF_CORRECTED;else if(config.odom_model_type == "omni-corrected")odom_model_type_ = ODOM_MODEL_OMNI_CORRECTED;if(config.min_particles > config.max_particles){ROS_WARN("You've set min_particles to be greater than max particles, this isn't allowed so they'll be set to be equal.");config.max_particles = config.min_particles;}min_particles_ = config.min_particles;max_particles_ = config.max_particles;alpha_slow_ = config.recovery_alpha_slow;alpha_fast_ = config.recovery_alpha_fast;tf_broadcast_ = config.tf_broadcast;do_beamskip_= config.do_beamskip; beam_skip_distance_ = config.beam_skip_distance; beam_skip_threshold_ = config.beam_skip_threshold; // Clear queued laser objects so that their parameters get updatedlasers_.clear();lasers_update_.clear();frame_to_laser_.clear();if( pf_ != NULL ){pf_free( pf_ );pf_ = NULL;} pf_ = pf_alloc(min_particles_, max_particles_,alpha_slow_, alpha_fast_,(pf_init_model_fn_t)AmclNode::uniformPoseGenerator,(void *)map_);pf_set_selective_resampling(pf_, selective_resampling_);pf_err_ = config.kld_err; pf_z_ = config.kld_z; pf_->pop_err = pf_err_;pf_->pop_z = pf_z_;// Initialize the filterpf_vector_t pf_init_pose_mean = pf_vector_zero();pf_init_pose_mean.v[0] = last_published_pose.pose.pose.position.x;pf_init_pose_mean.v[1] = last_published_pose.pose.pose.position.y;pf_init_pose_mean.v[2] = tf2::getYaw(last_published_pose.pose.pose.orientation);pf_matrix_t pf_init_pose_cov = pf_matrix_zero();pf_init_pose_cov.m[0][0] = last_published_pose.pose.covariance[6*0+0];pf_init_pose_cov.m[1][1] = last_published_pose.pose.covariance[6*1+1];pf_init_pose_cov.m[2][2] = last_published_pose.pose.covariance[6*5+5];pf_init(pf_, pf_init_pose_mean, pf_init_pose_cov);pf_init_ = false;// Instantiate the sensor objects// Odometrydelete odom_;odom_ = new AMCLOdom();ROS_ASSERT(odom_);odom_->SetModel( odom_model_type_, alpha1_, alpha2_, alpha3_, alpha4_, alpha5_ );// Laserdelete laser_;laser_ = new AMCLLaser(max_beams_, map_);ROS_ASSERT(laser_);if(laser_model_type_ == LASER_MODEL_BEAM)laser_->SetModelBeam(z_hit_, z_short_, z_max_, z_rand_,sigma_hit_, lambda_short_, 0.0);else if(laser_model_type_ == LASER_MODEL_LIKELIHOOD_FIELD_PROB){ROS_INFO("Initializing likelihood field model; this can take some time on large maps...");laser_->SetModelLikelihoodFieldProb(z_hit_, z_rand_, sigma_hit_,laser_likelihood_max_dist_, do_beamskip_, beam_skip_distance_, beam_skip_threshold_, beam_skip_error_threshold_);ROS_INFO("Done initializing likelihood field model with probabilities.");}else if(laser_model_type_ == LASER_MODEL_LIKELIHOOD_FIELD){ROS_INFO("Initializing likelihood field model; this can take some time on large maps...");laser_->SetModelLikelihoodField(z_hit_, z_rand_, sigma_hit_,laser_likelihood_max_dist_);ROS_INFO("Done initializing likelihood field model.");}odom_frame_id_ = stripSlash(config.odom_frame_id);base_frame_id_ = stripSlash(config.base_frame_id);global_frame_id_ = stripSlash(config.global_frame_id);delete laser_scan_filter_;laser_scan_filter_ = new tf2_ros::MessageFilter<sensor_msgs::LaserScan>(*laser_scan_sub_,*tf_,odom_frame_id_,100,nh_);laser_scan_filter_->registerCallback(boost::bind(&AmclNode::laserReceived,this, _1));initial_pose_sub_ = nh_.subscribe("initialpose", 2, &AmclNode::initialPoseReceived, this);//通過函數(shù)initialPoseReceived訂閱初始化的位姿 }void AmclNode::runFromBag(const std::string &in_bag_fn, bool trigger_global_localization) {//根據(jù)信息記錄包來運(yùn)行amclrosbag::Bag bag;bag.open(in_bag_fn, rosbag::bagmode::Read);std::vector<std::string> topics;topics.push_back(std::string("tf"));std::string scan_topic_name = "base_scan"; // TODO determine what topic this actually is from ROStopics.push_back(scan_topic_name);rosbag::View view(bag, rosbag::TopicQuery(topics));ros::Publisher laser_pub = nh_.advertise<sensor_msgs::LaserScan>(scan_topic_name, 100);ros::Publisher tf_pub = nh_.advertise<tf2_msgs::TFMessage>("/tf", 100);//發(fā)布tf// Sleep for a second to let all subscribers connectros::WallDuration(1.0).sleep();ros::WallTime start(ros::WallTime::now());// Wait for mapwhile (ros::ok()){{boost::recursive_mutex::scoped_lock cfl(configuration_mutex_);if (map_){ROS_INFO("Map is ready");break;}}ROS_INFO("Waiting for map...");ros::getGlobalCallbackQueue()->callAvailable(ros::WallDuration(1.0));}if (trigger_global_localization){std_srvs::Empty empty_srv;globalLocalizationCallback(empty_srv.request, empty_srv.response);}BOOST_FOREACH(rosbag::MessageInstance const msg, view){if (!ros::ok()){break;}// Process any ros messages or callbacks at this pointros::getGlobalCallbackQueue()->callAvailable(ros::WallDuration());tf2_msgs::TFMessage::ConstPtr tf_msg = msg.instantiate<tf2_msgs::TFMessage>();if (tf_msg != NULL){tf_pub.publish(msg);for (size_t ii=0; ii<tf_msg->transforms.size(); ++ii){tf_->setTransform(tf_msg->transforms[ii], "rosbag_authority");}continue;}sensor_msgs::LaserScan::ConstPtr base_scan = msg.instantiate<sensor_msgs::LaserScan>();if (base_scan != NULL){laser_pub.publish(msg);laser_scan_filter_->add(base_scan);if (bag_scan_period_ > ros::WallDuration(0)){bag_scan_period_.sleep();}continue;}ROS_WARN_STREAM("Unsupported message type" << msg.getTopic());}bag.close();double runtime = (ros::WallTime::now() - start).toSec();ROS_INFO("Bag complete, took %.1f seconds to process, shutting down", runtime);const geometry_msgs::Quaternion & q(last_published_pose.pose.pose.orientation);double yaw, pitch, roll;tf2::Matrix3x3(tf2::Quaternion(q.x, q.y, q.z, q.w)).getEulerYPR(yaw,pitch,roll);ROS_INFO("Final location %.3f, %.3f, %.3f with stamp=%f",last_published_pose.pose.pose.position.x,last_published_pose.pose.pose.position.y,yaw, last_published_pose.header.stamp.toSec());ros::shutdown(); }void AmclNode::savePoseToServer() {// We need to apply the last transform to the latest odom pose to get// the latest map pose to store. We'll take the covariance from// last_published_pose.tf2::Transform odom_pose_tf2;tf2::convert(latest_odom_pose_.pose, odom_pose_tf2);tf2::Transform map_pose = latest_tf_.inverse() * odom_pose_tf2;double yaw = tf2::getYaw(map_pose.getRotation());ROS_DEBUG("Saving pose to server. x: %.3f, y: %.3f", map_pose.getOrigin().x(), map_pose.getOrigin().y() );//amcl初始化的位置設(shè)置為了地圖的位置。private_nh_.setParam("initial_pose_x", map_pose.getOrigin().x());private_nh_.setParam("initial_pose_y", map_pose.getOrigin().y());private_nh_.setParam("initial_pose_a", yaw);private_nh_.setParam("initial_cov_xx", last_published_pose.pose.covariance[6*0+0]);private_nh_.setParam("initial_cov_yy", last_published_pose.pose.covariance[6*1+1]);private_nh_.setParam("initial_cov_aa", last_published_pose.pose.covariance[6*5+5]); }void AmclNode::updatePoseFromServer() {init_pose_[0] = 0.0;init_pose_[1] = 0.0;init_pose_[2] = 0.0;init_cov_[0] = 0.5 * 0.5;init_cov_[1] = 0.5 * 0.5;init_cov_[2] = (M_PI/12.0) * (M_PI/12.0);// Check for NAN on input from param server, #5239double tmp_pos;private_nh_.param("initial_pose_x", tmp_pos, init_pose_[0]);if(!std::isnan(tmp_pos))init_pose_[0] = tmp_pos;else ROS_WARN("ignoring NAN in initial pose X position");private_nh_.param("initial_pose_y", tmp_pos, init_pose_[1]);if(!std::isnan(tmp_pos))init_pose_[1] = tmp_pos;elseROS_WARN("ignoring NAN in initial pose Y position");private_nh_.param("initial_pose_a", tmp_pos, init_pose_[2]);if(!std::isnan(tmp_pos))init_pose_[2] = tmp_pos;elseROS_WARN("ignoring NAN in initial pose Yaw");private_nh_.param("initial_cov_xx", tmp_pos, init_cov_[0]);if(!std::isnan(tmp_pos))init_cov_[0] =tmp_pos;elseROS_WARN("ignoring NAN in initial covariance XX");private_nh_.param("initial_cov_yy", tmp_pos, init_cov_[1]);if(!std::isnan(tmp_pos))init_cov_[1] = tmp_pos;elseROS_WARN("ignoring NAN in initial covariance YY");private_nh_.param("initial_cov_aa", tmp_pos, init_cov_[2]);if(!std::isnan(tmp_pos))init_cov_[2] = tmp_pos;elseROS_WARN("ignoring NAN in initial covariance AA"); }void AmclNode::checkLaserReceived(const ros::TimerEvent& event) {ros::Duration d = ros::Time::now() - last_laser_received_ts_;if(d > laser_check_interval_){ROS_WARN("No laser scan received (and thus no pose updates have been published) for %f seconds. Verify that data is being published on the %s topic.",d.toSec(),ros::names::resolve(scan_topic_).c_str());} }void AmclNode::requestMap() {boost::recursive_mutex::scoped_lock ml(configuration_mutex_);// get map via RPCnav_msgs::GetMap::Request req;nav_msgs::GetMap::Response resp;ROS_INFO("Requesting the map...");while(!ros::service::call("static_map", req, resp)){ROS_WARN("Request for map failed; trying again...");ros::Duration d(0.5);d.sleep();}handleMapMessage( resp.map ); }void AmclNode::mapReceived(const nav_msgs::OccupancyGridConstPtr& msg) {if( first_map_only_ && first_map_received_ ) {return;}handleMapMessage( *msg );first_map_received_ = true; }void AmclNode::handleMapMessage(const nav_msgs::OccupancyGrid& msg) {boost::recursive_mutex::scoped_lock cfl(configuration_mutex_);ROS_INFO("Received a %d X %d map @ %.3f m/pix\n",msg.info.width,msg.info.height,msg.info.resolution);if(msg.header.frame_id != global_frame_id_)ROS_WARN("Frame_id of map received:'%s' doesn't match global_frame_id:'%s'. This could cause issues with reading published topics",msg.header.frame_id.c_str(),global_frame_id_.c_str());freeMapDependentMemory();// Clear queued laser objects because they hold pointers to the existing// map, #5202.lasers_.clear();lasers_update_.clear();frame_to_laser_.clear();map_ = convertMap(msg);#if NEW_UNIFORM_SAMPLING// Index of free spacefree_space_indices.resize(0);for(int i = 0; i < map_->size_x; i++)for(int j = 0; j < map_->size_y; j++)if(map_->cells[MAP_INDEX(map_,i,j)].occ_state == -1)free_space_indices.push_back(std::make_pair(i,j)); #endif// Create the particle filterpf_ = pf_alloc(min_particles_, max_particles_,alpha_slow_, alpha_fast_,(pf_init_model_fn_t)AmclNode::uniformPoseGenerator,(void *)map_);pf_set_selective_resampling(pf_, selective_resampling_);pf_->pop_err = pf_err_;pf_->pop_z = pf_z_;// Initialize the filterupdatePoseFromServer();pf_vector_t pf_init_pose_mean = pf_vector_zero();pf_init_pose_mean.v[0] = init_pose_[0];pf_init_pose_mean.v[1] = init_pose_[1];pf_init_pose_mean.v[2] = init_pose_[2];pf_matrix_t pf_init_pose_cov = pf_matrix_zero();pf_init_pose_cov.m[0][0] = init_cov_[0];pf_init_pose_cov.m[1][1] = init_cov_[1];pf_init_pose_cov.m[2][2] = init_cov_[2];pf_init(pf_, pf_init_pose_mean, pf_init_pose_cov);pf_init_ = false;// Instantiate the sensor objects// Odometrydelete odom_;odom_ = new AMCLOdom();ROS_ASSERT(odom_);odom_->SetModel( odom_model_type_, alpha1_, alpha2_, alpha3_, alpha4_, alpha5_ );// Laserdelete laser_;laser_ = new AMCLLaser(max_beams_, map_);ROS_ASSERT(laser_);if(laser_model_type_ == LASER_MODEL_BEAM)laser_->SetModelBeam(z_hit_, z_short_, z_max_, z_rand_,sigma_hit_, lambda_short_, 0.0);else if(laser_model_type_ == LASER_MODEL_LIKELIHOOD_FIELD_PROB){ROS_INFO("Initializing likelihood field model; this can take some time on large maps...");laser_->SetModelLikelihoodFieldProb(z_hit_, z_rand_, sigma_hit_,laser_likelihood_max_dist_, do_beamskip_, beam_skip_distance_, beam_skip_threshold_, beam_skip_error_threshold_);ROS_INFO("Done initializing likelihood field model.");}else{ROS_INFO("Initializing likelihood field model; this can take some time on large maps...");laser_->SetModelLikelihoodField(z_hit_, z_rand_, sigma_hit_,laser_likelihood_max_dist_);ROS_INFO("Done initializing likelihood field model.");}// In case the initial pose message arrived before the first map,// try to apply the initial pose now that the map has arrived.applyInitialPose();}void AmclNode::freeMapDependentMemory() {if( map_ != NULL ) {map_free( map_ );map_ = NULL;}if( pf_ != NULL ) {pf_free( pf_ );pf_ = NULL;}delete odom_;odom_ = NULL;delete laser_;laser_ = NULL; }/*** Convert an OccupancyGrid map message into the internal* representation. This allocates a map_t and returns it.*/ map_t* AmclNode::convertMap( const nav_msgs::OccupancyGrid& map_msg ) {map_t* map = map_alloc();ROS_ASSERT(map);map->size_x = map_msg.info.width;map->size_y = map_msg.info.height;map->scale = map_msg.info.resolution;map->origin_x = map_msg.info.origin.position.x + (map->size_x / 2) * map->scale;map->origin_y = map_msg.info.origin.position.y + (map->size_y / 2) * map->scale;// Convert to player formatmap->cells = (map_cell_t*)malloc(sizeof(map_cell_t)*map->size_x*map->size_y);ROS_ASSERT(map->cells);for(int i=0;i<map->size_x * map->size_y;i++){if(map_msg.data[i] == 0)map->cells[i].occ_state = -1;else if(map_msg.data[i] == 100)map->cells[i].occ_state = +1;elsemap->cells[i].occ_state = 0;}return map; }AmclNode::~AmclNode() {delete dsrv_;freeMapDependentMemory();delete laser_scan_filter_;delete laser_scan_sub_;// TODO: delete everything allocated in constructor }// 獲取機(jī)器人當(dāng)前的pose(x,y,方向) //注意,此處是通過tf來獲取的!!!!!! bool AmclNode::getOdomPose(geometry_msgs::PoseStamped& odom_pose,double& x, double& y, double& yaw,const ros::Time& t, const std::string& f) {//這里所謂的odom_pose其實(shí)就是odom->base_link //odom->base_link就是以輪式里程計(jì)來看base_link的運(yùn)動(dòng) //里程計(jì)開始計(jì)數(shù)的位置// Get the robot's posegeometry_msgs::PoseStamped ident;ident.header.frame_id = stripSlash(f);ident.header.stamp = t;tf2::toMsg(tf2::Transform::getIdentity(), ident.pose);try{this->tf_->transform(ident, odom_pose, odom_frame_id_);}catch(tf2::TransformException e){ROS_WARN("Failed to compute odom pose, skipping scan (%s)", e.what());return false;}x = odom_pose.pose.position.x;y = odom_pose.pose.position.y;yaw = tf2::getYaw(odom_pose.pose.orientation);return true; }pf_vector_t AmclNode::uniformPoseGenerator(void* arg) {map_t* map = (map_t*)arg; #if NEW_UNIFORM_SAMPLINGunsigned int rand_index = drand48() * free_space_indices.size();std::pair<int,int> free_point = free_space_indices[rand_index];pf_vector_t p;p.v[0] = MAP_WXGX(map, free_point.first);p.v[1] = MAP_WYGY(map, free_point.second);p.v[2] = drand48() * 2 * M_PI - M_PI; #elsedouble min_x, max_x, min_y, max_y;min_x = (map->size_x * map->scale)/2.0 - map->origin_x;max_x = (map->size_x * map->scale)/2.0 + map->origin_x;min_y = (map->size_y * map->scale)/2.0 - map->origin_y;max_y = (map->size_y * map->scale)/2.0 + map->origin_y;pf_vector_t p;ROS_DEBUG("Generating new uniform sample");for(;;){p.v[0] = min_x + drand48() * (max_x - min_x);p.v[1] = min_y + drand48() * (max_y - min_y);p.v[2] = drand48() * 2 * M_PI - M_PI;// Check that it's a free cellint i,j;i = MAP_GXWX(map, p.v[0]);j = MAP_GYWY(map, p.v[1]);if(MAP_VALID(map,i,j) && (map->cells[MAP_INDEX(map,i,j)].occ_state == -1))break;} #endifreturn p; }bool AmclNode::globalLocalizationCallback(std_srvs::Empty::Request& req,std_srvs::Empty::Response& res) {if( map_ == NULL ) {return true;}boost::recursive_mutex::scoped_lock gl(configuration_mutex_);ROS_INFO("Initializing with uniform distribution");pf_init_model(pf_, (pf_init_model_fn_t)AmclNode::uniformPoseGenerator,(void *)map_);ROS_INFO("Global initialisation done!");pf_init_ = false;return true; }// force nomotion updates (amcl updating without requiring motion) //無需運(yùn)動(dòng)即可更新amcl bool AmclNode::nomotionUpdateCallback(std_srvs::Empty::Request& req,std_srvs::Empty::Response& res) {m_force_update = true;//ROS_INFO("Requesting no-motion update");return true; }bool AmclNode::setMapCallback(nav_msgs::SetMap::Request& req,nav_msgs::SetMap::Response& res) {handleMapMessage(req.map);handleInitialPoseMessage(req.initial_pose);//對(duì)位置初始化res.success = true;return true; }//實(shí)現(xiàn)將激光雷達(dá)信息的tf轉(zhuǎn)換 //Amcl的業(yè)務(wù)邏輯總體就是在一個(gè)四五百行的巨大函數(shù)laserReceived中實(shí)現(xiàn)的 ///這個(gè)函數(shù)才是最重要的!!!!!! void AmclNode::laserReceived(const sensor_msgs::LaserScanConstPtr& laser_scan) {//函數(shù)以雷達(dá)掃描數(shù)據(jù)的指針為輸入?yún)?shù)std::string laser_scan_frame_id = stripSlash(laser_scan->header.frame_id);//從掃描數(shù)據(jù)中獲取雷達(dá)消息ID//同時(shí)用變量last_laser_received_ts_記錄下當(dāng)前的系統(tǒng)時(shí)間, 它用于判定是否長(zhǎng)時(shí)間未接收到雷達(dá)數(shù)據(jù)last_laser_received_ts_ = ros::Time::now();if( map_ == NULL ) {//如果沒有地圖,也會(huì)直接退出return;}boost::recursive_mutex::scoped_lock lr(configuration_mutex_);int laser_index = -1;//AmclNode借助一些vector和map的容器,來支持多傳感器// Do we have the base->base_laser Tx yet?通過在frame_to_laser_中查找雷達(dá)的坐標(biāo)ID//如果之前沒有收到該雷達(dá)的消息,將新建一個(gè)對(duì)象保存在lasers_中,// 并相應(yīng)的在lasers_update_中添加對(duì)應(yīng)更新狀態(tài), 建立映射關(guān)系。if(frame_to_laser_.find(laser_scan_frame_id) == frame_to_laser_.end()){ROS_DEBUG("Setting up laser %d (frame_id=%s)\n", (int)frame_to_laser_.size(), laser_scan_frame_id.c_str());lasers_.push_back(new AMCLLaser(*laser_));//記錄下當(dāng)前構(gòu)建的雷達(dá)對(duì)象lasers_update_.push_back(true);//標(biāo)記雷達(dá)的更新狀態(tài)laser_index = frame_to_laser_.size();geometry_msgs::PoseStamped ident;ident.header.frame_id = laser_scan_frame_id;ident.header.stamp = ros::Time();tf2::toMsg(tf2::Transform::getIdentity(), ident.pose);geometry_msgs::PoseStamped laser_pose;try{this->tf_->transform(ident, laser_pose, base_frame_id_);}catch(tf2::TransformException& e){ROS_ERROR("Couldn't transform from %s to %s, ""even though the message notifier is in use",laser_scan_frame_id.c_str(),base_frame_id_.c_str());return;}pf_vector_t laser_pose_v;laser_pose_v.v[0] = laser_pose.pose.position.x;laser_pose_v.v[1] = laser_pose.pose.position.y;// laser mounting angle gets computed later -> set to 0 here!laser_pose_v.v[2] = 0;lasers_[laser_index]->SetLaserPose(laser_pose_v);ROS_DEBUG("Received laser's pose wrt robot: %.3f %.3f %.3f",laser_pose_v.v[0],laser_pose_v.v[1],laser_pose_v.v[2]);frame_to_laser_[laser_scan_frame_id] = laser_index;} else {//直接通過frame_to_laser_獲取雷達(dá)對(duì)象的索引// we have the laser pose, retrieve laser indexlaser_index = frame_to_laser_[laser_scan_frame_id];//直接通過frame_to_laser_獲取雷達(dá)對(duì)象的索引}// Where was the robot when this scan was taken?pf_vector_t pose;//AmclNode通過成員函數(shù)getOdomPose獲取接收到雷達(dá)數(shù)據(jù)時(shí)刻的里程計(jì)位姿if(!getOdomPose(latest_odom_pose_, pose.v[0], pose.v[1], pose.v[2],laser_scan->header.stamp, base_frame_id_)){ROS_ERROR("Couldn't determine robot's pose associated with laser scan");return;}pf_vector_t delta = pf_vector_zero(); //如果里程計(jì)的數(shù)據(jù)顯示機(jī)器人已經(jīng)發(fā)生了明顯的位移或者旋轉(zhuǎn),就標(biāo)記所有的雷達(dá)更新標(biāo)記為true。 //pf_init_這個(gè)應(yīng)該只是初始化的標(biāo)志位if(pf_init_){// Compute change in pose//看看機(jī)器人是否已經(jīng)發(fā)生了明顯的位移或旋轉(zhuǎn)//delta = pf_vector_coord_sub(pose, pf_odom_pose_);delta.v[0] = pose.v[0] - pf_odom_pose_.v[0];delta.v[1] = pose.v[1] - pf_odom_pose_.v[1];delta.v[2] = angle_diff(pose.v[2], pf_odom_pose_.v[2]);// See if we should update the filterbool update = fabs(delta.v[0]) > d_thresh_ ||fabs(delta.v[1]) > d_thresh_ ||fabs(delta.v[2]) > a_thresh_;update = update || m_force_update;m_force_update=false;// Set the laser update flagsif(update)for(unsigned int i=0; i < lasers_update_.size(); i++)lasers_update_[i] = true;}bool force_publication = false;if(!pf_init_){// Pose at last filter updatepf_odom_pose_ = pose;// Filter is now initializedpf_init_ = true;// Should update sensor datafor(unsigned int i=0; i < lasers_update_.size(); i++)lasers_update_[i] = true;force_publication = true;resample_count_ = 0;}// If the robot has moved, update the filterelse if(pf_init_ && lasers_update_[laser_index]){//根據(jù)剛剛更新標(biāo)識(shí),我們調(diào)用里程計(jì)對(duì)象的UpdateAction接口完成運(yùn)動(dòng)更新。//這應(yīng)該是粒子濾波中,做狀態(tài)更新的//printf("pose\n");//pf_vector_fprintf(pose, stdout, "%.3f");AMCLOdomData odata;odata.pose = pose;// HACK// Modify the delta in the action data so the filter gets// updated correctlyodata.delta = delta;// Use the action data to update the filterodom_->UpdateAction(pf_, (AMCLSensorData*)&odata);// Pose at last filter update//this->pf_odom_pose = pose;}bool resampled = false;// If the robot has moved, update the filterif(lasers_update_[laser_index]){//接下來,我們根據(jù)激光的掃描數(shù)據(jù)更新濾波器AMCLLaserData ldata;//首先構(gòu)建AMCLLaserData對(duì)象ldata.sensor = lasers_[laser_index];//指定傳感器對(duì)象ldata.range_count = laser_scan->ranges.size();//指定傳感器量程//下面實(shí)現(xiàn)將接收到的傳感器數(shù)據(jù)拷貝到ldata對(duì)象中// To account for lasers that are mounted upside-down, we determine the// min, max, and increment angles of the laser in the base frame.//// Construct min and max angles of laser, in the base_link frame.tf2::Quaternion q;q.setRPY(0.0, 0.0, laser_scan->angle_min);geometry_msgs::QuaternionStamped min_q, inc_q;min_q.header.stamp = laser_scan->header.stamp;min_q.header.frame_id = stripSlash(laser_scan->header.frame_id);tf2::convert(q, min_q.quaternion);q.setRPY(0.0, 0.0, laser_scan->angle_min + laser_scan->angle_increment);inc_q.header = min_q.header;tf2::convert(q, inc_q.quaternion);try{tf_->transform(min_q, min_q, base_frame_id_);tf_->transform(inc_q, inc_q, base_frame_id_);}catch(tf2::TransformException& e){ROS_WARN("Unable to transform min/max laser angles into base frame: %s",e.what());return;}double angle_min = tf2::getYaw(min_q.quaternion);double angle_increment = tf2::getYaw(inc_q.quaternion) - angle_min;// wrapping angle to [-pi .. pi]angle_increment = fmod(angle_increment + 5*M_PI, 2*M_PI) - M_PI;ROS_DEBUG("Laser %d angles in base frame: min: %.3f inc: %.3f", laser_index, angle_min, angle_increment);// Apply range min/max thresholds, if the user supplied themif(laser_max_range_ > 0.0)ldata.range_max = std::min(laser_scan->range_max, (float)laser_max_range_);elseldata.range_max = laser_scan->range_max;double range_min;if(laser_min_range_ > 0.0)range_min = std::max(laser_scan->range_min, (float)laser_min_range_);elserange_min = laser_scan->range_min;// The AMCLLaserData destructor will free this memoryldata.ranges = new double[ldata.range_count][2];ROS_ASSERT(ldata.ranges);for(int i=0;i<ldata.range_count;i++){// amcl doesn't (yet) have a concept of min range. So we'll map short// readings to max range.if(laser_scan->ranges[i] <= range_min)ldata.ranges[i][0] = ldata.range_max;elseldata.ranges[i][0] = laser_scan->ranges[i];// Compute bearingldata.ranges[i][1] = angle_min +(i * angle_increment);}//最后通過激光傳感器對(duì)象的UpdateSensor接口完成粒子濾波器的測(cè)量更新。lasers_[laser_index]->UpdateSensor(pf_, (AMCLSensorData*)&ldata);//前面是狀態(tài)預(yù)測(cè),這里是測(cè)量更新lasers_update_[laser_index] = false;//然后清除lasers_update_的標(biāo)記pf_odom_pose_ = pose;//并更新濾波器的里程計(jì)位姿。(此處的粒子濾波的姿態(tài)是優(yōu)化后的)“pf_odom_pose_”// Resample the particles//根據(jù)重采樣計(jì)數(shù)器和設(shè)定的閾值,觸發(fā)對(duì)濾波器的重采樣if(!(++resample_count_ % resample_interval_)){pf_update_resample(pf_);resampled = true;}至此,amcl的主業(yè)務(wù)邏輯就基本結(jié)束了。//剩下的內(nèi)容就是:獲取當(dāng)前的粒子集合,從中提煉出機(jī)器人的位姿估計(jì),并將這些消息通過各個(gè)主題發(fā)布出去。pf_sample_set_t* set = pf_->sets + pf_->current_set;ROS_DEBUG("Num samples: %d\n", set->sample_count);// Publish the resulting cloud// TODO: set maximum rate for publishingif (!m_force_update){geometry_msgs::PoseArray cloud_msg;cloud_msg.header.stamp = ros::Time::now();cloud_msg.header.frame_id = global_frame_id_;cloud_msg.poses.resize(set->sample_count);for(int i=0;i<set->sample_count;i++){cloud_msg.poses[i].position.x = set->samples[i].pose.v[0];cloud_msg.poses[i].position.y = set->samples[i].pose.v[1];cloud_msg.poses[i].position.z = 0;tf2::Quaternion q;q.setRPY(0, 0, set->samples[i].pose.v[2]);tf2::convert(q, cloud_msg.poses[i].orientation);}particlecloud_pub_.publish(cloud_msg);//將粒子點(diǎn)群數(shù)據(jù)發(fā)布出去}}if(resampled || force_publication){// Read out the current hypotheses(假設(shè) )double max_weight = 0.0;int max_weight_hyp = -1;std::vector<amcl_hyp_t> hyps;hyps.resize(pf_->sets[pf_->current_set].cluster_count);for(int hyp_count = 0;hyp_count < pf_->sets[pf_->current_set].cluster_count; hyp_count++){double weight;pf_vector_t pose_mean;pf_matrix_t pose_cov;if (!pf_get_cluster_stats(pf_, hyp_count, &weight, &pose_mean, &pose_cov)){ROS_ERROR("Couldn't get stats on cluster %d", hyp_count);break;}hyps[hyp_count].weight = weight;hyps[hyp_count].pf_pose_mean = pose_mean;hyps[hyp_count].pf_pose_cov = pose_cov;if(hyps[hyp_count].weight > max_weight){max_weight = hyps[hyp_count].weight;max_weight_hyp = hyp_count;}}if(max_weight > 0.0){ROS_DEBUG("Max weight pose: %.3f %.3f %.3f",hyps[max_weight_hyp].pf_pose_mean.v[0],hyps[max_weight_hyp].pf_pose_mean.v[1],hyps[max_weight_hyp].pf_pose_mean.v[2]);/*puts("");pf_matrix_fprintf(hyps[max_weight_hyp].pf_pose_cov, stdout, "%6.3f");puts("");*/geometry_msgs::PoseWithCovarianceStamped p;//這個(gè)是amcl要發(fā)布的“amcl_pose”// Fill in the headerp.header.frame_id = global_frame_id_;p.header.stamp = laser_scan->header.stamp;// Copy in the posep.pose.pose.position.x = hyps[max_weight_hyp].pf_pose_mean.v[0];p.pose.pose.position.y = hyps[max_weight_hyp].pf_pose_mean.v[1];tf2::Quaternion q;//四元素q.setRPY(0, 0, hyps[max_weight_hyp].pf_pose_mean.v[2]);tf2::convert(q, p.pose.pose.orientation);// Copy in the covariance, converting from 3-D to 6-Dpf_sample_set_t* set = pf_->sets + pf_->current_set;for(int i=0; i<2; i++){for(int j=0; j<2; j++){// Report the overall filter covariance, rather than the// covariance for the highest-weight cluster//p.covariance[6*i+j] = hyps[max_weight_hyp].pf_pose_cov.m[i][j];p.pose.covariance[6*i+j] = set->cov.m[i][j];}}// Report the overall filter covariance, rather than the// covariance for the highest-weight cluster//p.covariance[6*5+5] = hyps[max_weight_hyp].pf_pose_cov.m[2][2];p.pose.covariance[6*5+5] = set->cov.m[2][2];/*printf("cov:\n");for(int i=0; i<6; i++){for(int j=0; j<6; j++)printf("%6.3f ", p.covariance[6*i+j]);puts("");}*/pose_pub_.publish(p);//將位姿發(fā)布出來last_published_pose = p;//AMCL算法是根據(jù)給定的地圖,結(jié)合粒子濾波獲取最佳定位點(diǎn)Mp,//這個(gè)定位點(diǎn)是相對(duì)于地圖map上的坐標(biāo),也就是base_link(也就是機(jī)器人的基坐標(biāo))相對(duì)map上的坐標(biāo)。// odom 的原點(diǎn)是機(jī)器人啟動(dòng)時(shí)刻的位置,它在map上的位置或轉(zhuǎn)換矩陣是未知的//AMCL可以根據(jù)最佳粒子的位置推算出 odom->map(就是說通過最佳粒子推算出機(jī)器人在地圖的中的位置)的tf轉(zhuǎn)換信息并發(fā)布到 tf主題上//因?yàn)閎ase_link->odom的tf轉(zhuǎn)換信息是每時(shí)每刻都在發(fā)布的,所以它是已知的//故此,此處有如下的關(guān)系://map->base_link(就是地圖中機(jī)器人的位置 是根據(jù)最佳粒子推算的)//base_link->odom(這是現(xiàn)實(shí)中機(jī)器人的位姿可以說是里程計(jì)的信息),由turtlebot3_core發(fā)布//故此有://機(jī)器人的里程計(jì)的信息 = 當(dāng)前地圖中的機(jī)器人的的位置 減去 地圖中機(jī)器人的起點(diǎn)位置// 轉(zhuǎn)為公式可以寫成 :map->odom = map->base_link - base_link->odom//或者寫為:base_link->odom = map->base_link - map->odom //所謂的base_link->odom是里程計(jì)開始計(jì)數(shù)的位置,用里程計(jì)來看base_link的運(yùn)動(dòng)//就是amcl中宣稱,它檢測(cè)的map->base_link是準(zhǔn)的//然后odom->base_link,這個(gè)里面其實(shí)是已經(jīng)有了odom的累計(jì)誤差了//那前者減去后者。那就是把累積誤差也減去了,所以得到的map->odom就準(zhǔn)確了//就是累計(jì)誤差已經(jīng)在odom->base_link中體現(xiàn)了,所以map->odom就沒有累積誤差ROS_DEBUG("New pose: %6.3f %6.3f %6.3f",hyps[max_weight_hyp].pf_pose_mean.v[0],hyps[max_weight_hyp].pf_pose_mean.v[1],hyps[max_weight_hyp].pf_pose_mean.v[2]);//hyps[max_weight_hyp].pf_pose_mean.v[0], [1], [2] 就代表了粒子濾波獲取最佳定位點(diǎn)Mp//也就是機(jī)器人的位姿那么位姿的格式是(x,y,theta)最后一個(gè)是yaw偏航角,// subtracting base to odom from map to base and send map to odom insteadgeometry_msgs::PoseStamped odom_to_map;//要發(fā)布的,從odom到map上的try{tf2::Quaternion q;//定義一個(gè)新的四元素q.setRPY(0, 0, hyps[max_weight_hyp].pf_pose_mean.v[2]);tf2::Transform tmp_tf(q, tf2::Vector3(hyps[max_weight_hyp].pf_pose_mean.v[0],hyps[max_weight_hyp].pf_pose_mean.v[1],0.0)); tmp_tf = 從map原點(diǎn)看base_link的位置 為yaw生成四元數(shù),最后的0.0是(x,y,z)的Z的值為0 因?yàn)檫@是在二維平面中g(shù)eometry_msgs::PoseStamped tmp_tf_stamped;tmp_tf_stamped.header.frame_id = base_frame_id_;tmp_tf_stamped.header.stamp = laser_scan->header.stamp;tf2::toMsg(tmp_tf.inverse(), tmp_tf_stamped.pose); tmp_tf.inverse()=最佳定位點(diǎn)Mp為坐標(biāo)的原點(diǎn),地圖map原點(diǎn)相對(duì)于Mp的位置,就是在base_link坐標(biāo)系下map原點(diǎn)的位置this->tf_->transform(tmp_tf_stamped, odom_to_map, odom_frame_id_);//進(jìn)行 base_link坐標(biāo)系下的點(diǎn)轉(zhuǎn)換到odom坐標(biāo)系,//也就是把map原點(diǎn)轉(zhuǎn)換到odom坐標(biāo)系下,等于從odom原點(diǎn)看map原點(diǎn)的位置。放到latest_tf_變量里面}catch(tf2::TransformException){ROS_DEBUG("Failed to subtract base to odom transform");return;}tf2::convert(odom_to_map.pose, latest_tf_);latest_tf_valid_ = true;if (tf_broadcast_ == true){// We want to send a transform that is good up until a// tolerance time so that odom can be usedros::Time transform_expiration = (laser_scan->header.stamp +transform_tolerance_);geometry_msgs::TransformStamped tmp_tf_stamped;tmp_tf_stamped.header.frame_id = global_frame_id_;tmp_tf_stamped.header.stamp = transform_expiration;tmp_tf_stamped.child_frame_id = odom_frame_id_;tf2::convert(latest_tf_.inverse(), tmp_tf_stamped.transform);this->tfb_->sendTransform(tmp_tf_stamped);sent_first_transform_ = true;}}else{ROS_ERROR("No pose!");}}else if(latest_tf_valid_){if (tf_broadcast_ == true){// Nothing changed, so we'll just republish the last transform, to keep// everybody happy.ros::Time transform_expiration = (laser_scan->header.stamp +transform_tolerance_);geometry_msgs::TransformStamped tmp_tf_stamped;tmp_tf_stamped.header.frame_id = global_frame_id_;tmp_tf_stamped.header.stamp = transform_expiration;tmp_tf_stamped.child_frame_id = odom_frame_id_;tf2::convert(latest_tf_.inverse(), tmp_tf_stamped.transform);this->tfb_->sendTransform(tmp_tf_stamped);}// Is it time to save our last pose to the param serverros::Time now = ros::Time::now();if((save_pose_period.toSec() > 0.0) &&(now - save_pose_last_time) >= save_pose_period){this->savePoseToServer();save_pose_last_time = now;}}diagnosic_updater_.update(); }//訂閱初始位姿,并調(diào)用initialPoseReceived函數(shù) //接收初始化的位姿數(shù)據(jù) void AmclNode::initialPoseReceived(const geometry_msgs::PoseWithCovarianceStampedConstPtr& msg) {handleInitialPoseMessage(*msg);//調(diào)用處理函數(shù) }//進(jìn)行的位姿初始化 //代碼部分實(shí)現(xiàn)將坐標(biāo)變換到map對(duì)應(yīng)的坐標(biāo)中的功能 //輸入?yún)?shù)為一個(gè)pose //http://docs.ros.org/en/melodic/api/geometry_msgs/html/msg/PoseWithCovarianceStamped.html void AmclNode::handleInitialPoseMessage(const geometry_msgs::PoseWithCovarianceStamped& msg) {//根據(jù)接收的初始位姿消息,在該位姿附近高斯采樣重新生成粒子 //只接受global_frame_id_(一般為map)的坐標(biāo),并重新生成粒子。在接收到的初始位姿附近采樣生成粒子集。boost::recursive_mutex::scoped_lock prl(configuration_mutex_);if(msg.header.frame_id == ""){// This should be removed at some pointROS_WARN("Received initial pose with empty frame_id. You should always supply a frame_id.");}// We only accept initial pose estimates in the global frame, #5148.else if(stripSlash(msg.header.frame_id) != global_frame_id_){ROS_WARN("Ignoring initial pose in frame \"%s\"; initial poses must be in the global frame, \"%s\"",stripSlash(msg.header.frame_id).c_str(),global_frame_id_.c_str());return;}// In case the client sent us a pose estimate in the past, integrate the// intervening odometric change.geometry_msgs::TransformStamped tx_odom;try{ros::Time now = ros::Time::now();// wait a little for the latest tf to become availabletx_odom = tf_->lookupTransform(base_frame_id_, msg.header.stamp,base_frame_id_, ros::Time::now(),odom_frame_id_, ros::Duration(0.5));}catch(tf2::TransformException e){// If we've never sent a transform, then this is normal, because the// global_frame_id_ frame doesn't exist. We only care about in-time// transformation for on-the-move pose-setting, so ignoring this// startup condition doesn't really cost us anything.if(sent_first_transform_)ROS_WARN("Failed to transform initial pose in time (%s)", e.what());tf2::convert(tf2::Transform::getIdentity(), tx_odom.transform);}tf2::Transform tx_odom_tf2;tf2::convert(tx_odom.transform, tx_odom_tf2);tf2::Transform pose_old, pose_new;tf2::convert(msg.pose.pose, pose_old);pose_new = pose_old * tx_odom_tf2;// Transform into the global frameROS_INFO("Setting pose (%.6f): %.3f %.3f %.3f",ros::Time::now().toSec(),pose_new.getOrigin().x(),pose_new.getOrigin().y(),tf2::getYaw(pose_new.getRotation()));// Re-initialize the filterpf_vector_t pf_init_pose_mean = pf_vector_zero();pf_init_pose_mean.v[0] = pose_new.getOrigin().x();pf_init_pose_mean.v[1] = pose_new.getOrigin().y();pf_init_pose_mean.v[2] = tf2::getYaw(pose_new.getRotation());pf_matrix_t pf_init_pose_cov = pf_matrix_zero();// Copy in the covariance, converting from 6-D to 3-Dfor(int i=0; i<2; i++){for(int j=0; j<2; j++){pf_init_pose_cov.m[i][j] = msg.pose.covariance[6*i+j];}}pf_init_pose_cov.m[2][2] = msg.pose.covariance[6*5+5];delete initial_pose_hyp_;initial_pose_hyp_ = new amcl_hyp_t();initial_pose_hyp_->pf_pose_mean = pf_init_pose_mean;initial_pose_hyp_->pf_pose_cov = pf_init_pose_cov;applyInitialPose(); }/*** If initial_pose_hyp_ and map_ are both non-null, apply the initial* pose to the particle filter state. initial_pose_hyp_ is deleted* and set to NULL after it is used.*/ void AmclNode::applyInitialPose() {boost::recursive_mutex::scoped_lock cfl(configuration_mutex_);if( initial_pose_hyp_ != NULL && map_ != NULL ) {pf_init(pf_, initial_pose_hyp_->pf_pose_mean, initial_pose_hyp_->pf_pose_cov);pf_init_ = false;delete initial_pose_hyp_;initial_pose_hyp_ = NULL;} }void AmclNode::standardDeviationDiagnostics(diagnostic_updater::DiagnosticStatusWrapper& diagnostic_status) {double std_x = sqrt(last_published_pose.pose.covariance[6*0+0]);double std_y = sqrt(last_published_pose.pose.covariance[6*1+1]);double std_yaw = sqrt(last_published_pose.pose.covariance[6*5+5]);diagnostic_status.add("std_x", std_x);diagnostic_status.add("std_y", std_y);diagnostic_status.add("std_yaw", std_yaw);diagnostic_status.add("std_warn_level_x", std_warn_level_x_);diagnostic_status.add("std_warn_level_y", std_warn_level_y_);diagnostic_status.add("std_warn_level_yaw", std_warn_level_yaw_);if (std_x > std_warn_level_x_ || std_y > std_warn_level_y_ || std_yaw > std_warn_level_yaw_){diagnostic_status.summary(diagnostic_msgs::DiagnosticStatus::WARN, "Too large");}else{diagnostic_status.summary(diagnostic_msgs::DiagnosticStatus::OK, "OK");} }

?

TF

在amcl中,tf同樣是非常重要的。

http://wiki.ros.org/tf/Overview

TF的數(shù)據(jù)類型

作為ROS中的特色,TF定義了它自己的數(shù)據(jù)類型

tf的消息類型為geometry_msgs/TransformStamped

std_mags/Header headeruint32 seqtime stampstring frame_id string child_frame_id geometry_msgs/Transform transformgeometry_msgs/Vector3 translationfloat64 xfloat64 yfloat64 zgeometry_msgs/Quaternion rotationfloat64 xfloat64 yflaot64 zfloat64 w

?

發(fā)布TF

為了在一個(gè)節(jié)點(diǎn)中發(fā)布tf,需要采用tf::TransformBroadcaster

首先需要構(gòu)建一個(gè)tf的發(fā)布者

tf::TransformBroadcaster();

然后發(fā)送

void sendTransform(const StampedTransform & transform); void sendTransform(const geometry_msgs::TransformStamped & transform);

例程:下面代碼將發(fā)布一個(gè)坐標(biāo)幀到tf2上,隨著turtles的運(yùn)動(dòng),發(fā)布其坐標(biāo)幀

1 #include <ros/ros.h>2 #include <tf2/LinearMath/Quaternion.h>3 #include <tf2_ros/transform_broadcaster.h> //發(fā)布pose到tf2的頭文件4 #include <geometry_msgs/TransformStamped.h>5 #include <turtlesim/Pose.h>6 7 std::string turtle_name;8 9 void poseCallback(const turtlesim::PoseConstPtr& msg){10 static tf2_ros::TransformBroadcaster br;//創(chuàng)建tf的發(fā)布者對(duì)象11 geometry_msgs::TransformStamped transformStamped; //要發(fā)布的信息 //關(guān)于這個(gè)消息類型可以參考:http://docs.ros.org/en/jade/api/geometry_msgs/html/msg/TransformStamped.html12 //往transformStamped里面加數(shù)據(jù)13 transformStamped.header.stamp = ros::Time::now();//需要給要發(fā)布的tf一個(gè)timestamp(時(shí)間戳)14 transformStamped.header.frame_id = "world";//設(shè)置父幀的名字15 transformStamped.child_frame_id = turtle_name;//設(shè)置子幀的名字(此處設(shè)置為自己)16 transformStamped.transform.translation.x = msg->x;17 transformStamped.transform.translation.y = msg->y;18 transformStamped.transform.translation.z = 0.0;19 tf2::Quaternion q;//定義四元數(shù)20 q.setRPY(0, 0, msg->theta);21 transformStamped.transform.rotation.x = q.x();22 transformStamped.transform.rotation.y = q.y();23 transformStamped.transform.rotation.z = q.z();24 transformStamped.transform.rotation.w = q.w();25 26 br.sendTransform(transformStamped);//然后發(fā)送出去(包含了位置和方向)27 }28 29 int main(int argc, char** argv){30 ros::init(argc, argv, "my_tf2_broadcaster");31 32 ros::NodeHandle private_node("~");33 if (! private_node.hasParam("turtle"))34 {35 if (argc != 2){ROS_ERROR("need turtle name as argument"); return -1;};36 turtle_name = argv[1];37 }38 else39 {40 private_node.getParam("turtle", turtle_name);41 }42 43 ros::NodeHandle node;44 ros::Subscriber sub = node.subscribe(turtle_name+"/pose", 10, &poseCallback);45 46 ros::spin();47 return 0;48 };

參考:http://wiki.ros.org/tf2/Tutorials/Writing%20a%20tf2%20broadcaster%20%28C%2B%2B%29#CA-655844a832971675d8cca09ac2412a96f7a2c2c6_1

訂閱TF

例程:

#include <ros/ros.h> #include <tf2_ros/transform_listener.h>//transform_listener就是接收tf的接收器 #include <geometry_msgs/TransformStamped.h> #include <geometry_msgs/Twist.h> #include <turtlesim/Spawn.h>int main(int argc, char** argv){ros::init(argc, argv, "my_tf2_listener");//初始化ROS節(jié)點(diǎn)ros::NodeHandle node;//創(chuàng)建節(jié)點(diǎn)句柄ros::service::waitForService("spawn");ros::ServiceClient spawner =node.serviceClient<turtlesim::Spawn>("spawn");turtlesim::Spawn turtle;turtle.request.x = 4;turtle.request.y = 2;turtle.request.theta = 0;turtle.request.name = "turtle2";spawner.call(turtle);// 創(chuàng)建發(fā)布turtle2速度控制指令的發(fā)布者ros::Publisher turtle_vel =node.advertise<geometry_msgs::Twist>("turtle2/cmd_vel", 10);//發(fā)布者tf2_ros::Buffer tfBuffer;tf2_ros::TransformListener tfListener(tfBuffer);//創(chuàng)建TransformListener的對(duì)象 //一旦tf的接收者創(chuàng)建成功,就會(huì)開始訂閱tfros::Rate rate(10.0);while (node.ok()){geometry_msgs::TransformStamped transformStamped;//接收的tf樣本try{獲取turtle1與turtle2坐標(biāo)系之間的tf數(shù)據(jù)transformStamped = tfBuffer.lookupTransform("turtle2", "turtle1",ros::Time(0));}catch (tf2::TransformException &ex) {ROS_WARN("%s",ex.what());ros::Duration(1.0).sleep();continue;}geometry_msgs::Twist vel_msg;// 根據(jù)turtle1與turtle2坐標(biāo)系之間的位置關(guān)系,發(fā)布turtle2的速度控制指令vel_msg.angular.z = 4.0 * atan2(transformStamped.transform.translation.y,transformStamped.transform.translation.x);vel_msg.linear.x = 0.5 * sqrt(pow(transformStamped.transform.translation.x, 2) +pow(transformStamped.transform.translation.y, 2));turtle_vel.publish(vel_msg);rate.sleep();}return 0; };

參考:http://wiki.ros.org/tf2/Tutorials/Writing%20a%20tf2%20listener%20%28C%2B%2B%29

關(guān)于TF的一些測(cè)試

我們先來看看正常開啟機(jī)器人的tf樹和開啟amcl后的tf樹有什么區(qū)別

正常開啟機(jī)器人的時(shí)候,由turtlebot3_core發(fā)布了一個(gè)odom到base——footprint的tf

若由amcl進(jìn)行定位。可以看到由amcl發(fā)布了一個(gè)由map到odom的

若由robot_localization進(jìn)行定位

若同時(shí)發(fā)布兩個(gè)回怎么樣呢?

結(jié)果發(fā)現(xiàn)顯示的東東跟上面一樣。

單純開啟gmapping的時(shí)候的tf

?

參考資料

http://wiki.ros.org/tf/Overview

https://www.guyuehome.com/279

https://gaoyichao.com/Xiaotu/?book=turtlebot&title=amcl%E7%9A%84%E6%80%BB%E4%BD%93%E4%B8%9A%E5%8A%A1%E9%80%BB%E8%BE%91

https://www.cnblogs.com/li-yao7758258/p/7672521.html

?

?

?

?

?

?

?

?

?

總結(jié)

以上是生活随笔為你收集整理的ROS学习笔记之——amcl源码的解读的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。

亚洲区欧美区综合区自拍区 | 亚洲aⅴ无码成人网站国产app | 国产成人无码区免费内射一片色欲 | 国产艳妇av在线观看果冻传媒 | 精品久久久久久人妻无码中文字幕 | 网友自拍区视频精品 | 日本xxxx色视频在线观看免费 | а√天堂www在线天堂小说 | 国产办公室秘书无码精品99 | 特黄特色大片免费播放器图片 | 麻豆蜜桃av蜜臀av色欲av | 国产精品久久福利网站 | 人人澡人人妻人人爽人人蜜桃 | 少妇无码av无码专区在线观看 | 18无码粉嫩小泬无套在线观看 | 亚洲熟女一区二区三区 | 国产疯狂伦交大片 | 日日摸夜夜摸狠狠摸婷婷 | 性欧美牲交xxxxx视频 | 婷婷综合久久中文字幕蜜桃三电影 | 亚洲综合色区中文字幕 | 西西人体www44rt大胆高清 | 女人高潮内射99精品 | 色综合视频一区二区三区 | 窝窝午夜理论片影院 | 97色伦图片97综合影院 | 精品熟女少妇av免费观看 | 秋霞成人午夜鲁丝一区二区三区 | 国产亚洲日韩欧美另类第八页 | 一本久道高清无码视频 | 欧美人与牲动交xxxx | 免费观看的无遮挡av | 日本成熟视频免费视频 | 丝袜美腿亚洲一区二区 | 激情国产av做激情国产爱 | 国产美女精品一区二区三区 | 欧美熟妇另类久久久久久不卡 | 久久这里只有精品视频9 | 鲁一鲁av2019在线 | 久久综合久久自在自线精品自 | 成人性做爰aaa片免费看不忠 | 一本无码人妻在中文字幕免费 | 国产农村妇女高潮大叫 | 久久久www成人免费毛片 | 男人的天堂av网站 | 欧美丰满少妇xxxx性 | 成人性做爰aaa片免费看不忠 | 亚洲日韩av一区二区三区中文 | 少妇无码吹潮 | 精品厕所偷拍各类美女tp嘘嘘 | 久久精品国产99精品亚洲 | 国产做国产爱免费视频 | 少妇无码av无码专区在线观看 | 自拍偷自拍亚洲精品10p | 荫蒂被男人添的好舒服爽免费视频 | 狠狠躁日日躁夜夜躁2020 | 亚洲综合无码一区二区三区 | 成年女人永久免费看片 | 国产无遮挡又黄又爽免费视频 | 鲁一鲁av2019在线 | 爱做久久久久久 | 狠狠色色综合网站 | 丰满岳乱妇在线观看中字无码 | 国产偷国产偷精品高清尤物 | 国产精品无码一区二区桃花视频 | 亚洲欧洲无卡二区视頻 | 最新版天堂资源中文官网 | 国产成人午夜福利在线播放 | 在线看片无码永久免费视频 | 欧美怡红院免费全部视频 | 欧美激情内射喷水高潮 | 成人一在线视频日韩国产 | 欧美亚洲日韩国产人成在线播放 | 亚洲男女内射在线播放 | 国产精品亚洲一区二区三区喷水 | 玩弄人妻少妇500系列视频 | 亚洲中文字幕av在天堂 | 国产农村妇女aaaaa视频 撕开奶罩揉吮奶头视频 | 国内揄拍国内精品人妻 | 亚洲精品国偷拍自产在线观看蜜桃 | 国产午夜亚洲精品不卡下载 | 欧美激情一区二区三区成人 | 免费网站看v片在线18禁无码 | 欧美老熟妇乱xxxxx | 成年美女黄网站色大免费视频 | 国产人妻人伦精品 | 精品久久8x国产免费观看 | 国产福利视频一区二区 | 内射老妇bbwx0c0ck | 男人和女人高潮免费网站 | aⅴ在线视频男人的天堂 | 日日碰狠狠躁久久躁蜜桃 | 午夜性刺激在线视频免费 | 亚洲熟妇色xxxxx欧美老妇y | 日本一卡二卡不卡视频查询 | 国产av久久久久精东av | 人妻体内射精一区二区三四 | 国内揄拍国内精品少妇国语 | 精品久久久久久人妻无码中文字幕 | 欧美日韩亚洲国产精品 | 日本熟妇人妻xxxxx人hd | 国产网红无码精品视频 | 久久亚洲中文字幕无码 | 国产综合色产在线精品 | 欧美大屁股xxxxhd黑色 | 国产精品二区一区二区aⅴ污介绍 | 国产精品久久精品三级 | 亚洲自偷自偷在线制服 | 动漫av网站免费观看 | 亚洲乱亚洲乱妇50p | 亚洲综合精品香蕉久久网 | 成人试看120秒体验区 | 国产精品亚洲综合色区韩国 | 国产suv精品一区二区五 | 国内精品一区二区三区不卡 | 亚洲精品一区二区三区四区五区 | 中文字幕乱码中文乱码51精品 | 一本精品99久久精品77 | 波多野结衣乳巨码无在线观看 | 日韩无码专区 | 久久久精品国产sm最大网站 | 成人一区二区免费视频 | 精品夜夜澡人妻无码av蜜桃 | 亚洲а∨天堂久久精品2021 | 熟女俱乐部五十路六十路av | 男女下面进入的视频免费午夜 | 99精品久久毛片a片 | 乱人伦人妻中文字幕无码 | 国内少妇偷人精品视频免费 | 日本一卡2卡3卡四卡精品网站 | 欧美丰满老熟妇xxxxx性 | 玩弄人妻少妇500系列视频 | 四虎国产精品免费久久 | 露脸叫床粗话东北少妇 | 国产精品无码久久av | 国产精品人人爽人人做我的可爱 | 亚洲天堂2017无码中文 | 国产精品无码久久av | 欧美日本免费一区二区三区 | 波多野结衣高清一区二区三区 | 日本免费一区二区三区最新 | 国产成人一区二区三区别 | 天堂久久天堂av色综合 | 红桃av一区二区三区在线无码av | 色欲av亚洲一区无码少妇 | 蜜桃视频韩日免费播放 | 欧美午夜特黄aaaaaa片 | 狠狠色丁香久久婷婷综合五月 | 久热国产vs视频在线观看 | 免费播放一区二区三区 | 国模大胆一区二区三区 | 国产亚洲精品久久久久久国模美 | 亚洲中文字幕在线无码一区二区 | 九九热爱视频精品 | 国产精品久久国产精品99 | 久久国产精品偷任你爽任你 | 国产成人精品久久亚洲高清不卡 | 亚洲精品一区二区三区四区五区 | 欧美成人家庭影院 | 精品国产麻豆免费人成网站 | 成人一在线视频日韩国产 | 漂亮人妻洗澡被公强 日日躁 | 亚洲の无码国产の无码步美 | 久久久www成人免费毛片 | 久久精品视频在线看15 | 精品一区二区三区无码免费视频 | 精品久久久久久人妻无码中文字幕 | 亚洲色欲久久久综合网东京热 | 在线成人www免费观看视频 | 亚洲乱码日产精品bd | 成人三级无码视频在线观看 | 久久婷婷五月综合色国产香蕉 | 亚洲高清偷拍一区二区三区 | 无遮无挡爽爽免费视频 | 国产精品亚洲一区二区三区喷水 | 老头边吃奶边弄进去呻吟 | 九九热爱视频精品 | 久久精品女人的天堂av | 九九在线中文字幕无码 | 东北女人啪啪对白 | 欧美日韩亚洲国产精品 | 欧美老人巨大xxxx做受 | 扒开双腿疯狂进出爽爽爽视频 | 麻豆果冻传媒2021精品传媒一区下载 | 性欧美大战久久久久久久 | 国产无遮挡又黄又爽又色 | 伊人久久大香线焦av综合影院 | 国产成人精品无码播放 | 漂亮人妻洗澡被公强 日日躁 | 国产午夜视频在线观看 | 狠狠躁日日躁夜夜躁2020 | 亚洲精品一区二区三区在线观看 | 久久综合久久自在自线精品自 | а天堂中文在线官网 | 成人试看120秒体验区 | 99riav国产精品视频 | 亚拍精品一区二区三区探花 | 国产精品欧美成人 | 67194成是人免费无码 | 天下第一社区视频www日本 | 丰满少妇女裸体bbw | 国产99久久精品一区二区 | 野外少妇愉情中文字幕 | ass日本丰满熟妇pics | 无码成人精品区在线观看 | 色一情一乱一伦一区二区三欧美 | 精品日本一区二区三区在线观看 | 国产电影无码午夜在线播放 | 亚洲精品午夜无码电影网 | 欧美丰满老熟妇xxxxx性 | 特级做a爰片毛片免费69 | 国内少妇偷人精品视频免费 | 伊人久久大香线焦av综合影院 | 欧美成人高清在线播放 | 亚洲国产精品毛片av不卡在线 | 欧美亚洲日韩国产人成在线播放 | 2019nv天堂香蕉在线观看 | 少女韩国电视剧在线观看完整 | 漂亮人妻洗澡被公强 日日躁 | 曰本女人与公拘交酡免费视频 | 日本熟妇人妻xxxxx人hd | 国产熟妇另类久久久久 | 精品人妻中文字幕有码在线 | 日本在线高清不卡免费播放 | 小鲜肉自慰网站xnxx | 亚洲欧美中文字幕5发布 | 免费国产成人高清在线观看网站 | 人人妻人人澡人人爽欧美一区 | 日韩人妻无码中文字幕视频 | 亚洲国产精品成人久久蜜臀 | 国产成人无码av片在线观看不卡 | 少妇无套内谢久久久久 | 亚洲精品一区二区三区四区五区 | 欧美丰满熟妇xxxx | 中文字幕乱码人妻二区三区 | 免费无码的av片在线观看 | 麻花豆传媒剧国产免费mv在线 | 久9re热视频这里只有精品 | 欧美人与禽zoz0性伦交 | 亚洲国产精品久久久天堂 | 亚洲精品综合一区二区三区在线 | 亚洲精品国产精品乱码不卡 | 免费国产成人高清在线观看网站 | 亚洲色www成人永久网址 | 精品久久综合1区2区3区激情 | 久久国产自偷自偷免费一区调 | 欧美zoozzooz性欧美 | 亚洲va欧美va天堂v国产综合 | 色老头在线一区二区三区 | 搡女人真爽免费视频大全 | 欧美性生交活xxxxxdddd | 99国产欧美久久久精品 | 日韩成人一区二区三区在线观看 | 初尝人妻少妇中文字幕 | 狂野欧美性猛xxxx乱大交 | 国内精品久久久久久中文字幕 | 中文久久乱码一区二区 | 99久久人妻精品免费二区 | 色偷偷人人澡人人爽人人模 | 亚洲另类伦春色综合小说 | 国产精品国产自线拍免费软件 | 国产后入清纯学生妹 | 黑森林福利视频导航 | 精品国产精品久久一区免费式 | 乱中年女人伦av三区 | 5858s亚洲色大成网站www | 国产美女极度色诱视频www | 日日橹狠狠爱欧美视频 | 18黄暴禁片在线观看 | 成人精品一区二区三区中文字幕 | 波多野结衣一区二区三区av免费 | 丁香花在线影院观看在线播放 | 成年美女黄网站色大免费全看 | 精品无人区无码乱码毛片国产 | 在线观看国产午夜福利片 | 国产又爽又猛又粗的视频a片 | 亚洲男人av天堂午夜在 | 无码人妻丰满熟妇区毛片18 | 一区二区三区高清视频一 | 欧美人与禽zoz0性伦交 | 久久亚洲a片com人成 | 香蕉久久久久久av成人 | 未满小14洗澡无码视频网站 | 久久精品国产一区二区三区 | 天堂久久天堂av色综合 | 亚洲精品国产精品乱码不卡 | 亚洲伊人久久精品影院 | 日产精品高潮呻吟av久久 | 成人免费视频一区二区 | 97资源共享在线视频 | 人人妻人人澡人人爽精品欧美 | 亚洲欧洲日本无在线码 | 亚洲一区二区三区偷拍女厕 | 少妇人妻大乳在线视频 | 麻豆果冻传媒2021精品传媒一区下载 | 亚洲一区二区三区国产精华液 | 乱人伦人妻中文字幕无码久久网 | 巨爆乳无码视频在线观看 | 久久精品国产精品国产精品污 | 国模大胆一区二区三区 | 国产精品理论片在线观看 | 成人欧美一区二区三区黑人 | 伊人久久大香线焦av综合影院 | 精品国产一区av天美传媒 | 欧美丰满熟妇xxxx性ppx人交 | 人妻中文无码久热丝袜 | 成人性做爰aaa片免费看不忠 | 欧美黑人巨大xxxxx | 亚洲欧美国产精品久久 | 国产97在线 | 亚洲 | 麻豆精产国品 | 亚洲色大成网站www | 樱花草在线社区www | 亚洲欧美精品aaaaaa片 | 亚洲精品综合五月久久小说 | √天堂中文官网8在线 | 粉嫩少妇内射浓精videos | 国内老熟妇对白xxxxhd | 国产精品无码久久av | 亚洲区小说区激情区图片区 | 色爱情人网站 | 免费网站看v片在线18禁无码 | 成 人 网 站国产免费观看 | 人妻少妇精品无码专区二区 | 国产亚洲精品久久久久久久久动漫 | 亚洲国产欧美日韩精品一区二区三区 | 最新国产麻豆aⅴ精品无码 | 国产另类ts人妖一区二区 | 色爱情人网站 | 亚洲精品中文字幕乱码 | 成人免费视频视频在线观看 免费 | а√天堂www在线天堂小说 | 国产精品欧美成人 | 精品无码成人片一区二区98 | 亚洲日本va午夜在线电影 | 漂亮人妻洗澡被公强 日日躁 | 妺妺窝人体色www婷婷 | 国产成人精品必看 | aⅴ亚洲 日韩 色 图网站 播放 | 水蜜桃av无码 | 国产无av码在线观看 | 99久久久无码国产精品免费 | 中文字幕亚洲情99在线 | 国产极品美女高潮无套在线观看 | 色婷婷av一区二区三区之红樱桃 | 波多野结衣高清一区二区三区 | 国产无遮挡又黄又爽免费视频 | 高潮毛片无遮挡高清免费视频 | 少妇的肉体aa片免费 | a片免费视频在线观看 | 性啪啪chinese东北女人 | 欧美老妇交乱视频在线观看 | 永久免费观看国产裸体美女 | 亚洲熟妇色xxxxx欧美老妇y | 中文亚洲成a人片在线观看 | 国产舌乚八伦偷品w中 | 大乳丰满人妻中文字幕日本 | 少妇愉情理伦片bd | 丝袜足控一区二区三区 | 2019nv天堂香蕉在线观看 | 真人与拘做受免费视频一 | 玩弄中年熟妇正在播放 | 丰满人妻一区二区三区免费视频 | 无码乱肉视频免费大全合集 | 麻花豆传媒剧国产免费mv在线 | 亚洲综合无码一区二区三区 | 四虎影视成人永久免费观看视频 | 国产精品久久久久久亚洲影视内衣 | 亚洲日韩av一区二区三区四区 | 亚洲区小说区激情区图片区 | 欧美成人高清在线播放 | 成人无码影片精品久久久 | 久久久久久国产精品无码下载 | 蜜桃视频韩日免费播放 | 婷婷色婷婷开心五月四房播播 | 久久国产精品偷任你爽任你 | 成人影院yy111111在线观看 | 中文字幕+乱码+中文字幕一区 | 搡女人真爽免费视频大全 | 国产乱人无码伦av在线a | 亚洲日韩精品欧美一区二区 | 国产精品成人av在线观看 | 久久久精品国产sm最大网站 | 亚洲精品久久久久中文第一幕 | 18精品久久久无码午夜福利 | 国精产品一区二区三区 | 亚洲爆乳精品无码一区二区三区 | 无码国产激情在线观看 | 亚洲国产精品久久久天堂 | 午夜时刻免费入口 | 四虎国产精品一区二区 | 黑人粗大猛烈进出高潮视频 | 熟妇人妻无乱码中文字幕 | 麻豆精产国品 | 亚洲综合伊人久久大杳蕉 | 在线 国产 欧美 亚洲 天堂 | 波多野42部无码喷潮在线 | 亚欧洲精品在线视频免费观看 | 精品国产一区二区三区av 性色 | 老头边吃奶边弄进去呻吟 | 日韩亚洲欧美中文高清在线 | 亚洲精品国产a久久久久久 | 成人综合网亚洲伊人 | 久久无码中文字幕免费影院蜜桃 | 欧美丰满少妇xxxx性 | 玩弄少妇高潮ⅹxxxyw | 4hu四虎永久在线观看 | 装睡被陌生人摸出水好爽 | 亚洲中文字幕av在天堂 | 亚洲中文字幕无码中字 | 国产成人一区二区三区别 | 无套内谢的新婚少妇国语播放 | 久久亚洲精品中文字幕无男同 | 国产精品久久久久久无码 | 欧美zoozzooz性欧美 | 国产精品自产拍在线观看 | 日本熟妇大屁股人妻 | 久久久久久久女国产乱让韩 | 一本大道伊人av久久综合 | 一本色道婷婷久久欧美 | 成人精品视频一区二区三区尤物 | 色综合久久中文娱乐网 | 蜜桃视频插满18在线观看 | 亚洲色偷偷偷综合网 | 啦啦啦www在线观看免费视频 | 色噜噜亚洲男人的天堂 | 亚洲成av人在线观看网址 | 丰满少妇女裸体bbw | 99国产精品白浆在线观看免费 | 少妇厨房愉情理9仑片视频 | 少妇久久久久久人妻无码 | аⅴ资源天堂资源库在线 | 人妻有码中文字幕在线 | 亚洲精品一区二区三区四区五区 | 国产一区二区三区四区五区加勒比 | 中文字幕乱妇无码av在线 | 国产内射老熟女aaaa | 伊人久久大香线蕉午夜 | 国产精品久久久久影院嫩草 | 精品久久8x国产免费观看 | 国产激情综合五月久久 | 高潮喷水的毛片 | 久久久成人毛片无码 | 对白脏话肉麻粗话av | 日韩精品无码免费一区二区三区 | 18无码粉嫩小泬无套在线观看 | 午夜成人1000部免费视频 | 色欲av亚洲一区无码少妇 | 亚洲综合无码一区二区三区 | 久久久久久久久蜜桃 | 亚洲综合无码久久精品综合 | 欧美丰满熟妇xxxx性ppx人交 | 亚洲 日韩 欧美 成人 在线观看 | 久久国产精品精品国产色婷婷 | 男人扒开女人内裤强吻桶进去 | 精品熟女少妇av免费观看 | 一区二区传媒有限公司 | 天天躁夜夜躁狠狠是什么心态 | 国产精品内射视频免费 | 国产内射爽爽大片视频社区在线 | 国产精品久久久午夜夜伦鲁鲁 | 久久久久亚洲精品中文字幕 | a片在线免费观看 | 久久久久成人精品免费播放动漫 | 久久国语露脸国产精品电影 | 欧美性色19p | 成人性做爰aaa片免费看不忠 | 国产综合色产在线精品 | 97夜夜澡人人爽人人喊中国片 | av在线亚洲欧洲日产一区二区 | 国内综合精品午夜久久资源 | 四虎永久在线精品免费网址 | 久久99国产综合精品 | 国产精品久久久久无码av色戒 | 99国产精品白浆在线观看免费 | 久久精品国产亚洲精品 | 亚洲熟妇色xxxxx欧美老妇y | 一区二区三区高清视频一 | 综合网日日天干夜夜久久 | a片在线免费观看 | 欧美成人午夜精品久久久 | 亚洲中文无码av永久不收费 | 亚洲国产精品毛片av不卡在线 | 午夜无码区在线观看 | 2019午夜福利不卡片在线 | 亚洲精品一区二区三区四区五区 | 亚洲爆乳无码专区 | 亚洲精品国产精品乱码不卡 | 欧美高清在线精品一区 | 精品久久久久香蕉网 | 蜜桃视频韩日免费播放 | 日韩欧美成人免费观看 | 国产在线精品一区二区三区直播 | 国产成人无码av片在线观看不卡 | 免费观看的无遮挡av | 久久精品国产亚洲精品 | 99久久久无码国产aaa精品 | 午夜丰满少妇性开放视频 | 久久久久亚洲精品中文字幕 | 无码人妻av免费一区二区三区 | 婷婷五月综合缴情在线视频 | 色婷婷久久一区二区三区麻豆 | 久久精品人人做人人综合 | 一本色道久久综合狠狠躁 | 初尝人妻少妇中文字幕 | 国产成人av免费观看 | 久久97精品久久久久久久不卡 | 国产精品沙发午睡系列 | 精品久久综合1区2区3区激情 | 国产精品美女久久久网av | 天堂а√在线地址中文在线 | 成人精品天堂一区二区三区 | 精品久久久无码人妻字幂 | 日本一本二本三区免费 | 亚洲国精产品一二二线 | 国产精品久久久久9999小说 | 亚洲欧美日韩国产精品一区二区 | 亚洲日韩av一区二区三区中文 | 中文字幕av无码一区二区三区电影 | 国产国语老龄妇女a片 | 国产成人综合在线女婷五月99播放 | 天天躁日日躁狠狠躁免费麻豆 | 午夜福利一区二区三区在线观看 | 日韩av无码中文无码电影 | 国产精品人人爽人人做我的可爱 | 男女猛烈xx00免费视频试看 | 免费观看的无遮挡av | 久久97精品久久久久久久不卡 | 99久久99久久免费精品蜜桃 | 搡女人真爽免费视频大全 | 99久久精品日本一区二区免费 | 鲁一鲁av2019在线 | 亚洲精品一区二区三区在线 | yw尤物av无码国产在线观看 | 理论片87福利理论电影 | 久久无码中文字幕免费影院蜜桃 | 日本熟妇大屁股人妻 | 国产成人人人97超碰超爽8 | 俺去俺来也在线www色官网 | 久久久中文久久久无码 | 5858s亚洲色大成网站www | 奇米影视888欧美在线观看 | a在线亚洲男人的天堂 | 国产亚洲精品久久久久久久 | 久久精品99久久香蕉国产色戒 | 久久久av男人的天堂 | 亚洲日韩乱码中文无码蜜桃臀网站 | 国产超碰人人爽人人做人人添 | 婷婷五月综合激情中文字幕 | 久久亚洲中文字幕精品一区 | 少妇被粗大的猛进出69影院 | 国产成人精品无码播放 | 日日天干夜夜狠狠爱 | 麻豆国产97在线 | 欧洲 | 日本大乳高潮视频在线观看 | 亚洲精品久久久久中文第一幕 | 巨爆乳无码视频在线观看 | 亚洲另类伦春色综合小说 | 无码精品人妻一区二区三区av | 欧美野外疯狂做受xxxx高潮 | 99久久久无码国产aaa精品 | 99久久精品国产一区二区蜜芽 | 国产绳艺sm调教室论坛 | 一本色道婷婷久久欧美 | 少妇被粗大的猛进出69影院 | 亚洲理论电影在线观看 | 曰本女人与公拘交酡免费视频 | 久久久久亚洲精品中文字幕 | 亚洲国产欧美在线成人 | 国产香蕉尹人综合在线观看 | 永久黄网站色视频免费直播 | 97夜夜澡人人双人人人喊 | 妺妺窝人体色www在线小说 | 全黄性性激高免费视频 | 熟妇人妻激情偷爽文 | 久久人人97超碰a片精品 | 麻豆成人精品国产免费 | 性做久久久久久久免费看 | 爆乳一区二区三区无码 | 天堂亚洲免费视频 | 国产成人无码a区在线观看视频app | 欧美肥老太牲交大战 | 久久亚洲精品中文字幕无男同 | 日日碰狠狠丁香久燥 | 亚洲区欧美区综合区自拍区 | 欧美熟妇另类久久久久久不卡 | 亚洲经典千人经典日产 | 极品嫩模高潮叫床 | 久久午夜无码鲁丝片 | 久久国产精品_国产精品 | 无码免费一区二区三区 | 精品国产乱码久久久久乱码 | 亚洲中文无码av永久不收费 | 精品无码av一区二区三区 | 中文字幕av无码一区二区三区电影 | 丰满人妻被黑人猛烈进入 | 欧美人与善在线com | 国产精品鲁鲁鲁 | 国产热a欧美热a在线视频 | 麻豆av传媒蜜桃天美传媒 | 日本饥渴人妻欲求不满 | 欧美丰满少妇xxxx性 | 天天爽夜夜爽夜夜爽 | 香蕉久久久久久av成人 | 久久精品中文字幕一区 | 亚洲人成人无码网www国产 | 国产欧美精品一区二区三区 | 双乳奶水饱满少妇呻吟 | 成 人影片 免费观看 | 欧美日韩色另类综合 | 久久久婷婷五月亚洲97号色 | 粉嫩少妇内射浓精videos | 无码成人精品区在线观看 | 亚洲中文字幕无码一久久区 | 亚洲人成网站在线播放942 | 99视频精品全部免费免费观看 | 牲欲强的熟妇农村老妇女 | av无码久久久久不卡免费网站 | 国产亚洲美女精品久久久2020 | 久久精品成人欧美大片 | 欧美日本免费一区二区三区 | 日本精品高清一区二区 | 国产人妻大战黑人第1集 | 无码精品国产va在线观看dvd | 色五月五月丁香亚洲综合网 | 国模大胆一区二区三区 | 欧美一区二区三区视频在线观看 | 97精品人妻一区二区三区香蕉 | 97色伦图片97综合影院 | 久久久精品欧美一区二区免费 | 久久99精品久久久久婷婷 | 日韩欧美成人免费观看 | 亚洲s色大片在线观看 | 丰满肥臀大屁股熟妇激情视频 | 天堂а√在线地址中文在线 | 国产成人精品优优av | 中文字幕乱码中文乱码51精品 | 午夜精品久久久久久久 | 图片区 小说区 区 亚洲五月 | 娇妻被黑人粗大高潮白浆 | 日韩亚洲欧美精品综合 | 蜜桃av抽搐高潮一区二区 | 76少妇精品导航 | 天下第一社区视频www日本 | 久久熟妇人妻午夜寂寞影院 | 少妇邻居内射在线 | 狠狠亚洲超碰狼人久久 | 亚洲一区二区三区无码久久 | 精品久久综合1区2区3区激情 | 国产性生交xxxxx无码 | 97精品国产97久久久久久免费 | 国内少妇偷人精品视频免费 | 人妻少妇精品无码专区动漫 | 国产sm调教视频在线观看 | 亚洲国产精品美女久久久久 | 无码av岛国片在线播放 | 午夜福利电影 | 99精品视频在线观看免费 | 亚洲国产精品无码久久久久高潮 | 国产亚洲精品久久久久久久 | 亚洲无人区一区二区三区 | 国产性生交xxxxx无码 | 88国产精品欧美一区二区三区 | 亚洲爆乳精品无码一区二区三区 | 国产亲子乱弄免费视频 | 精品人妻中文字幕有码在线 | 在线观看国产一区二区三区 | 精品久久久久久亚洲精品 | 国产精品久久久久久久9999 | 精品国产福利一区二区 | 久久久国产一区二区三区 | 未满成年国产在线观看 | 亚洲精品美女久久久久久久 | www国产亚洲精品久久网站 | 亚洲色大成网站www国产 | 成 人 网 站国产免费观看 | 久久人人爽人人爽人人片av高清 | 久久人人爽人人人人片 | 久久亚洲精品中文字幕无男同 | 国产色xx群视频射精 | 日本肉体xxxx裸交 | 熟女少妇在线视频播放 | 国产精品无码永久免费888 | 东京热男人av天堂 | 国产一区二区三区精品视频 | 人人澡人人透人人爽 | 国产内射爽爽大片视频社区在线 | 四虎影视成人永久免费观看视频 | 男人和女人高潮免费网站 | 牲欲强的熟妇农村老妇女 | 亚洲人亚洲人成电影网站色 | 久久久久久久久888 | 午夜福利一区二区三区在线观看 | 97夜夜澡人人双人人人喊 | 成人精品天堂一区二区三区 | 欧美亚洲国产一区二区三区 | 亚洲国产精华液网站w | 成熟女人特级毛片www免费 | 免费无码午夜福利片69 | 久久久久久九九精品久 | 蜜桃av蜜臀av色欲av麻 999久久久国产精品消防器材 | 亚洲欧洲中文日韩av乱码 | 国产欧美熟妇另类久久久 | 国产午夜精品一区二区三区嫩草 | 亚洲日韩乱码中文无码蜜桃臀网站 | 国产电影无码午夜在线播放 | 波多野结衣一区二区三区av免费 | 亚洲日本va午夜在线电影 | 亚洲乱码日产精品bd | 国产精品久久精品三级 | 国产小呦泬泬99精品 | 亚洲爆乳无码专区 | 天堂无码人妻精品一区二区三区 | 无码精品人妻一区二区三区av | 偷窥日本少妇撒尿chinese | 蜜桃无码一区二区三区 | 亚洲 另类 在线 欧美 制服 | 久久久久久九九精品久 | 精品国产国产综合精品 | 中国女人内谢69xxxx | 成人性做爰aaa片免费看不忠 | 日本大乳高潮视频在线观看 | 日韩精品无码一区二区中文字幕 | 亚洲日韩中文字幕在线播放 | 真人与拘做受免费视频一 | 精品亚洲韩国一区二区三区 | 97精品人妻一区二区三区香蕉 | 偷窥日本少妇撒尿chinese | 香蕉久久久久久av成人 | 小sao货水好多真紧h无码视频 | 好屌草这里只有精品 | 亚洲中文字幕无码一久久区 | 欧美乱妇无乱码大黄a片 | 成人aaa片一区国产精品 | 在线播放免费人成毛片乱码 | 色综合久久久久综合一本到桃花网 | 亚洲a无码综合a国产av中文 | 亚洲热妇无码av在线播放 | 麻豆人妻少妇精品无码专区 | 亚洲人成影院在线观看 | 日韩人妻无码一区二区三区久久99 | 国产在线一区二区三区四区五区 | 国内精品久久久久久中文字幕 | 欧美人与禽zoz0性伦交 | 国产成人综合色在线观看网站 | 色婷婷综合激情综在线播放 | 日本www一道久久久免费榴莲 | 精品人妻中文字幕有码在线 | 国内揄拍国内精品少妇国语 | 乌克兰少妇xxxx做受 | 国产人妻大战黑人第1集 | 人妻少妇精品无码专区动漫 | 18禁黄网站男男禁片免费观看 | 国产精品内射视频免费 | 精品熟女少妇av免费观看 | 国产精品毛多多水多 | www一区二区www免费 | 国产熟妇另类久久久久 | 亚洲人成网站色7799 | 成人毛片一区二区 | 亚洲熟女一区二区三区 | 四虎永久在线精品免费网址 | 亚洲 激情 小说 另类 欧美 | 成人精品天堂一区二区三区 | 亚洲人成无码网www | 麻豆国产人妻欲求不满谁演的 | 国产高潮视频在线观看 | 成 人影片 免费观看 | 久久国产精品萌白酱免费 | 亚洲色在线无码国产精品不卡 | 性做久久久久久久久 | 大肉大捧一进一出好爽视频 | 日韩亚洲欧美精品综合 | 在线天堂新版最新版在线8 | 99久久久无码国产精品免费 | 熟妇激情内射com | 亚洲熟妇自偷自拍另类 | 九月婷婷人人澡人人添人人爽 | 麻豆果冻传媒2021精品传媒一区下载 | 亚洲日韩中文字幕在线播放 | 女人被男人躁得好爽免费视频 | 最近免费中文字幕中文高清百度 | 国产精品99爱免费视频 | 人人妻在人人 | 久久精品中文闷骚内射 | 欧美日韩一区二区免费视频 | 国产在线aaa片一区二区99 | 色 综合 欧美 亚洲 国产 | 夜夜躁日日躁狠狠久久av | 青青草原综合久久大伊人精品 | 久久久久久国产精品无码下载 | 无码免费一区二区三区 | 无码人妻少妇伦在线电影 | 天天综合网天天综合色 | 亚洲国产成人av在线观看 | 亚洲人成影院在线观看 | 麻豆国产人妻欲求不满谁演的 | 77777熟女视频在线观看 а天堂中文在线官网 | 亚洲中文字幕在线无码一区二区 | 人妻熟女一区 | 人人妻人人澡人人爽人人精品浪潮 | 国产精品久久久久久无码 | 精品国产麻豆免费人成网站 | 内射白嫩少妇超碰 | 午夜无码人妻av大片色欲 | 国产乱人伦偷精品视频 | 国产舌乚八伦偷品w中 | 亚洲精品国产精品乱码不卡 | 久久久精品国产sm最大网站 | 99riav国产精品视频 | 99国产欧美久久久精品 | 美女扒开屁股让男人桶 | 在线观看欧美一区二区三区 | 亚洲日韩av一区二区三区四区 | 丰满少妇人妻久久久久久 | 一区二区三区乱码在线 | 欧洲 | 色婷婷av一区二区三区之红樱桃 | 成人女人看片免费视频放人 | 亚洲热妇无码av在线播放 | 中文字幕乱码人妻无码久久 | 国产欧美熟妇另类久久久 | 久久人人爽人人爽人人片ⅴ | 一二三四社区在线中文视频 | 性生交大片免费看l | 国产av一区二区三区最新精品 | 国产片av国语在线观看 | 狠狠综合久久久久综合网 | 欧美老人巨大xxxx做受 | 国产成人精品无码播放 | 扒开双腿疯狂进出爽爽爽视频 | 999久久久国产精品消防器材 | 99久久99久久免费精品蜜桃 | 伊人久久大香线蕉亚洲 | 久久久成人毛片无码 | 99er热精品视频 | 狠狠色噜噜狠狠狠7777奇米 | 熟妇激情内射com | 动漫av网站免费观看 | 中文字幕乱妇无码av在线 | 亚洲 日韩 欧美 成人 在线观看 | 天天综合网天天综合色 | 亚洲一区二区三区播放 | 野狼第一精品社区 | 成人免费视频视频在线观看 免费 | 久青草影院在线观看国产 | 欧美性色19p | 无码人妻少妇伦在线电影 | 无码精品人妻一区二区三区av | 久久久久人妻一区精品色欧美 | 亚洲精品欧美二区三区中文字幕 | 亚洲欧洲日本综合aⅴ在线 | 日韩少妇白浆无码系列 | 精品国产乱码久久久久乱码 | 国产人妖乱国产精品人妖 | 国产av人人夜夜澡人人爽麻豆 | 国产在线无码精品电影网 | 亚洲日韩一区二区三区 | 人妻互换免费中文字幕 | 六月丁香婷婷色狠狠久久 | 国产午夜福利100集发布 | 高潮毛片无遮挡高清免费 | 76少妇精品导航 | 一本一道久久综合久久 | 亚洲精品一区二区三区在线观看 | 色窝窝无码一区二区三区色欲 | 国产性生大片免费观看性 | 男女超爽视频免费播放 | 国产精品成人av在线观看 | 特级做a爰片毛片免费69 | 中文字幕人妻无码一区二区三区 | 精品国产av色一区二区深夜久久 | 鲁大师影院在线观看 | 两性色午夜视频免费播放 | 亚洲欧洲无卡二区视頻 | 天天综合网天天综合色 | 成人免费无码大片a毛片 | 精品人妻人人做人人爽 | 精品国产成人一区二区三区 | 又大又紧又粉嫩18p少妇 | 国内丰满熟女出轨videos | 十八禁视频网站在线观看 | 欧美激情内射喷水高潮 | 无码精品国产va在线观看dvd | 国产亚av手机在线观看 | 性色av无码免费一区二区三区 | 亚洲精品国产第一综合99久久 | 东京热一精品无码av | 亚洲乱亚洲乱妇50p | 综合激情五月综合激情五月激情1 | 狠狠亚洲超碰狼人久久 | 欧美日韩一区二区三区自拍 | 国产精品久久久午夜夜伦鲁鲁 | av人摸人人人澡人人超碰下载 | 免费乱码人妻系列无码专区 | 日本精品少妇一区二区三区 | 人人妻人人澡人人爽人人精品浪潮 | 欧美黑人乱大交 | aⅴ亚洲 日韩 色 图网站 播放 | 又色又爽又黄的美女裸体网站 | 亚洲熟悉妇女xxx妇女av | 欧美性生交活xxxxxdddd | 两性色午夜视频免费播放 | 麻豆蜜桃av蜜臀av色欲av | 日韩精品a片一区二区三区妖精 | 国产特级毛片aaaaaaa高清 | 国产一区二区不卡老阿姨 | 国产色精品久久人妻 | 无码一区二区三区在线 | 精品久久久久久亚洲精品 | 好男人www社区 | 国产莉萝无码av在线播放 | 老子影院午夜伦不卡 | 真人与拘做受免费视频 | 国产精品亚洲五月天高清 | 国产绳艺sm调教室论坛 | 国产精品亚洲专区无码不卡 | 久久人人爽人人爽人人片ⅴ | 日本一区二区三区免费播放 | 亚洲欧洲日本综合aⅴ在线 | 亚洲欧美日韩综合久久久 | 国产成人av免费观看 | 亚洲中文无码av永久不收费 | 色婷婷综合激情综在线播放 | 18禁止看的免费污网站 | 无码一区二区三区在线 | 人人澡人人透人人爽 | 国产精品久久久久久久9999 | 亚洲色大成网站www | 免费看男女做好爽好硬视频 | 国产一精品一av一免费 | 久激情内射婷内射蜜桃人妖 | 国产成人无码av在线影院 | 熟妇人妻无码xxx视频 | 国产成人无码av一区二区 | 天天拍夜夜添久久精品大 | 免费播放一区二区三区 | 中文精品无码中文字幕无码专区 | 无码纯肉视频在线观看 | 欧美人与善在线com | 377p欧洲日本亚洲大胆 | 欧美性生交活xxxxxdddd | 香港三级日本三级妇三级 | 少妇性荡欲午夜性开放视频剧场 | 欧美 丝袜 自拍 制服 另类 | 狠狠色欧美亚洲狠狠色www | 久久久久免费看成人影片 | 东北女人啪啪对白 | 台湾无码一区二区 | 亚洲中文字幕无码一久久区 | 高清国产亚洲精品自在久久 | 国产真实夫妇视频 | 国产亚洲精品久久久久久大师 | 又大又黄又粗又爽的免费视频 | 377p欧洲日本亚洲大胆 | 少妇太爽了在线观看 | 丰满妇女强制高潮18xxxx | 日本xxxx色视频在线观看免费 | 国产精品国产三级国产专播 | 亚洲欧洲无卡二区视頻 | 国产免费观看黄av片 | 国产精品久久精品三级 | 露脸叫床粗话东北少妇 | 国产乱子伦视频在线播放 | 人妻无码αv中文字幕久久琪琪布 | 国产精品成人av在线观看 | 国语精品一区二区三区 | 久久久无码中文字幕久... | 欧美熟妇另类久久久久久不卡 | 97se亚洲精品一区 | 爽爽影院免费观看 | 精品久久久无码人妻字幂 | 天天躁日日躁狠狠躁免费麻豆 | 99视频精品全部免费免费观看 | 亚洲精品中文字幕久久久久 | 亚洲一区二区三区国产精华液 | 性色欲情网站iwww九文堂 | 亚洲小说春色综合另类 | yw尤物av无码国产在线观看 | 久久97精品久久久久久久不卡 | 少妇高潮喷潮久久久影院 | 精品国产青草久久久久福利 | 风流少妇按摩来高潮 | 亚洲精品综合五月久久小说 | 任你躁在线精品免费 | 欧美肥老太牲交大战 | 欧美大屁股xxxxhd黑色 | 久久精品女人的天堂av | 无遮无挡爽爽免费视频 | 亚洲aⅴ无码成人网站国产app | 日本精品久久久久中文字幕 | 国产九九九九九九九a片 | 国产内射老熟女aaaa | 欧美成人高清在线播放 | 久久久久se色偷偷亚洲精品av | 中文字幕无码av激情不卡 | 国产人成高清在线视频99最全资源 | 久久人妻内射无码一区三区 | 国产精品久久久久久无码 | 亚洲中文字幕在线观看 | 日本一区二区三区免费播放 | 精品久久久无码中文字幕 | 麻豆蜜桃av蜜臀av色欲av | 国产精品国产三级国产专播 | 真人与拘做受免费视频 | 精品人妻人人做人人爽 | 性欧美牲交在线视频 | 呦交小u女精品视频 | 成人精品视频一区二区三区尤物 | 久久久久久久人妻无码中文字幕爆 | 亚洲人成人无码网www国产 | 国产精品美女久久久 | 国产精品爱久久久久久久 | 久久亚洲精品中文字幕无男同 | 亚洲乱码中文字幕在线 | 国产性生大片免费观看性 | 黑人巨大精品欧美黑寡妇 | 亚洲性无码av中文字幕 | 欧美兽交xxxx×视频 | 性生交片免费无码看人 | 免费无码午夜福利片69 | 无码任你躁久久久久久久 | 国色天香社区在线视频 | 欧美色就是色 | 日欧一片内射va在线影院 | 婷婷色婷婷开心五月四房播播 | 国产av无码专区亚洲a∨毛片 | 国产xxx69麻豆国语对白 | 一个人看的视频www在线 | 激情综合激情五月俺也去 | 国产亚洲美女精品久久久2020 | 精品成人av一区二区三区 | 国产女主播喷水视频在线观看 | 国产成人午夜福利在线播放 | 国内揄拍国内精品少妇国语 | 午夜肉伦伦影院 | 精品国产国产综合精品 | 中文字幕无线码 | 蜜桃无码一区二区三区 | 在线a亚洲视频播放在线观看 | 国产莉萝无码av在线播放 | 台湾无码一区二区 | 99re在线播放 | 荫蒂添的好舒服视频囗交 | 国产无遮挡又黄又爽免费视频 | 日韩人妻少妇一区二区三区 | 欧美日韩视频无码一区二区三 | 久久亚洲精品中文字幕无男同 | 水蜜桃亚洲一二三四在线 | 国产激情艳情在线看视频 | 久久久久免费看成人影片 | 一本大道伊人av久久综合 | 亚欧洲精品在线视频免费观看 | 久久人妻内射无码一区三区 | 成人免费无码大片a毛片 | 男女下面进入的视频免费午夜 | 中国大陆精品视频xxxx | 国产综合在线观看 | 国产精品久久国产精品99 | 国产精品人人妻人人爽 | 激情亚洲一区国产精品 | 一个人免费观看的www视频 | 内射巨臀欧美在线视频 | 国产精品99爱免费视频 | 一个人看的视频www在线 | 一本精品99久久精品77 | 日本大香伊一区二区三区 | 俺去俺来也在线www色官网 | 激情人妻另类人妻伦 | 中文字幕av日韩精品一区二区 | 野狼第一精品社区 | yw尤物av无码国产在线观看 | 中文字幕精品av一区二区五区 | 午夜丰满少妇性开放视频 | 亚洲国产精品无码一区二区三区 | 玩弄中年熟妇正在播放 | 樱花草在线社区www | 人妻无码久久精品人妻 | 野狼第一精品社区 | 成人片黄网站色大片免费观看 | 99re在线播放 | 日本肉体xxxx裸交 | 亚洲男女内射在线播放 | 国产97色在线 | 免 | 一本大道伊人av久久综合 | 人妻天天爽夜夜爽一区二区 | 成人免费视频在线观看 | 国产精品亚洲五月天高清 | 国产精品va在线观看无码 | 永久免费观看美女裸体的网站 | 国产乱人伦av在线无码 | 狠狠色色综合网站 | 粉嫩少妇内射浓精videos | 青春草在线视频免费观看 | 撕开奶罩揉吮奶头视频 | 久久久久久久久蜜桃 | 东京热无码av男人的天堂 | 麻豆国产人妻欲求不满 | 天天爽夜夜爽夜夜爽 | 99国产精品白浆在线观看免费 | 久久国产精品萌白酱免费 | 无码帝国www无码专区色综合 | 无码av免费一区二区三区试看 | 天天燥日日燥 | 免费视频欧美无人区码 | 欧美黑人性暴力猛交喷水 | 国产人妻久久精品二区三区老狼 | 亚洲自偷精品视频自拍 | 中国女人内谢69xxxxxa片 | 夫妻免费无码v看片 | 日韩无码专区 | 国产精品亚洲lv粉色 | 亚洲日韩av一区二区三区四区 | 国内揄拍国内精品人妻 | 狠狠色丁香久久婷婷综合五月 | 亚洲一区二区三区偷拍女厕 | 欧美黑人乱大交 | 2020久久超碰国产精品最新 | 午夜精品久久久久久久 | 欧美变态另类xxxx | 青春草在线视频免费观看 | 丁香啪啪综合成人亚洲 | 欧美人与物videos另类 | 麻豆果冻传媒2021精品传媒一区下载 | 欧美日本精品一区二区三区 | 亚洲中文字幕无码中文字在线 | 野外少妇愉情中文字幕 | 激情五月综合色婷婷一区二区 | 欧美熟妇另类久久久久久多毛 | 国产明星裸体无码xxxx视频 | 综合人妻久久一区二区精品 | 波多野结衣乳巨码无在线观看 | 精品少妇爆乳无码av无码专区 | 色婷婷综合激情综在线播放 | 偷窥日本少妇撒尿chinese | 在线精品国产一区二区三区 | 十八禁真人啪啪免费网站 | 一二三四社区在线中文视频 | 在教室伦流澡到高潮hnp视频 | 日韩视频 中文字幕 视频一区 | 亚洲欧美中文字幕5发布 | 国产国语老龄妇女a片 | 性欧美疯狂xxxxbbbb | 亚洲七七久久桃花影院 | 国产精品.xx视频.xxtv | 亚洲欧洲日本无在线码 | 强开小婷嫩苞又嫩又紧视频 | 又色又爽又黄的美女裸体网站 | 亚洲欧美精品伊人久久 | 亚洲理论电影在线观看 | 中文久久乱码一区二区 | 国产成人久久精品流白浆 | 一区二区三区高清视频一 | 国产麻豆精品一区二区三区v视界 | 51国偷自产一区二区三区 | www国产亚洲精品久久网站 | 亚洲男人av天堂午夜在 | 粗大的内捧猛烈进出视频 | 老子影院午夜精品无码 | 日本免费一区二区三区最新 | 国产乱人无码伦av在线a | 亚洲码国产精品高潮在线 | 午夜无码人妻av大片色欲 | 亚洲乱码中文字幕在线 | 又紧又大又爽精品一区二区 | 熟女体下毛毛黑森林 | 亚洲国精产品一二二线 | 国产九九九九九九九a片 | 国产香蕉尹人综合在线观看 | 妺妺窝人体色www婷婷 | 国产乱人无码伦av在线a | 免费国产黄网站在线观看 | 丰满人妻一区二区三区免费视频 | 精品夜夜澡人妻无码av蜜桃 | 久久精品人妻少妇一区二区三区 | 对白脏话肉麻粗话av | 国产人妻大战黑人第1集 | 福利一区二区三区视频在线观看 | 少妇人妻大乳在线视频 | 99精品久久毛片a片 | 欧美熟妇另类久久久久久不卡 | 久久国产精品精品国产色婷婷 | 曰韩无码二三区中文字幕 | 色窝窝无码一区二区三区色欲 | 人人妻人人澡人人爽欧美精品 | 双乳奶水饱满少妇呻吟 | 图片区 小说区 区 亚洲五月 | 俺去俺来也在线www色官网 | 无码免费一区二区三区 | 亚洲中文无码av永久不收费 | 一本久道久久综合狠狠爱 | 亚洲精品午夜无码电影网 | 人人妻人人藻人人爽欧美一区 | 亚洲精品成人福利网站 | 久久综合网欧美色妞网 | 国产猛烈高潮尖叫视频免费 | 国产亲子乱弄免费视频 | 久9re热视频这里只有精品 | 欧美刺激性大交 | 67194成是人免费无码 | 国产极品美女高潮无套在线观看 | 国产人妻大战黑人第1集 | 乱人伦人妻中文字幕无码久久网 | 无码人妻精品一区二区三区不卡 | 亚洲国产欧美日韩精品一区二区三区 | 国产成人人人97超碰超爽8 | 玩弄人妻少妇500系列视频 | 亚洲区欧美区综合区自拍区 | 国产日产欧产精品精品app | 男女性色大片免费网站 | 波多野结衣一区二区三区av免费 | 国产成人一区二区三区在线观看 | 国产成人精品优优av | 国产手机在线αⅴ片无码观看 | 欧美阿v高清资源不卡在线播放 | 免费国产成人高清在线观看网站 | 国产精品久久国产精品99 | 午夜精品久久久内射近拍高清 | 国内少妇偷人精品视频免费 | 久久亚洲中文字幕精品一区 | 在线亚洲高清揄拍自拍一品区 | 精品无人国产偷自产在线 | 免费国产成人高清在线观看网站 | 国产精品亚洲lv粉色 | 久久精品99久久香蕉国产色戒 | 日本精品少妇一区二区三区 | 久久99久久99精品中文字幕 | 色欲久久久天天天综合网精品 | 亚洲人成影院在线无码按摩店 | 3d动漫精品啪啪一区二区中 | 一个人看的www免费视频在线观看 | 午夜熟女插插xx免费视频 | 亚洲精品中文字幕乱码 | 国产无av码在线观看 | 无码任你躁久久久久久久 | 丰满肥臀大屁股熟妇激情视频 | 极品尤物被啪到呻吟喷水 | 国产成人一区二区三区在线观看 | 亚洲精品一区国产 | 天堂无码人妻精品一区二区三区 | 亚洲欧美色中文字幕在线 | 丰满人妻一区二区三区免费视频 | 精品无码av一区二区三区 | 国内精品一区二区三区不卡 | 一本久久a久久精品vr综合 | 久久综合激激的五月天 | 久久久成人毛片无码 | 午夜男女很黄的视频 | 国产性生大片免费观看性 | 国内精品久久毛片一区二区 | 一本一道久久综合久久 | 最近免费中文字幕中文高清百度 | 亚洲a无码综合a国产av中文 | 久久久久成人片免费观看蜜芽 | 国产精品久久久av久久久 | 又大又硬又黄的免费视频 | 久久人人爽人人爽人人片ⅴ | 亚洲狠狠婷婷综合久久 | 久久精品国产一区二区三区肥胖 | 国产精品久久久久9999小说 | 伊人久久婷婷五月综合97色 | 国产另类ts人妖一区二区 | 麻豆国产人妻欲求不满谁演的 | 男女猛烈xx00免费视频试看 | 国内揄拍国内精品人妻 | 色一情一乱一伦一区二区三欧美 | 亚洲综合无码久久精品综合 | 亚洲 另类 在线 欧美 制服 | 国产激情无码一区二区app | 国产色精品久久人妻 | 一本久久a久久精品亚洲 | ass日本丰满熟妇pics | 久久zyz资源站无码中文动漫 | 国产一区二区三区影院 | 精品欧洲av无码一区二区三区 | 欧美精品无码一区二区三区 | 99久久精品午夜一区二区 | 少妇性荡欲午夜性开放视频剧场 | 中文字幕无码日韩欧毛 | 又粗又大又硬又长又爽 | 欧美黑人乱大交 | 久久无码人妻影院 | 日韩欧美中文字幕公布 | 人人澡人人妻人人爽人人蜜桃 | 精品亚洲成av人在线观看 | 久久久精品成人免费观看 | 熟妇女人妻丰满少妇中文字幕 | 亚洲成av人片天堂网无码】 | 亚洲欧洲日本综合aⅴ在线 | 久久亚洲中文字幕无码 | 亚洲男人av天堂午夜在 | 国产成人精品视频ⅴa片软件竹菊 | 亚洲自偷精品视频自拍 | 国产精品亚洲а∨无码播放麻豆 | 国产女主播喷水视频在线观看 | 国产真实乱对白精彩久久 | 精品无码一区二区三区爱欲 | 免费视频欧美无人区码 | 亚洲国产精品毛片av不卡在线 | 无码人妻丰满熟妇区毛片18 | 日日橹狠狠爱欧美视频 | 国产亚洲欧美日韩亚洲中文色 | 精品国产精品久久一区免费式 | 国产内射爽爽大片视频社区在线 | 国产精品视频免费播放 | 亚洲国产欧美在线成人 | 正在播放老肥熟妇露脸 | 天堂а√在线地址中文在线 | 免费人成网站视频在线观看 | 国产区女主播在线观看 | 免费无码av一区二区 | 国产精品美女久久久久av爽李琼 | 精品aⅴ一区二区三区 | 色综合久久久无码网中文 | 欧美性生交活xxxxxdddd | 99久久久国产精品无码免费 | 国产av无码专区亚洲a∨毛片 | 久久久久亚洲精品中文字幕 | 人人妻在人人 | 国色天香社区在线视频 | 任你躁国产自任一区二区三区 | 99久久精品日本一区二区免费 | 麻豆国产丝袜白领秘书在线观看 | 永久黄网站色视频免费直播 | 人人妻人人澡人人爽欧美精品 | 久久久久久久女国产乱让韩 | 蜜臀aⅴ国产精品久久久国产老师 | 亚洲成a人片在线观看无码 | 国产精品第一区揄拍无码 | 人人妻在人人 | 久久久久久久久蜜桃 | 日本熟妇乱子伦xxxx | 天堂在线观看www | 377p欧洲日本亚洲大胆 | 又大又硬又黄的免费视频 | 亚洲熟悉妇女xxx妇女av | 国产无套内射久久久国产 | 精品无人区无码乱码毛片国产 | 任你躁在线精品免费 | 少妇邻居内射在线 | 麻豆国产97在线 | 欧洲 | 蜜臀aⅴ国产精品久久久国产老师 | 波多野42部无码喷潮在线 | 成人精品一区二区三区中文字幕 | 精品成人av一区二区三区 | 日韩av无码一区二区三区 | 欧洲精品码一区二区三区免费看 | 鲁大师影院在线观看 | 精品偷自拍另类在线观看 | 国产精品99爱免费视频 | 粗大的内捧猛烈进出视频 | √天堂资源地址中文在线 | 在线天堂新版最新版在线8 | 乱人伦人妻中文字幕无码久久网 | 国产超碰人人爽人人做人人添 | 日日噜噜噜噜夜夜爽亚洲精品 | 久激情内射婷内射蜜桃人妖 | 亚洲一区av无码专区在线观看 | 国产成人亚洲综合无码 | 亚洲爆乳精品无码一区二区三区 | 日韩成人一区二区三区在线观看 | 又大又黄又粗又爽的免费视频 | 国产超级va在线观看视频 | 欧美一区二区三区视频在线观看 | 久久精品无码一区二区三区 | 亚洲狠狠婷婷综合久久 | 婷婷六月久久综合丁香 | 免费看男女做好爽好硬视频 | 国产亚洲视频中文字幕97精品 | 日本丰满熟妇videos | 成人无码视频在线观看网站 | 久久99精品国产麻豆 | 国产suv精品一区二区五 | 性欧美videos高清精品 | 好男人www社区 | 国产乱人伦偷精品视频 | 久久国产精品偷任你爽任你 | 亚洲午夜久久久影院 | 国产精品亚洲五月天高清 | 国产在线精品一区二区高清不卡 | 国产精品美女久久久久av爽李琼 | 丰满诱人的人妻3 | 国产无套内射久久久国产 | 一本色道久久综合亚洲精品不卡 | 亚洲热妇无码av在线播放 | 欧洲熟妇精品视频 | 中国女人内谢69xxxxxa片 | 亚洲日韩精品欧美一区二区 | 国产日产欧产精品精品app | 成年女人永久免费看片 | 99精品无人区乱码1区2区3区 | 无码人妻精品一区二区三区不卡 | 在线成人www免费观看视频 | 99精品视频在线观看免费 | 无码av岛国片在线播放 | 欧美丰满熟妇xxxx性ppx人交 | 亚洲精品久久久久中文第一幕 | 亚洲综合无码久久精品综合 | 亚洲精品国产品国语在线观看 | 亚洲の无码国产の无码步美 | 人妻少妇精品无码专区动漫 | 亚洲国产精品美女久久久久 | 亚洲熟女一区二区三区 | 99久久婷婷国产综合精品青草免费 | 国产偷抇久久精品a片69 | 久久精品国产99精品亚洲 | 亚洲成a人片在线观看无码3d | 欧美日韩久久久精品a片 | 日本一卡2卡3卡4卡无卡免费网站 国产一区二区三区影院 | 中国大陆精品视频xxxx | av香港经典三级级 在线 | 三上悠亚人妻中文字幕在线 | 日韩成人一区二区三区在线观看 | 久久99精品久久久久婷婷 | 亚洲人成网站在线播放942 | 377p欧洲日本亚洲大胆 | 性欧美熟妇videofreesex | 高清国产亚洲精品自在久久 | 东京热无码av男人的天堂 | 欧美一区二区三区 | 性色欲网站人妻丰满中文久久不卡 | 欧美性生交xxxxx久久久 | 欧美高清在线精品一区 | 国产艳妇av在线观看果冻传媒 | 中文字幕无码日韩欧毛 | 全黄性性激高免费视频 | 成人aaa片一区国产精品 | 久久天天躁夜夜躁狠狠 | 粉嫩少妇内射浓精videos | 夜先锋av资源网站 | 18禁止看的免费污网站 | 玩弄少妇高潮ⅹxxxyw | 无码国产色欲xxxxx视频 | 色欲综合久久中文字幕网 | 亚洲人亚洲人成电影网站色 | 乱中年女人伦av三区 | 国产性生大片免费观看性 | 在线播放免费人成毛片乱码 | 国产免费久久久久久无码 | 性色av无码免费一区二区三区 | 性生交大片免费看l | 久久久成人毛片无码 | 亚洲一区av无码专区在线观看 | 亚洲熟妇色xxxxx欧美老妇y | 免费国产黄网站在线观看 | 欧美人与牲动交xxxx | 美女极度色诱视频国产 | 国产成人无码区免费内射一片色欲 | 丝袜足控一区二区三区 | 欧美丰满熟妇xxxx | 亚洲成av人片在线观看无码不卡 | 久久久久se色偷偷亚洲精品av | 日本大香伊一区二区三区 | 大胆欧美熟妇xx | 国产亚洲精品久久久久久国模美 | 久久久精品成人免费观看 | 亚洲爆乳无码专区 | 秋霞成人午夜鲁丝一区二区三区 | 精品国产aⅴ无码一区二区 | 亚洲国产高清在线观看视频 | 中文字幕 人妻熟女 | 久久综合香蕉国产蜜臀av | 亚洲日韩一区二区三区 | 又黄又爽又色的视频 | 男女猛烈xx00免费视频试看 | 国产精品美女久久久 | 久久精品女人的天堂av | 国产 浪潮av性色四虎 | 中文字幕色婷婷在线视频 | 中文字幕人妻丝袜二区 | 亚洲第一无码av无码专区 | 精品国产一区二区三区av 性色 | 麻豆国产97在线 | 欧洲 | 精品人妻人人做人人爽夜夜爽 | 国产精品a成v人在线播放 | 秋霞成人午夜鲁丝一区二区三区 | 色 综合 欧美 亚洲 国产 | 久久成人a毛片免费观看网站 | 久久国产劲爆∧v内射 | 国产成人综合在线女婷五月99播放 | 欧美阿v高清资源不卡在线播放 | 久久久精品欧美一区二区免费 | 美女扒开屁股让男人桶 | 日韩少妇白浆无码系列 | 亚洲综合无码一区二区三区 | 国产麻豆精品一区二区三区v视界 | 日韩精品a片一区二区三区妖精 | 少妇人妻偷人精品无码视频 | 成年美女黄网站色大免费视频 | 中文字幕乱码人妻二区三区 | 亚洲国产精品久久人人爱 | 国产亚洲日韩欧美另类第八页 | 国产亚洲精品精品国产亚洲综合 | 亚洲精品一区国产 | 久久精品国产日本波多野结衣 | 亚洲成a人片在线观看无码 | 国产香蕉尹人视频在线 | 丰满人妻精品国产99aⅴ | 国产色xx群视频射精 | 麻花豆传媒剧国产免费mv在线 | 精品国产一区av天美传媒 | 久久99精品国产麻豆蜜芽 | 久久精品中文字幕一区 | 国内揄拍国内精品少妇国语 | 人人爽人人澡人人人妻 | 狠狠色噜噜狠狠狠狠7777米奇 | 国产欧美精品一区二区三区 | 国产又粗又硬又大爽黄老大爷视 | 成人精品天堂一区二区三区 | 国内老熟妇对白xxxxhd | 强辱丰满人妻hd中文字幕 | 精品久久久久久亚洲精品 | 亚洲欧洲中文日韩av乱码 | аⅴ资源天堂资源库在线 | aa片在线观看视频在线播放 | 久久久久99精品国产片 | 日本精品久久久久中文字幕 | 夜夜躁日日躁狠狠久久av | 国产精品二区一区二区aⅴ污介绍 | 牲欲强的熟妇农村老妇女 | 鲁大师影院在线观看 | 国产无套粉嫩白浆在线 | 亚洲精品一区二区三区婷婷月 | 精品国产aⅴ无码一区二区 | 在线成人www免费观看视频 | 中文无码成人免费视频在线观看 | 久久无码中文字幕免费影院蜜桃 | 在线欧美精品一区二区三区 | 亚洲国精产品一二二线 | 精品厕所偷拍各类美女tp嘘嘘 | 欧美丰满老熟妇xxxxx性 | 在线播放无码字幕亚洲 | 99国产欧美久久久精品 | 国产精品自产拍在线观看 | 丝袜 中出 制服 人妻 美腿 | 天天av天天av天天透 | 精品一二三区久久aaa片 | 亚洲大尺度无码无码专区 | 久久精品中文闷骚内射 | 99久久久国产精品无码免费 | 国产精品无套呻吟在线 | 国产网红无码精品视频 | 青青青爽视频在线观看 | 国模大胆一区二区三区 | 欧美喷潮久久久xxxxx | 国产精品自产拍在线观看 | 丰满人妻被黑人猛烈进入 | 国产精品a成v人在线播放 | 美女毛片一区二区三区四区 | 国产高清av在线播放 | 国产熟女一区二区三区四区五区 | 日日躁夜夜躁狠狠躁 | 久久久久久久久蜜桃 | а√资源新版在线天堂 | 国产人妻人伦精品1国产丝袜 | 一区二区传媒有限公司 | 丰腴饱满的极品熟妇 | 18精品久久久无码午夜福利 | 国产精品a成v人在线播放 | 色婷婷av一区二区三区之红樱桃 | 少妇一晚三次一区二区三区 | av香港经典三级级 在线 | 性做久久久久久久免费看 | 亚洲国产精品成人久久蜜臀 | 国产免费久久久久久无码 | 少妇无码av无码专区在线观看 | 亚洲va中文字幕无码久久不卡 | 亚洲精品鲁一鲁一区二区三区 | 中文无码成人免费视频在线观看 | 精品久久久久久亚洲精品 | 成人影院yy111111在线观看 | 小sao货水好多真紧h无码视频 | 久久 国产 尿 小便 嘘嘘 | 欧美人妻一区二区三区 | www国产亚洲精品久久网站 | 色欲久久久天天天综合网精品 | 欧洲精品码一区二区三区免费看 | 中文字幕日韩精品一区二区三区 | 成人试看120秒体验区 | 国产无av码在线观看 | 香港三级日本三级妇三级 | 午夜无码人妻av大片色欲 | 欧美人与善在线com | 曰韩无码二三区中文字幕 | 国产精品久久久久久亚洲影视内衣 | 人妻aⅴ无码一区二区三区 | 少女韩国电视剧在线观看完整 | 妺妺窝人体色www婷婷 | 女人被爽到呻吟gif动态图视看 | 麻豆国产人妻欲求不满谁演的 | 丰满人妻一区二区三区免费视频 | 久久天天躁狠狠躁夜夜免费观看 | 人人爽人人澡人人人妻 | 欧美一区二区三区 | 少妇邻居内射在线 | 999久久久国产精品消防器材 | 97久久超碰中文字幕 | 在线播放无码字幕亚洲 | 又色又爽又黄的美女裸体网站 | 荫蒂被男人添的好舒服爽免费视频 | 人人妻人人藻人人爽欧美一区 | 图片小说视频一区二区 | 亚洲午夜久久久影院 | 在线播放免费人成毛片乱码 | 樱花草在线社区www | 久久久中文久久久无码 | 影音先锋中文字幕无码 | 18禁黄网站男男禁片免费观看 | 欧美日韩一区二区免费视频 | 荫蒂被男人添的好舒服爽免费视频 | 少妇性俱乐部纵欲狂欢电影 | 人妻夜夜爽天天爽三区 | 天天av天天av天天透 | 日本一区二区三区免费高清 | √天堂资源地址中文在线 | 少妇性俱乐部纵欲狂欢电影 | 少妇被黑人到高潮喷出白浆 | 无人区乱码一区二区三区 | 在线观看欧美一区二区三区 | 亚洲欧美国产精品专区久久 | 国产成人精品无码播放 | 免费国产成人高清在线观看网站 | 99国产精品白浆在线观看免费 | 人妻中文无码久热丝袜 | 99久久人妻精品免费一区 | 人人澡人人妻人人爽人人蜜桃 | 国产熟妇高潮叫床视频播放 | 久久aⅴ免费观看 | 日本熟妇浓毛 | 天天拍夜夜添久久精品大 | 3d动漫精品啪啪一区二区中 | 精品亚洲成av人在线观看 | 国産精品久久久久久久 | 亚洲精品一区二区三区在线观看 | 黑人大群体交免费视频 | 性开放的女人aaa片 | 鲁大师影院在线观看 | 国产人妻精品一区二区三区 | 久久午夜无码鲁丝片午夜精品 | 国产真实夫妇视频 | 亚洲综合色区中文字幕 | 初尝人妻少妇中文字幕 | 亚洲自偷自拍另类第1页 | 国产精品福利视频导航 | 亚洲日韩乱码中文无码蜜桃臀网站 | www一区二区www免费 | 51国偷自产一区二区三区 | 对白脏话肉麻粗话av | 亚洲一区二区三区偷拍女厕 | 男女超爽视频免费播放 | 草草网站影院白丝内射 | 国产精品.xx视频.xxtv | 国产又爽又猛又粗的视频a片 | 美女扒开屁股让男人桶 | 日本大香伊一区二区三区 | 丰满诱人的人妻3 | 免费看男女做好爽好硬视频 | 欧美成人家庭影院 | 亚洲精品一区二区三区大桥未久 | 99久久人妻精品免费二区 | 色五月五月丁香亚洲综合网 | 精品无码国产一区二区三区av | 最新国产乱人伦偷精品免费网站 | 精品国精品国产自在久国产87 | 久久精品国产一区二区三区 | a片在线免费观看 | 国产精品久久久久无码av色戒 | 欧美乱妇无乱码大黄a片 | 呦交小u女精品视频 | 日本又色又爽又黄的a片18禁 | yw尤物av无码国产在线观看 | 女高中生第一次破苞av | 大色综合色综合网站 | 国产在线精品一区二区高清不卡 | 国产精品无套呻吟在线 | 国产深夜福利视频在线 | 欧美 丝袜 自拍 制服 另类 |