Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1766 1453 (N08W63 L=311)  1437    306        8 (1707/1369-1841/1518) -45.0 W58-W70 N01-N11
      2  125 1170 (N31E63 L=184)   615    132        1 ( 100/1145- 150/1190) -36.9 E66-E60 N29-N32
      3  708 1013 (N10E18 L=229)   253     85        1 ( 660/ 987- 742/1028) -11.3 E21-E16 N08-N11
      4  878  904 (N00E11 L=236)   636    181        5 ( 792/ 878- 956/ 936) -24.0 E15-E06 S02-N03
      5  279  525 (S09E68 L=179)  3714    383        5 ( 167/ 420- 344/ 638) -56.3 E82-E61 S16-S01
      6  567  484 (S17E43 L=204)   286    124        0 ( 524/ 462- 598/ 517)   0.0 E46-E40 S18-S15
      7  549  437 (S20E47 L=201)   301     49        1 ( 523/ 415- 571/ 453) -90.0 E48-E44 S21-S18
  

    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