Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  849 1661 (N26E31 L=216)   379     86        0 ( 805/1634- 888/1680)   0.0 E34-E27 N24-N27
      2  336 1398 (S02E56 L=191)   458     72        1 ( 309/1371- 365/1429)  56.3 E58-E52 S03-N00
      3  655 1358 (N03E30 L=217)   271     59        1 ( 636/1331- 677/1387) -36.9 E32-E28 N01-N04
      4 1638 1070 (N14W35 L=282)  2366    488        8 (1580/ 863-1678/1334) -55.8 W25-W42 N03-N28
      5  150 1232 (S15E70 L=177)   762    108        2 ( 126/1201- 172/1268)  76.0 E72-E66 S16-S12
      6 1246 1166 (N07W08 L=255)   659    197        1 (1173/1096-1342/1238) -13.0 W02-W15 N05-N09
      7  251 1085 (S22E53 L=194)   381    177        1 ( 205/1060- 301/1121)  -5.7 E58-E48 S22-S20
      8  578 1094 (S14E28 L=219)   609    189        2 ( 504/1061- 664/1121)  10.3 E32-E21 S16-S11
      9  423 1069 (S19E38 L=209)   245    106        0 ( 384/1039- 468/1092)  45.0 E40-E34 S21-S16
     10  180  809 (S40E58 L=190)  1403    383        1 ( 125/ 751- 244/ 851)  28.3 E63-E50 S44-S36
     11 1044  769 (S20W09 L=256)   255     90        1 (1002/ 757-1098/ 786)  18.4 W05-W11 S20-S18
     12  537  597 (S45E21 L=226)   576    149        0 ( 473/ 580- 594/ 624)  24.4 E27-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