Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  855 1660 (N26E30 L=216)   316     81        1 ( 819/1633- 892/1681)  -9.5 E33-E27 N24-N27
      2  339 1395 (S02E55 L=192)   380     58        1 ( 313/1372- 365/1429)  56.3 E57-E52 S03-N00
      3 1642 1067 (N14W35 L=282)  2010    502        7 (1582/ 863-1679/1334) -55.8 W25-W42 N03-N28
      4  151 1229 (S15E70 L=177)   676     89        1 ( 128/1197- 174/1261) -76.0 E72-E66 S16-S12
      5 1253 1161 (N07W09 L=256)   551    190        1 (1179/1093-1346/1230) -14.0 W03-W15 N05-N09
      6  303 1110 (S19E49 L=198)   335    112        1 ( 270/1090- 329/1143)   0.0 E51-E46 S20-S17
      7  577 1093 (S14E28 L=219)   507    179        0 ( 496/1069- 649/1115)  21.8 E32-E22 S17-S11
      8  181  806 (S40E57 L=189)  1374    391        0 ( 124/ 750- 246/ 851)  28.3 E63-E50 S44-S36
      9  547  595 (S45E20 L=227)   478    127        0 ( 496/ 576- 596/ 623)  32.0 E24-E16 S47-S41
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y