matlab提取FY4数据CLM

clc;
close all;
load 'G:\theone\GeoRef';
load 'G:\theone\latlon_china_4000_grid';
load 'G:\theone\latlon_china_4000_FY';

FilePath='Z:\龚雨婷\FY'; %待处理数据所存放的文件夹
files=dir(fullfile(FilePath,'*.hdf'));
% files=dir(fullfile(FilePath,'*.NC'));
FilesCount=length(files);
% [a,R]=geotiffread('C:\Users\Administrator\Desktop\data_fusion\30x30\2017_111_new\gapfill1+xiaozu\clm_0429_0300.tif');
% info=geotiffinfo('C:\Users\Administrator\Desktop\data_fusion\30x30\2017_111_new\gapfill1+xiaozu\clm_0429_0300.tif');
for i=1:FilesCount
filename = fullfile(FilePath,files(i).name);

datasetname = '/CLM';      %待提取的字段
    CLM = h5read(filename,datasetname);
    CLM = uint16(CLM);
    CLM=CLM'; % 0表示确定是云
    China_lst = CLM(185:935,555:2180);
%     China_cloud
%     sizeOfA=size(China_cloud);
%     R=sizeOfA(1);% 矩阵有多少行
%     L=sizeOfA(2);% 矩阵有多少列
%     data=zeros(R,L);
%     for k=1:R
%        for j=1:L
%            if China_cloud(k,j)>=3 %取绝对晴空 3
%                data(k,j)=1;
%            else              %都是有云的情况
%                data(k,j)=0;
%            end
%        end
%     end
%     res=data.*double(China_lst); %提取

China_lst = double(reshape(China_lst,751*1626,1));
Ref_values = griddata(FY_lon4,FY_lat4,China_lst,grid4_lon,grid4_lat,'nearest');% 1417500*1 'linear' 基于三角形的线性插补法(缺省),'cubic' 基于三角形的三次插补法,'nearest' 最近邻插补法
Ref_values = reshape(Ref_values,900,1575);
Ref_values(Ref_values==0)=255;
Ref_values(Ref_values~=255)=0;
geotiffwrite(['G:\theone\yangzriver\2019\FY\geotif\',files(i).name,'.tif'],Ref_values,GeoRef); %后缀必须是tif
% geotiffwrite(fullfile(output_dir,files(i).name,'.tif'),China_lst,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);

end

matlab提取FY4数据CLM相关推荐

利用matlab对rosbag数据,通过matlab提取rosbag数据
通过matlab提取rosbag数据代码如下: 最近在做AGV导航实验,导航小车是基于ROS实现的,判断小车导航是否运行稳定,可以直接通过rviz直观看出来.但是如何通过数据定性地分析数据,这就需要 ...
使用Matlab提取ADC数据占空比变化的波形
提取变化占空比的规则如下: 1. 如果占空比发生 >=5% 的变化,则认为占空比发生变化: 2. 正常的1个周期100个,其中连续 <0.1V 变化的点约为40个左右: 3. 那么,40* ...
Matlab提取矩阵数据
首先,Matlab中,()的地位跟Python中[]的地位有点类似. 特别是在提取矩阵数据中的情况在演示之前,先建立矩阵,并看一下 M = [1,2; 3,4; 5,6] 输出的结果是: M =1 ...
matlab提取矩阵数据：行、列、块
#2020/4/17 @matlab 矩阵行列和矩阵块的提取: 先随意设定一个6维的随机数矩阵取矩阵a第2行 a(2,:) 取矩阵a第3列 a(:,3) 取矩阵a第2至第3列取矩阵a第2至第3行 ...
使用MATLAB快速提取图片数据
今天介绍一下如何使用MATLAB提取图片数据.在MATLAB中其实也有很多好用的附加功能可免费使用,而今天介绍的这个就是小编在附加功能中发现所发现,一起来了解下吧! 1.GRABIT下载 MATLAB ...
matlab提取数据的一部分,matlab如何提取数组中的满足一定范围的一段数据
给定一个数组,如何让matlab生成一个这个数组中的一个随机数? x=[102030];x(randi(length(x)));其中randi(length(x))生成从1~(x的长度)这几个自然数中 ...
matlab输出多个参数到表格固定列,matlab将excle的列提取出来?如何使用Matlab提取Excel中指定列的数据？...
MATLAB 在txt文件中怎么提取数据,并将那三列数据转化为EXCLE格式? 在matlab中 file import data number of text header lines 改成6 下一 ...
matlab提取数据的一部分,matlab处理excel数据【怎么用MATLAB从excel中提取部分数据】...
已经把excel的一列导入到matlab的工作区,怎么计算这一列数的均值.标准差.最小值.最大值以及中位数? 把这个EXCEL放到matlab目录下的work里,假设名为1.xls,然后在matlab ...
使用Matlab提取ADC采样数据中的噪声
正常波形与噪声波形的区别正常方波信号的特征:一个周期内的数据按时间单调变化,则仅有2种情况:① 数据在上升或下降边沿时,其两边的数据则是一边大一边小(如下图中C和D点):② 数据处于平稳时,则两边大 ...

matlab提取FY4数据CLM

matlab提取FY4数据CLM相关推荐

最新文章

热门文章