Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1545 1685 (N51W24 L=171)   644     50        1 (1520/1650-1563/1715)  74.1 W21-W27 N47-N54
      2  589 1509 (N12E38 L=108)   345     66        0 ( 561/1477- 616/1537)  76.0 E39-E37 N10-N14
      3 1617 1477 (N36W30 L=177)   340     72        0 (1585/1448-1647/1500) -21.8 W27-W33 N34-N36
      4 1551 1411 (N29W25 L=171)   382    112        0 (1514/1357-1587/1459) -33.7 W21-W28 N26-N31
      5 1032  607 (S31W11 L=157)   539    243        0 ( 991/ 547-1114/ 692) -21.8 W07-W17 S34-S25
      6 1380  623 (S21W32 L=179)   539    103        3 (1351/ 574-1410/ 674) -81.9 W29-W33 S24-S17
      7  905  272 (S58W18 L=165)  1009    117        2 ( 841/ 247- 957/ 315)  24.4 W13-W24 S62-S53
  

    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