Changeset 27529

Timestamp:

01/19/23 05:01:08 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: simplified epsg2proj

File:

• issm/trunk-jpl/src/m/coordsystems/epsg2proj.m (modified) (1 diff)

issm/trunk-jpl/src/m/coordsystems/epsg2proj.m

@@ -10 +10 @@
 %      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
+%Call PROJ library
+[status, string]=system(['projinfo -o PROJ -q epsg:' num2str(epsg)]);
 
-        version_major=str2num(r);
+%Check status
+if status~=0; error(string); end
 
-        [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);
+%remove trailing  blanks
+string = deblank(string);