Defining the central longitude in Arctic borders in the Arctic Mapping Tool

Question

Ashfaq Ahmed 2023-12-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2056674-defining-the-central-longitude-in-arctic-borders-in-the-arctic-mapping-tool

评论： Ashfaq Ahmed 2023-12-5

Hi! I am curretnly using ArcticMappingToolbox to plot the Arctic borders. But automatically it uses 45°W as the central longitude (based on the NSIDC stereographic system). Like in this figure -

But I want the 0* to be the central longitude like this -

Can anyone please help me with either a) how can I rotate the arcticborders.m function? or b) create a stereographic grid that has 0* as central longitude?

4 个评论
显示 2更早的评论隐藏 2更早的评论

Ashfaq Ahmed 2023-12-5

@Dyuman Joshi I meant, as you can see from the 2nd plot, the figure it 45° offset. I want to do that. I want to plot like Figure 2 in my figure 1.

Dyuman Joshi 2023-12-5

You mean the plot, not the function, ok.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Chad Greene 2023-12-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2056674-defining-the-central-longitude-in-arctic-borders-in-the-arctic-mapping-tool#answer_1366144

在 MATLAB Online 中打开

Try this. Into the Command Window type

open arcticborders

Then find the line where it loads the border data. Try adding this right after that line:

load('arcticborderdata.mat'); 
for k = 1:177
    % unproject: 
    [lat,lon] = psn2ll(x{k},y{k});
    % reproject: 
    [x{k},y{k}] = ll2psn(lat,lon,'meridian',0); 
end

2 个评论
显示无隐藏无

Ashfaq Ahmed 2023-12-5

在 MATLAB Online 中打开

Ok, I think I got it, it worked

%% Plot
for k = 1:177
    % unproject: 
    [lat,lon] = psn2ll(x{k},y{k});
    % reproject: 
    [x{k},y{k}] = ll2psn(lat,lon,'meridian',0); 
end
hold on
if faceplot
   n = 1; 
   for k = 1:177
      xk = x{k}; 
      yk = y{k}; 
       
      nanz = [0,find(isnan(xk))];
      for kk = 1:length(nanz)-1
         h(n) = patch(xk(nanz(kk)+1:nanz(kk+1)-1),yk(nanz(kk)+1:nanz(kk+1)-1),.5*[1 1 1],varargin{:}); 
         n=n+1; 
      end
   end
else
      h = plot(cell2mat(x),cell2mat(y),varargin{:}); 
end