Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  933 1840 (N42E40 L=134)   956    168        1 ( 861/1802-1003/1879)   0.0 E50-E32 N40-N43
      2  455 1652 (N16E61 L=114)  2113    337        2 ( 331/1611- 576/1693)   8.1 E77-E48 N11-N18
      3 1493 1544 (N40W14 L=189)   437    131        1 (1476/1493-1522/1610) -53.1 W11-W17 N36-N44
      4 1546 1579 (N45W19 L=193)   313     89        1 (1517/1544-1576/1606) -18.4 W15-W21 N43-N46
      5  492 1415 (N02E44 L=131)   302     76        0 ( 460/1389- 530/1432)  45.0 E45-E41 N00-N03
      6  814 1253 (N01E18 L=157)   780    231        0 ( 756/1151- 889/1340) -26.6 E23-E10 S03-N04
      7 1714 1114 (N20W40 L=215)  1126    154        1 (1650/1063-1764/1165)  -4.8 W33-W45 N18-N21
      8  310 1124 (S19E50 L=125)   428     66        0 ( 291/1084- 333/1158)  68.2 E50-E47 S21-S16
      9  215 1125 (S20E60 L=114)   354     26        1 ( 198/1100- 237/1147) -71.6 E62-E57 S21-S18
     10  156 1058 (S25E68 L=107)  2711    226        4 ( 119/ 977- 193/1143)  51.3 E75-E61 S30-S20
     11  885 1087 (S07E09 L=165)   326     66        1 ( 856/1049- 916/1119)  78.7 E10-E08 S09-S04
     12 1074  698 (S23W13 L=187)   669    180        0 (1000/ 646-1155/ 733)  60.9 W10-W16 S27-S18
     13 1354  644 (S17W31 L=205)  1201    287        0 (1212/ 614-1492/ 669)  34.5 W23-W39 S22-S11
     14  822  557 (S39W02 L=177)   527    135        1 ( 795/ 497- 859/ 614) -69.4 W00-W04 S42-S34
     15 1122  467 (S34W26 L=200)   245     92        0 (1081/ 446-1158/ 491)   0.0 W22-W29 S35-S33
     16 1576  459 (S19W57 L=231)   852    185        0 (1520/ 428-1642/ 489)  49.4 W52-W62 S22-S15
     17  825  440 (S46W08 L=182)   315     58        1 ( 796/ 418- 857/ 467)  -8.1 W04-W11 S47-S45
  

    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