matlab 读写shp

### 如何在 MATLAB 中读取和写入 Shapefile 文件

MATLAB 提供了内置函数来处理地理空间数据，包括读取和写入 Shapefile 文件。这些功能主要通过 Mapping Toolbox 实现。

#### 使用 `shaperead` 函数读取 Shapefile 数据
为了读取 Shapefile 文件中的几何对象及其属性，可以使用 `shaperead` 函数[^1]：

S = shaperead('filename.shp');

此命令会返回一个结构体数组 S，其中每个元素代表 shapefile 中的一个要素，并包含其几何坐标和其他属性字段的信息。

#### 使用 `mapwrite` 和 `shapeinfo` 函数创建新的 Shapefiles 或更新现有文件
对于写操作，可以通过组合使用 `mapwrite` 和 `shapeinfo` 来实现。先定义要保存的数据格式，再调用 `mapwrite` 将数据写出到指定路径下的新 shapefile 文件中：

% 定义输出文件名及路径
outputFile = 'new_shape_file.shp';

% 创建一个新的 shapefile 并设置投影信息 (如果适用)
if exist(outputFile,'file')
    delete(outputFile);
end
fid = mapwrite(outputFile);

% 添加图层并填充数据...
for i=1:length(S)
    % 假设我们有一个名为 "geometry" 的字段存储着 geometry 对象，
    % 可能还有其他一些附加属性列。
    info = shapeinfo;
    info.Geometry = S(i).Geometry;  % 设置几何形状
    
    % 如果有额外的属性，则继续添加它们
    setfield(info, 'Attribute', value); 

    fwrite(fid, info);
end

fclose(fid);

上述代码片段展示了如何遍历已有的 shape 结构体并将之逐条记录至目标文件内；需要注意的是实际应用时应根据具体需求调整字段名称与值。