Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1131 1586 (N40W04 L=340)   438     59        2 (1084/1566-1167/1604)  15.9 W00-W07 N38-N41
      2  656 1476 (N28E30 L=306)   725    217        0 ( 576/1435- 767/1515) -10.0 E37-E20 N26-N30
      3 1645 1388 (N29W45 L= 21)  1190    246        0 (1568/1301-1713/1478) -47.7 W40-W50 N23-N34
      4 1427 1421 (N30W27 L=  3)   343     89        0 (1388/1388-1459/1444) -21.8 W23-W29 N27-N31
      5  464 1328 (N17E41 L=295)  1421    318        0 ( 371/1246- 567/1433) -22.2 E53-E31 N12-N23
      6  946 1307 (N19E07 L=329)   725    251        1 ( 847/1238-1033/1372) -26.6 E14-W00 N15-N22
      7  769  730 (S18E15 L=321)   292     72        1 ( 729/ 708- 800/ 750)  11.3 E17-E12 S19-S16
      8 1430  691 (S16W29 L=  5)   371     81        1 (1407/ 648-1462/ 724)  68.2 W27-W31 S19-S13
      9  288  543 (S35E65 L=271)   737     85        0 ( 246/ 524- 327/ 570)   9.9 E79-E56 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