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

function string = epsg2proj(epsg)
%FUNCTION 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) " | grep PROJ.4 | tr -d '\n' | sed 's/PROJ.4 : //'"]);

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

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

        string=(r);