在我们阅读文献的时候,经常出现如下类型公式:
上述公式出自 Chen J L, Wilson C R, Seo K W. Optimized smoothing of Gravity Recovery and Climate Experiment (GRACE) time‐variable gravity observations[J]. Journal of Geophysical Research: Solid Earth, 2006, 111(B6). 中的公式(6)在JGR: Solid Earth本篇文章网站显示所有引用这篇文章的学者论文
大多数文献中对于该公式的详细解释和具体操作是缺乏的,这里我想试着根据自己的理解大致实现这一公式。如有问题欢迎各位友友大胆批评指正!
1、数据准备
(1)陆地和海洋的掩膜数据(UTCSR网站:GRACE/GRACE-FO - Gravity Recovery and Climate Experiment (utexas.edu))
(2)GLDAS_NOAH025_M.A.021数据从2002年4月至2009年12月,缺失月份和GRACE对应上
去除2004.01-2009.12的均值,代码中“gldas02_09.mat”即为该格网数据
(3)GRACE CSR RL06 2002.04-2009.12去平均,不经过任何滤波的质量数据(cs_mss02_09.mat),经gmt_cs2grid即可转为格网EWH
2、代码展示
clearvars -except gldas02_09;
addpath E:\Data\CSR\result\RL06;
load E:\Code\GLDAS\GLDAS2006m025\gldas02_09.mat
%%如何计算误差项,是除去季节项和趋势项之后剩余的部分么?
load cs_mss02_09.mat;cs_mass(:,:)=cs_mss02_09(17,:,:);%以2003年11月份数据为例
grid=gmt_cs2grid(cs_mass,0,0.25);
file_land='E:\Data\CSR\Mascon\data\CSR_GRACE_GRACE-FO_RL06_Mascons_v02_LandMask.nc';
file_ocean='E:\Data\CSR\Mascon\data\CSR_GRACE_GRACE-FO_RL06_Mascons_v02_OceanMask.nc';
% info = ncinfo('E:\Data\CSR\Mascon\data\CSR_GRACE_GRACE-FO_RL06_Mascons_v02_OceanMask.nc');
lon=ncread(file_land,'lon');lat=ncread(file_land,'lat');
L_val=ncread(file_land,'LO_val');L_msk=flipud(L_val');
O_val=ncread(file_ocean,'LO_val');O_msk=flipud(O_val');
grid_weight=zeros(720,1440);
for j=1:720
for i=1:1440
grid_weight(j,i)=cos(lat(j)*pi/180);
end
end
land_weight=grid_weight.*L_msk;ocean_weight=grid_weight.*O_msk;
l_weight_sum=sum(sum(land_weight));o_weight_sum=sum(sum(ocean_weight));
temp1(:,:) = grid(:,:).*100;temp2(:,:) = gldas02_09(:,:,17)*0.1;
RMS_L=sqrt(sum(sum(((temp1-temp2).*L_msk).^2))./sum(sum(L_msk)));
RMS_O=sqrt(sum(sum(((temp1-temp2).*O_msk).^2))./sum(sum(O_msk)));
SNR=RMS_L./RMS_O;
SNR_log=10*log10(SNR);
3、结果与讨论
计算得到GRACE无约束解的RMS比率为1.5755;后期实验高斯300滤波的RMS比率为3.5005、Swenson滤波为2.4583。结果符合我们的认知,经过滤波处理后RMS比率有所提升,说明滤波处理能提高信噪比。但是上述代码仍有缺陷,没有完全与RMS比率公式匹配上,缺乏Err项的计算。对于Err误差项如何计算,我目前还未找到合适的方法,如果有明白其中细节的,欢迎私信交流(q:2629512206)
如果对你有帮助,多多转发多多集赞哦。Thanks!
参考文献:
[1] Chen J L, Wilson C R, Seo K W. Optimized smoothing of Gravity Recovery and Climate Experiment (GRACE) time‐variable gravity observations[J]. Journal of Geophysical Research: Solid Earth, 2006, 111(B6)
[2] 魏征强. GRACE时变重力数据球谐域滤波方法与质量分布反演应用研究[D]. 中国矿业大学, 2023. DOI:10.27623/d.cnki.gzkyu.2023.000516