Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1006 1791 (N41E25 L=152)   345     45        1 ( 978/1775-1036/1810)  56.3 E27-E22 N39-N43
      2  618 1588 (N17E40 L=136)   430    128        1 ( 568/1559- 658/1629)  51.3 E43-E38 N15-N20
      3  724 1572 (N18E32 L=144)   286     61        0 ( 703/1537- 745/1604) -45.0 E34-E29 N16-N20
      4 1751 1372 (N33W43 L=220)   460     59        1 (1726/1343-1776/1397) -36.9 W40-W46 N31-N34
      5  913 1324 (N07E13 L=163)   210    101        0 ( 892/1289- 947/1350)  63.4 E14-E11 N05-N09
      6 1221 1207 (N08W07 L=184)   862    274        0 (1122/1166-1353/1234)  30.3 W02-W15 N02-N10
      7 1743 1138 (N18W43 L=220)   443     95        1 (1711/1102-1777/1172) -23.2 W40-W47 N16-N19
      8  633 1025 (S15E23 L=153)   364    160        1 ( 601/ 961- 672/1100)  77.5 E24-E21 S19-S10
      9  231  952 (S27E56 L=120)  1631    239        4 ( 184/ 863- 279/1060)  55.3 E63-E51 S32-S19
     10 1749  980 (N09W47 L=224)   500     70        0 (1712/ 959-1786/1000)  15.9 W43-W50 N07-N10
     11 1475  581 (S20W41 L=217)   309     53        1 (1444/ 562-1507/ 596)  11.3 W38-W43 S21-S18
  

    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