CTD数据处理及可视化.doc

第7周上机操作
一、大家以前接触过CTD数据,关于CTD数据的导出各位同学已有基础,不知道的同学可以向同学学一下,这里我们用CTD导出的ASCII码数据(后缀名为*.cnv),进行读取等简单的处理画图。
要求:
1 用fopen、fgetl、fscanf以及fclose命令实现数据的读入压力、温度及盐度;
2、首先看文件的格式,弄清楚:数据从多少行开始,多少行结束、有多少列以及数据的各列含义;
3. 实现操作为
fid=fopen('v','r');
文件数据读取操作
for i=1:96
fgetl(fid); %把文件的说明部分略过
end
data=scanf(fid,'%f',[12 inf]);% 把数据读入到矩阵,注意为什么行要用12?
data=data';
%%% 大家注意 name 0 v数据的第一列
d_pressure=data(:,2); % 压力(等价于深度),也即说明文件的 name 1 对应于第2列
T=data(:,3); %温度, 也即说明文件的 name 2 对应于第3列
S=data(:,8); %盐度, 也即说明文件的 name 7 对应于第8列
%绘出温度-深度图、盐度深度图以及T-S散点图
subplot(2,2,1)
绘制点线图,大家认真体会
plot(t,dpth,'--k');
set(gca,'ydir','reverse','fontsize',15);
title('Depth-Temp Diagram')
xlabel('Temperature (^0C)');
ylabel('Depth (m)');
ylim([0 100])
subplot(2,2,2)
plot(s,dpth,'--b')
set(gca,'ydir','reverse','fontsize',15);
title('Depth-Salinity Diagram')
xlabel('Salinity (psu)');
ylabel('Depth (m)');
ylim([0 100])
subplot(2,2,3)
scatter(s,t)
4 要求对其余几个数据做同样操作,绘出相应的图给我看。记得用figure打开新的图形窗口
二、 HDF数据的加载:
什么是HDF?
HDF是一个能够自我描述、多目标、用于科学数据存储和分发的数据格式(http://baike./view/)。。可用命令 hdftool以及hdfread
读写操作():
首先用hdftool看文件的格式(或用hdfview),知道第一维、第二维内容,
维数“(主要查看:Northernmost Latitude; Southernmost Longitude;Western Longitude;以及Latitude step以及Longitude step);理解数据是由北向南,由东向西排列,这很重要。
用hdftool浏览: 1)在命令行敲入 hdftool回车
2)当出现HDF Import Tool界面时,在File