source: issm/trunk-jpl/src/m/coordsystems/epsg2proj.m@ 25003

function string = epsg2proj(epsg)
%use 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 '
%

        [s,r]=system(['gdalsrsinfo  epsg:' num2str(epsg)  ' | grep PROJ.4 | sed "s/PROJ.4 : //"']);
        r=r(2:end-2);

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