cycling/shape2struct.m

   1 function [ perim_struct ] = shape2struct( prefix )
   2 %[ perim_struct ] = shape2struct( prefix )
   3 % inputs:
   4 %   prefix - string, path to directory with shapefiles
   5 % output:
   6 %   perim_struct - structure with perimeter information
   7
   8 % to do:
   9 % put names in sthe struct file as in UT-SLD-HU2S Patch Springs 8-20-2013
  10 % 2150
  11
  12 d=dir([prefix,'*.shp'])
  13 d={d.name};
  14 if(isempty(d)), error(['No files found for ',prefix]),end
  15
  16 % order the files in time
  17 nfiles=length(d);
  18 t=zeros(1,nfiles);
  19 date_format = 'yyyymmdd_HHMM';
  20 %date_format = 'mm-dd-yyyy_HHMM';
  21 %perim times are local, need to convert to UTC
  22 zone_shift = -8;
  23
  24
  25 for i=1:nfiles
  26     f{i}=[prefix,d{i}];
  27     date = d{i}(end-21:end-9);
  28     %date = d{i}(end-18:end-4);
  29     %t(i) = datenum(date,date_format)+zone_shift/24;
  30     t(i) = datenum(date,date_format);
  31 end
  32 [t,i]=sort(t);
  33 p.file={d{i}};
  34 p.time=t;
  35
  36 %combine into struct
  37 for i = 1:length(p.file)
  38     perim_struct(i) = shaperead(f{i});
  39     %perim_struct(i).Lon = perim_struct(i).X
  40     %perim_struct(i).Lat = perim_struct(i).Y
  41 end
  42
  43 for i = 1:length(p.file)
  44     %perim_struct(i) = shaperead(f{i})
  45     perim_struct(i).Lon = perim_struct(i).X;
  46     perim_struct(i).Lat = perim_struct(i).Y;
  47     perim_struct(i).Name = replace(p.file{i},'_',' ');
  48     perim_struct(i).Name = perim_struct(i).Name(1:end-9);
  49 end
  50
  51 % perim_struct.Lon = perim_struct.X
  52 % perim_struct.Lat = perim_struct.Y
  53 end
  54