Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  811 1805 (N46E30 L=260)  1085     86        3 ( 771/1777- 868/1836) -17.1 E36-E23 N43-N48
      2 1173 1698 (N41W03 L=293)   419    103        0 (1127/1676-1225/1717)  15.9 W00-W07 N38-N42
      3 1628 1616 (N42W45 L=335)   719    123        1 (1579/1590-1684/1638)   0.0 W38-W51 N39-N43
      4  486 1474 (N18E41 L=249)  1486    276        2 ( 401/1383- 588/1579) -26.6 E53-E31 N12-N25
      5 1087 1512 (N26E01 L=289)   303    157        1 (1054/1482-1122/1555) -59.0 E03-W01 N23-N28
      6 1390 1252 (N12W20 L=310)   503    275        0 (1323/1172-1471/1372) -45.0 W15-W26 N08-N19
      7 1740 1276 (N19W47 L=337)   739    179        0 (1670/1234-1819/1321)  -8.1 W41-W55 N17-N20
      8 1698 1149 (N10W43 L=333)   318     45        0 (1672/1128-1726/1171) -14.0 W41-W45 N09-N11
      9  693  723 (S26E19 L=271)  1224    377        1 ( 597/ 572- 809/ 869)  54.5 E27-E11 S37-S16
     10 1790  594 (S21W67 L=357)  1149     76        2 (1754/ 556-1833/ 623)   0.0 W61-W75 S22-S19
  

    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