source: issm/branches/trunk-larour-NatGeoScience2016/src/m/exp/exp2shp.m@ 21759

function exp2shp(expfilename,shpfilename,geometry)
%SHPWRITE - write a shape file from a contour structure
%
%   Usage:
%      exp2shp(expfilename,shpfilename,geometry)
%
%   Example:
%      exp2shp('domainoutline.exp','domainoutline.shp')
%      exp2shp('domainoutline.exp','domainoutline.shp','Polygon')
%      exp2shp('massfluxgate.exp','massfluxgate.shp','Line')
%
%   See also SHPREAD, SHPWRITE, SHP2EXP

%check file extensions
[pathstr,name,ext] = fileparts(shpfilename);
if ~strcmp(ext,'.shp'),
        error(['Shapefile ' shpfilename ' does not have an extension .shp']);
end

[pathstr,name,ext] = fileparts(expfilename);
if ~strcmp(ext,'.exp'),
        error(['Exp file ' expfilename ' does not have an extension .exp']);
end

shp=expread(expfilename);

%initialize number of profile
count=0;

contours=struct([]);
for i=1:length(shp),
        if nargin < 3
                if length(shp(1).x) == 1
                        geometry = 'Point';
                elseif length(shp(1).x) < 3
                        geometry = 'Line';
                else 
                        if (shp(1).x(end)==shp(1).x(1) && shp(1).y(end)==shp(1).y(1)),
                                geometry = 'Polygon';
                        else
                                geometry = 'Line';
                        end
                end
        end
        contours(i).Geometry=geometry;
        contours(i).id=i;
        contours(i).X=shp(i).x;
        contours(i).Y=shp(i).y;
end
        
shapewrite(contours,shpfilename);