plot以船形状画点

葫芦爷儿

群里有人问到了这样的船舶如何画，见图1。我是这样做，1 得到图像的外形，写入ship函数 2 获取目标位置信息、航首向信息 3 通过旋转矩阵对ship进行变形 4画图
其中3的步骤是：

• 利用M矩阵平移图形中心到原点
• 利用M1矩阵根据heading进行旋转
• 利用M2矩阵平移至坐标点
  

得到结果如图2.
参考网址：https://blog.csdn.net/zhang11wu4/article/details/49761121
                https://blog.csdn.net/Bryan_QAQ/article/details/78805201


参考网址中有欧拉公式；在ship.m有说明；附件ship测试.zip 里有函数、测试数据、测试脚本。如有错误，感谢指正～

调用方法：

1. index =100;
   
2. sta =250;
   
3. longtitude_use = longitude(sta:end-index);
   
4. latitude_use = latitude(sta:end-index);
   
5. heading_use = heading(sta:end-index);
   
6. %plot(longtitude_use,latitude_use,'r')
   
7. hold on
   
8. for index  = 1:1:length(heading_use)
   
9.     ship(heading_use(index),longtitude_use(index),latitude_use(index), 0.00003)
   
10. end
    
11. axis equal

复制代码

图1 目标效果

图2 结果

1. % 本程序画个小船…
   
2. % 输入数据 ：longitude，latitude，theta， scale
   
3. %                     分别是经度、纬度、船首向(deg)，船相对大小
   
4. % 思路：1 船的形状 2 先平移动到原点 3 旋转矩阵 4 平移到实际位置
   
5. %
   
6. % 调用方式：
   
7. % for index  = 1:5:length(heading)
   
8. %     ship(heading(index),longtitude(index),latitude(index), 0.00003)
   
9. % end
   
10. % 参考网址：https://blog.csdn.net/zhang11wu4/article/details/49761121
    
11. %                      https://blog.csdn.net/Bryan_QAQ/article/details/78805201
    
12. % 作者：葫芦爷儿  
    
13. % 联系方式 719087413@qq.com
    
14. 
    
15. function ship(theta, longitude, latitude,scale)
    
16. % 船的形状
    
17. shipy = -scale.*([1, 3, 5, 7, 15, 15, 7, 5, 3, 1]-8);
    
18. shipx = -scale.*([8, 9, 10, 11, 11, 5, 5, 6, 7, 8]-8);
    
19. % 船的中心
    
20. center_point = -scale.* [8,8];
    
21. % heading转换弧度
    
22. theta = pi*theta/180;
    
23. %%  平移到原点，因为shipy-8，已经平移了
    
24. M=[1    0     0;
    
25.        0    1     0;
    
26.        0    0     1;];
    
27. R(1,:)=shipx;
    
28. R(2,:)=shipy;
    
29. R(3,:)=1;
    
30. R1=M*R;
    
31. %% 旋转矩阵
    
32. M1=[cos(theta)   sin(theta)     0;
    
33.     -sin(theta)   cos(theta)    0;
    
34.         0           0         1];
    
35. R2=M1*R1;
    
36. %% 平移到实际位置
    
37. M3=[1    0    longitude;
    
38.        0    1     latitude;
    
39.        0    0     1;];
    
40. R3=M3*R2;
    
41. plot(longitude, latitude,'ko','MarkerFaceColor','k','MarkerSize',2)
    
42. plot(R3(1,:), R3(2,:),'k')

复制代码

sreveryone

感谢分享十分给力

hunterren

非常给力，感谢

zoeyyuanyuan

没看懂这是干嘛用的