In [107]:

    

import numpy as np
from glob import glob


    

In [108]:

    

files  = glob('*.vec')


    

In [109]:

    

d = np.loadtxt(files[0],skiprows=1,delimiter=',',dtype={'names': ('x', 'y', 'u', 'v'),'formats': ('f4', 'f4', 'f4', 'f4')},usecols=(0,1,2,3))


    

In [110]:

    

d


    

    
Out[110]:





array([(0.31248000264167786, -0.31248000264167786, 0.0, -0.0),
       (0.6249600052833557, -0.31248000264167786, 0.0, -0.0),
       (0.9374399781227112, -0.31248000264167786, 0.0, -0.0), ...,
       (19.061279296875, -19.68623924255371, 0.0, -0.0),
       (19.373760223388672, -19.68623924255371, 0.0, -0.0),
       (19.68623924255371, -19.68623924255371, 0.0, -0.0)], 
      dtype=[('x', '<f4'), ('y', '<f4'), ('u', '<f4'), ('v', '<f4')])

In [111]:

    

d['u']


    

    
Out[111]:





array([ 0.,  0.,  0., ...,  0.,  0.,  0.], dtype=float32)

In [112]:

    

quiver(d['x'],d['y'],d['u'],d['v'])


    

    
Out[112]:





<matplotlib.quiver.Quiver at 0x96ca310>






    

In [113]:

    

fid = open('Run000001.T000.D000.P000.H001.L.vec','r')


    

In [114]:

    

header = fid.readline()


    

In [115]:

    

fid.close()


    

In [116]:

    

header


    

    
Out[116]:





'TITLE="E:\\2CM_FP500_5%G_68K\\C001H001S0015CC\\Soapfilmone\\Analysis\\Run000001.T000.D000.P000.H001.L.vec" VARIABLES="X mm", "Y mm", "U m/s", "V m/s", "CHC", DATASETAUXDATA Application="PIV" DATASETAUXDATA SourceImageWidth="1024" DATASETAUXDATA SourceImageHeight="1024" DATASETAUXDATA MicrometersPerPixelX="19.530001" DATASETAUXDATA MicrometersPerPixelY="19.530001" DATASETAUXDATA LengthUnit="mm" DATASETAUXDATA OriginInImageX="0.000000" DATASETAUXDATA OriginInImageY="0.000000" DATASETAUXDATA MicrosecondsPerDeltaT="2000.000000" DATASETAUXDATA TimeUnit="ms" DATASETAUXDATA SecondaryPeakNumber="0" DATASETAUXDATA DewarpedImageSource="0" ZONE I=63, J=63, F=POINT\r\n'

In [128]:

    

header_list = header.replace(',',' ').replace('=',' ').replace('"',' ').split()


    

In [133]:

    

header_list.index('VARIABLES')


    

    
Out[133]:





2

In [157]:

    

header_list[3:12][::2]


    

    
Out[157]:





['X', 'Y', 'U', 'V', 'CHC']

In [158]:

    

header_list[4:12][::2]


    

    
Out[158]:





['mm', 'mm', 'm/s', 'm/s']

In [ ]: