MATLAB——文件读写(2)


一、importdata函数

1、

txt 如图,提取经纬度。

程序如下

clear all test=importdata('经纬度.txt'); [r,c]=size(test.data);%row行,column列 a=max(size(test.data));%行 l=length(test.data); for i=1:r lon(i)=test.data(i,8);%经度 lat(i)=test.data(i,7); end plot(lon,lat); xlabel('经度','fontsize',13,'fontweight','bold'); ylabel('纬度');
title('轨迹')  grid on

      

可以看出importdata将读取的数据生成了一个struct结构体,结构体将txt中的数据分类。所用的数据如下

循环调用经纬度可以画出如图所示的图像。

 2、

clear all a=importdata('road.png') image(a)

二、excel写入

clear all filename = 'testdata.xls'; A = [12.7 5.02 -98 63.9 0 -.2 56]; xlswrite(filename,A)

 

clear all filename = 'testdata.xls'; A = [12.7,5.02,-98;63.9,0.2,56]; xlswrite(filename,A)

 

clear all filename = 'testdata.xlsx'; A = {'Time','Temperature'; 12,98; 13,99; 14,97}; sheet = 2;%第几个表格 xlRange = 'E1';%第一个位置 xlswrite(filename,A,sheet,xlRange)

 

三、Excel读取

1、

数据格式是文本类型

clear all a=xlsread('气象生态环境监测系统报表0102.xls'); [num,txt,raw] = xlsread('气象生态环境监测系统报表0102.xls')

可以看出返回的及时元胞数组,a中是空的,因为Excel中没有数据类型,只有文本,对比可以看出下面的num中也是空的。所以计算或者画图时需要转换文本格式。

2、

clear all a=xlsread('气象生态环境监测系统报表0102.xls'); [num,txt,raw] = xlsread('气象生态环境监测系统报表0102.xls'); date=txt(4:end,1);%读取日期 tem=txt(4:end,7);%读取温度 date(end-1:end)=[]%去除无效数据

 日期670和671是没有数据的,通过[]去除。温度同理。再有,通过end,可以任意读取某一行到最后的行和列的内容。

3、剔除异常值

clear all a=xlsread('0102.xls'); [num,txt,raw] = xlsread('0102.xls'); date=txt(4:end,1);%读取日期 tem=txt(4:end,7);%读取温度 date(end-1:end)=[];%去除日期空数据 tem(end-1:end)=[];%去除温度空数据 n=datenum(date,'yyyy-mm-dd HH:MM');%日期数字 t=str2num(char(tem));%将cell类型,转换成double类型 s=length(t); w=1+0.4*log(s-2); a=mean(t);%均值 e=std(t);%标准偏差 yichang=abs(t-a)>w*e; %判断语句,肖威勒方法,返回0或者1. t(yichang)=a;  %将异常值重新赋值 figure subplot(211),plot(n,t); subplot(212),plot(n,t); datetick('x','HH:MM');%更改坐标形式 grid on

 

 

四、多文件读取

 1、

clear all dirs=dir('C:\Users\Administrator\Desktop\test\*.xls')

生成这样一个结构体,只需要将名字提取出来,然后循环提取就可以处理这个文件夹下的所有的xls文件了。注意正则表达式的应用。

clear all dirs=dir('C:\Users\Administrator\Desktop\test\*.xls'); l=length(dirs); for i =1:l names{i}=dirs(i).name; [num{i},txt{i},raw{i}] = xlsread(names{i}); date{i}=[txt{i}(4:end,1)];%读取日期 tem{i}=[txt{i}(4:end,7)]; date{i}(end-1:end)=[]; tem{i}(end-1:end)=[]; d{i}=datenum(date{i},'yyyy-mm-dd HH:MM'); t{i}=str2num(char(tem{i})); figure plot(d{i},t{i}); datetick('x','HH:MM'); grid on %filename=[names{i} '.jpg']; saveas(gcf,[names{i} '.jpg']);%注意空格,否则出错 end

