source: issm/trunk-jpl/src/m/shp/Shp2Exp.m@ 16852

function Shp2Exp(expfilename,shapefilename)
%SHP2EXP- transform shape file to Argus .exp file
%
%   Usage:
%      Shp2Exp(expfilename,shapefilename);
%
%   Example:
%      Shp2Exp('Domain.exp','Domain.shp');
%
%   See also EXPMASTER, EXPDOC

        if ~exist(shapefilename,'file'),
                error(['Shapefile ' shapefilename ' does not exist']);
        end
        shp=shaperead(shapefilename);

        exp=struct([]);
        for i=1:length(shp),
                if strcmpi(shp(i).Geometry,'Polygon'),
                        x=shp(i).X; y=shp(i).Y;
                        ids=find(isnan(x));
                        x(ids)=[]; y(ids)=[];

                        exp(end+1).x=x;
                        exp(end).y=y;
                        exp(end).nods=length(x);
                        exp(end).density=1;
                        exp(end).closed=1;
                        exp(end).name=num2str(shp(i).id);
                elseif strcmpi(shp(i).Geometry,'Point'),
                                x=shp(i).X; y=shp(i).Y;

                                exp(end+1).x=x;
                                exp(end).y=y;
                                exp(end).nods=length(x);
                                exp(end).density=1;
                                exp(end).closed=1;
                                %exp(end).name=num2str(shp(i).id);
                end
        end

        expwrite(exp,expfilename);