etopo地形数据转换为xyz格式

 clear all;
ncload（‘etopo2.nc‘）;
sum=0;
for i=8377:9421    
a=lon（i1）;
for j=2728:3931
b=lat（j1）;
c=topo（ji）;
d=[a b c];
aaa=size（d）;
pp（sum+1:sum+aaa:）=d;
sum=sum+aaa;
end 
end
save etopo.dat -ascii pp
b=find（pp（:1）>111.48 & pp（:1）<112.26355 &pp（:2）>21.085217 &pp（:2）<24.14）;
pp（b:）=[];
c=find（pp（:1）>112.26355 & pp（:1）<113.0711608 &pp（:2）>21.496067 &pp（:2）<24.14）;
pp（c:）=[];
d=find（pp（:1）>113.0711608 & pp（:1）<114.4086373&pp（:2）>21.76037 &pp（:2）<24.14）;