Changeset 21349


Ignore:
Timestamp:
11/07/16 09:58:46 (8 years ago)
Author:
dfelikson
Message:

BUG: changed so that exps from line and point shapefiles are not closed

Location:
issm/trunk-jpl/src
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/shp/shp2exp.m

    r17542 r21349  
    3333                        expstruct(end).nods=length(x);
    3434                        expstruct(end).density=1;
    35                         expstruct(end).closed=1;
     35                        expstruct(end).closed=0;
    3636                        %exp(end).name=num2str(shp(i).id);
    3737                elseif strcmpi(shp(i).Geometry,'Line'),
    3838                        x=shp(i).X; y=shp(i).Y;
    39                         x(end)=x(1); y(end)=y(1);
    4039                        expstruct(end+1).x=x;
    4140                        expstruct(end).y=y;
    4241                        expstruct(end).nods=length(x);
    4342                        expstruct(end).density=1;
    44                         expstruct(end).closed=1;
     43                        expstruct(end).closed=0;
    4544                end
    4645        end

  • issm/trunk-jpl/src/m/shp/shp2exp.py

    r19423 r21349  
    2626
    2727        shp=shapefile.Reader(shapefilename)
    28         expdict=dict(closed=1,density=1)
     28        expdict=dict(density=1)
    2929
    3030        x=[]
    … …  
    3333                geom=shp.shapes()[i].shapeType
    3434                if geom==5: # polygon
     35         expdict['closed']=1
    3536                        tmpx=[p[0] for p in shp.shapes()[i].points]
    3637                        tmpy=[q[1] for q in shp.shapes()[i].points]
    … …  
    3839                        y.append(tmpy)
    3940                elif geom==3: # line
     41         expdict['closed']=0
    4042                        tmpx=[p[0] for p in shp.shapes()[i].points]
    4143                        tmpy=[q[1] for q in shp.shapes()[i].points]
    … …  
    4345                        y.append(tmpy)
    4446                elif geom==1: # point
     47         expdict['closed']=0
    4548                        x.append(shp.shapes()[i].points[0][0])
    4649                        y.append(shp.shapes()[i].points[0][1])

  • issm/trunk-jpl/src/py3/shp/shp2exp.py

    r19895 r21349  
    2626
    2727        shp=shapefile.Reader(shapefilename)
    28         expdict=dict(closed=1,density=1)
     28        expdict=dict(density=1)
    2929
    3030        x=[]
    … …  
    3333                geom=shp.shapes()[i].shapeType
    3434                if geom==5: # polygon
     35         expdict['closed']=1
    3536                        tmpx=[p[0] for p in shp.shapes()[i].points]
    3637                        tmpy=[q[1] for q in shp.shapes()[i].points]
    … …  
    3839                        y.append(tmpy)
    3940                elif geom==3: # line
     41         expdict['closed']=0
    4042                        tmpx=[p[0] for p in shp.shapes()[i].points]
    4143                        tmpy=[q[1] for q in shp.shapes()[i].points]
    … …  
    4345                        y.append(tmpy)
    4446                elif geom==1: # point
     47         expdict['closed']=0
    4548                        x.append(shp.shapes()[i].points[0][0])
    4649                        y.append(shp.shapes()[i].points[0][1])
Note: See TracChangeset for help on using the changeset viewer.