Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  648 1427 (N15E32 L=286)  3340    957       15 ( 509/1120- 829/1710)  12.0 E56-E13 N00-N30
      2 1438 1631 (N47W18 L=336)   693    186        2 (1385/1592-1498/1673) -10.3 W12-W24 N44-N50
      3 1583 1510 (N42W31 L=350)   436     81        1 (1552/1487-1614/1535)  26.6 W28-W34 N39-N43
      4  539 1011 (S12E30 L=289)   363    104        3 ( 502/ 962- 563/1047) -68.2 E32-E28 S14-S09
      5 1630 1018 (N10W39 L=357)   347     93        0 (1592/ 993-1668/1039)  -8.1 W35-W42 N09-N11
      6  342  941 (S20E44 L=274)   909    213        2 ( 266/ 922- 414/ 966)  13.0 E50-E37 S21-S17
      7 1706  625 (S11W57 L= 15)  1109    278        3 (1634/ 596-1797/ 654)  12.5 W49-W67 S13-S08
      8 1078  546 (S29W14 L=333)   595    118        2 (1017/ 514-1131/ 584)  54.5 W10-W16 S32-S25
      9  765  501 (S40E06 L=312)  1342    387        4 ( 583/ 479- 921/ 535)  17.1 E21-W05 S45-S34
  

    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