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

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-2);
        end

        string=(r);