Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  577 1677 (N31E51 L=  7)   832    185        1 ( 560/1582- 607/1745) -47.5 E56-E43 N24-N36
      2 1520 1507 (N40W28 L= 86)   311     50        0 (1503/1478-1540/1529)  59.0 W26-W29 N37-N42
      3  861 1495 (N25E20 L= 38)   245     45        0 ( 843/1477- 886/1516) -14.0 E22-E18 N23-N25
      4  750 1451 (N20E26 L= 32)   540     93        2 ( 715/1423- 806/1479)   8.1 E29-E22 N17-N21
      5 1337 1421 (N29W13 L= 71)   313     89        0 (1289/1406-1385/1436)   8.1 W09-W16 N28-N30
      6 1648 1263 (N25W39 L= 97)   267     39        1 (1630/1241-1663/1287) -63.4 W37-W40 N24-N26
      7  173 1243 (S03E71 L=348)   588     54        1 ( 154/1209- 190/1270) -20.6 E76-E67 S04-S01
      8 1684 1165 (N20W42 L=101)   747     97        2 (1640/1133-1724/1216)  -7.1 W38-W46 N17-N23
      9  412  871 (S21E39 L= 20)  3328    661        7 ( 196/ 711- 636/1108)  36.0 E68-E24 S35-S03
     10 1741 1067 (N15W48 L=107)   557     82        2 (1698/1041-1771/1091)  23.2 W44-W51 N13-N16
     11  151  988 (S19E69 L=349)  3307    308        6 ( 109/ 891- 199/1067)  50.7 E80-E61 S24-S13
     12 1680  948 (N07W45 L=103)   298     59        0 (1661/ 923-1704/ 975) -21.8 W42-W47 N05-N08
     13 1492  890 (S00W31 L= 90)   598    104        1 (1439/ 863-1539/ 917)  18.4 W28-W34 S02-N01
     14 1148  699 (S18W14 L= 72)   882    362        4 ( 986/ 656-1312/ 727)  27.8 W04-W23 S23-S13
     15 1694  584 (S15W58 L=117)  2049    325        3 (1597/ 491-1792/ 708)  59.0 W53-W64 S21-S05
     16  661  539 (S39E16 L= 42)   365     39        2 ( 646/ 511- 680/ 567) -63.4 E17-E14 S40-S36
     17 1461  396 (S31W49 L=107)   934    112        2 (1400/ 371-1506/ 436)  29.7 W45-W52 S34-S28
     18 1277  379 (S37W35 L= 93)   339     59        1 (1251/ 358-1301/ 400) -21.8 W32-W37 S37-S35
  

    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