Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  816 1463 (N38E01 L=159)   912    228        4 ( 723/1397- 905/1504)   0.0 E08-W06 N35-N39
      2 1896 1031 (S19W61 L=220)  1234    134        2 (1878/ 968-1925/1091)  50.2 W57-W66 S23-S15
      3  594 1063 (N20E25 L=135)   208     88        0 ( 549/1054- 636/1076) -18.4 E27-E21 N18-N20
      4 1635  892 (S18W33 L=193)   967    297        2 (1541/ 776-1706/1001)  25.2 W24-W41 S21-S13
      5 1399  886 (S11W18 L=177)   221     72        1 (1379/ 840-1416/ 916)  63.4 W15-W19 S12-S08
      6  322  795 (N10E52 L=108)   307     50        0 ( 303/ 769- 341/ 825) -63.4 E52-E50 N07-N11
      7  861  619 (S12E20 L=140)  2024    368        5 ( 793/ 437- 926/ 798)  68.6 E24-E15 S24-S01
  

    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