MATLAB TIFF转Shape、TIFF和Shape的读写

• 读TIFF

[A,R] = geotiffread(‘Filename.tif’);
info = geotiffinfo(‘Filename.tif’);

• 写TIFF

geotiffwrite(‘Filename’,A,R,‘GeoKeyDirectoryTag’,info.GeoTIFFTags.GeoKeyDirectoryTag)

• 读Shape

A = shaperead(‘Filename.shp’);

• 写Shape

shapewrite(Map,‘Filename’); % 这里的Map为结构体，且需带有数据类型（点线面）
% 若要将数组A写为Shape
A = num2cell(A);
fields = {‘X’,‘Y’,‘Z’};
A = cell2struct(A,fields,dim); % fields为N×1时，dim=1；当fields为1×M时，dim=2
[A.Geometry] = deal(‘Point’); % 以点为例
shapewrite(A,‘Filename’);

• Tiff2Shp

[A,R,~] = geotiffread(‘Filename.tif’);
[I,J] = ind2sub(size(A),find(A>0)); % pixel value >0
Map_XY