Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  687 1824 (N40E56 L= 59)  2167    150        1 ( 648/1763- 733/1890) -14.8 E79-E45 N35-N45
      2  927 1652 (N33E24 L= 92)   362    110        0 ( 902/1611- 967/1699)  90.0 E24-E21 N29-N37
      3  382 1555 (N15E60 L= 56)  3813    436       12 ( 265/1446- 516/1647)  48.8 E67-E50 N06-N23
      4  659 1570 (N21E38 L= 78)   257     55        0 ( 638/1549- 694/1590) -21.8 E39-E34 N20-N22
      5  254 1302 (S03E60 L= 56)   774    162        1 ( 209/1244- 284/1374)  90.0 E63-E57 S06-N01
      6 1054 1054 (S00W01 L=117)  1242    430        3 (1010/ 841-1145/1212) -54.2 E04-W10 S10-N08
      7 1826 1070 (N19W55 L=170)  1882    274        1 (1743/ 956-1888/1175) -25.3 W45-W64 N13-N23
      8  536 1097 (S09E31 L= 85)   439    148        0 ( 477/1074- 596/1146)  45.0 E34-E27 S11-S04
      9 1009  799 (S16W04 L=120)   195     87        0 ( 974/ 775-1054/ 820)  36.9 W02-W06 S17-S13
     10 1616  592 (S13W50 L=166)  1517    264        2 (1496/ 541-1733/ 639)  34.5 W42-W58 S18-S07
     11 1602  396 (S25W63 L=179)  2339    181        5 (1514/ 365-1699/ 426)  12.8 W53-W75 S27-S21
     12  960  356 (S46W16 L=132)   713     88        2 ( 918/ 330-1008/ 402)  14.0 W12-W20 S48-S42
  

    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