Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  751 1864 (N38E56 L=136)   723     64        3 ( 723/1838- 791/1885)   0.0 E61-E51 N35-N40
      2  688 1813 (N32E54 L=138)   529     51        0 ( 669/1785- 709/1844) -56.3 E57-E51 N29-N35
      3  585 1733 (N24E54 L=137)   385     53        1 ( 560/1713- 604/1758) -76.0 E56-E52 N22-N26
      4  832 1622 (N23E30 L=162)   667    237        3 ( 788/1509- 892/1714) -36.3 E36-E21 N17-N28
      5  456 1581 (N12E54 L=137)  1556    180        6 ( 378/1547- 539/1625)  37.9 E61-E47 N08-N15
      6 1135 1523 (N25E08 L=183)   316    125        0 (1118/1461-1154/1577) -60.3 E10-E05 N22-N29
      7  381 1481 (N04E55 L=136)  1133    158        4 ( 350/1401- 406/1540) -41.2 E59-E50 N00-N07
      8 1225 1263 (N13W04 L=196)   393    191        0 (1195/1180-1241/1341) -60.9 W00-W06 N07-N16
      9  170 1281 (S11E70 L=122)   667     86        2 ( 145/1253- 196/1331)  90.0 E73-E65 S12-S08
     10  438 1067 (S18E37 L=155)  1317    329        2 ( 302/ 954- 531/1164)  50.5 E46-E31 S27-S10
     11  178 1086 (S22E63 L=128)   473     89        1 ( 154/1049- 205/1129) -33.7 E66-E58 S24-S19
     12  549  835 (S29E24 L=167)  2374    620        6 ( 276/ 745- 836/ 915)  29.7 E49-E06 S41-S17
     13 1518  756 (S06W36 L=228)   218     59        0 (1498/ 729-1537/ 789) -26.6 W34-W38 S07-S04
     14 1200  614 (S23W23 L=214)   333    117        1 (1145/ 593-1259/ 636)  45.0 W20-W25 S25-S20
     15 1695  570 (S10W59 L=250)   360     83        0 (1683/ 519-1715/ 597) -45.0 W56-W61 S12-S07
     16 1636  525 (S14W56 L=248)   632     92        2 (1608/ 480-1664/ 574) -71.6 W54-W58 S16-S10
     17  575  468 (S53E16 L=176)   419    121        0 ( 542/ 428- 621/ 520) -29.1 E19-E10 S54-S49
     18 1652  439 (S17W65 L=256)   453     43        1 (1629/ 418-1678/ 454)  76.0 W61-W67 S19-S15
     19 1030  286 (S49W32 L=223)  1750    342        3 ( 872/ 218-1152/ 380)   0.0 W14-W47 S51-S45
  

    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