Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1271 1817 (N57E13 L=353)  1518    168        5 (1186/1775-1351/1850)   9.0 E21-E02 N53-N61
      2  652 1667 (N24E45 L=320)  2630    410        6 ( 607/1512- 725/1835) -32.5 E68-E32 N14-N35
      3 1683 1616 (N57W43 L= 49)   898     80        1 (1660/1568-1700/1653)  48.0 W37-W55 N52-N63
      4 1738 1219 (N27W43 L= 48)   336     77        0 (1706/1183-1762/1242) -21.8 W39-W45 N25-N28
      5  350 1007 (S22E44 L=322)   368     90        1 ( 310/ 990- 395/1024)   7.1 E47-E39 S23-S20
      6  543  863 (S26E26 L=340)   608    177        2 ( 492/ 788- 573/ 939) -41.2 E30-E22 S29-S22
      7 1898  896 (N14W68 L= 74)  1585    118        5 (1865/ 851-1929/ 934)  36.9 W63-W76 N10-N16
      8  175  894 (S32E64 L=301)   638     88        1 ( 145/ 865- 208/ 928)  24.4 E71-E58 S34-S29
      9  704  714 (S30E11 L=354)   254    115        0 ( 645/ 701- 751/ 736)  35.5 E15-E08 S32-S27
     10  754  627 (S35E06 L=360)   487    128        0 ( 687/ 605- 806/ 645)  12.5 E10-E01 S35-S33
     11 1394  572 (S19W36 L= 42)   372    120        0 (1334/ 557-1452/ 593)  35.5 W32-W39 S21-S16
     12 1421  272 (S35W62 L= 68)  1242    135        3 (1367/ 247-1476/ 305)   7.6 W53-W69 S37-S32
  

    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