批量读取不规则数组的shp文件经纬度

卑微江小澄

S1 = shaperead('F:\毕设\shp\中国区域1~5级水系图\我国二、三级河流.shp');%读入shp文件

c=1;
for i = 1:2143 %2143为S1中总区域个数
    a = size(S1(i).X,2); %S1中每个区域含坐标个数
    b = S1(i).X;%X为S1中储存经度的变量，Y同理
    d = S1(i).Y;
    k = 1;
    for j = c:c+a-1
      lon(j,1) = b(1,k);%所有区域的经度
      lat(j,1) = d(1,k);
      k = k+1;
    end
    c = c+a;
end