Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  592 1532 (N34E28 L=130)   362     60        1 ( 571/1495- 610/1561) -68.2 E29-E25 N31-N36
      2  878 1523 (N31E06 L=151)   307    109        0 ( 833/1497- 922/1553) -39.8 E08-E02 N28-N33
      3  363 1396 (N27E46 L=112)   671    155        1 ( 335/1335- 394/1476) -48.0 E51-E41 N22-N32
      4  658 1395 (N24E21 L=136)   553     98        0 ( 617/1358- 712/1432) -35.5 E24-E17 N21-N26
      5  246 1285 (N21E56 L=102)  1618    183        9 ( 189/1226- 293/1346) -34.7 E65-E50 N16-N25
      6 1396 1077 (S01W23 L=180)   533     82        1 (1363/1046-1434/1118) -26.6 W20-W25 S02-N01
      7 1858  979 (S10W60 L=217)  1390    222        3 (1767/ 924-1933/1015)   5.4 W50-W71 S13-S07
      8 1798  860 (S16W54 L=211)   621     95        3 (1769/ 822-1823/ 906)  78.7 W51-W57 S18-S13
      9 1857  858 (S17W61 L=219)   399     66        1 (1831/ 823-1874/ 883)  90.0 W58-W63 S18-S15
     10 1340  798 (S17W18 L=175)   238     95        0 (1304/ 766-1382/ 824)  21.8 W15-W20 S19-S15
     11  701  720 (S17E23 L=134)   219     63        0 ( 684/ 691- 735/ 741)  33.7 E24-E20 S19-S15
     12 1127  505 (S35W03 L=160)  3161    938        5 ( 852/ 365-1511/ 739)  13.2 E17-W30 S44-S22
     13  269  701 (S15E58 L= 99)   453     65        1 ( 247/ 673- 300/ 722)  15.9 E61-E54 S16-S13
     14 1729  636 (S30W52 L=209)   884    167        1 (1689/ 604-1782/ 670) -13.0 W46-W59 S32-S27
     15  461  621 (S21E43 L=114)   255     45        2 ( 438/ 597- 482/ 649)  45.0 E44-E41 S22-S19
     16  357  556 (S24E55 L=102)   715    103        1 ( 328/ 513- 395/ 600) -71.6 E57-E52 S27-S21
     17  745  549 (S28E23 L=134)   363     80        0 ( 709/ 527- 794/ 566)   0.0 E26-E19 S29-S27
     18 1390  246 (S60W33 L=191)  1222    125        0 (1337/ 224-1453/ 284) -10.8 W24-W46 S62-S55
  

    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