Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1518 1687 (N28W56 L= 79)  2578    325        5 (1389/1640-1680/1730) -21.8 W44-W70 N21-N31
      2 1630 1514 (N16W53 L= 76)   477    127        1 (1573/1476-1682/1544) -50.2 W49-W56 N12-N18
      3  233 1223 (N32E54 L=329)  1532    276        3 ( 143/1110- 318/1293) -30.0 E70-E44 N23-N38
      4 1868 1173 (S09W61 L= 84)  1070    175        2 (1831/1104-1899/1230)  60.3 W56-W65 S12-S05
      5 1802 1157 (S08W54 L= 77)   885    134        2 (1776/1095-1826/1228)  74.1 W51-W56 S11-S04
      6 1703 1180 (S04W45 L= 69)   415    105        2 (1673/1133-1742/1215)  23.2 W42-W49 S06-S02
      7 1409 1128 (N00W24 L= 47)   342    149        0 (1355/1089-1469/1149)  -8.1 W20-W27 S02-N01
      8  253 1049 (N21E52 L=331)   427     84        1 ( 218/1019- 292/1076) -23.2 E55-E47 N18-N22
      9 1106  932 (S04W02 L= 25)  1894    537        4 ( 933/ 751-1250/1058) -61.2 E05-W06 S16-N05
     10 1436  981 (S09W22 L= 45)   219     46        1 (1419/ 956-1456/1006)  45.0 W21-W24 S10-S06
     11  539  965 (N10E30 L=353)   551    168        2 ( 451/ 946- 620/ 990) -24.4 E35-E24 N07-N13
     12 1588  735 (S27W29 L= 52)  1388    246        6 (1506/ 689-1688/ 815) -26.6 W23-W37 S30-S20
     13 1329  472 (S37W06 L= 29)  3512    686        9 (1168/ 270-1426/ 764)  53.4 E12-W14 S49-S18
     14  687  666 (S09E28 L=355)   228     30        0 ( 666/ 643- 708/ 689) -71.6 E29-E26 S10-S07
     15  748  496 (S21E30 L=353)  1417    292        2 ( 664/ 427- 811/ 636)  38.2 E36-E22 S24-S13
     16  632  601 (S12E34 L=349)   433     82        2 ( 592/ 572- 678/ 628)   0.0 E36-E29 S13-S10
     17  540  404 (S22E50 L=333)  1702    271        3 ( 448/ 333- 634/ 471)   9.1 E64-E39 S25-S18
  

    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