Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  898 1761 (N40E33 L=124)   292     35        1 ( 882/1736- 916/1782) -63.4 E34-E31 N37-N41
      2  998 1677 (N36E21 L=136)  1919    248        6 ( 916/1570-1064/1765) -27.8 E30-E11 N30-N41
      3  549 1714 (N28E56 L=100)   407     68        0 ( 527/1681- 571/1744) -21.8 E61-E51 N25-N29
      4 1612 1597 (N50W33 L=190)  1361    167        3 (1552/1554-1691/1637)  -3.4 W25-W43 N47-N52
      5  431 1319 (N01E44 L=113)  1522    220        4 ( 392/1223- 485/1405) -45.0 E48-E38 S04-N05
      6 1301 1324 (N21W09 L=166)   515    173        0 (1232/1260-1354/1379) -21.8 W03-W13 N18-N22
      7  556  959 (S17E27 L=130)   237     27        1 ( 538/ 934- 577/ 982) -71.6 E28-E25 S18-S15
      8 1052  823 (S14W06 L=163)   990    316        1 ( 914/ 753-1190/ 870)  42.7 W00-W13 S20-S08
      9 1506  792 (S05W34 L=191)   317    103        0 (1471/ 744-1526/ 833) -38.7 W30-W36 S07-S02
     10 1180  596 (S24W21 L=177)  1120    248        2 (1074/ 576-1315/ 634)  18.4 W13-W28 S26-S20
     11  470  565 (S45E31 L=126)   369     58        1 ( 441/ 544- 496/ 589)  38.7 E33-E27 S46-S42
  

    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