Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  887 1726 (N32E31 L=140)  1571    354        1 ( 807/1616- 975/1783)  85.9 E36-E22 N22-N36
      2 1290 1671 (N40E02 L=169)   707    191        1 (1209/1614-1336/1727) -20.0 E09-W02 N36-N42
      3  373 1663 (N17E71 L=100)   792    108        3 ( 340/1614- 394/1701) -60.3 E76-E65 N13-N20
      4  486 1479 (N07E46 L=125)   781    140        3 ( 415/1447- 533/1527)  45.0 E50-E42 N03-N10
      5 1804 1407 (N41W50 L=220)   587     98        2 (1773/1359-1829/1444) -29.7 W45-W52 N37-N43
      6 1510 1202 (N16W24 L=194)   214     94        0 (1496/1157-1528/1246)  90.0 W23-W24 N12-N18
      7  878 1127 (S05E11 L=160)   229     89        0 ( 864/1082- 892/1170) -78.7 E12-E09 S07-S02
      8  315  893 (S31E46 L=124)  1329    345        1 ( 230/ 780- 384/1045)  48.4 E57-E41 S39-S20
      9  855  862 (S20E06 L=165)   423    151        0 ( 807/ 795- 889/ 907) -29.7 E09-E02 S23-S18
     10  695  840 (S26E15 L=156)   234     63        1 ( 669/ 808- 717/ 869) -36.9 E17-E13 S27-S24
     11  455  763 (S36E32 L=138)   676    220        0 ( 423/ 683- 512/ 851) -45.0 E35-E26 S40-S30
     12 1599  726 (S07W43 L=214)   527    100        2 (1576/ 694-1645/ 770) -68.2 W40-W46 S09-S04
     13 1454  450 (S26W47 L=218)   732    175        0 (1370/ 422-1532/ 486)  41.6 W42-W51 S30-S21
     14 1107  338 (S45W32 L=203)   476    112        0 (1083/ 292-1139/ 388) -54.5 W30-W36 S47-S40
     15  688  304 (S64E01 L=170)  1827    226        5 ( 607/ 269- 798/ 365)   0.0 E13-W13 S68-S60
     16 1569  335 (S28W69 L=239)  1682    158        3 (1486/ 304-1631/ 363)  12.3 W58-W82 S31-S25
     17 1173  280 (S46W42 L=213)   295     49        0 (1158/ 252-1187/ 303) -45.0 W38-W44 S48-S44
     18 1262  230 (S46W55 L=226)  2283    164        6 (1184/ 191-1343/ 286)  20.6 W45-W66 S49-S42
  

    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