Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1127 1409 (N26W17 L=  6)  1087    417        4 ( 932/1377-1301/1466) -15.6 W04-W29 N22-N30
      2  577 1373 (N39E20 L=328)   267     42        1 ( 552/1354- 605/1389) -26.6 E22-E18 N37-N40
      3 1591 1317 (N07W42 L= 31)   452     58        0 (1566/1282-1615/1350)  76.0 W40-W44 N05-N09
      4 1501 1205 (N04W32 L= 21)   399     91        1 (1450/1187-1539/1231)   0.0 W28-W35 N02-N05
      5  456 1121 (N26E34 L=315)  1380    258        4 ( 330/1084- 579/1157) -22.8 E44-E25 N21-N29
      6  174  954 (N21E65 L=284)  1600    192        2 (  89/ 917- 252/ 984)  -6.3 E82-E54 N17-N23
      7 1213  805 (S10W05 L=353)   909    234        0 (1142/ 718-1275/ 893)  28.6 E01-W10 S13-S05
      8 1528  830 (S18W24 L= 13)   270     98        0 (1492/ 805-1576/ 856)   0.0 W21-W28 S18-S16
      9  381  782 (N08E47 L=302)   435    111        1 ( 355/ 728- 413/ 841) -90.0 E48-E45 N04-N11
     10  697  582 (S10E31 L=317)   433    178        2 ( 660/ 517- 747/ 656) -90.0 E32-E29 S15-S04
     11  961  577 (S17E16 L=333)   443    129        0 ( 921/ 524-1000/ 633)  26.6 E19-E11 S19-S14
     12  755  578 (S12E28 L=321)   215     36        0 ( 732/ 551- 777/ 598)  26.6 E29-E25 S12-S10
     13 1710  569 (S42W37 L= 26)   630     80        1 (1677/ 544-1748/ 596)   8.1 W33-W42 S43-S39
     14  957  425 (S26E21 L=327)   269     59        1 ( 936/ 396- 975/ 457)  45.0 E23-E18 S28-S24
     15  484  345 (S21E61 L=288)  1073    135        1 ( 469/ 280- 502/ 414)  34.7 E69-E54 S25-S16
     16 1036  240 (S43E27 L=321)   698    112        1 (1021/ 188-1054/ 302)  48.0 E31-E22 S47-S37
     17 1322  225 (S56E07 L=341)  1483    221        3 (1231/ 189-1410/ 269) -23.6 E16-W03 S60-S50
     18  980  181 (S46E37 L=312)   548     66        1 ( 952/ 157-1008/ 203)  71.6 E39-E32 S49-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