Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  737 1630 (N33E33 L= 17)  1819    487        7 ( 608/1532- 911/1724) -11.3 E47-E17 N28-N37
      2 1173 1592 (N38E00 L= 50)   302     72        0 (1130/1573-1211/1605) -11.3 E03-W02 N37-N39
      3  739 1480 (N22E27 L= 23)   299     81        0 ( 700/1454- 770/1498)   0.0 E29-E24 N21-N23
      4 1426 1467 (N34W20 L= 71)   297     81        1 (1384/1438-1466/1484)  -9.5 W17-W23 N32-N35
      5 1796 1455 (N41W56 L=107)   857    256        1 (1754/1438-1840/1466)  15.3 W51-W62 N39-N42
      6 1157 1370 (N23W02 L= 53)   209     64        1 (1129/1336-1190/1394) -21.8 W00-W05 N21-N23
      7  265 1182 (S03E56 L=354)   327     88        0 ( 236/1167- 302/1195)  26.6 E59-E52 S04-S01
      8  956  527 (S33W05 L= 55)   299     58        0 ( 930/ 505- 989/ 556) -21.8 W02-W07 S33-S30
  

    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