豆粑粑 matlab 极坐标下的散点分页

meatball1982

在极坐标下，把两组散点用等高线的方式表达出来。
因为散点基本在一起，分不清具体的方向。

图1 是两组数的极点坐标下的位置。
图2 是连线的版本。
图3 是出现次数的统计(平滑后的)。两个不同的colormap代表了不同的数据的分页，散点根据colormap的最低color着色。

这时的统计和帖子
https://www.ilovematlab.cn/thread-539582-1-1.html
思路一致。
数据和程序在附件中。









主程序：

1. clear all
   
2. clc
   
3. close all
   
4. 
   
5. 
   
6. dat = load('../file_data/dat.txt');
   
7. 
   
8. %
   
9. R=dat(:,1);
   
10. % Rmax = max(R);
    
11. R1=R;
    
12. R2=R;
    
13. the1 = dat(:,2)*pi/180;
    
14. the2 = dat(:,3)*pi/180;
    
15. 
    
16. 
    
17. 
    
18. h=figure(1)
    
19. set(h, 'Position', [1000, 100, 400, 400]);
    
20. h1 = polar(the1,R1);
    
21. set(h1,'color','b','linewidth',0.5)
    
22. hold on
    
23. h2=polar(the2,R2);
    
24. set(h2,'color','r','linewidth',0.5)
    
25. 
    
26. 
    
27. h=figure(2)
    
28. set(h, 'Position', [1000, 100, 400, 400]);
    
29. h1 = polar(the1,R1,'o');
    
30. set(h1,'markersize',0.8,'markerfacecolor','b','markeredgecolor','none')
    
31. hold on
    
32. h2=polar(the2,R2,'o');
    
33. set(h2,'markersize',0.8,'markerfacecolor','r','markeredgecolor','none')
    
34. 
    
35. 
    
36. 
    
37. figure(3)
    
38. pl_lev =[0.07 0.4 0.7];
    
39. [h]=fun_polar_2group(R1,the1,R2,the2,pl_lev)
    
40. 
    
41. 
    
42. 
    

复制代码

画图函数：

1. function [h]=fun_polar_2group(R1,the1,R2,the2,pl_lev)
   
2. % [h]=fun_polar_2group(R1,the1,R2,the2,pl_lev)
   
3. % inputs
   
4. % R1:     rad of 1st group data
   
5. % the1:   theta of 1st group data
   
6. % R2:     rad of 2nd group data
   
7. % the2:   theta of 2nd group data
   
8. % pl_lev : the contour line value of plot each group
   
9. % outputs
   
10. % h:      handle of the figure
    
11. % first typed by mm
    
12. % contact me : meatball1982@163.com
    
13. 
    
14. %% main
    
15. % cart cor
    
16. x1 = R1.*cos(the1);
    
17. y1 = R1.*sin(the1);
    
18. x2 = R2.*cos(the2);
    
19. y2 = R2.*sin(the2);
    
20. 
    
21. 
    
22. 
    
23. Rmax = max(max(R1),max(R2));
    
24. 
    
25. % grid number
    
26. n_x = 200;
    
27. n_y = 199;
    
28. 
    
29. x_lin = linspace(-1.1*Rmax,1.1*Rmax,n_x);
    
30. y_lin = linspace(-1.1*Rmax,1.1*Rmax,n_y);
    
31. 
    
32. [X_mat,Y_mat]=meshgrid(x_lin,y_lin);
    
33. Z1_mat =zeros(size(X_mat));
    
34. Z2_mat =zeros(size(X_mat));
    
35. 
    
36. 
    
37. % static numb in each grid
    
38. for i = 1:(n_x-1)
    
39.     for j = 1:(n_y-1)
    
40.         ind = x1 > x_lin(i) & x1<=x_lin(i+1) & ...
    
41.               y1 > y_lin(j) & y1<=y_lin(j+1);
    
42.           Z1_mat(j,i)=sum(ind);
    
43.     end
    
44. end
    
45. 
    
46. 
    
47. for i = 1:(n_x-1)
    
48.     for j = 1:(n_y-1)
    
49.         ind = x2 > x_lin(i) & x2<=x_lin(i+1) & ...
    
50.               y2 > y_lin(j) & y2<=y_lin(j+1);
    
51.           Z2_mat(j,i)=sum(ind);
    
52.     end
    
53. end
    
54. 
    
55. % smooth the data, change 9 if necessary.
    
56. Z1_mat_sm = smooth2(Z1_mat,9,9);
    
57. Z2_mat_sm = smooth2(Z2_mat,9,9);
    
58. 
    
59. % colormap
    
60. tm = jet(120);
    
61. col_h1 = flipud(tm(10:50,:));
    
62. tm = hot(100);
    
63. col_h2 = flipud(tm(20:70,:));
    
64. 
    
65. 
    
66. % plot data 1
    
67. ax1 = axes;
    
68. axis off
    
69. axis equal
    
70. hidden off
    
71. h1 = polar(ax1,the1,R1,'o');
    
72. set(h1,'markersize',1,'markeredgecolor','none','markerfacecolor',col_h1(1,:));
    
73. hold on
    
74. [~,hContour1]=contour(ax1,X_mat,Y_mat,Z1_mat_sm,pl_lev,'linewidth',2);
    
75. 
    
76. % plot data 2
    
77. ax2= axes;
    
78. hold on
    
79. plot(x2,y2,'o','markeredgecolor','none','markerfacecolor',col_h2(1,:),'markersize',1);
    
80. axis off
    
81. axis equal
    
82. hidden off
    
83. [~,hContour2]=contour(ax2,X_mat,Y_mat,Z2_mat_sm,pl_lev,'linewidth',2);
    
84. 
    
85. % link two axes
    
86. linkaxes([ax1,ax2]);
    
87. ax2.Visible = 'off';
    
88. ax2.XTick=[];
    
89. ax2.YTick =[];
    
90. 
    
91. colormap(ax1,col_h1);
    
92. colormap(ax2,col_h2);
    
93. 
    
94. % plot position
    
95. set([ax1,ax2],'Position',[0.05 .11 0.685 0.815],'fontsize',12);
    
96. cb1 = colorbar(ax1,'Position',[.75 .11 .0675 .815],'fontsize',12);
    
97. cb2 = colorbar(ax2,'Position',[.88 .11 .0675 .815],'fontsize',12);
    
98. axis([-1.2* Rmax 1.2* Rmax -1.2 * Rmax 1.2* Rmax])
    
99. 
    
100. % output
     
101. h=gcf;

复制代码