MATLAB循环读取文件

一夜荒原

我选择某几个省份 想通过for循环进行读取shp文件，然后是先将地名放到cell里边 然后进行循环读取，但在括号里，一直显示错误？请问大佬们该如何解决 谢谢

灰色节能君

sta=p{i}%这里用大括号
infoL=shapeinfo(['chinamap\province\',sta,'.shp']);%字符串可以这样连

一夜荒原

只能输出第一个省份的

一夜荒原

灰色节能君 发表于 2020-3-12 15:14
sta=p{i}%这里用大括号
infoL=shapeinfo(['chinamap\province\',sta,'.shp']);%字符串可以这样连

只能输出第一个省份的
警告: Unsupported shape type PolyLineZ (type code = 13).
> In openShapeFiles>readHeaderTypeCode (line 147)
  In openShapeFiles (line 23)
  In shapeinfo (line 73)
  In NE_China (line 8)
错误使用 openShapeFiles>readHeaderTypeCode (line 145)
Unsupported shape type PolyLineZ (type code = 13).

出错 openShapeFiles (line 23)
headerTypeCode = readHeaderTypeCode(shpFileId,callingFcn);

出错 shaperead (line 210)
    = openShapeFiles(filename,'shaperead');

出错 NE_China (line 9)
    proL=shaperead(['chinamap\province\',p{i},'.shp']);

一夜荒原

灰色节能君 发表于 2020-3-12 15:14
sta=p{i}%这里用大括号
infoL=shapeinfo(['chinamap\province\',sta,'.shp']);%字符串可以这样连

。。。吉林省的文件有问题

一夜荒原

灰色节能君 发表于 2020-3-12 15:14
sta=p{i}%这里用大括号
infoL=shapeinfo(['chinamap\province\',sta,'.shp']);%字符串可以这样连

我发现一个问题  ，，只要是接壤的省份就读不出来 不接壤的就可以  比如可以读出来 辽宁和黑龙江 。。这是为什么？

灰色节能君

一夜荒原 发表于 2020-3-12 16:53
我发现一个问题  ，，只要是接壤的省份就读不出来 不接壤的就可以  比如可以读出来 辽宁和黑龙江 。。这 ...

不太清楚是什么问题，不过我用mmap画中国省份是这样画的，你参考一下吧
infoL=shapeinfo('D:\matlab123\bin\bou2_4l.shp');
infoP=shapeinfo('D:\matlab123\bin\bou2_4p.shp');
ChinaL=shaperead('D:\matlab123\bin\bou2_4l.shp');
ChinaP=shaperead('D:\matlab123\bin\bou2_4p.shp');
bou2_4lx=[ChinaL(:).X];
bou2_4ly=[ChinaL(:).Y];
bou2_4px=[ChinaP(:).X];
bou2_4py=[ChinaP(:).Y];
m_proj('miller','lon',LONLIMS,'lat',LATLIMS);
m_plot(bou2_4lx,bou2_4ly,'k');

一夜荒原

灰色节能君 发表于 2020-3-12 17:02
不太清楚是什么问题，不过我用mmap画中国省份是这样画的，你参考一下吧
infoL=shapeinfo('D:\matlab123\bi ...

我知道这个 不过我只需要一部分 ，我发现可能是因为我分省的shp文件有问题。

灰色节能君

刚刚看到你在找各省的边界文件，不知道解决了没有，我试了一下，单独画个别省也是可以画的啊
p=cell(1,3);
p{1}='Jiangsu_province';
p{2}='Zhejiang_province';
p{3}='Anhui_Province';
for i=1:3
infoL=shapeinfo(['D:\',p{i},'.shp']);
infoP=shapeinfo(['D:\',p{i},'.shp']);
ChinaL=shaperead(['D:\',p{i},'.shp']);
ChinaP=shaperead(['D:\',p{i},'.shp']);
bou2_4lx=[ChinaL(:).X];
bou2_4ly=[ChinaL(:).Y];
bou2_4px=[ChinaP(:).X];
bou2_4py=[ChinaP(:).Y];
m_plot(bou2_4lx,bou2_4ly,'k');
hold on
end

一夜荒原

灰色节能君 发表于 2020-3-12 17:54
刚刚看到你在找各省的边界文件，不知道解决了没有，我试了一下，单独画个别省也是可以画的啊
p=cell(1,3); ...

没解决 你这个是只用那个bou_l2文件画的吗