注意保存时字符串和图片符号之间的空格,否则出现如下错误。

 

 

 

 

 

 

 

 

 

 

2、

>> dirs(1) ans = name: '气象生态环境监测系统报表0101.xls' date: '02-一月-2017 09:04:38' bytes: 436151 isdir: 0 datenum: 7.3670e+05 >> dirs{1} 元胞内容引用自非元胞数组对象。

需要根据对方的类型选择大括号和小括号。

 

>> a=txt(1) a = {556x36 cell}

>> a=txt{1};

可以看出用{}比()提前一步获取cell中的信息。

a=txt{1} b=a(:,1)

b(1:3)=[]

>> b{end-1:end}=[]
需要大括号或点索引表达式中的一个输出,但结果有 2 个。

 无效的空赋值。’

 

a=txt(1)

b=a{:,1}

 

殊途同归。

 3、保存形式

所以能用下标访问

 五、mat 文件生成读写

 

clear all sheng=shaperead('maps/bou2_4p.shp', 'UseGeoCoords', true) geoname=[sheng.NAME]

生成了一个字符串 

clear all

sheng=shaperead('maps/bou2_4p.shp', 'UseGeoCoords', true)
geoname=sheng.NAME

只有黑龙江省

 

clear all

sheng=shaperead('maps/bou2_4p.shp', 'UseGeoCoords', true)
geoname={sheng.NAME}

 

生成元胞数组

 

clear all p=shaperead('maps/bou2_4p.dbf'); l=length(p); for i=1:l a{i}=p(i).NAME; %读取省份,注意用元胞数组保存,否则保存不了 end a=a'; %转置
 save mydata a % 保存为dat文件

双击,载入mydata.mat 文件

六、时间异常处理

Excel如下

clear;clc;close all
[Das,Xls]=xlsread('中心80.xlsx');  % 依次得到.xls的数据,文本矩阵,存放在IMs元胞数组里
Das(1,1)
Xls(1,1)

 

接下来将时间格式转换下

datestr(Xls{1,1},29)

出现这种情况的原因是还有星期二文本,所以失败,正确程序如下

datestr(Xls{1,1}(1:8),29)%只读取星期二之前的数据

 

将时分秒,从double类型,变为时分秒类型 21:26:05

Das(1,1)
datestr(Das(1,1),'HH:MM:SS')

出现差错

s=mod(Das(1,1),100)
m=(mod(Das(1,1),10000)-s)/100
h=floor(Das(1,1)/10000)

分别提取出时分秒

t=hours(h)+minutes(m)+seconds(s)
datestr(t,'HH:MM:SS')

还原出时分秒

七、图片读取保存

(1)读取

clear all
[filename,pathname]=uigetfile({'*.jpg';'*.bmp';'*.gif'},'选择图片');
if isequal(filename,0)
    disp('Users Selected Canceled');
else
str=[pathname filename];
im = imread(str);
imshow(im);
end

(2)保存

clear all
x=0:0.01:2*pi;
plot(x,sin(x));

[filename,pathname]=uiputfile({'*.bmp';},'保存图片');%路径和图片名
if ~isequal(filename,0)
    str = [pathname filename];%路径名
    data= getframe(gcf);%图片内容数据
    imwrite(data.cdata,str,'bmp');%数据,文件名,图片格式
%     saveas(gcf,str,'bmp');%两种方式都可以保存图片,gcf返回当前图片句柄
    close(gcf);
 else
    disp('保存失败');
end

(3)imwrite函数

clear all
A = rand(50);
imwrite(A,'C:\Users\ning wang\Desktop\123.jpg')

 


免责声明!

本站转载的文章为个人学习借鉴使用,本站对版权不负任何法律责任。如果侵犯了您的隐私权益,请联系本站邮箱yoyou2525@163.com删除。



 
粤ICP备18138465号  © 2018-2025 CODEPRJ.COM