color_ratio_frequency.rar

  • keyan1
    了解作者
  • matlab
    开发工具
  • 1KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 3
    下载次数
  • 2019-02-25 21:36
    上传日期
用于所下载的calipso数据中每年气溶胶色比的频率统计
color_ratio_frequency.rar
  • color_ratio_frequency.m
    4.3KB
内容介绍
clear all %清除工作区所有变量 clc %清除命令窗口命令记录 path_data='G:\calipso\level2\alay\2016\1\'; % 将字符串(数据所在文件夹路径)赋值给path_data files0=dir(path_data); % dir函数,得到path_data路径下,所有文件和文件夹的名称,赋值给变量files0 files=files0(arrayfun(@(x) ~strcmp(x.name(1),'.'),files0)); % 将首字符为“.”的文件或文件夹剔除,注:首字符为“.”的为隐藏文件(夹) files=files(arrayfun(@(x) strcmp(x.name(1:3),'CAL'),files)); % 获得前三个字符为CAL的所有文件名称 nfiles=length(files); % 文件个数 i=0; for a=0:5 for b=0:1 for c=0:3 for d=0:3 for e=0:3 i=i+1; allvalue(i)=a*8192+b*4096+c*128+d*32+e*8; % Feature_Classification_Flags 中可以判断为某一类气溶胶(需加上特定数值)所有可能数值 end end end end end count1=0; % 计数 count2=0; count3=0; count4=0; count5=0; count6=0; count7=0; for k=1:nfiles % 循环,每个数据文件都统计 dirname=files(k).name; %第k个数据文件的文件名 Latitude = hdfread([path_data,dirname], '/Latitude'); %读取数据文件中的纬度变量 Longitude = hdfread([path_data,dirname], '/Longitude'); % 读取数据文件中的经度变量 Feature_Classification_Flags = hdfread([path_data,dirname], '/Feature_Classification_Flags'); % 读取 Feature Classification Flags 数据,可根据此数据判断气溶胶类型 Integrated_Particulate_Depolarization_Ratio = hdfread([path_data,dirname], '/Integrated_Particulate_Depolarization_Ratio'); % 读取层积分退偏比数据 % 类似,若读取色比数据,则用以下命令: % Integrated_Particulate_Color_Ratio = hdfread([path_data,dirname], '/Integrated_Particulate_Color_Ratio'); % 注:退偏比取值范围一般为0-1.0之间,色比取值范围要更大,若要统计色比数据,需合稍加修改程序 Lat=Latitude(:,2); % 中心纬度 Lon=Longitude(:,2); % 中心经度 flagx=find(Lon>70 & Lon <105); % 找出经度 >70° 且 <105°的数据 % 根据实际需要修改 flagy=find(Lat>27 & Lat <40); % 找出纬度 >27° 且 <40°的数据 % 根据实际需要修改 domain=intersect(flagx,flagy); % intersect函数,找出两个数组的交集,这里即找出 70°<Lon<105° & 27°<Lat<40°的数据 FCF=Feature_Classification_Flags(domain,:); % 得到满足经纬度范围条件的 Feature Classification Flags 数据 PDR=Integrated_Particulate_Depolarization_Ratio(domain,:); % 得到满足经纬度范围条件的退偏比数据 flag=find(ismember(FCF,allvalue+1027)==1); % 找出气溶胶类型为沙尘的数据记录 % allvalue+1027 即为沙尘气溶胶的所有可能的Feature_Classification_Flags值 % allvalue+515 则气溶胶类型为 clean marine % allvalue+1539 则气溶胶类型为 polluted continental % allvalue+2051 则气溶胶类型为 clean continental % allvalue+2563 则气溶胶类型为 polluted dust % allvalue+3075 则气溶胶类型为 smoke % 因此可以用此命令统计不同类型气溶胶的频率分布 % flag=find(ismember(FCF,allvalue+ X )==1)得到某气溶胶的数据记录 % count=length(flag); 即可得到记录次数 % ismember函数,用法: C=ismember(A,B); 判断A数组中的数是否等于B数组中的任意一个数,返回一个和A等大小的逻辑数组C % 例如:A=[1,2,3;4,5,6]; B=[2,8;6,9] % C=ismember(A,B); 则 C=[0,1,0;0,0,1]; class1=find(PDR>0 & PDR<=0.1); % 找出退偏比 >0 且 ≤0.1 的退偏比的数据记录 class2=find(PDR>0.1 & PDR<=0.2); % 找出退偏比 >0.1 且 ≤0.2 的退偏比的数据记录 class3=find(PDR>0.2 & PDR<=0.3); % 找出退偏比 >0.2 且 ≤0.3 的退偏比的数据记录 class4=find(PDR>0.3 & PDR<=0.4); % 找出退偏比 >0.3 且 ≤0.4 的退偏比的数据记录 class5=find(PDR>0.4 & PDR<=0.5); % 找出退偏比 >0.4 且 ≤0.5 的退偏比的数据记录 class6=find(PDR>0.5 & PDR<=0.6); % 找出退偏比 >0.5 且 ≤0.6 的退偏比的数据记录 class7=find(PDR>0.6 & PDR <=1.0); % 找出退偏比 >0.6 且 ≤1.0 的退偏比的数据记录 count1=count1+length(intersect(flag,class1)); % 求和,退偏比范围在0-0.1的数据记录次数,下面类似 count2=count2+length(intersect(flag,class2)); count3=count3+length(intersect(flag,class3)); count4=count4+length(intersect(flag,class4)); count5=count5+length(intersect(flag,class5)); count6=count6+length(intersect(flag,class6)); count7=count7+length(intersect(flag,class7)); end total=count1+count2+count3+count4+count5+count6+count7; % 数据记录总数 f1=count1/total; % 沙尘气溶胶退偏比范围为0-0.1的出现频率 f2=count2/total; f3=count3/total; f4=count4/total; f5=count5/total; f6=count6/total; f7=count7/total;
评论
    相关推荐
    • depolarization_ratio_frequency.rar
      用于所下载的calipso数据中每年的气溶胶退偏比频率统计
    • GaussDB_100_1.0.1-DATABASE-REDHAT-64bit.tar.gz
      guassdb100在redhat上安装包,单机部署的包,安装步骤请看我的文中介绍,经过大量实验搭建总结出来的文档
    • SIM800C_MQTT.rar
      使用SIM800C模块,使用MQTT协议,连接中国移动onenet平台,能实现数据的订阅、发布、存储等
    • 卷积神经网络
      这是卷积神经网络的一个实际用例,已经调试好了,能够在matlab上成功运行,适合从事卷积神经网络(CNN)研究的人员学习使用。
    • android从bootloader到launcher启动流程整理
      讲述android 开机流程 从boot rom---bootloader---init--zygote---systemserver---ams 并附上自己整理的每个流程的流程图 ,清晰熟悉android 启动流程
    • 分数阶混沌系统:分数阶混沌系统的数值解。-matlab开发
      该工具箱包含可用于模拟一些著名的分数阶混沌系统的函数,例如: - 陈的系统, - Arneodo的系统, - Genesio-Tesi 的系统, - 洛伦兹系统, - 牛顿-莱普尼克系统, - 罗斯勒的系统, - Lotka-Volterra系统, - 达芬的系统, - 范德波尔的振荡器, -伏打的系统- 陆氏系统, - 刘的系统, - Chua的系统, - 金融系统, - 3 细胞 CNN。 这些函数以数值方式计算描述混沌系统的分数阶非线性微分方程的解。 每个函数返回总模拟时间的状态轨迹(吸引器)。 更多详情请看书: Ivo Petras,分数阶非线性系统:建模、分析和仿真,Springer,系列:非线性物理科学,2011,ISBN 978-3-642-18100-9。 http://www.springer.com/engineering/control/book/978-3-
    • matlabpam代码-dsp-library:用于光通信的MATLABDSP函数库
      matlab pam代码 DSP库 该项目包含几种不同的功能,可将DSP算法应用于光通信。 这些功能适用于相干和非相干(PAM,DMT)光通信。 大多数功能是相互独立的。 因此,该代码中的功能可以轻松使用,并与其他DSP功能结合使用。 用法 功能列表以及简短说明在文件中。 输入和输出参数的描述在每个函数的标题中。 讯号 通常,输入(和输出)信号在第一维度上具有时间(例如,列向量),而第二维度用于一次管理多个信号(例如,不同的极化,不同的参数等)。 之所以选择这种约定,是因为MATLAB通过将列保留在内存的连续部分中来存储矩阵,因此,这种约定比其他方法(行向量)要快。 参数 大多数功能使用参数结构作为输入参数。 函数中使用的参数的描述通常在标头中,而coherent-dsp函数的默认参数在文件中。 参考 执照 此代码在下发布。
    • 有关多目标跟踪的PHD滤波的一些资料
      一些关于多目标跟踪的新的文献,主要是有关概率假设密度(PHD)的
    • 基于python开发的全国新工商采集工具 v1.2版本
      这个软件是通过scrapy爬虫框架结合代理IP池再加上request模拟请求技术以及验证码识别技术,可以做到日更新采集全国新工商信息。采集的数据自动存储在mysql数据库表里,可下载全量1.8亿多企业工商基本信息和36维度的详细信息.支持sql和excel导出数据包格式。
    • Aerosim Blockset
      The AeroSim aeronautical simulation blockset provides a complete set of tools for the rapid development of nonlinear 6-degree of freedom aircraft dynamic models. In addition to the basic aircraft dynamics blocks, the library also includes complete aircraft models which can be customized through parameter files.