将地球按经纬度划分成360*180的网格,现在每个交点的高程给出了,是一个360*180的矩阵,如何做高程分布图

将地球按经纬度划分成360*180的网格,现在每个交点的高程给出了,是一个360*180的矩阵,如何做高程分布图
我编了程序,但是无法实现：
>> h=xlsread('topogrd1.xls'); %读取高程数据
H=reshape(h,180,360); %将高程数据转换成360*180的矩阵
zz=9.8+9.5.*tanh((H+1200)./1632.5); %将高程数据做一个代数变换
xx=1:360; %生成经线
yy=1:180; %生成纬线
[x,y,z]=meshgrid(xx,yy,zz); %栅格化
surf(X,Y,zz)
Maximum variable size allowed by the program is exceeded.
有什么更好的实现方法吗?