Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1205 1807 (N52W12 L= 89)   390     49        1 (1177/1782-1234/1828) -18.4 W09-W15 N49-N53
      2 1473 1655 (N40W33 L=110)   317     53        0 (1453/1628-1489/1684)  45.0 W30-W35 N37-N42
      3 1554 1482 (N28W36 L=112)   289     58        1 (1531/1459-1581/1502) -33.7 W34-W37 N26-N28
      4  288 1362 (N15E55 L= 21)   689    108        0 ( 258/1310- 328/1417)  74.1 E57-E52 N11-N18
      5  501 1269 (N09E35 L= 42)   258     75        0 ( 465/1254- 547/1287)   0.0 E37-E31 N08-N10
      6  265 1170 (N03E53 L= 23)   626    199        0 ( 241/1078- 295/1257) -68.2 E55-E49 S02-N08
      7 1848 1102 (N06W59 L=135)   709    104        3 (1816/1054-1872/1144)  78.7 W55-W62 N03-N08
      8 1139  867 (S13W08 L= 84)   323    206        1 (1055/ 808-1201/ 936)  45.0 W02-W11 S16-S08
      9 1002  881 (S13E01 L= 76)   186     97        0 ( 956/ 864-1051/ 893)  21.8 E03-W02 S14-S11
     10  211  845 (S16E60 L= 16)   592    112        0 ( 156/ 813- 263/ 872)   7.6 E69-E54 S17-S14
     11 1178  748 (S20W11 L= 87)   610    156        2 (1119/ 696-1221/ 816)  76.0 W07-W13 S23-S15
     12  180  753 (S22E68 L=  9)  1073     74        2 ( 137/ 711- 206/ 789) -45.0 E75-E64 S24-S19
     13  391  765 (S21E43 L= 33)   253     30        1 ( 371/ 741- 410/ 784)  63.4 E45-E41 S22-S20
     14  910  735 (S22E06 L= 70)   450    187        1 ( 832/ 705- 992/ 771) -15.3 E11-W00 S23-S20
     15 1650  719 (S19W45 L=122)   498     89        0 (1611/ 692-1686/ 752)  45.0 W41-W47 S20-S16
     16 1314  520 (S34W24 L=100)   532     92        0 (1276/ 493-1353/ 552) -23.2 W20-W27 S35-S32
     17 1011  162 (S68W08 L= 84)   486     57        0 ( 986/ 139-1040/ 188) -26.6 W03-W13 S71-S64
  

    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