Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  792 1671 (N48W01 L=220)   323     55        0 ( 752/1653- 817/1687) -36.9 E02-W03 N46-N49
      2 1600 1624 (N23W58 L=277)   933    131        5 (1549/1579-1656/1673) -50.2 W55-W62 N19-N27
      3  597 1627 (N51E18 L=201)   461     81        0 ( 557/1601- 627/1661)  23.2 E21-E14 N49-N53
      4 1257 1606 (N31W31 L=250)  2049    422        2 (1067/1557-1434/1644) -26.6 W17-W41 N24-N35
      5 1638 1534 (N17W55 L=274)  1601    283        6 (1518/1478-1747/1583) -11.3 W45-W70 N14-N21
      6 1681 1332 (N05W49 L=268)  1013    173        3 (1601/1282-1735/1369) -29.1 W43-W54 N02-N08
      7  188 1256 (N34E61 L=158)   473     73        1 ( 156/1222- 218/1281) -21.8 E66-E55 N31-N35
      8  613 1171 (N20E21 L=198)  1386    426        4 ( 423/1103- 774/1234) -34.3 E34-E12 N12-N27
      9  171 1115 (N25E61 L=158)  1152    156        4 ( 106/1079- 244/1150)  -2.9 E72-E52 N23-N27
     10  321  963 (N14E46 L=173)   427    163        1 ( 291/ 913- 354/1057) -83.7 E48-E43 N10-N19
     11 1435  839 (S16W20 L=239)   504    224        3 (1370/ 796-1498/ 890)   6.3 W15-W24 S19-S14
     12  882  714 (S11E15 L=204)   424    158        0 ( 860/ 656- 924/ 771)  60.3 E16-E11 S15-S08
     13 1879  743 (S35W62 L=281)  1086    116        0 (1822/ 713-1938/ 761)  -5.0 W52-W77 S36-S33
     14 1494  640 (S30W21 L=240)   922    239        0 (1386/ 581-1564/ 712)  11.3 W12-W27 S32-S26
     15  860  441 (S28E25 L=194)  1819    452        4 ( 676/ 364-1064/ 553)   5.4 E40-E08 S29-S24
     16  459  434 (S20E54 L=165)   377     52        0 ( 433/ 407- 481/ 454) -63.4 E55-E52 S21-S17
     17 1234  340 (S45E04 L=215)  1794    357        2 (1068/ 294-1371/ 436)   2.2 E17-W09 S47-S40
  

    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