Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1343 1745 (N54W35 L=210)   366     46        0 (1325/1718-1364/1766)  29.7 W32-W39 N51-N55
      2 1395 1742 (N53W42 L=217)   432     43        0 (1374/1718-1421/1758)  20.6 W38-W46 N51-N54
      3  972 1729 (N53E06 L=170)   295     39        0 ( 954/1707- 990/1752)  59.0 E07-E03 N50-N55
      4  850 1731 (N53E18 L=157)   323     49        0 ( 822/1716- 881/1744)  21.8 E20-E15 N51-N54
      5  683 1549 (N37E28 L=148)   659    162        3 ( 618/1492- 745/1580) -24.8 E34-E21 N33-N39
      6 1730 1548 (N36W69 L=244)  1054     84        0 (1693/1526-1776/1563)  -3.2 W62-W83 N34-N36
      7  855 1497 (N34E13 L=163)   495    145        1 ( 800/1470- 928/1518)  18.4 E16-E07 N32-N35
      8 1053 1460 (N31W02 L=177)   290     63        0 (1021/1444-1090/1476) -26.6 W00-W04 N30-N32
      9  228 1268 (N17E63 L=112)   477     77        1 ( 197/1234- 243/1304) -38.7 E67-E60 N14-N19
     10  354 1274 (N18E49 L=126)   447    104        0 ( 319/1240- 404/1302)  23.2 E51-E44 N15-N19
     11  369 1143 (N10E46 L=130)   321     40        1 ( 346/1120- 395/1164) -26.6 E47-E43 N08-N11
     12 1215  672 (S19W13 L=188)   354    160        0 (1136/ 655-1289/ 692)   5.7 W07-W17 S19-S17
     13  333  656 (S21E53 L=122)   353     94        1 ( 308/ 636- 358/ 690)  33.7 E56-E49 S22-S18
     14  658  488 (S32E28 L=148)  3130    510        6 ( 499/ 346- 827/ 667)  33.3 E51-E12 S44-S18
     15  372  523 (S31E55 L=121)  1337    226        0 ( 341/ 425- 423/ 618)  31.6 E70-E44 S39-S23
     16 1508  506 (S31W37 L=213)   314     69        0 (1476/ 487-1542/ 528)  21.8 W34-W39 S32-S29
     17 1395  418 (S37W30 L=205)   422     63        0 (1367/ 389-1417/ 452)  68.2 W28-W32 S39-S34
     18 1424  321 (S46W39 L=214)   483     71        1 (1394/ 297-1440/ 358)  90.0 W36-W42 S48-S42
     19 1480  299 (S49W48 L=224)   507     35        1 (1455/ 282-1504/ 322)   0.0 W44-W53 S50-S46
  

    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