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=:r
    lon(i)=test.data(i,);%经度
    lat(i)=test.data(i,);
end
plot(lon,lat);
xlabel('经度','fontsize',,'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 - 63.9  -. ];
xlswrite(filename,A)

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

clear all
filename = 'testdata.xlsx';
A = {'Time','Temperature'; ,; ,; ,};
sheet = ;%第几个表格
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(:end,);%读取日期
tem=txt(:end,);%读取温度
date(end-:end)=[]%去除无效数据

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

3、剔除异常值

clear all
a=xlsread('0102.xls');
[num,txt,raw] = xlsread('0102.xls');
date=txt(:end,);%读取日期
tem=txt(:end,);%读取温度
date(end-:end)=[];%去除日期空数据
tem(end-:end)=[];%去除温度空数据

n=datenum(date,'yyyy-mm-dd HH:MM');%日期数字
t=str2num(char(tem));%将cell类型，转换成double类型

s=length(t);
w=+0.4*log(s-);
a=mean(t);%均值
e=std(t);%标准偏差
yichang=abs(t-a)>w*e; %判断语句,肖威勒方法,返回0或者1.
t(yichang)=a;  %将异常值重新赋值

figure
subplot(),plot(n,t);
subplot(),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 =:l
    names{i}=dirs(i).name;
    [num{i},txt{i},raw{i}] = xlsread(names{i});
    date{i}=[txt{i}(:end,)];%读取日期
    tem{i}=[txt{i}(:end,)];
    date{i}(end-:end)=[];
    tem{i}(end-: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()

ans = 

       name: '气象生态环境监测系统报表0101.xls'
       date: '02-一月-2017 09:04:38'
      bytes:
      isdir:
    datenum: 7.3670e+05

>> dirs{}
元胞内容引用自非元胞数组对象。

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

>> a=txt()

a = 

    {556x36 cell}

>> a=txt{};

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

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

b(:)=[]

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

无效的空赋值。’

a=txt()

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=: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')