Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  769 1534 (N16E29 L=321)   491     70        1 ( 749/1495- 790/1573) -68.2 E30-E26 N13-N18
      2 1514 1327 (N23W22 L= 12)  1062    334        2 (1453/1179-1569/1457) -53.7 W15-W28 N15-N30
      3  651 1020 (S16E22 L=328)   368    156        0 ( 585/ 975- 717/1056)  45.0 E25-E18 S20-S12
      4 1022  751 (S23W07 L=357)   267     70        0 ( 987/ 729-1058/ 771)  53.1 W05-W09 S24-S20
      5  728  700 (S33E10 L=340)   271     66        0 ( 698/ 667- 761/ 729) -21.8 E12-E07 S34-S31
      6 1518  551 (S19W46 L= 35)   346    131        0 (1465/ 526-1554/ 596)  63.4 W42-W47 S21-S15
      7  461  427 (S61E36 L=314)   695     94        0 ( 434/ 396- 498/ 471) -71.6 E41-E30 S63-S56
      8  630  360 (S61E12 L=338)   354     51        2 ( 603/ 339- 656/ 390) -15.9 E15-E07 S62-S58
  

    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