Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  729 1657 (N41E31 L=290)  4021    563        6 ( 488/1525- 931/1771) -11.0 E73-E11 N32-N47
      2 1138 1539 (N36W04 L=325)  1520    321        1 (1000/1495-1295/1565)   5.2 E06-W16 N32-N38
      3  414 1444 (N23E49 L=272)   797    108        1 ( 392/1381- 450/1503) -58.0 E53-E45 N18-N26
      4  243 1418 (N19E68 L=253)  1160    104        3 ( 222/1365- 271/1457)  90.0 E72-E64 N15-N21
      5 1111 1344 (N22W03 L=325)   438    184        4 (1071/1281-1147/1442) -63.4 W00-W05 N18-N28
      6 1215 1360 (N24W10 L=331)   429    138        0 (1200/1296-1230/1424)  90.0 W08-W11 N19-N28
      7 1462 1321 (N23W27 L=349)   264     38        0 (1435/1301-1487/1342) -26.6 W25-W29 N22-N24
      8  336 1139 (N03E48 L=273)  2200    345        4 ( 301/ 974- 377/1292) -86.8 E50-E45 S06-N12
      9 1838 1240 (N20W62 L= 24)   967     99        1 (1803/1203-1875/1276) -45.0 W58-W67 N17-N22
     10 1623 1246 (N19W40 L=  1)   300     81        0 (1580/1235-1669/1256)   0.0 W36-W43 N18-N20
     11 1897 1161 (N15W70 L= 31)  1312    110        1 (1855/1127-1929/1203)  21.8 W63-W78 N13-N17
     12 1903  958 (N03W70 L= 31)   655     81        1 (1872/ 922-1934/ 981)  21.8 W65-W76 N00-N04
     13 1585  920 (S01W37 L=358)   228     59        0 (1567/ 902-1612/ 948) -14.0 W35-W39 S01-N00
     14  559  769 (S17E29 L=293)   238     81        0 ( 518/ 754- 594/ 789) -11.3 E31-E26 S18-S16
     15 1329  658 (S19W23 L=344)  1146    301        2 (1212/ 618-1471/ 723)  26.6 W15-W31 S21-S13
     16 1022  521 (S30W04 L=325)   704    139        1 ( 956/ 496-1091/ 544)  14.0 W00-W08 S32-S28
     17  662  376 (S45E25 L=296)   325     25        1 ( 642/ 360- 685/ 398) -33.7 E27-E23 S46-S43
  

    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