Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  537 1605 (N15E49 L=262)   382     81        2 ( 523/1569- 554/1637) -59.0 E51-E47 N12-N17
      2  324 1538 (N07E65 L=246)  1814    170        4 ( 274/1470- 370/1614)  90.0 E69-E63 N02-N12
      3  775 1567 (N19E31 L=280)   264     91        0 ( 738/1541- 810/1598)  -8.1 E34-E27 N17-N20
      4  953 1420 (N15E15 L=296)   305    100        0 ( 925/1376- 975/1469) -80.5 E16-E14 N11-N17
      5 1530 1077 (N12W28 L=339)   391     97        1 (1505/1030-1567/1125) -23.2 W24-W31 N10-N14
      6 1649 1093 (N16W36 L=347)   400    118        0 (1616/1055-1687/1121)  63.4 W33-W38 N13-N19
      7 1687  890 (N07W44 L=355)  2373    428        5 (1608/ 767-1732/ 984) -31.6 W36-W51 N01-N11
      8 1845  827 (N09W61 L= 13)   709    142        1 (1818/ 774-1868/ 873) -26.6 W56-W66 N06-N10
      9  549  703 (S37E22 L=290)   518    151        1 ( 491/ 670- 618/ 731)  37.9 E26-E17 S40-S33
     10 1678  438 (S17W66 L= 18)  1244    306        2 (1651/ 384-1697/ 505) -18.4 W61-W73 S18-S13
  

    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