Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  864 1857 (N46E48 L= 77)  2312    255        5 ( 790/1775- 964/1924) -10.1 E75-E30 N42-N51
      2  714 1678 (N28E41 L= 83)   306     75        1 ( 687/1644- 746/1707) -18.4 E45-E36 N26-N29
      3  783 1516 (N19E28 L= 96)   341    153        0 ( 732/1480- 846/1568) -15.3 E33-E22 N18-N21
      4 1379 1196 (N15W16 L=141)   447    160        4 (1338/1152-1409/1248)  90.0 W13-W18 N12-N18
      5  707 1138 (S05E21 L=103)   269     81        2 ( 673/1115- 749/1157)  36.9 E23-E18 S06-S02
      6  598  903 (S21E23 L=101)   500    213        0 ( 515/ 828- 669/ 953)  54.0 E27-E19 S27-S16
      7  147  859 (S34E72 L= 53)   660     62        0 ( 126/ 828- 172/ 882)  17.1 E79-E65 S35-S31
      8 1669  732 (S03W49 L=173)   728    149        1 (1633/ 672-1689/ 810) -66.8 W45-W51 S06-N01
      9 1433  533 (S21W41 L=165)   363     94        2 (1382/ 507-1470/ 551)   7.1 W36-W44 S21-S19
     10  698  277 (S64E00 L=124)   402     59        1 ( 674/ 251- 723/ 306) -33.7 E04-W04 S66-S61
  

    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