Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1766 1490 (N27W63 L= 14)   636    227        0 (1741/1453-1807/1509)  12.5 W59-W68 N25-N28
      2 1441 1398 (N25W31 L=342)   268     83        1 (1400/1384-1476/1419)  -9.5 W28-W34 N24-N26
      3 1631 1300 (N17W44 L=355)   984    231        4 (1567/1246-1693/1350) -54.5 W38-W48 N13-N20
      4 1852 1296 (N14W67 L= 18)   846    109        1 (1827/1247-1870/1339) -80.5 W64-W70 N11-N17
      5 1886  955 (S07W67 L= 17)  1131    270        2 (1836/ 913-1937/ 980) -16.4 W59-W76 S10-S04
      6 1565  755 (S16W35 L=345)   893    309        1 (1460/ 644-1658/ 894)  47.0 W27-W42 S22-S07
      7  180  829 (S06E67 L=243)   772    138        1 ( 160/ 798- 210/ 884) -51.3 E70-E62 S08-S02
      8 1656  563 (S30W45 L=356)  1431    260        5 (1558/ 529-1726/ 628)  21.0 W37-W51 S32-S26
      9  479  534 (S25E43 L=268)   739    174        4 ( 416/ 495- 548/ 583) -33.7 E47-E37 S27-S21
  

    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