Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  359 1586 (N27E62 L=146)   782    145        2 ( 332/1552- 384/1633) -80.5 E64-E58 N24-N30
      2 1518 1413 (N31W30 L=238)   725    192        1 (1434/1362-1595/1451)  -8.7 W23-W36 N28-N33
      3  552 1343 (N14E35 L=173)   458    109        0 ( 541/1285- 568/1396) -66.8 E37-E34 N10-N17
      4 1305 1333 (N23W14 L=223)   292     92        1 (1285/1288-1329/1382) -63.4 W13-W16 N20-N26
      5 1134 1296 (N19W04 L=212)   218     49        0 (1107/1274-1164/1319) -26.6 W01-W05 N17-N19
      6 1853  960 (N06W63 L=271)  4388    593       11 (1803/ 756-1914/1197) -81.3 W58-W69 S05-N21
      7 1817 1090 (N14W57 L=265)  1087    155        4 (1788/1018-1839/1156)  83.7 W53-W59 N09-N18
      8  194  890 (S18E61 L=147)  3193    552       10 ( 132/ 801- 257/ 981)  79.7 E69-E53 S23-S12
      9  351  930 (S13E45 L=163)   338     64        2 ( 329/ 905- 374/ 967) -56.3 E46-E42 S15-S11
     10  994  864 (S09W00 L=208)   209     78        0 ( 977/ 827-1007/ 907)  90.0 W00-W00 S11-S06
     11  317  794 (S23E48 L=160)   479    118        2 ( 283/ 769- 358/ 820)  23.2 E50-E43 S24-S20
     12 1609  695 (S12W44 L=252)  1410    306        4 (1498/ 603-1701/ 784)  52.4 W37-W51 S19-S06
     13  841  723 (S20E08 L=200)   374    134        0 ( 791/ 698- 914/ 758)  -7.1 E11-E03 S21-S18
     14 1425  736 (S12W29 L=238)   427    133        0 (1362/ 722-1483/ 748)  14.0 W25-W33 S13-S10
     15 1261  633 (S21W21 L=229)   277     79        0 (1216/ 607-1298/ 659)  45.0 W18-W22 S22-S18
     16 1634  630 (S16W48 L=256)   645     89        1 (1606/ 593-1670/ 655)  63.4 W45-W51 S18-S14
     17 1599  531 (S23W49 L=257)   474    136        2 (1563/ 495-1640/ 567)   0.0 W45-W53 S25-S20
     18  608  530 (S37E24 L=184)   409     98        0 ( 568/ 515- 644/ 558)  26.6 E27-E21 S38-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