Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1266 1820 (N57E13 L=353)  1585    171        3 (1182/1777-1343/1849)   9.5 E22-E03 N53-N61
      2  641 1673 (N24E47 L=320)  2944    407        7 ( 598/1537- 715/1842) -24.8 E72-E33 N15-N35
      3 1678 1620 (N57W42 L= 49)   854     77        1 (1656/1572-1696/1655)  38.2 W37-W54 N52-N63
      4 1728 1187 (N25W42 L= 49)   649    138        1 (1696/1127-1749/1250) -54.5 W38-W44 N21-N28
      5  534  863 (S26E26 L=340)   685    184        2 ( 483/ 789- 564/ 937) -41.2 E31-E23 S29-S22
      6 1900  900 (N14W69 L= 75)  1346     81        3 (1879/ 861-1927/ 936)  50.2 W64-W74 N11-N17
      7  921  755 (S22W01 L=  8)   237     59        0 ( 905/ 721- 942/ 786) -45.0 W00-W03 S23-S20
      8  689  718 (S31E12 L=354)   219     95        0 ( 643/ 706- 730/ 739)  38.7 E15-E10 S32-S28
      9 1064  642 (S25W13 L= 20)   368    126        0 (1041/ 592-1101/ 700) -35.5 W10-W17 S26-S21
     10  759  629 (S34E05 L=  1)   340     82        0 ( 721/ 613- 802/ 649)   9.5 E08-E02 S35-S33
     11 1381  577 (S19W35 L= 42)   390    121        0 (1323/ 562-1440/ 595)  35.5 W31-W38 S21-S16
     12 1426  274 (S35W62 L= 69)  1134    110        2 (1383/ 252-1474/ 306)  39.8 W56-W68 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