matlab读取.bmp文件显示图形？

sunminmin

用readbmp.m能读取dx.bmp显示图形，为什么我把I=imread('dx.bmp'); %读入一幅图片改成I=imread('map.bmp'); 把gao=[0.1:0.1:0.7];[x1,y1]=meshgrid(1:0.1:5)做相应修改，还是不能读取map.bmp显示图形。提示??? Index exceeds matrix dimensions. Error in ==> dizhenyuce at 22   BW=II(:,:,1)+II(:,:,2)+II(:,:,3);

下面为readbmp.m文件：

1. I=imread('dx.bmp'); %读入一幅图片
   
2. %I=I(33:375,74:508,:); %须事先人工确定行标和列标
   
3. 
   
4. %提取坐标框内部分
   
5. I1=I(:,:,1); %提取红色灰度矩阵
   
6. pixel_black=(I1==0);
   
7. number_col=sum(pixel_black);
   
8. number_row=sum(pixel_black');
   
9. j1=find(number_col==max(number_col));
   
10. number_col(j1)=-1;
    
11. j2=find(number_col==max(number_col));
    
12. 
    
13. i1=find(number_row==max(number_row));
    
14. number_row(i1)=-1;
    
15. i2=find(number_row==max(number_row));
    
16. I=I(min(i1,i2):max(i1,i2),min(j1,j2):max(j1,j2),:);
    
17. 
    
18. %提取各等高线上颜色灰度值
    
19. II=I;
    
20. II(II==255)=0;
    
21. II(II>0)=255;
    
22. BW=II(:,:,1)+II(:,:,2)+II(:,:,3);
    
23. [row_id,col_id]=find(BW);
    
24. up_row_id=min(row_id);
    
25. up_col_id=col_id(round(mean(find(row_id==up_row_id))));
    
26. left_col_id=min(col_id);
    
27. left_row_id=row_id(round(mean(find(col_id==left_col_id))));
    
28. pixel_BW=BW(up_row_id:left_row_id,up_col_id);
    
29. row_id=find(pixel_BW)-1+up_row_id;
    
30. yanse(:,1:3)=I(row_id,up_col_id,:);
    
31. %yanse=[0 0 143;
    
32. % 0 47 255
    
33. % 0 223 255
    
34. % 143 255 111
    
35. % 255 207 0
    
36. % 255 31 0
    
37. % 127 0 0]; %须事先人工确定各等高线颜色灰度值，用到了getpts函数
    
38. 
    
39. [m,n]=size(I(:,:,1));
    
40. BW=zeros([m,n]);
    
41. x_xishu=4/(n-1);
    
42. y_xishu=4/(m-1);
    
43. gao=[0.1:0.1:0.7];
    
44. 
    
45. %提取各等高线上点的三维坐标
    
46. xyz=[];
    
47. for i=1:7
    
48. BW1=BW;
    
49. BW1(I(:,:,1)==yanse(i,1)&I(:,:,2)==yanse(i,2)&I(:,:,3)==yanse(i,3))=1;
    
50. [y,x]=find(BW1);
    
51. xyz=[xyz;1+(x-1)*x_xishu 1+(y-1)*y_xishu gao(i)*ones(size(x))];
    
52. end
    
53. xyz=xyz(1:3:end,:);
    
54. xyz=[xyz;1 1 0;1 2 0;1 3 0;1 4 0;1 5 0;2 1 0;3 1 0;4 1 0;5 1 0;5 2 0;5 3 0;5 4 0;5 5 0;2 5 0;3 5 0;4 5 0];
    
55. x=xyz(:,1);
    
56. y=xyz(:,2);
    
57. z=xyz(:,3);
    
58. 
    
59. %产生网格数据，进行散乱节点插值拟合
    
60. [x1,y1]=meshgrid(1:0.1:5)
    
61. z1=griddata(x,y,z,x1,y1,'v4');
    
62. 
    
63. [n,m]=size(x1);
    
64. 
    
65. xyz1=[reshape(x1,n*m,1) reshape(y1,n*m,1) reshape(z1,n*m,1)]
    
66. 
    
67. 
    
68. %画三维面图
    
69. figure
    
70. surf(x1,y1,z1)
    
71. shading flat; %各小曲面之间不要网格
    
72. zlim([0 0.8])
    
73. xlabel('X(单位：千米)')
    
74. ylabel('Y(单位：千米)')
    
75. zlabel('Z(单位：千米)')
    
76. 
    
77. 
    

复制代码

怎样能读取任何.bmp文件显示图形，做到不需要每次读一个.bmp文件就要重新写一个代码？

ChaChing

回复 1 # sunminmin 的帖子

Ref 12F, 常见的程序出错问题整理 http://forum.vibunion.com/thread-46001-1-1.html
(from http://forum.vibunion.com/home-s ... -blog-id-18250.html)

sunminmin

回复 2 # ChaChing 的帖子

谢谢！

1. a = [1];
   
2. b = [1];
   
3. c = [1];
   
4. plot3(a,b,c,'s')

复制代码

用上面的代码能生成一个点，得到用一个小方格表示的，怎样写才能做到用黑点表示该点？

ChaChing

回复 3 # sunminmin 的帖子

plot3(a,b,c,'k.','MarkerSize',12)

sunminmin

回复 4 # ChaChing 的帖子

怎样用一条直线把两个点（1，2,3）（5,2,8）连起来？

ChaChing

回复 5 # sunminmin 的帖子

我都已上高速公路跑了几百公里回工作地了, 楼主还在!?:@)
plot3([1,5],[2,2],[3,8])
或
line([1,5],[2,2],[3,8]); view(3)

meiyongyuandeze

回复 6 # ChaChing 的帖子

几百公里去工作地，好辛苦！

ChaChing

回复 1 # sunminmin 的帖子

clc; clear; I=imread('dx.bmp'); I2=imread('map.bmp');
whos看看两者维数的不同(421x561x3, 698x1280)!

ChaChing

回复 7 # meiyongyuandeze 的帖子

工作地与老家较远些! 不过相较内地应该还好:@)

