Changeset 27038


Ignore:
Timestamp:
06/01/22 05:59:13 (3 years ago)
Author:
Mathieu Morlighem
Message:

CHG: added Bamber projection

Location:
issm/trunk-jpl/src/m/coordsystems
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/coordsystems/ll2xy.m

    r26958 r27038  
    4444        disp('Calling CoordTransform instead, make sure to change your MATLAB script');
    4545        if sgn==+1
    46                 assert(delta==45); assert(slat ==70);
    47                 [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3413');
    48                 return;
     46                if delta==45 && slat==70      %BedMachine
     47                        [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3413'); return;
     48                elseif delta==0 && slat==75   %IBCAO
     49                        [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3996'); return;
     50                elseif delta==39 && slat==71  %Bamber
     51                        Bamber_proj = '+proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs';
     52                        [x y]=CoordTransform(lat,lon,'EPSG:4326',Bamber_proj); return;
     53                else
     54                        error('not supported yet');
     55                end
    4956        else
    5057                assert(delta==0); assert(slat ==71);

  • issm/trunk-jpl/src/m/coordsystems/xy2ll.m

    r26958 r27038  
    4040        disp('Calling CoordTransform instead, make sure to change your MATLAB script');
    4141        if sgn==+1
    42                 assert(delta==45); assert(slat ==70);
    43                 [lat lon]=CoordTransform(x, y,'EPSG:3413','EPSG:4326');
    44                 return;
     42                if delta==45 && slat==70      %BedMachine
     43                        [lat lon]=CoordTransform(x, y,'EPSG:3413','EPSG:4326'); return;
     44                elseif delta==0 && slat==75   %IBCAO
     45                        [lat lon]=CoordTransform(x, y,'EPSG:3996','EPSG:4326'); return;
     46                elseif delta==39 && slat==71  %Bamber
     47                        Bamber_proj = '+proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs';
     48                        [lat lon]=CoordTransform(x, y, Bamber_proj, 'EPSG:4326'); return;
     49                else
     50                        error('not supported yet');
     51                end
    4552        else
    4653                assert(delta==0); assert(slat ==71);
Note: See TracChangeset for help on using the changeset viewer.