Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  582 1805 (N45E54 L=264)  1230    106        3 ( 527/1774- 634/1831)  -9.5 E68-E44 N42-N47
      2  835 1735 (N40E23 L=295)   333    109        1 ( 791/1718- 890/1748)  -6.3 E27-E18 N38-N41
      3 1277 1708 (N43W13 L=331)   461     98        2 (1236/1681-1328/1736) -15.9 W09-W16 N41-N45
      4 1667 1564 (N37W47 L=  6)  1479    170        2 (1604/1506-1723/1616) -37.9 W41-W51 N33-N41
      5  921 1347 (N13E09 L=309)   787    323        1 ( 869/1245-1005/1493) -54.5 E13-E03 N07-N22
      6  221 1432 (N16E68 L=251)   636     97        2 ( 192/1405- 253/1459) -12.1 E75-E61 N13-N17
      7 1300 1062 (S00W16 L=334)   982    250        3 (1232/ 977-1411/1140)  45.0 W12-W22 S05-N05
      8 1692 1011 (N01W44 L=  2)   414     78        1 (1648/ 995-1729/1032)   8.1 W40-W47 N00-N02
      9 1602  964 (S03W37 L=356)   317     63        0 (1579/ 927-1626/ 993)  63.4 W35-W38 S05-S01
     10 1028  954 (S09W01 L=319)   197     94        0 (1000/ 922-1069/ 988)  53.1 W00-W03 S11-S06
     11 1542  686 (S20W38 L=357)   295    127        0 (1480/ 671-1592/ 714)  26.6 W33-W41 S21-S17
     12  898  552 (S36E04 L=314)  1614    371        2 ( 749/ 475-1019/ 682)  43.3 E15-W03 S42-S25
     13  324  596 (S36E55 L=263)   895    175        0 ( 308/ 512- 339/ 681)  47.7 E59-E49 S41-S30
     14  505  609 (S35E35 L=283)   466     65        1 ( 467/ 592- 551/ 626)   8.1 E38-E31 S35-S33
     15  945  389 (S47W01 L=319)   304     80        1 ( 905/ 374- 984/ 405)   0.0 E02-W04 S48-S45
     16 1486  319 (S44W55 L= 13)   449     56        0 (1450/ 302-1518/ 337)  31.0 W51-W58 S46-S42
  

    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