Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1547 1532 (N18W49 L= 56)  1135    224        1 (1452/1507-1639/1557) -30.3 W42-W55 N13-N21
      2 1817 1414 (N02W68 L= 75)   953    254        2 (1778/1379-1863/1447)  -3.2 W62-W80 N01-N03
      3 1616 1170 (S03W39 L= 46)   788    175        3 (1582/1086-1663/1258)  77.5 W36-W42 S08-N01
      4 1782 1102 (S12W50 L= 57)   382     93        1 (1749/1063-1802/1139)  53.1 W46-W52 S14-S10
      5  160  967 (N22E63 L=304)   822    144        3 ( 127/ 940- 190/ 997)   7.1 E66-E58 N20-N24
      6  388  960 (N18E41 L=326)   261    108        1 ( 353/ 942- 427/ 987)  -9.5 E43-E37 N16-N20
      7 1539  745 (S25W23 L= 30)   241     76        0 (1526/ 713-1550/ 790) -90.0 W21-W24 S26-S21
      8 1692  705 (S33W34 L= 41)   477    123        0 (1661/ 645-1723/ 772)  45.0 W31-W37 S35-S29
      9  210  716 (N07E66 L=300)  1242    217        1 ( 163/ 680- 264/ 750) -26.6 E71-E61 N04-N09
     10  428  363 (S18E63 L=303)   572    107        1 ( 410/ 331- 452/ 381)  63.4 E66-E60 S20-S16
  

    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