Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  725 1697 (N26E38 L=119)   678    168        3 ( 663/1652- 795/1743) -12.1 E45-E31 N24-N28
      2  381 1640 (N16E63 L= 94)  1211    236        4 ( 348/1589- 415/1714) -63.4 E66-E59 N12-N22
      3  463 1482 (N07E47 L=110)  1002    360        3 ( 341/1417- 591/1540)  29.4 E55-E37 N01-N10
      4  934 1448 (N15E16 L=141)   251     79        0 ( 894/1420- 972/1477)  -9.5 E18-E12 N13-N15
      5  199  739 (S41E57 L=100)   439     87        1 ( 174/ 725- 228/ 758)   9.5 E59-E53 S42-S39
      6 1766  626 (S08W62 L=219)   904    181        3 (1713/ 592-1823/ 668)  36.9 W57-W67 S11-S04
      7 1600  353 (S27W63 L=220)   595     86        0 (1559/ 338-1631/ 373)  38.7 W59-W65 S29-S25
      8  860  210 (S65W24 L=181)   625     77        3 ( 831/ 187- 894/ 240)  15.9 W18-W27 S67-S62
  

    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