Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1597 1544 (N19W56 L=232)  2561    497        2 (1437/1405-1777/1624) -51.0 W45-W66 N06-N27
      2 1194 1413 (N24W21 L=198)   378    117        0 (1131/1389-1243/1435) -45.0 W18-W24 N22-N27
      3 1616 1415 (N12W49 L=226)   311     83        0 (1567/1404-1655/1425) -23.2 W45-W52 N10-N13
      4  692 1344 (N34E13 L=164)   455    215        0 ( 614/1281- 772/1399) -48.0 E18-E08 N28-N39
      5  234 1028 (N24E57 L=119)  1329    260        2 ( 196/ 942- 298/1159) -61.9 E62-E51 N17-N32
      6  339 1039 (N23E46 L=131)   283     80        0 ( 318/1004- 366/1076) -51.3 E47-E43 N20-N25
      7 1832  999 (S18W55 L=231)   484     78        0 (1810/ 960-1865/1032)  14.0 W51-W59 S19-S16
      8 1532  961 (S10W28 L=205)   367     82        0 (1510/ 921-1553/1006) -90.0 W27-W29 S12-S06
      9  226  823 (N12E62 L=115)  1017    197        1 ( 167/ 747- 263/ 898) -48.4 E68-E59 N07-N16
     10  188  747 (N07E71 L=106)  1649    156        4 ( 156/ 669- 216/ 809) -83.7 E76-E65 N02-N11
     11 1706  551 (S44W40 L=216)   507     55        0 (1684/ 534-1733/ 574)   0.0 W36-W43 S45-S40
  

    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