Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  327 1437 (N29E58 L=302)   698    174        0 ( 272/1382- 366/1470)  45.0 E62-E52 N25-N31
      2 1321 1369 (N25W21 L= 21)   411    181        1 (1246/1314-1379/1402) -36.9 W15-W24 N21-N27
      3 1525 1338 (N22W36 L= 36)   298     66        1 (1493/1309-1544/1368)  45.0 W33-W37 N20-N24
      4  286 1328 (N22E58 L=302)   496     91        1 ( 264/1296- 323/1367) -26.6 E62-E53 N20-N24
      5 1722 1344 (N22W54 L= 55)   657    117        1 (1673/1313-1770/1366) -16.7 W49-W59 N19-N23
      6  782 1303 (N21E16 L=344)   242     97        0 ( 753/1269- 810/1345)  59.0 E17-E14 N19-N24
      7  344 1111 (N09E48 L=313)   308     86        0 ( 327/1071- 361/1146)  90.0 E49-E45 N06-N10
      8 1749  866 (S08W52 L= 52)   486     41        1 (1721/ 845-1772/ 889)  26.6 W48-W54 S09-S06
      9 1824  794 (S13W62 L= 62)  1515    227        3 (1791/ 709-1858/ 888)  84.8 W56-W67 S18-S07
     10 1753  774 (S14W53 L= 54)   692     57        1 (1718/ 753-1792/ 793)   0.0 W49-W57 S15-S12
     11 1653  711 (S18W45 L= 45)   385    167        0 (1609/ 657-1684/ 774)  63.4 W41-W46 S21-S13
     12 1419  699 (S17W26 L= 26)   567    136        3 (1356/ 669-1490/ 741)  18.4 W21-W30 S19-S14
     13  476  670 (S19E39 L=321)   414    106        0 ( 458/ 631- 490/ 725)  56.3 E41-E37 S21-S15
     14  276  688 (S18E59 L=302)   393     29        0 ( 260/ 666- 294/ 710)  56.3 E61-E56 S19-S16
     15 1722  644 (S23W54 L= 54)   481     82        2 (1688/ 617-1749/ 676)  38.7 W50-W56 S24-S20
     16 1778  610 (S26W63 L= 63)  1212    150        2 (1718/ 580-1846/ 641)   2.9 W54-W75 S27-S23
     17  834  456 (S34E15 L=345)  2741    360        7 ( 727/ 388- 978/ 572) -18.4 E22-E04 S39-S25
     18  578  293 (S49E49 L=312)   591     60        1 ( 564/ 258- 599/ 323)  30.3 E56-E44 S52-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