Resolution Reduction (Image Processing)

// Function for resolution reduction (useful for transforming/handling large images)
// Original code written by Peter Rehbein, Goethe-University, Frankfurt, Germany
Function Shrink(image, shrinkfactor)
 
Wave image
Variable shrinkfactor 
 
if(ShrinkFactor<=0)
    DoAlert /T="Alert" 0,"ShrinkFactor must be greater than 0!"
return -1
endif
 
Duplicate /O image adjusted_image
print "ShrinkFactor:"+num2str(shrinkfactor)
print "====="
Variable height = dimsize(adjusted_image,0)
Variable width = dimsize(adjusted_image,1)
 
print "Height (original):"+num2str(height)
print "Width (original):"+num2str(width)
print "====="
// Check whether image dimensions are dividable by shrink factor:
 
Variable check_height
Variable check_width
 
Variable new_height
Variable new_width
 
Variable u,v
Variable d,e
 
//Loop while height is not dividable by shrinkfactor and enlarge picture appropriately
//(New rows are filled with copied adjacent data)
check_height = mod(height,shrinkfactor)
 
if(check_height!=0)
print "Height not dividable by "+num2str(shrinkfactor)+"!"
    print "Starting correction!"
    print "====="
for(u=0;check_height!=0;u+=1)
    print "Actual Height: "+num2str(height)
    new_height = height+1
    Redimension/N=(new_height,-1) adjusted_image
        for(e=0;e<width;e+=1)
            //print "e: "+num2str(e)
            adjusted_image [new_height-1] [e] = adjusted_image [new_height-2] [e] 
        endfor
    height = new_height
    print "New Height: "+num2str(new_height)
    check_height = mod(height,shrinkfactor)
endfor
endif
// end check
    
//Loop while width is not dividable by shrinkfactor and enlarge picture appropriately
//(New columns are filled with copied adjacent data)
check_width = mod(width,shrinkfactor)
if(check_width!=0)
print "Width not dividable by "+num2str(shrinkfactor)+"!"
    print "Starting correction!"
    print "====="
for(v=0;check_width!=0;v+=1)
    print "Actual Width: "+num2str(width)
    new_width = width+1
    Redimension/N=(-1,new_width) adjusted_image
        for(d=0;d<height;d+=1)
            //print "d: "+num2str(d)
            adjusted_image [d] [new_width-1] = adjusted_image [d] [new_width-2] 
        endfor
    width = new_width
    print "New Width: "+num2str(new_width)
    check_width = mod(width,shrinkfactor)
endfor
endif
// end check
 
Variable sImage_Height = Height/ShrinkFactor
Variable sImage_Width = Width/ShrinkFactor
print "====="
print "SmallImage Height :"+num2str(sImage_Height)
print "SmallImage Width: "+num2str(sImage_Width)
 
Make /O /N=(sImage_Height, sImage_Width) sImage
 
Variable i,j,k,l
 
Variable avg
 
//print "====Report===="
for(i=0;i<sImage_Height;i+=1) //correct
    for(j=0;j<sImage_Width;j+=1) //correct      
    
        k=i*ShrinkFactor
        l=j*ShrinkFactor
                
        //print "===="
        //print "k:"+num2str(k)
        //print "l:"+num2str(l)
                    
        Wave adjusted_image
        Duplicate/O/R=[k,k+ShrinkFactor-1][l,l+ShrinkFactor-1] adjusted_image,region // Extract region
        WaveStats /Q Region // Determine average pixel value
        avg = V_avg 
        //print "i: "+num2str(i)
        //print "j: "+num2str(j)
        //print "Average (region): "+num2str(avg)
        //print "i: "+num2str(i)
        //print "j: "+num2str(j)
        //print "===="
        sImage [i] [j] = avg // write average pixel value into smaller image    
    endfor
endfor
 
End