Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1744 1539 (N26W63 L=  4)  1233    136        2 (1691/1482-1802/1571) -18.4 W57-W69 N22-N28
      2 1181  832 (S14W08 L=309)   555    214        0 (1114/ 748-1238/ 915) -55.0 W04-W11 S19-S09
      3 1856  876 (S15W61 L=  2)   414     49        0 (1837/ 846-1876/ 907) -53.1 W58-W64 S16-S12
      4 1632  842 (S16W39 L=340)   604    157        0 (1574/ 778-1688/ 899)  41.2 W35-W43 S19-S12
      5  634  702 (S19E28 L=273)   714    284        1 ( 577/ 590- 705/ 800)  42.7 E34-E21 S25-S13
      6  518  747 (S15E35 L=266)   357     54        2 ( 486/ 725- 551/ 772)   0.0 E37-E32 S16-S13
      7  245  627 (S20E65 L=236)   654     56        1 ( 223/ 591- 271/ 659) -90.0 E67-E61 S22-S18
      8  982  544 (S31E06 L=294)   308     51        0 ( 956/ 523-1011/ 562) -18.4 E07-E04 S32-S30
      9  821  289 (S49E26 L=274)   314     46        1 ( 791/ 268- 849/ 305)   9.5 E29-E22 S51-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