1、准备60阶的球谐系数格网(以GRACE CSR RL06 Level-2数据为例),有两种方式得到球谐系数:
①直接从机构数据读取,经过一阶项、二阶项和取平均等预处理操作后的61*121的球谐系数sc矩阵;
②例如经过高斯300km之后的等效水柱高,再利用gmt_grid2cs函数转化为球谐系数
2、利用pcolor函数绘制球谐系数格网图
下面展示第二种方式绘制格网图的代码:
clearvars -except cs300;
close all;
load E:\Data\result\grid\grid300_CSR.mat;
num=60;lmax=60;
cs=gmt_grid2cs(grid300_CSR(:,:,1:num),lmax);
for ii=1:num
temp(:,:)=cs(ii,:,:);
cs300(ii,:,:)=gmt_mc2gc(temp);
end
% save cs300.mat cs300
sc=gmt_cs2sc(reshape(cs300(8,:,:),lmax+1,lmax+1));sc(sc==0)=NaN;
m1=-60:60;l1=0:60;
[m,l]=meshgrid(m1,l1);
figure(1)
pcolor(m,l,sc./1e-12);shading flat;
caxis([-10,10]);colormap('jet');grid off;
xlabel('Order(m)','Fontname','Time New Roman',"FontSize",20);
ylabel('Degree(l)','Fontname','Time New Roman',"FontSize",20);
set(gca, 'YDir', 'reverse')
xticks([-60,-30,0,30,60]);yticks([0,30,60]);
title('Spherical harmonic coefficient','FontName','Time New Roman','Fontsize',20,'fontweight','bold')
h=colorbar('v','FontSize',15,'fontweight','bold');
set(get(h,'title'),'string','$10\times^{-12}$', 'Interpreter', 'latex','FontSize',10,'fontweight','bold');