Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  929 1877 (N70W32 L=184)   462     53        0 ( 912/1850- 952/1895)  14.0 W24-W41 N67-N72
      2 1521 1523 (N25W47 L=199)  1664    268        2 (1418/1432-1603/1606) -75.1 W43-W50 N17-N32
      3  391 1432 (N44E45 L=106)   406     52        0 ( 368/1405- 417/1452)  33.7 E48-E41 N41-N44
      4 1147 1299 (N21W13 L=164)   334    155        0 (1103/1244-1202/1350) -63.4 W11-W15 N16-N24
      5 1560  998 (S06W33 L=184)   458    109        1 (1505/ 970-1601/1026) -33.7 W29-W35 S07-S02
      6 1651  929 (S12W38 L=190)   265     52        0 (1637/ 897-1663/ 962)  56.3 W36-W39 S13-S10
      7  608  907 (N07E28 L=123)   275     43        1 ( 583/ 890- 634/ 931)  14.0 E29-E25 N05-N08
      8  274  701 (S01E61 L= 91)   916     95        3 ( 247/ 650- 300/ 739) -90.0 E63-E56 S04-N01
      9 1687  682 (S29W40 L=191)   989    134        2 (1633/ 632-1735/ 735)  12.5 W35-W45 S32-S26
     10  268  552 (S11E71 L= 80)   837     75        2 ( 241/ 520- 290/ 578) -63.4 E74-E67 S13-S09
     11  450  511 (S15E51 L=100)   575     90        3 ( 433/ 468- 475/ 560)  40.6 E54-E47 S18-S12
  

    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