function [VMPmax, VMPmin] = vmpGetMaxMin(path, vmp);

% syntax: [VMPmax, VMPmin] = vmpGetMaxMin(path, vmp)
% 
% this func returns the minimum and maximum of the vmp map
% in the file [path vmp '.vmp']
% 
% This function was written by:
% Hagar Gelbard
% Rafi Malach's Lab
% Weizmann Institute of Science
% Rehovot, Israel
% hagar.gelbard@weizmann.ac.il


% open the vmp file for reading
vmp_fid = fopen([path vmp '.vmp'],'r');

ver = fread(vmp_fid,1,'short');
NrOfMaps = fread(vmp_fid,1,'short');
MapType = fread(vmp_fid,1,'short');
NrOfLags = fread(vmp_fid,1,'short');
ClusterSize = fread(vmp_fid,1,'short');
CritThreshold = fread(vmp_fid,1,'float');
MaxThreshold = fread(vmp_fid,1,'float');
df1 = fread(vmp_fid,1,'short');
df2 = fread(vmp_fid,1,'short');
% read map name
a = 'a';
while a ~= 0
  a = fread(vmp_fid,1,'char');
end
DimX = fread(vmp_fid,1,'short');
DimY = fread(vmp_fid,1,'short');
DimZ = fread(vmp_fid,1,'short');
XStart = fread(vmp_fid,1,'short');
XEnd = fread(vmp_fid,1,'short');
YStart = fread(vmp_fid,1,'short');
YEnd = fread(vmp_fid,1,'short');
ZStart = fread(vmp_fid,1,'short');
ZEnd = fread(vmp_fid,1,'short');
resolution = fread(vmp_fid,1,'short');

cols = (XEnd-XStart+1)/resolution
rows = (YEnd-YStart+1)/resolution
slices = (ZEnd-ZStart+1)/resolution


VMPmat = fread(vmp_fid,cols*rows*slices,'float32');
VMPmax = max(max(max(VMPmat)));
VMPmin = min(min(min(VMPmat)));

fclose('all');