豆粑粑 matlab 在圆形的区域内contourf，再挖掉一部分

meatball1982

别人的问题，在一个圆里面画contourf的图，再挖去一块。

http://www.ilovematlab.cn/forum. ... p;extra=#pid3221245

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. [x,y,z]=peaks(40);
   
6. r_big=linspace(0,3,400);
   
7. the=linspace(0,2*pi,720);
   
8. the=[2*pi,the];
   
9. 
   
10. [R_big,The]=meshgrid(r_big,the);
    
11. X_cir_big=R_big.*cos(The);
    
12. Y_cir_big=R_big.*sin(The);
    
13. 
    
14. Z_cir_big=griddata(x(:),y(:),z(:),X_cir_big,Y_cir_big);
    
15. 
    
16. Cir_big_x =r_big(end)*cos(the);
    
17. Cir_big_y =r_big(end)*sin(the);
    
18. 
    
19. 
    
20. r_sma=0.9;
    
21. cen_x=-1.7;
    
22. cen_y=0.3;
    
23. 
    
24. Cir_sma_x =r_sma*cos(the)+cen_x;
    
25. Cir_sma_y =r_sma*sin(the)+cen_y;
    
26. 
    
27. 
    
28. 
    
29. ind=((X_cir_big-cen_x).^2+(Y_cir_big-cen_y).^2)<r_sma^2;
    
30. 
    
31. % X_cir_big(ind)=nan;
    
32. % Y_cir_big(ind)=nan;
    
33. Z_cir_big(ind)=nan;
    
34. hold on
    
35. % surf(X_cir_big,Y_cir_big,Z_cir_big,'edgecolor','none');
    
36. [C,h]=contourf(X_cir_big,Y_cir_big,Z_cir_big,20,'linecolor','k','linestyle','--');
    
37. clabel(C,'fontSize',8)
    
38. colorbar
    
39. plot(Cir_big_x,Cir_big_y,'k','linewidth',1)
    
40. plot(Cir_sma_x,Cir_sma_y,'k','linewidth',1)
    
41. axis equal
    
42. axis off
    
43. colormap(jet)
    
44. h=gcf;
    
45. fi_na=['./fig_cont_2D'];
    
46. fun_work_li_035_myfig_out(h,fi_na,3);

复制代码