Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1373  810 (S17W21 L=311)   947    370        1 (1272/ 705-1454/ 956) -53.7 W14-W26 S23-S07
      2 1762  853 (S16W51 L=341)   615    139        1 (1707/ 800-1813/ 898)  29.1 W46-W55 S18-S13
      3  445  828 (S10E39 L=251)   259     55        0 ( 418/ 809- 475/ 852) -26.6 E40-E36 S11-S08
      4  704  772 (S15E22 L=268)   441    141        2 ( 645/ 743- 758/ 801)  14.0 E25-E17 S17-S13
      5  767  700 (S20E18 L=271)   478    227        2 ( 725/ 624- 812/ 788)  55.0 E21-E14 S24-S14
      6  196  682 (S17E69 L=221)  1107    139        4 ( 178/ 617- 223/ 746)  63.4 E75-E64 S20-S12
      7  335  653 (S20E53 L=237)   478     86        0 ( 315/ 608- 353/ 687)  78.7 E54-E50 S22-S17
      8 1159  563 (S31W06 L=296)   442    100        1 (1114/ 540-1209/ 584)  -9.5 W03-W09 S32-S29
      9  647  288 (S48E44 L=246)   337     45        1 ( 627/ 266- 666/ 318)  68.2 E46-E41 S50-S45
     10  923  292 (S50E16 L=274)   335     82        0 ( 899/ 274- 963/ 317)   0.0 E18-E11 S52-S48
  

    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