Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  971 1749 (N49W10 L=164)   331     34        0 ( 953/1722- 992/1769)  36.9 W08-W12 N46-N50
      2  849 1565 (N36E04 L=151)   288     81        0 ( 806/1552- 894/1577) -26.6 E06-W00 N34-N37
      3  733 1537 (N36E12 L=142)   325     86        0 ( 701/1497- 768/1574)  39.8 E15-E09 N33-N38
      4 1301  725 (S22W13 L=167)  1178    390        0 (1128/ 667-1434/ 829)   5.4 W01-W22 S24-S16
      5  401  640 (S14E48 L=106)   700     87        1 ( 363/ 597- 428/ 687)  45.0 E51-E44 S16-S11
      6  829  584 (S24E19 L=135)   516    183        0 ( 802/ 496- 871/ 661)  45.0 E23-E14 S28-S19
      7  215  552 (S17E77 L= 77)  1562     85        0 ( 193/ 507- 234/ 591)  39.8 E82-E73 S19-S14
      8  481  413 (S29E53 L=101)   811    114        0 ( 456/ 365- 505/ 468)  69.4 E56-E50 S33-S25
      9  570  372 (S34E48 L=106)  1047    147        0 ( 507/ 335- 614/ 427)  14.9 E55-E40 S35-S30
  

    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