Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1616 1656 (N27W68 L= 61)  1467    234        2 (1511/1599-1712/1698) -40.6 W60-W77 N20-N32
      2 1508 1546 (N24W48 L= 41)   588    117        0 (1484/1489-1526/1604)  54.5 W45-W51 N20-N27
      3 1165 1536 (N34W24 L= 16)   761    159        0 (1090/1490-1231/1571) -53.1 W19-W26 N29-N37
      4 1192 1250 (N16W15 L=  7)   608    247        1 (1090/1198-1285/1304) -48.0 W10-W19 N11-N21
      5  782 1173 (N21E11 L=341)   258     65        1 ( 753/1141- 808/1200)   0.0 E13-E09 N19-N22
      6  471  997 (N17E36 L=316)   798    213        0 ( 455/ 899- 499/1092)  80.5 E37-E33 N10-N22
      7 1269  953 (S02W13 L=  5)   383    181        1 (1214/ 880-1315/1012)  29.1 W07-W16 S05-N00
      8  168  821 (N10E71 L=281)   833     87        1 ( 125/ 796- 203/ 852)  10.0 E82-E64 N07-N11
      9 1202  646 (S19W02 L=354)   580    224        1 (1123/ 604-1245/ 735)  45.0 E03-W05 S21-S13
     10 1421  642 (S25W16 L=  8)  1121    299        1 (1307/ 574-1557/ 716)   6.0 W07-W26 S26-S23
     11  472  574 (S09E47 L=305)   497     74        1 ( 436/ 545- 504/ 604)   8.1 E51-E44 S10-S07
     12  618  357 (S26E47 L=305)   766    128        1 ( 584/ 310- 638/ 426)  38.7 E52-E41 S29-S21
  

    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