sunminmin

回复 6 # ChaChing 的帖子

呵呵，:@)我是想实现下面这个用三维坐标（x,y,z）画折线的功能

1. 1.        Engine *ep;
   
2. 2.       
   
3. 3.        if (!(ep = engOpen(NULL))) //打开Matlab引擎
   
4. 4.       
   
5. 5.        {
   
6. 6.       
   
7. 7.        MessageBox ("Can't start MATLAB engine", MB_OK);
   
8. 8.       
   
9. 9.        exit(-1);
   
10. 10.        }
    
11. 11.        mxArray *xx=mxCreateDoubleMatrix(1,6,mxREAL);
    
12. 12.        mxArray *yy=mxCreateDoubleMatrix(1,6,mxREAL);
    
13. 13.        mxArray *zz=mxCreateDoubleMatrix(1,6,mxREAL);
    
14. 14.       
    
15. 15.        double x[] ={0,1,2,3,4,5}, y[] = {0,1,4,9,16,25}, z[] = {1,2,3,4,5,6};
    
16. 16.       
    
17. 17.        engPutVariable(ep,"xx",xx);
    
18. 18.        engPutVariable(ep,"yy",yy);
    
19. 19.        engPutVariable(ep,"zz",zz);
    
20. 20.       
    
21. 21.        memcpy(mxGetPr(xx),x,6*sizeof(double));
    
22. 22.        memcpy(mxGetPr(yy),y,6*sizeof(double));
    
23. 23.        memcpy(mxGetPr(zz),z,6*sizeof(double));
    
24. 24.       
    
25. 25.        engEvalString(ep,"绘制6个三维点连成的线");？？？
    
26. 26.        接下来用循环如何实现6点相连接？？？
    
27. 
    

复制代码

就是利用matlab引擎在三维空间将这些坐标（x,y,z）所表示的点用黑点表示出来,前后两个点之间用直线连接，（0,0,1）连接（1,1,2）连接（2,4,3）连接（3,9,4）连接（4,16,5）连接（5,25,6），形成一条折线。

sunminmin

回复 10 # sunminmin 的帖子

VC++与matlab混合编程利用matlab引擎做的想要实现在三维空间画折线的功能，上面的代码接下来该如何处理??

smtmobly

plot3不行?

meiyongyuandeze

一直很纳闷，楼主费劲周折在matlab中画出曲面后再求出曲面点，在用这些点画直线，好像程序都是在matlab中编写的！