matlab程序，size(A)返回矩阵维度信息。
nx=195;ny=153;sta=1551;nz=2;
output = zeros(sta,4);
for k=1:sta
IOUT = 0;
JOUT = 0;
dis=1000000.0;
for j=1:ny
for i=1:nx
tmp_dis=sqrt((XLAT_M(i,j)-lonlat(k,2))^2+(XLONG_M(i,j)-lonlat(k,1))^2);
if (tmp_dis <= dis)
dis = tmp_dis;
IOUT = i;
JOUT = j;
end
end
end
output(k,1)=IOUT;
output(k,2)=JOUT;
output(k,3)=XLAT_M(IOUT,JOUT);
output(k,4)=XLONG_M(IOUT,JOUT);
end
网友评论