function [xmin xmax ymin ymax] = compute_and_display_shot_count_histogram(shotspermeter,histogram_breaks, histogram_max)
%given a grid of lidar shots per square meter, the number of total
%histogram breaks, and the maximum shot count to show in the histogram
%range, compute and plot a histogram of points per square meter
%
%Ramon Arrowsmith, May 2010

mn=size(shotspermeter);
len=mn(1)*mn(2);
shotsforhist=reshape(shotspermeter,len,1);
hist(shotsforhist,histogram_breaks)
xlabel('Points per sq. meter')
ylabel('Frequency')
axis([0 histogram_max -inf inf])
mode(shotsforhist)