859 lines
36 KiB
Python
859 lines
36 KiB
Python
import math
|
||
import time
|
||
import cv2
|
||
|
||
import sys
|
||
import os
|
||
|
||
sys.path.append(os.path.dirname(os.path.dirname(os.path.abspath(__file__))))
|
||
from utils.detect_track import detect_horizontal_track_edge
|
||
from base_move.turn_degree import turn_degree
|
||
|
||
def align_to_horizontal_line(ctrl, msg, observe=False, max_attempts=3):
|
||
"""
|
||
控制机器人旋转到横向线水平的位置
|
||
|
||
参数:
|
||
ctrl: Robot_Ctrl 对象，包含里程计信息
|
||
msg: robot_control_cmd_lcmt 对象，用于发送命令
|
||
observe: 是否输出中间状态信息和可视化结果，默认为False
|
||
max_attempts: 最大尝试次数，默认为3次
|
||
|
||
返回:
|
||
bool: 是否成功校准
|
||
"""
|
||
attempts = 0
|
||
aligned = False
|
||
image = ctrl.image_processor.get_current_image()
|
||
|
||
while attempts < max_attempts and not aligned:
|
||
print(f"尝试次数: {attempts+1}/{max_attempts}")
|
||
# 检测横向线边缘
|
||
edge_point, edge_info = detect_horizontal_track_edge(ctrl.image_processor.get_current_image(), observe=observe, delay=1000 if observe else 0)
|
||
|
||
if edge_point is None or edge_info is None:
|
||
print("未检测到横向线，无法进行校准")
|
||
return False
|
||
|
||
# 获取检测到的斜率和其他信息
|
||
slope = edge_info["slope"]
|
||
is_horizontal = edge_info["is_horizontal"]
|
||
|
||
if observe:
|
||
print(f"检测到横向线，斜率: {slope:.6f}")
|
||
print(f"是否足够水平: {is_horizontal}")
|
||
|
||
# 如果已经水平，则无需旋转
|
||
if is_horizontal:
|
||
print("横向线已经水平，无需校准")
|
||
return True
|
||
|
||
# 计算需要旋转的角度
|
||
# 斜率 = tan(θ)，因此 θ = arctan(斜率)
|
||
angle_rad = math.atan(slope)
|
||
angle_deg = math.degrees(angle_rad)
|
||
|
||
# 调整角度方向
|
||
# 正的斜率意味着线条从左到右上升，需要逆时针旋转校正
|
||
# 负的斜率意味着线条从左到右下降，需要顺时针旋转校正
|
||
# 注意旋转方向: 顺时针为负角度，逆时针为正角度
|
||
angle_to_rotate = -angle_deg # 取负值使旋转方向正确
|
||
|
||
if observe:
|
||
print(f"需要旋转的角度: {angle_to_rotate:.2f}度")
|
||
|
||
# 执行旋转
|
||
# 如果角度很小，增加一个小的偏移以确保旋转足够
|
||
if abs(angle_to_rotate) < 3.0:
|
||
angle_to_rotate *= 1.5 # 对小角度进行放大以确保效果
|
||
|
||
# 限制旋转角度，避免过度旋转
|
||
angle_to_rotate = max(-30, min(30, angle_to_rotate))
|
||
|
||
# 使用turn_degree函数执行旋转
|
||
turn_success = turn_degree(ctrl, msg, angle_to_rotate, absolute=False)
|
||
|
||
if observe:
|
||
print(f"旋转结果: {'成功' if turn_success else '失败'}")
|
||
|
||
# 增加尝试次数
|
||
attempts += 1
|
||
|
||
# 在旋转后重新获取图像，这里需要调用获取图像的函数
|
||
# 代码中没有提供获取实时图像的方法，假设每次外部会更新image参数
|
||
|
||
# 检查是否已经对齐
|
||
# 对于实际应用，应该在旋转后重新捕获图像并检测横向线
|
||
# 这里简单地根据旋转是否成功和旋转角度是否足够小来判断
|
||
if turn_success and abs(angle_to_rotate) < 5.0:
|
||
aligned = True
|
||
|
||
return aligned
|
||
|
||
def calculate_distance_to_line(edge_info, camera_height, camera_tilt_angle_deg=0, observe=False):
|
||
"""
|
||
根据相机参数和图像中横线位置计算相机到横线的实际距离
|
||
|
||
几何模型说明:
|
||
1. 相机位于高度camera_height处，向下倾斜camera_tilt_angle_deg度
|
||
2. 图像底部对应相机视场的下边缘，横线在图像中的位置通过像素坐标确定
|
||
3. 计算相机视线到横线的角度，然后使用三角函数计算实际距离
|
||
|
||
参数:
|
||
edge_info: 边缘信息字典，包含distance_to_bottom等信息
|
||
camera_height: 相机高度(米)
|
||
camera_tilt_angle_deg: 相机向下倾斜的角度（度）
|
||
observe: 是否打印中间计算值
|
||
|
||
返回:
|
||
float: 到横向线的X轴水平距离(米)
|
||
"""
|
||
if edge_info is None or "distance_to_bottom" not in edge_info:
|
||
return None
|
||
|
||
# 1. 获取图像中交点到底部的距离(像素)
|
||
distance_in_pixels = edge_info["distance_to_bottom"]
|
||
if observe:
|
||
print(f"图像中交点到底部的像素距离: {distance_in_pixels}")
|
||
|
||
# 2. 获取相机参数
|
||
horizontal_fov_rad = 1.46608 # 水平视场角（弧度）约84度
|
||
image_height_px = 1080 # 图像高度（像素）
|
||
image_width_px = 1920 # 图像宽度（像素）
|
||
|
||
# 3. 计算垂直视场角
|
||
aspect_ratio = image_width_px / image_height_px # 宽高比
|
||
vertical_fov_rad = horizontal_fov_rad / aspect_ratio # 垂直视场角(弧度)
|
||
vertical_fov_deg = math.degrees(vertical_fov_rad) # 垂直视场角(度)
|
||
|
||
if observe:
|
||
print(f"相机参数: 水平FOV={math.degrees(horizontal_fov_rad):.1f}°, 垂直FOV={vertical_fov_deg:.1f}°")
|
||
print(f"图像尺寸: {image_width_px}x{image_height_px}, 宽高比: {aspect_ratio:.2f}")
|
||
|
||
# 4. 直接计算视线角度
|
||
|
||
# 计算图像底部到相机视场中心的角度
|
||
half_vfov_rad = vertical_fov_rad / 2
|
||
|
||
# 计算图像底部到横线的角度比例
|
||
# 比例 = 底部到横线的像素距离 / 图像总高度
|
||
pixel_ratio = distance_in_pixels / image_height_px
|
||
|
||
# 计算从图像底部到横线的角度
|
||
bottom_to_line_angle_rad = pixel_ratio * vertical_fov_rad
|
||
|
||
# 计算从相机视场中心到横线的角度
|
||
# 负值表示横线在视场中心以下，正值表示在中心以上
|
||
center_to_line_angle_rad = bottom_to_line_angle_rad - half_vfov_rad
|
||
|
||
# 考虑相机倾斜角度
|
||
# 相机向下倾斜为正值，此时视场中心相对水平线向下
|
||
camera_tilt_rad = math.radians(camera_tilt_angle_deg)
|
||
|
||
# 计算横线相对于水平面的视线角度
|
||
# 负值表示视线向下看到横线，正值表示视线向上看到横线
|
||
view_angle_rad = center_to_line_angle_rad - camera_tilt_rad
|
||
|
||
if observe:
|
||
print(f"视场角度关系:")
|
||
print(f" - 图像底部到横线角度: {math.degrees(bottom_to_line_angle_rad):.2f}°")
|
||
print(f" - 视场中心到横线角度: {math.degrees(center_to_line_angle_rad):.2f}°")
|
||
print(f" - 相机倾斜角度: {camera_tilt_angle_deg}°")
|
||
print(f" - 最终视线角度: {math.degrees(view_angle_rad):.2f}° ({'向下' if view_angle_rad < 0 else '向上'})")
|
||
|
||
# 5. 防止除零错误或异常值
|
||
# 确保视线角度不接近于0(水平视线无法确定地面交点)
|
||
min_angle_rad = 0.01 # 约0.57度
|
||
if abs(view_angle_rad) < min_angle_rad:
|
||
if observe:
|
||
print(f"视线角度过小({math.degrees(view_angle_rad):.2f}°)，使用最小角度: {math.degrees(min_angle_rad):.2f}°")
|
||
view_angle_rad = -min_angle_rad # 设为向下的最小角度
|
||
|
||
# 6. 计算水平距离
|
||
# 仅当视线向下时计算地面距离
|
||
if view_angle_rad < 0: # 视线向下
|
||
# 基本几何关系: 水平距离 = 高度 / tan(视线向下的角度)
|
||
# 注意角度为负，所以需要取负
|
||
ground_distance = camera_height / math.tan(-view_angle_rad)
|
||
|
||
if observe:
|
||
print(f"计算公式: 距离 = 相机高度({camera_height}米) / tan(|视线角度|({abs(math.degrees(view_angle_rad)):.2f}°))")
|
||
print(f"计算结果: 距离 = {ground_distance:.3f}米")
|
||
else: # 视线平行或向上，无法确定地面交点
|
||
if observe:
|
||
print(f"视线向上或水平，无法计算地面距离")
|
||
return 0.5 # 返回一个默认值
|
||
|
||
# 7. 应用校正和限制
|
||
# 可选的校正因子(通过实验校准)
|
||
correction_factor = 1.0
|
||
distance = ground_distance * correction_factor
|
||
|
||
# 设置合理的范围限制
|
||
min_distance = 0.1 # 最小距离(米)
|
||
|
||
# 限制结果在合理范围内
|
||
final_distance = max(min_distance, distance)
|
||
|
||
if observe and final_distance != distance:
|
||
print(f"应用范围限制: 原始距离 {distance:.3f}米 -> 最终距离 {final_distance:.3f}米")
|
||
elif observe:
|
||
print(f"最终距离: {final_distance:.3f}米")
|
||
|
||
return final_distance
|
||
|
||
def move_to_hori_line(ctrl, msg, target_distance=0.5, observe=False):
|
||
"""
|
||
控制机器人校准并移动到横向线前的指定距离
|
||
|
||
参数:
|
||
ctrl: Robot_Ctrl 对象，包含里程计信息
|
||
msg: robot_control_cmd_lcmt 对象，用于发送命令
|
||
target_distance: 目标位置与横向线的距离(米)，默认为0.5米
|
||
observe: 是否输出中间状态信息和可视化结果，默认为False
|
||
|
||
返回:
|
||
bool: 是否成功到达目标位置
|
||
"""
|
||
# 首先校准到水平
|
||
print("校准到横向线水平")
|
||
aligned = align_to_horizontal_line(ctrl, msg, observe=observe)
|
||
|
||
if not aligned:
|
||
print("无法校准到横向线水平，停止移动")
|
||
return False
|
||
|
||
# 检测横向线
|
||
image = ctrl.image_processor.get_current_image()
|
||
edge_point, edge_info = detect_horizontal_track_edge(image, observe=observe)
|
||
|
||
if edge_point is None or edge_info is None:
|
||
print("无法检测到横向线，停止移动")
|
||
return False
|
||
|
||
# 获取相机高度
|
||
camera_height = 0.355 # 单位: 米 # INFO from TF-tree
|
||
|
||
# 计算当前距离
|
||
current_distance = calculate_distance_to_line(edge_info, camera_height, observe=observe)
|
||
|
||
if current_distance is None:
|
||
print("无法计算到横向线的距离，停止移动")
|
||
return False
|
||
|
||
if observe:
|
||
print(f"当前距离: {current_distance:.3f}米, 目标距离: {target_distance:.3f}米")
|
||
|
||
# 计算需要移动的距离
|
||
distance_to_move = current_distance - target_distance
|
||
|
||
if abs(distance_to_move) < 0.05: # 如果已经很接近目标距离
|
||
print("已经达到目标距离，无需移动")
|
||
return True
|
||
|
||
# 设置移动命令
|
||
msg.mode = 11 # Locomotion模式
|
||
msg.gait_id = 26 # 自变频步态
|
||
|
||
# 移动方向设置
|
||
forward = distance_to_move > 0 # 判断是前进还是后退
|
||
|
||
# 设置移动速度
|
||
# 根据需要移动的距离动态调整移动速度
|
||
if abs(distance_to_move) > 1.0:
|
||
move_speed = 1.0 # 距离较远时用最大速度
|
||
elif abs(distance_to_move) > 0.5:
|
||
move_speed = 0.6 # 中等距离用中等速度
|
||
elif abs(distance_to_move) > 0.2:
|
||
move_speed = 0.3 # 较近距离用较慢速度
|
||
else:
|
||
move_speed = 0.15 # 非常接近时用最慢速度
|
||
if forward:
|
||
msg.vel_des = [move_speed, 0, 0] # 设置前进速度
|
||
else:
|
||
msg.vel_des = [-move_speed, 0, 0] # 设置后退速度
|
||
|
||
# 获取起始位置
|
||
start_position = list(ctrl.odo_msg.xyz) # 转换为列表，因为xyz是元组
|
||
|
||
if observe:
|
||
print(f"起始位置: {start_position}")
|
||
print(f"开始移动: {'前进' if forward else '后退'} {abs(distance_to_move):.3f}米")
|
||
# 在起点放置绿色标记
|
||
if hasattr(ctrl, 'place_marker'):
|
||
ctrl.place_marker(start_position[0], start_position[1], start_position[2] if len(start_position) > 2 else 0.0, 'green', observe=True)
|
||
|
||
# 估算移动时间，但实际上会通过里程计控制
|
||
move_time = abs(distance_to_move) / move_speed
|
||
msg.duration = 0 # wait next cmd
|
||
msg.step_height = [0.06, 0.06] # 抬腿高度
|
||
msg.life_count += 1
|
||
|
||
# 发送命令
|
||
ctrl.Send_cmd(msg)
|
||
|
||
# 使用里程计进行实时监控移动距离
|
||
distance_moved = 0
|
||
start_time = time.time()
|
||
timeout = move_time + 2 # 增加超时时间为预计移动时间加2秒
|
||
|
||
# 监控移动距离，增加移动距离的系数
|
||
while distance_moved < abs(distance_to_move) * 1.05 and time.time() - start_time < timeout:
|
||
# 计算已移动距离
|
||
current_position = ctrl.odo_msg.xyz
|
||
dx = current_position[0] - start_position[0]
|
||
dy = current_position[1] - start_position[1]
|
||
distance_moved = math.sqrt(dx*dx + dy*dy)
|
||
|
||
if observe and time.time() % 0.5 < 0.02: # 每0.5秒左右打印一次
|
||
print(f"已移动: {distance_moved:.3f}米, 目标: {abs(distance_to_move):.3f}米")
|
||
|
||
time.sleep(0.05) # 小间隔检查位置
|
||
|
||
# 使用平滑停止
|
||
# ctrl.base_msg.stop_smooth()
|
||
ctrl.base_msg.stop()
|
||
|
||
if observe:
|
||
print(f"移动完成，平稳停止，通过里程计计算的移动距离: {distance_moved:.3f}米")
|
||
# 在终点放置红色标记
|
||
end_position = list(ctrl.odo_msg.xyz)
|
||
if hasattr(ctrl, 'place_marker'):
|
||
ctrl.place_marker(end_position[0], end_position[1], end_position[2] if len(end_position) > 2 else 0.0, 'red', observe=True)
|
||
|
||
# 放宽判断标准
|
||
return abs(distance_moved - abs(distance_to_move)) < 0.15 # 如果误差小于15厘米，则认为成功
|
||
|
||
def arc_turn_around_hori_line(ctrl, msg, angle_deg=90, left=True, target_distance=0.2,
|
||
pass_align=False,
|
||
scan_qrcode=False, qr_check_interval=0.3,
|
||
observe=False,
|
||
bad_big_angle_corret=False):
|
||
"""
|
||
对准前方横线，然后以计算出来的距离为半径，做一个向左或向右的圆弧旋转。
|
||
参数:
|
||
ctrl: Robot_Ctrl对象
|
||
msg: robot_control_cmd_lcmt对象
|
||
target_distance: 横线前目标保持距离，默认0.5米
|
||
angle_deg: 旋转角度，支持90或180度
|
||
left: True为左转，False为右转
|
||
pass_align: 是否跳过对准步骤，默认为False
|
||
observe: 是否打印调试信息
|
||
scan_qrcode: 是否在旋转过程中扫描QR码，默认为False
|
||
qr_check_interval: QR码检查间隔时间(秒)，默认为0.3秒
|
||
|
||
bad_big_angle_corret: 是否对大角度错误进行修正，默认为False # TODO clear
|
||
返回:
|
||
bool或元组: 如果scan_qrcode为False，返回bool表示是否成功完成操作；
|
||
如果scan_qrcode为True，返回(bool, str)元组，表示(是否成功完成操作, QR码扫描结果)
|
||
"""
|
||
# 启动异步QR码扫描（如果需要）
|
||
qr_result = None
|
||
if scan_qrcode:
|
||
# 确保image_processor存在
|
||
if not hasattr(ctrl, 'image_processor') or ctrl.image_processor is None:
|
||
print("警告: 无法启用QR码扫描，image_processor不存在")
|
||
scan_qrcode = False
|
||
else:
|
||
# 启动异步扫描
|
||
try:
|
||
ctrl.image_processor.start_async_scan(interval=0.2)
|
||
print("已启动异步QR码扫描")
|
||
except Exception as e:
|
||
print(f"启动QR码扫描失败: {e}")
|
||
scan_qrcode = False
|
||
|
||
# 1. 对准横线
|
||
if not pass_align:
|
||
print("校准到横向线水平")
|
||
aligned = align_to_horizontal_line(ctrl, msg, observe=observe)
|
||
if not aligned:
|
||
print("无法校准到横向线水平，停止动作")
|
||
if scan_qrcode:
|
||
ctrl.image_processor.stop_async_scan()
|
||
return False, None
|
||
return False
|
||
else:
|
||
print("跳过对准步骤")
|
||
|
||
# 2. 检测横线并计算距离
|
||
image = ctrl.image_processor.get_current_image()
|
||
edge_point, edge_info = detect_horizontal_track_edge(image, observe=observe)
|
||
if edge_point is None or edge_info is None:
|
||
print("无法检测到横向线，停止动作")
|
||
if scan_qrcode:
|
||
ctrl.image_processor.stop_async_scan()
|
||
return False, qr_result
|
||
return False
|
||
|
||
camera_height = 0.355 # 单位: 米
|
||
r = calculate_distance_to_line(edge_info, camera_height, observe=observe)
|
||
|
||
# 减去目标距离
|
||
r -= target_distance
|
||
|
||
if r is None:
|
||
print("无法计算到横向线的距离，停止动作")
|
||
if scan_qrcode:
|
||
ctrl.image_processor.stop_async_scan()
|
||
return False, qr_result
|
||
return False
|
||
|
||
if observe:
|
||
print(f"当前距离: {r:.3f}米")
|
||
|
||
# 3. 计算圆弧运动参数
|
||
# 根据角度大小调整时间补偿系数
|
||
if angle_deg <= 70:
|
||
time_compensation = 0.78 # 对70度及以下角度使用更小的补偿系数
|
||
elif angle_deg <= 90:
|
||
time_compensation = 0.85 # 90度左右使用稍大的系数
|
||
else:
|
||
# 对180度旋转，增加补偿系数，使旋转时间更长
|
||
time_compensation = 1.4 # 增加40%的时间
|
||
|
||
# 调整实际目标角度，针对不同角度应用不同的缩放比例
|
||
actual_angle_deg = angle_deg
|
||
if angle_deg <= 70:
|
||
actual_angle_deg = angle_deg * 0.85 # 70度及以下角度减少到85%
|
||
elif angle_deg <= 90:
|
||
actual_angle_deg = angle_deg * 0.9 # 90度减少到90%
|
||
elif angle_deg >= 170: # 处理大角度旋转，特别是180度情况
|
||
actual_angle_deg = angle_deg # 对于180度旋转，使用完整角度，不进行缩放
|
||
|
||
if observe and actual_angle_deg != angle_deg:
|
||
print(f"应用角度补偿: 目标{angle_deg}度 -> 实际指令{actual_angle_deg:.1f}度")
|
||
|
||
angle_rad = math.radians(actual_angle_deg)
|
||
|
||
# 设定角速度（rad/s），可根据实际调整
|
||
# 减小基础角速度，增加精度
|
||
# if angle_deg >= 170:
|
||
# base_w = 0.5 # 大角度旋转时使用更小的基础角速度，提高精度
|
||
# else:
|
||
base_w = 0.8 # 原来是0.6
|
||
|
||
w = base_w if left else -base_w # 左转为正，右转为负
|
||
|
||
# 确保r有一个最小值，避免速度过低
|
||
min_radius = 0.2 # 最小半径，单位：米
|
||
effective_r = max(r, min_radius)
|
||
|
||
# 线速度，必须与角速度和半径的关系一致：v = w * r
|
||
# 这样才能保证走圆弧
|
||
v = abs(w * effective_r) # 线速度，正负号与角速度方向一致
|
||
|
||
if observe:
|
||
if effective_r != r:
|
||
print(f"注意: 实际半径 {r:.3f}米 太小，使用最小半径 {effective_r:.3f}米 计算速度")
|
||
|
||
t = abs(angle_rad / w) # 运动时间，取绝对值保证为正
|
||
|
||
# 应用时间补偿系数
|
||
t *= time_compensation
|
||
|
||
if observe:
|
||
print(f"圆弧半径: {effective_r:.3f}米, 角速度: {w:.3f}rad/s, 线速度: {v:.3f}m/s")
|
||
print(f"理论运动时间: {abs(angle_rad / w):.2f}s, 应用补偿系数{time_compensation}后: {t:.2f}s")
|
||
|
||
# 4. 发送圆弧运动指令
|
||
msg.mode = 11
|
||
msg.gait_id = 26
|
||
msg.vel_des = [v, 0, w] # [前进速度, 侧向速度, 角速度]
|
||
msg.duration = 0 # wait next
|
||
msg.step_height = [0.06, 0.06]
|
||
msg.life_count += 1
|
||
|
||
ctrl.Send_cmd(msg)
|
||
|
||
# 记录起始角度和位置
|
||
start_yaw = ctrl.odo_msg.rpy[2]
|
||
start_position = list(ctrl.odo_msg.xyz)
|
||
|
||
# 在起点放置蓝色标记
|
||
# TEST 获取当前姿态
|
||
yaw = ctrl.odo_msg.rpy[2]
|
||
ctrl.place_marker(start_position[0], start_position[1],
|
||
start_position[2] if len(start_position) > 2 else 0.0,
|
||
'blue', observe=True)
|
||
print(f"在起点位置放置蓝色标记: {start_position}")
|
||
|
||
# 计算圆心位置 - 相对于机器人前进方向的左侧或右侧
|
||
# 圆心与机器人当前位置的距离就是圆弧半径
|
||
center_x = start_position[0] - effective_r * math.sin(yaw) * (1 if left else -1)
|
||
center_y = start_position[1] + effective_r * math.cos(yaw) * (1 if left else -1)
|
||
center_z = start_position[2] if len(start_position) > 2 else 0.0
|
||
|
||
# 在圆心放置绿色标记
|
||
if observe and hasattr(ctrl, 'place_marker'):
|
||
ctrl.place_marker(center_x, center_y, center_z, 'green', observe=True)
|
||
print(f"在旋转圆心放置绿色标记: [{center_x:.3f}, {center_y:.3f}, {center_z:.3f}]")
|
||
|
||
# 计算目标点的位置 - 旋转后的位置
|
||
# 使用几何关系：目标点是从起点绕圆心旋转angle_rad弧度后的位置
|
||
target_angle = yaw + (angle_rad if left else -angle_rad)
|
||
# 计算从圆心到目标点的向量
|
||
target_x = center_x + effective_r * math.sin(target_angle) * (1 if left else -1)
|
||
target_y = center_y - effective_r * math.cos(target_angle) * (1 if left else -1)
|
||
target_z = start_position[2] if len(start_position) > 2 else 0.0
|
||
|
||
# 在目标点放置红色标记
|
||
if observe and hasattr(ctrl, 'place_marker'):
|
||
ctrl.place_marker(target_x, target_y, target_z, 'red', observe=True)
|
||
print(f"在目标位置放置红色标记: [{target_x:.3f}, {target_y:.3f}, {target_z:.3f}]")
|
||
|
||
# 计算圆弧长度
|
||
arc_length = abs(angle_rad * effective_r)
|
||
|
||
# 进入循环，实时判断
|
||
distance_moved = 0
|
||
angle_turned = 0
|
||
start_time = time.time()
|
||
timeout = t + 2 # 给个超时保护
|
||
last_qr_check_time = 0
|
||
|
||
# 初始速度
|
||
current_v = v
|
||
current_w = w
|
||
|
||
# 对小角度提前开始减速
|
||
if angle_deg <= 70:
|
||
slowdown_threshold = 0.65 # 当达到目标角度的65%时开始减速
|
||
elif angle_deg >= 170: # 针对大角度(如180度)采用更晚的减速时机
|
||
slowdown_threshold = 0.7 # 当达到目标角度的70%时才开始减速
|
||
else:
|
||
slowdown_threshold = 0.75 # 对于大角度，75%时开始减速（原来是80%）
|
||
|
||
# 根据角度大小调整旋转停止阈值
|
||
if angle_deg >= 170:
|
||
rotation_completion_threshold = 0.85 # 大角度旋转提前停止，避免惯性导致过度旋转
|
||
else:
|
||
rotation_completion_threshold = 0.95 # 默认旋转到95%时停止
|
||
|
||
# 监控旋转进度并自行控制减速
|
||
# 对于180度旋转使用特殊处理
|
||
if angle_deg >= 170:
|
||
# 简化旋转逻辑，直接持续旋转一段时间
|
||
rotation_time = t # 使用计算的理论时间
|
||
if observe:
|
||
print(f"大角度旋转({angle_deg}度)，将持续旋转 {rotation_time:.2f} 秒")
|
||
|
||
# 发送固定时间的旋转命令
|
||
start_time = time.time()
|
||
|
||
# 前段使用全速，但减少前段旋转时间比例，减少过度旋转
|
||
first_segment_time = rotation_time * 0.65 # 原来是0.7，减少到0.65
|
||
elapsed_time = 0
|
||
|
||
while elapsed_time < first_segment_time and time.time() - start_time < rotation_time + 2: # 加2秒保护
|
||
elapsed_time = time.time() - start_time
|
||
|
||
# 计算已旋转角度（仅用于打印）
|
||
if observe and time.time() % 0.5 < 0.02:
|
||
current_yaw = ctrl.odo_msg.rpy[2]
|
||
angle_turned = current_yaw - start_yaw
|
||
# 角度归一化处理
|
||
while angle_turned > math.pi:
|
||
angle_turned -= 2 * math.pi
|
||
while angle_turned < -math.pi:
|
||
angle_turned += 2 * math.pi
|
||
print(f"全速阶段 - 已旋转: {math.degrees(angle_turned):.2f}度, 已用时间: {elapsed_time:.2f}s/{rotation_time:.2f}s")
|
||
|
||
# 检查QR码扫描结果（如果启用）
|
||
if scan_qrcode:
|
||
current_time = time.time()
|
||
if current_time - last_qr_check_time >= qr_check_interval:
|
||
qr_data, scan_time = ctrl.image_processor.get_last_qr_result()
|
||
if qr_data and scan_time > start_time:
|
||
qr_result = qr_data
|
||
print(f"🔍 在旋转过程中扫描到QR码: {qr_data}")
|
||
last_qr_check_time = current_time
|
||
|
||
time.sleep(0.05)
|
||
|
||
# 后段减速，减少减速幅度，保持更高速度
|
||
if observe:
|
||
print(f"进入减速阶段")
|
||
|
||
# 减速到70%，保持前进速度和角速度的比例关系
|
||
reduced_w = w * 0.7 # 原来是0.5
|
||
reduced_v = abs(reduced_w * effective_r) # 维持圆弧关系
|
||
|
||
msg.mode = 11
|
||
msg.gait_id = 26
|
||
msg.vel_des = [reduced_v, 0, reduced_w] # 维持圆弧关系
|
||
msg.duration = 0
|
||
msg.step_height = [0.06, 0.06]
|
||
msg.life_count += 1
|
||
ctrl.Send_cmd(msg)
|
||
|
||
# 继续旋转直到总时间结束，减少总旋转时间
|
||
total_rotation_time = rotation_time * 1.15 # 原来是1.2，减少到1.15
|
||
while time.time() - start_time < total_rotation_time:
|
||
# 计算已旋转角度（仅用于打印）
|
||
if observe and time.time() % 0.5 < 0.02:
|
||
current_yaw = ctrl.odo_msg.rpy[2]
|
||
angle_turned = current_yaw - start_yaw
|
||
# 角度归一化处理
|
||
while angle_turned > math.pi:
|
||
angle_turned -= 2 * math.pi
|
||
while angle_turned < -math.pi:
|
||
angle_turned += 2 * math.pi
|
||
print(f"减速阶段 - 已旋转: {math.degrees(angle_turned):.2f}度, 已用时间: {time.time() - start_time:.2f}s/{total_rotation_time:.2f}s")
|
||
|
||
# 检查QR码扫描结果（如果启用）
|
||
if scan_qrcode:
|
||
current_time = time.time()
|
||
if current_time - last_qr_check_time >= qr_check_interval:
|
||
qr_data, scan_time = ctrl.image_processor.get_last_qr_result()
|
||
if qr_data and scan_time > start_time:
|
||
qr_result = qr_data
|
||
print(f"🔍 在旋转过程中扫描到QR码: {qr_data}")
|
||
last_qr_check_time = current_time
|
||
|
||
time.sleep(0.05)
|
||
|
||
# 停止
|
||
ctrl.base_msg.stop()
|
||
else:
|
||
# 非180度的常规旋转逻辑保持不变
|
||
while abs(angle_turned) < abs(angle_rad) * rotation_completion_threshold and time.time() - start_time < timeout:
|
||
# 计算已移动距离
|
||
current_position = ctrl.odo_msg.xyz
|
||
dx = current_position[0] - start_position[0]
|
||
dy = current_position[1] - start_position[1]
|
||
distance_moved = math.sqrt(dx*dx + dy*dy)
|
||
|
||
# 计算已旋转角度
|
||
current_yaw = ctrl.odo_msg.rpy[2]
|
||
angle_turned = current_yaw - start_yaw
|
||
# 角度归一化处理
|
||
while angle_turned > math.pi:
|
||
angle_turned -= 2 * math.pi
|
||
while angle_turned < -math.pi:
|
||
angle_turned += 2 * math.pi
|
||
|
||
# 检查QR码扫描结果（如果启用）
|
||
if scan_qrcode:
|
||
current_time = time.time()
|
||
if current_time - last_qr_check_time >= qr_check_interval:
|
||
qr_data, scan_time = ctrl.image_processor.get_last_qr_result()
|
||
if qr_data and scan_time > start_time:
|
||
qr_result = qr_data
|
||
print(f"🔍 在旋转过程中扫描到QR码: {qr_data}")
|
||
last_qr_check_time = current_time
|
||
|
||
# 计算角度完成比例
|
||
completion_ratio = abs(angle_turned) / abs(angle_rad)
|
||
|
||
# 当完成一定比例时开始减速
|
||
if completion_ratio > slowdown_threshold:
|
||
# 实现多阶段减速，特别针对大角度旋转
|
||
if angle_deg >= 170:
|
||
if completion_ratio < 0.7: # 第一阶段：60%-70%
|
||
speed_factor = 0.65 # 原来是0.7，减少到0.65
|
||
elif completion_ratio < 0.8: # 第二阶段：70%-80%
|
||
speed_factor = 0.35 # 原来是0.4，减少到0.35
|
||
else: # 第三阶段：80%以上
|
||
speed_factor = 0.15 # 原来是0.2，减少到0.15
|
||
else:
|
||
# 标准减速逻辑保持不变
|
||
# 剩余比例作为速度系数
|
||
speed_factor = 1.0 - (completion_ratio - slowdown_threshold) / (1.0 - slowdown_threshold)
|
||
|
||
# 根据角度调整最小速度因子
|
||
min_speed_factor = 0.15 # 默认最小速度因子
|
||
# 确保速度系数不会太小，但可以更慢一些
|
||
speed_factor = max(min_speed_factor, speed_factor)
|
||
|
||
# 应用减速
|
||
current_v = v * speed_factor
|
||
current_w = w * speed_factor
|
||
|
||
# 更新速度命令
|
||
msg.mode = 11
|
||
msg.gait_id = 26
|
||
msg.vel_des = [current_v, 0, current_w]
|
||
msg.duration = 200 # 短时间更新
|
||
msg.step_height = [0.06, 0.06]
|
||
msg.life_count += 1
|
||
ctrl.Send_cmd(msg)
|
||
|
||
# 针对180度旋转，在减速命令后稍作等待，以确保减速效果
|
||
if angle_deg >= 170:
|
||
time.sleep(0.03) # 短暂延迟，让减速命令更有效
|
||
|
||
if observe and time.time() % 0.5 < 0.02: # 每0.5秒左右打印一次
|
||
print(f"减速阶段 - 已旋转: {math.degrees(angle_turned):.2f}度, 速度因子: {speed_factor:.2f}, 当前角速度: {current_w:.3f}rad/s")
|
||
else:
|
||
if observe and time.time() % 0.5 < 0.02: # 每0.5秒左右打印一次
|
||
print(f"已移动: {distance_moved:.3f}米, 已旋转: {math.degrees(angle_turned):.2f}度")
|
||
|
||
time.sleep(0.05)
|
||
|
||
if observe:
|
||
print(f"旋转过程结束，当前速度: [{current_v:.3f}, 0, {current_w:.3f}]")
|
||
print(f'{math.degrees(ctrl.odo_msg.rpy[2]):.2f}')
|
||
ctrl.base_msg.stop()
|
||
if observe:
|
||
print(f'{math.degrees(ctrl.odo_msg.rpy[2]):.2f}')
|
||
|
||
# 记录停止前的角度
|
||
pre_stop_yaw = ctrl.odo_msg.rpy[2]
|
||
|
||
# 等待机器人完全停止
|
||
time.sleep(0.5)
|
||
|
||
# 停下来后的最终角度
|
||
final_yaw = ctrl.odo_msg.rpy[2]
|
||
final_angle_turned = final_yaw - start_yaw
|
||
|
||
# 记录实际结束位置
|
||
end_position = list(ctrl.odo_msg.xyz)
|
||
|
||
# 在实际结束位置放置黄色标记，与预计的目标位置进行对比
|
||
if observe and hasattr(ctrl, 'place_marker'):
|
||
ctrl.place_marker(end_position[0], end_position[1],
|
||
end_position[2] if len(end_position) > 2 else 0.0,
|
||
'yellow', observe=True)
|
||
print(f"在实际结束位置放置黄色标记: {end_position}")
|
||
|
||
# 计算实际位置与预期位置的偏差
|
||
position_error = math.sqrt((end_position[0] - target_x)**2 +
|
||
(end_position[1] - target_y)**2)
|
||
print(f"位置误差: {position_error:.3f}米")
|
||
|
||
# 角度归一化处理
|
||
while final_angle_turned > math.pi:
|
||
final_angle_turned -= 2 * math.pi
|
||
while final_angle_turned < -math.pi:
|
||
final_angle_turned += 2 * math.pi
|
||
|
||
final_angle_deg = math.degrees(final_angle_turned)
|
||
# 这里使用原始目标角度进行比较
|
||
if not bad_big_angle_corret:
|
||
target_angle_deg = angle_deg if left else -angle_deg
|
||
else:
|
||
target_angle_deg = -angle_deg if left else angle_deg
|
||
|
||
if observe:
|
||
# 输出停止过程中角度变化
|
||
stop_angle_diff = math.degrees(final_yaw - pre_stop_yaw)
|
||
while stop_angle_diff > 180:
|
||
stop_angle_diff -= 360
|
||
while stop_angle_diff < -180:
|
||
stop_angle_diff += 360
|
||
print(f"停止过程中旋转了: {stop_angle_diff:.2f}度")
|
||
print(f"旋转完成，目标角度: {target_angle_deg:.2f}度, 实际角度: {final_angle_deg:.2f}度")
|
||
|
||
# 检查是否需要微调
|
||
angle_error = abs(final_angle_deg - target_angle_deg)
|
||
|
||
# 针对不同角度使用不同的微调阈值
|
||
adjustment_threshold = 3.0 # 默认微调阈值
|
||
if angle_deg >= 170:
|
||
adjustment_threshold = 10.0 # 大角度旋转使用更大的微调阈值，因为初步旋转已经确保大致方向正确
|
||
|
||
if angle_error > adjustment_threshold: # 如果误差超过阈值
|
||
if observe:
|
||
print(f"角度误差: {angle_error:.2f}度，进行微调")
|
||
|
||
# 计算需要微调的角度
|
||
adjust_angle = target_angle_deg - final_angle_deg
|
||
|
||
# 增加补偿系数，因为实际旋转常会超出理论计算值
|
||
if angle_deg <= 70:
|
||
compensation_factor = 0.55 # 比之前更保守，从0.65降至0.55
|
||
elif angle_deg >= 170: # 单独处理180度附近的大角度
|
||
# 动态调整补偿因子，更保守处理
|
||
if angle_error > 60: # 如果误差超过60度
|
||
compensation_factor = 0.4 # 从0.5降至0.4
|
||
elif angle_error > 30: # 如果误差超过30度
|
||
compensation_factor = 0.35 # 从0.4降至0.35
|
||
else:
|
||
compensation_factor = 0.25 # 从0.3降至0.25，更保守
|
||
else:
|
||
compensation_factor = 0.6 # 从0.7降至0.6
|
||
|
||
# 对于特别大的误差进行更保守的单次校准
|
||
if angle_deg >= 170 and angle_error > 30:
|
||
# 大误差直接一步校准，使用更保守的补偿系数
|
||
adjusted_compensation = compensation_factor * 0.85 # 在已有补偿基础上再降低15%
|
||
adjust_angle = adjust_angle * adjusted_compensation
|
||
|
||
if observe:
|
||
print(f"误差较大({angle_error:.2f}度)，使用更保守的补偿系数{adjusted_compensation:.2f}进行一次性校准: {adjust_angle:.2f}度")
|
||
|
||
# 使用turn_degree函数进行微调
|
||
turn_success = turn_degree(ctrl, msg, adjust_angle, absolute=False)
|
||
|
||
# 等待稳定
|
||
time.sleep(1.0)
|
||
|
||
# 再次检查剩余误差(仅用于显示)
|
||
if observe:
|
||
current_yaw = ctrl.odo_msg.rpy[2]
|
||
current_angle_turned = current_yaw - start_yaw
|
||
# 角度归一化处理
|
||
while current_angle_turned > math.pi:
|
||
current_angle_turned -= 2 * math.pi
|
||
while current_angle_turned < -math.pi:
|
||
current_angle_turned += 2 * math.pi
|
||
|
||
current_angle_deg = math.degrees(current_angle_turned)
|
||
remaining_error = target_angle_deg - current_angle_deg
|
||
print(f"校准后的角度: {current_angle_deg:.2f}度, 剩余误差: {remaining_error:.2f}度")
|
||
else:
|
||
adjust_angle *= compensation_factor
|
||
|
||
if observe:
|
||
print(f"应用补偿系数{compensation_factor}后的微调角度: {adjust_angle:.2f}度")
|
||
|
||
# 使用turn_degree函数进行微调
|
||
turn_success = turn_degree(ctrl, msg, adjust_angle, absolute=False)
|
||
|
||
# 等待稳定(仅用于显示)
|
||
if observe:
|
||
time.sleep(0.8)
|
||
current_yaw = ctrl.odo_msg.rpy[2]
|
||
current_angle_turned = current_yaw - start_yaw
|
||
# 角度归一化处理
|
||
while current_angle_turned > math.pi:
|
||
current_angle_turned -= 2 * math.pi
|
||
while current_angle_turned < -math.pi:
|
||
current_angle_turned += 2 * math.pi
|
||
|
||
current_angle_deg = math.degrees(current_angle_turned)
|
||
remaining_error = target_angle_deg - current_angle_deg
|
||
print(f"校准后的角度: {current_angle_deg:.2f}度, 剩余误差: {remaining_error:.2f}度")
|
||
|
||
if observe:
|
||
print(f"微调结果: {'成功' if turn_success else '失败'}")
|
||
|
||
# 移动完成后再检查一次QR码扫描结果
|
||
if scan_qrcode:
|
||
qr_data, scan_time = ctrl.image_processor.get_last_qr_result()
|
||
if qr_data and (qr_result is None or scan_time > last_qr_check_time):
|
||
qr_result = qr_data
|
||
print(f"🔍 旋转完成后最终扫描到QR码: {qr_data}")
|
||
|
||
# 停止异步扫描
|
||
ctrl.image_processor.stop_async_scan()
|
||
|
||
if observe:
|
||
print("圆弧转弯完成")
|
||
|
||
# 根据scan_qrcode参数返回不同格式的结果
|
||
if scan_qrcode:
|
||
return True, qr_result
|
||
else:
|
||
return True
|
||
|
||
# 用法示例
|
||
if __name__ == "__main__":
|
||
move_to_hori_line(None, None, observe=True)
|
||
# 90度左转
|
||
arc_turn_around_hori_line(None, None, angle_deg=90, left=True, observe=True)
|
||
# 180度右转
|
||
arc_turn_around_hori_line(None, None, angle_deg=180, left=False, observe=True)
|