popup menu help

Imagestats is better adapted to your problem than Wavestats. You can directly enter the section coordinates like this
ImageStats/GS={xmin, xmax, ymin, ymax} m_interpolatedimage
There is no need to copy a subset of the image to a temporary wave. Note that the results are NOT printed in the history window as is the case for WaveStats. Rather, they are returned as global variables (if you execute ImageStats from the command line) or as local variables (if called from a function) named V_avg, V_min, etc.

If you want to "automate" this, I guess there are two ways to go about it. I have written some quick sample code below.

1) If you want to operate on many images and if the section coordinates are always the same, preparing a table with those coordinates (as you did) seems the best solution. You can then write a function that operates on those coordinates and puts the ImageStats results in a wave.

I have called my sample function "ImageStatsFromWave". In my code, the ImageStats results are appended as new rows to your "SectionX" waves. If you prepare a two-dimensional waves (say, "AllSections") by executing
Make/O/N=(4,5) AllSections; Edit AllSections
on the command line, and by copy+pasting the contents of your table-of-1D-waves into the new table, you'll have to call the function only once:
ImageStatsFromWave(AllSections, m_interpolatedimage)
Otherwise, you'll have to call it for each of your 1D waves "Section1", "Section2", etc.

2) If you want to be more "interactive", you might want to mark the region of interest on your image with the mouse and then select a context menu item (by clicking in the Marquee rectangle) that lets you directly print the ImageStats results to the history window. This can be done using Igor's "GetMarquee" function. The code below consists of two parts, the function ImageStatsFromMarquee() and the menu definition Menu "GraphMarquee" (see Igor's help files for more on that).

// version 1: sections (xmin..ymax) are in a normal (1D) wave or in a 2D wave
//  the stats results are appended as new rows
function ImageStatsFromWave(sections, imageWave)
wave sections, imageWave
    Redimension/N=(8,-1) sections
    variable k, kmax = DimSize(sections,1)
    if( kmax == 0 ) // 1D wave
        kmax = 1    
    endif
    for( k=0; k < kmax; k+=1 )
        ImageStats/GS={sections[0][k], sections[1][k], sections[2][k], sections[3][k]} imageWave
        sections[4][k] = V_avg
        sections[5][k] = V_min
        sections[6][k] = V_max
        sections[7][k] = V_sdev
    endfor
end
 
// version 2: section is taken from a graph marquee
//  the stats results are printed into the history window
function ImageStatsFromMarquee()
    wave imageWave = $(StringFromList(0, ImageNameList("", ";")))   // get topmost image
    if( !WaveExists(imageWave) )
        print "ImageStatsFromMarquee works only on images!"
        return -1
    endif
    string axes = ReplaceString(";",RemoveEnding(AxisList("")),",") // get axes' names as comma-separated list
    Execute/Q/Z "GetMarquee/K "+axes    // have to execute since axis names (not strings) are expected
    if( V_flag )
        print "ImageStatsFromMarquee error: no marquee?"    // should never happen
        return -2
    endif
    NVAR V_left, V_top, V_right, V_bottom
    ImageStats/GS={V_left, V_right, V_top, V_bottom} imageWave
    if( V_flag )
        print "ImageStatsFromMarquee error: illegal coordinates?"   // may happen if marquee larger than image
        return -3
    endif
    string format = "%s(%g,%g)(%g,%g): avg: %g; min: %g; max: %g; sdev: %g\r"
    printf format, NameOfWave(imageWave), V_left, V_right, V_top, V_bottom, V_avg, V_min, V_max, V_sdev
end
 
Menu "GraphMarquee"
    "Print Local Image Statistics", /Q, ImageStatsFromMarquee()
End

If you copy the text above into your Procedure window and compile that, you can experiment with the sample code...

HTH,
Wolfgang

July 17, 2010 at 12:11 pm - Permalink