Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  732 1616 (N20E36 L=134)   902    235        3 ( 677/1553- 810/1669) -15.9 E42-E28 N18-N22
      2  714 1572 (N17E35 L=135)   277     53        1 ( 690/1554- 740/1594) -26.6 E36-E32 N15-N18
      3  313 1568 (N08E70 L=101)   599     84        0 ( 286/1553- 336/1588)  21.8 E72-E66 N07-N09
      4 1815 1161 (N27W48 L=218)   911    168        1 (1770/1092-1857/1236) -26.6 W42-W53 N23-N30
      5 1035  814 (S17W07 L=177)   710    179        1 ( 942/ 763-1123/ 871)  61.4 W02-W09 S22-S11
      6 1732  738 (N00W54 L=224)  1091    213        3 (1700/ 636-1759/ 831) -59.0 W50-W57 S05-N05
      7 1240  623 (S21W25 L=195)   704    250        2 (1144/ 569-1341/ 697)   0.0 W16-W32 S22-S19
      8 1696  589 (S08W58 L=228)   614     83        1 (1659/ 562-1725/ 617)  18.4 W54-W60 S09-S06
      9 1469  539 (S19W43 L=213)   260     80        3 (1448/ 514-1486/ 586) -78.7 W41-W45 S20-S15
     10  720  506 (S45E03 L=167)  1411    355        5 ( 556/ 474- 878/ 538)  25.6 E17-W08 S52-S40
     11 1433  343 (S30W52 L=223)  1228    473        1 (1291/ 299-1571/ 414)  13.5 W38-W63 S33-S25
  

    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