关于Matlab在雷达数据处理中的一些经验分享

置顶珞瑜·

已于 2022-03-18 19:45:27 修改

阅读量7.1k

点赞数 20

分类专栏： # Matlab基础内容文章标签： matlab

于 2020-05-27 20:01:57 首次发布

本文链接：https://blog.csdn.net/weixin_45770896/article/details/106388969

版权

Matlab基础内容专栏收录该内容

110 篇文章 80 订阅

订阅专栏

1.滤波：

2.交换行与行之间的数据：

3.matlab矩阵运算很快捷，用一下整体乘法，提高计算效率。

4.[SNR,I] = min(Data_Ze,[],1,'omitnan');

10.matlab求和取均值以及计算标准差时

11.Matlab 简化运算

12.完成仰角-方位-距离三维元素位置的转换

13.shading阴影函数

14.Matlab中计算程序运行时间的方法

15.数据计算均值、标准偏差时，边缘数据处理不到

1.滤波：

函数medfilt1,smooth,filter,无法解决已存在的数据NAN，污染周边数据。

2.交换行与行之间的数据：

aa([1,3],:) = aa([3,1],:)。注释：交换二维数组第1行与第3行之间的数据。

3.matlab矩阵运算很快捷，用一下整体乘法，提高计算效率。

相比for循环，进入每个数据去加减乘除，换个角度去思考，用矩阵整体运算，例如：Arry1 .* Arry2;

4.[SNR,I] = min(Data_Ze,[],1,'omitnan');

%排除每一列中的NAN，求每一列中的最小值，SNR为每列的最小值，I为每列最小值的下标。

5.ncdisp()：

可以充分地体现nc文件的解码内容。

6.命令行窗口显示信息：

X = ['仰角',num2str(elevtation_number), '下的ML总点数：']; % 命令窗口要显示的信息

disp(X);

7.弹出提示框显示信息

msgbox('中间数据图像显示完毕！','温馨提示','modal');

errordlg('参数超出范围或者参数输入不足，请重新配置！', 'Warning');

8.在Figure图中标记数据

str=[ num2str(r') 'km']; % 显示的文本

text(x,y,cellstr(str)); % (x,y)要显示的位置

9.标记Figure图中线条内容

l = legend('温度廓线','露点温度廓线'); % 标注线条代表内容

title(l,'线条表示内容'); % 线条标题

10.matlab求和取均值以及计算标准差时

数据有nan,可以使用下面三个特定的函数（nansum，nanmean，nanstd）

11.Matlab 简化运算

ZH_Data1=mod((ZH_Data0+33)*2,256)*0.5-33; % 重新解码校正

ZH_Data1(ZH_Data1==-33)=NaN; % 将数值为-33的数据标记为NAN无效数据

ZH_Data=ZH_Data1(1:1192,:,1); % 截取距离库1:1192的数据

12.完成仰角-方位-距离三维元素位置的转换

permute(ZH_data,[3,2,1]);

13.shading阴影函数

控制曲面和图形对象的颜色着色，即用来处理色彩效果的，包括以下三种形式:

shading faceted：默认模式，在曲面或图形对象上叠加黑色的网格线；

shading flat：是在shading faceted的基础上去掉图上的网格线；

shading interp：对曲面或图形对象的颜色着色进行色彩的插值处理，使色彩平滑过渡；

14.Matlab中计算程序运行时间的方法

方法一：

tic; Function(); toc;

注：toc计算的是与最后一次运行的tic之间的时间。

方法二：

T1=clock；Function(); T2=clock; Etime(T2,T1);

方法三：

T0=cputime; Function(); T1=cputime-T0;

15.数据计算均值、标准偏差时，边缘数据处理不到

在计算前，定义一个nan数组,即可解决边缘数据在赋值的问题；

Hei = nan(elevtation_number,radialNum,maxBinNum);

16.matalb中的希腊字体标识

17.色标右侧文字位置控制

hh=colorbar('location','EastOutside');hold on; % 色标右侧文字

set(hh,'ytick',1.5:1:length(Z_index)+1); % 1.5是初始位置，1为步长

18.选择文件并存储文件路径

[Filename, Pathname]  = uigetfile('*.*', '选择所用的雷达数据');   % 获取文件名和路径
filePath= [Pathname  Filename];                                  % 存储文件路径

Pathname ：文件上层路径

Filename ：文件名

19.读数据时，如何确定数据尾的标志？

FID = fopen(FILENAME)；
ST = feof(FID) ；

feof 解释：如果ST = feof（FID）的文件结束，指示符返回1。

returns 1 if the end-of-file indicator for the file with file identifier FID has been set, and 0 otherwise.
The end-of-file indicator is set when a read operation on the file associated with the FID attempts to read past the end of the file.

20.fseek在读取数据时的用法

fseek Set file position indicator. 
    
STATUS = fseek(FID, OFFSET, ORIGIN) repositions the file position indicator in the file associated with the given FID.  fseek sets the position indicator to the byte with the specified OFFSET relative to ORIGIN.
 
FID is an integer file identifier obtained from FOPEN.
 
OFFSET values are interpreted as follows:
   >= 0    Move position indicator OFFSET bytes after ORIGIN.
   < 0    Move position indicator OFFSET bytes before ORIGIN.
 
ORIGIN values are interpreted as follows:
    'bof' or -1   Beginning of file
    'cof' or  0   Current position in file
    'eof' or  1   End of file
 
STATUS is 0 on success and -1 on failure.  If an error occurs, use FERROR to get more information.

翻译如下：

fseek：查找设置文件位置指示器。

STATUS = fseek（FID，OFFSET，ORIGIN）在与给定FID关联的文件中重新定位文件位置指示符。 
fseek将位置指示器设置为相对于ORIGIN具有指定OFFSET的字节。
 
FID是从FOPEN获得的整数文件标识符。
 
OFFSET值的解释如下：
    > = 0，在ORIGIN之后移动位置指示器OFFSET字节。
    < 0，在ORIGIN之前移动位置指示器OFFSET字节。
 
ORIGIN值的解释如下：
    'bof'或-1，文件开头
    'cof'或0，文件中的当前位置
    'eof'或1，文件结尾
 
STATUS成功时为0，失败时为-1。 如果发生错误，请使用FERROR获取更多信息。