function string = epsg2proj(epsg)
%EPSG2PROJ - uses gdalsrsinfo to provide PROJ.4 compatible string from 
%EPSG code
%
%   Usage:
%      proj4string = epsg2proj(4326);
%
%   Example:
%      proj4string = epsg2proj(4326);
%      return proj4string='+proj=longlat +datum=wgs84 +no_defs'
%
	% First, get GDAL version
	[s,r]=system(['gdalsrsinfo --version | awk ''{print $2}'' | cut -d ''.'' -f1']);

	if s~=0, 
		error(r);
	end

	version_major=str2num(r);

	[s,r]=system(['gdalsrsinfo epsg:' num2str(epsg) ' | command grep PROJ.4 | tr -d ''\n'' | sed ''s/PROJ.4 : //''']);

	if s~=0, 
		error(r);
	end

	if version_major==1,
		r=r(2:end-2);
	end

	string=(r);
