Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1550 1733 (N30W64 L= 73)  1363    152        1 (1501/1684-1614/1765)  -6.0 W56-W75 N27-N31
      2 1810 1190 (S06W56 L= 65)  1232    235        4 (1777/1099-1846/1279) -85.2 W52-W58 S12-N00
      3 1902 1220 (S07W68 L= 77)  1084    107        4 (1872/1169-1928/1271) -90.0 W64-W73 S10-S04
      4  206 1145 (N27E56 L=313)   394     61        1 ( 177/1127- 234/1171)  15.9 E60-E52 N26-N28
      5 1332 1015 (S04W17 L= 26)  1865    523        3 (1167/ 840-1453/1154) -61.2 W09-W21 S16-N05
      6  424 1117 (N21E36 L=334)   234     40        1 ( 396/1097- 452/1137) -36.9 E37-E33 N19-N22
      7 1577 1120 (S04W35 L= 44)   258     36        1 (1551/1101-1605/1137) -33.7 W33-W36 S05-S02
      8 1679 1016 (S13W40 L= 49)   990    286        3 (1604/ 938-1731/1097)  22.6 W33-W45 S16-S09
      9 1610 1058 (S09W36 L= 45)   342     67        1 (1588/1018-1630/1092)  76.0 W34-W37 S10-S06
     10  749 1044 (N10E15 L=354)   525    152        0 ( 669/1025- 825/1070) -29.1 E19-E10 N07-N13
     11 1687  821 (S24W38 L= 47)   508     92        1 (1667/ 783-1716/ 884) -80.5 W36-W40 S26-S20
     12 1777  792 (S29W46 L= 55)   905    116        0 (1715/ 771-1834/ 817)  -4.8 W40-W53 S30-S27
     13 1566  581 (S36W26 L= 35)  3005    437        8 (1483/ 443-1642/ 778)  77.2 W18-W32 S46-S23
     14  890  733 (S10E14 L=355)   313     75        1 ( 860/ 691- 915/ 759) -76.0 E15-E12 S12-S08
     15  943  711 (S13E12 L=358)   295     68        1 ( 914/ 690- 975/ 736) -33.7 E13-E09 S13-S11
     16  765  475 (S23E30 L=339)  3273    812        9 ( 514/ 246-1015/ 696)  15.9 E64-E08 S33-S14
     17  488  584 (S10E44 L=325)   460    140        0 ( 421/ 564- 545/ 619)  -5.2 E50-E38 S11-S09
     18 1424  480 (S40W13 L= 22)   293     54        1 (1396/ 455-1445/ 504)  45.0 W11-W15 S41-S37
  

    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