Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1055 1390 (N27W12 L=  2)   507    166        1 ( 968/1374-1125/1407) -24.0 W06-W16 N24-N29
      2  567 1368 (N38E21 L=328)   235     47        0 ( 541/1347- 594/1383) -26.6 E23-E19 N37-N40
      3 1583 1305 (N07W41 L= 31)   453     61        2 (1561/1274-1610/1337)  76.0 W38-W43 N04-N08
      4 1493 1200 (N04W32 L= 21)   405    100        1 (1442/1174-1540/1225)  -8.1 W28-W35 N02-N05
      5  444 1118 (N26E35 L=315)  1443    265        6 ( 326/1078- 567/1157) -22.8 E45-E26 N21-N30
      6  187  943 (N20E63 L=286)  2170    250        2 (  89/ 894- 285/ 981)  -7.6 E83-E52 N16-N22
      7 1203  804 (S10W04 L=354)   903    233        0 (1135/ 725-1270/ 907)  32.5 E01-W10 S12-S05
      8  376  772 (N07E47 L=302)   611    160        0 ( 342/ 699- 416/ 839) -77.5 E48-E45 N02-N11
      9  946  571 (S17E17 L=333)   466    108        1 ( 906/ 523- 986/ 626)  29.7 E20-E12 S18-S14
     10  684  575 (S10E32 L=317)   244     65        1 ( 667/ 532- 697/ 603)  63.4 E33-E30 S12-S08
     11  744  574 (S12E28 L=321)   228     39        0 ( 721/ 548- 766/ 594)  26.6 E30-E26 S12-S10
     12 1699  564 (S42W36 L= 26)   516     58        0 (1668/ 543-1732/ 584) -20.6 W32-W40 S43-S39
     13  782  513 (S16E28 L=321)   329     99        1 ( 762/ 459- 813/ 544)  33.7 E31-E25 S19-S14
     14 1623  493 (S44W28 L= 17)   407     49        0 (1587/ 478-1654/ 508)  -8.1 W24-W31 S46-S42
     15  477  349 (S21E61 L=288)  1189    109        3 ( 455/ 285- 499/ 413)  32.7 E70-E55 S25-S16
     16 1023  250 (S41E27 L=322)   737    138        1 ( 998/ 194-1047/ 335)  55.3 E31-E22 S46-S33
     17 1316  221 (S57E08 L=342)  1528    211        0 (1226/ 194-1405/ 263) -23.6 E16-W02 S60-S50
     18  971  179 (S46E37 L=312)   503     27        1 ( 946/ 157- 994/ 201)  68.2 E40-E34 S48-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