时空谱分析

木市街8号

有谁会时空谱分析啊？求帮助，感激不尽

kongfeng0824

可以尝试小波分析看看。

木市街8号

kongfeng0824 发表于 2014-3-7 11:20

                               
登录/注册后可看大图

可以尝试小波分析看看。

哦，谢谢。主要想用到三维的时空谱分析（主要通过fft变换），获得纬向波数、经向波数和频率，小波应该实现不了吧

kongfeng0824

木市街8号 发表于 2014-3-13 09:13

                               
登录/注册后可看大图

哦，谢谢。主要想用到三维的时空谱分析（主要通过fft变换），获得纬向波数、经向波数和频率，小波应该实现 ...

我没用过快速傅里叶，用过小波分析

木市街8号

kongfeng0824 发表于 2014-3-13 09:28

                               
登录/注册后可看大图

我没用过快速傅里叶，用过小波分析

那小波分析可以处理这种三维的，获得这三个参数吗？请问

木市街8号

木市街8号 发表于 2014-3-13 10:26

                               
登录/注册后可看大图

那小波分析可以处理这种三维的，获得这三个参数吗？请问

能否发个例子给我分享下，谢谢了

kongfeng0824

木市街8号 发表于 2014-3-13 10:27

                               
登录/注册后可看大图

能否发个例子给我分享下，谢谢了

clc,clear; %清屏，清空
ncid = netcdf.open('D:\precip.mon.1981-2010.ltm.nc','NOWRITE'); %打开nc文件
ncdisp('precip.mon.1981-2010.ltm.nc'); %在命令窗中显示nc文件的变量
% [ndims,nvars,ngatts,unlimdimid]=netcdf.inq(ncid)
PrecipData  = ncread('precip.mon.1981-2010.ltm.nc','precip'); %读入变量precip
TimeData  = ncread('precip.mon.1981-2010.ltm.nc','time'); %读入变量time
LonData  = ncread('precip.mon.1981-2010.ltm.nc','lon'); %读入变量lon
LatData  = ncread('precip.mon.1981-2010.ltm.nc','lat'); %读入变量lat
Valid_yr_countData  = ncread('precip.mon.1981-2010.ltm.nc','valid_yr_count'); %读入变量validprecip_yr_count
ClimatologyData  = ncread('precip.mon.1981-2010.ltm.nc','climatology_bounds'); %读入变量climatology_bounds



[X, Y] = meshgrid(LatData,LonData);
contourf(Y-180.0,X,PrecipData(:,:,9)); %画9月份等值面图
shading flat; %去掉等值线
colorbar('SouthOutside','Position',[0.142,0.03,0.75,0.04]); %添加颜色条 [左右,上下,长，宽]


load coast %加载全球海岸线，但是不显示出来
geoshow(lat,long,'DisplayType','line','Color','White'); %显示出海岸线，lat和long是coast的属性
hold on;


set(gca,'LineWidth',1,'FontSize',10,'Ylim',[-90,90],'Xlim',[-180,180],'Position',[0.142,0.105,0.75,0.84]...
   ,'XTick',[-180:60:180],'XTicklabel',{'-180W','-120W','-60W','0','60E','120E','180E'},'XGrid','on'...
   ,'YTick',[-90:30:90],'YTicklabel',{'-90S','-60S','-30S','0','30N','60N','90N'}, 'YGrid','on');  %添加经纬度信息
xlabel({'经度'}); %添加x轴名称
ylabel({'纬度'}); %添加y轴名称
title(['\fontsize{16}precip.mon.1981-2010.ltm']);  %添加地图题目
annotation('textbox',[0.838711566617862,0.045,0.0431918008784773,0.0579710144927536],'String',{'Unit(mm/d)'},'EdgeColor','none'); %添加文本框
hold off;

木市街8号

kongfeng0824 发表于 2014-3-13 11:12

                               
登录/注册后可看大图

clc,clear; %清屏，清空
ncid = netcdf.open('D:\precip.mon.1981-2010.ltm.nc','NOWRITE'); %打开nc文件 ...

谢谢，可是这里面好像没用到具体的小波分析啊？

kongfeng0824

木市街8号 发表于 2014-3-13 11:37

                               
登录/注册后可看大图

谢谢，可是这里面好像没用到具体的小波分析啊？

这个只是维度的处理。没有小波分析。ctw是内置的小波分析

maiyatang_xihui

画图画图？？