Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1270 1818 (N57E13 L=353)  1488    177        3 (1185/1778-1349/1845)  12.5 E21-E03 N53-N61
      2  648 1668 (N24E45 L=321)  2593    385        4 ( 611/1535- 716/1835) -27.2 E68-E33 N15-N35
      3 1680 1618 (N57W43 L= 49)   837     80        2 (1658/1569-1697/1653)  35.5 W37-W54 N52-N63
      4  538  866 (S26E26 L=340)   639    169        2 ( 487/ 792- 569/ 943) -33.7 E31-E22 S29-S22
      5 1889  882 (N13W67 L= 74)  1653    137        5 (1836/ 841-1928/ 933)  31.6 W60-W74 N08-N16
      6  167  888 (S33E66 L=300)   674     98        1 ( 142/ 845- 194/ 920)  28.6 E72-E61 S36-S30
      7  927  752 (S22W02 L=  8)   195     70        0 ( 912/ 719- 947/ 784) -45.0 W00-W03 S23-S20
      8  698  716 (S30E12 L=354)   213     99        0 ( 650/ 704- 740/ 737)  38.7 E14-E09 S32-S28
      9  759  628 (S34E05 L=  1)   338     62        0 ( 726/ 609- 795/ 646)   0.0 E07-E02 S35-S33
     10 1389  574 (S19W36 L= 42)   369    118        0 (1330/ 560-1447/ 594)  39.8 W32-W39 S21-S16
     11 1423  273 (S35W62 L= 68)  1282    117        3 (1369/ 249-1476/ 304)  20.0 W54-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