Changeset 25035

Timestamp:

06/12/20 20:03:09 (5 years ago)

Author:

jdquinn

Message:

CHG: Saving changes for MATLAB -> Python API for solid earth

Location:

issm/trunk-jpl/src/m

Files:

• classes/sealevelmodel.py (added)
• coordsystems/epsg2proj.m (modified) (1 diff)
• coordsystems/epsg2proj.py (added)
• coordsystems/gdaltransform.m (modified) (1 diff)
• coordsystems/gdaltransform.py (added)
• mesh/TwoDToThreeD.py (added)
• mesh/meshintersect3d.py (added)

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

@@ -1 +1 @@
 function string = epsg2proj(epsg)
-%use gdalsrsinfo to provide Proj.4 compatible string from EPSG code
+%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 '
+%   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 : //"']);
+        [s,r]=system(['gdalsrsinfo epsg:' num2str(epsg) ' | grep PROJ.4 | sed "s/PROJ.4 : //"']);
 
         if s~=0, 

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

@@ -2 +2 @@
 %GDALTRANSFORM - switch from one projection system to another 
 %
-%   Usage:
-%      [x,y] = gdaltransform(x1,y1,epsg_in, epsg_out);
+%       Usage:
+%               [x,y] = gdaltransform(x1,y1,epsg_in, epsg_out);
 %
-%   Example: 
-%      [x,y] = gdaltransform(md.mesh.long,md.mesh.lat,'EPSG:4326','EPSG:3031')
+%       Example:
+%               [x,y] = gdaltransform(md.mesh.long,md.mesh.lat,'EPSG:4326','EPSG:3031')
 %
-%   For reference: 
-%      EPSG: 4326 (lat,long)
-%      EPSG: 3411 (Greenland,  UPS 45W, 70N)
-%      EPSG: 3031 (Antarctica, UPS 0E,  71S)
+%       For reference:
+%               EPSG: 4326 (lat,long)
+%               EPSG: 3411 (Greenland,  UPS 45W, 70N)
+%               EPSG: 3031 (Antarctica, UPS 0E,  71S)
 %
-%      ll2xy default projection Antarctica:
-%        +proj=stere +lat_0=-90 +lat_ts=-71 +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378273 +b=6356889.448564109 +units=m +no_defs
-%      ll2xy default projection Greenland:
-%        +proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +a=6378273 +b=6356889.448564109 +units=m +no_defs
-%      Bamber's Greeland projection
-%        +proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39 +k=1 +x_0=0 +y_0=0 +a=6378273 +b=6356889.448564109 +units=m +no_defs
+%               ll2xy default projection Antarctica:
+%                       +proj=stere +lat_0=-90 +lat_ts=-71 +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378273 +b=6356889.448564109 +units=m +no_defs
+%               ll2xy default projection Greenland:
+%                       +proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +a=6378273 +b=6356889.448564109 +units=m +no_defs
+%               Bamber's Greeland projection
+%                       +proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39 +k=1 +x_0=0 +y_0=0 +a=6378273 +b=6356889.448564109 +units=m +no_defs
 %
-%      To get proj.4 string from EPSG, use gdalsrsinfo. Example:    gdalsrsinfo "EPSG:4326" | grep "PROJ.4"
+%               To get proj.4 string from EPSG, use gdalsrsinfo. Example:
+%
+%                       gdalsrsinfo epsg:4326 | grep "PROJ.4" | sed "s/PROJ.4 : //"
 
         %give ourselves unique file names