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

조회 수: 2 (최근 30일)
Ashfaq Ahmed 2023년 12월 5일
댓글: 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.

댓글을 달려면 로그인하십시오.

채택된 답변

Chad Greene 2023년 12월 5일
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:
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일
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
Ashfaq Ahmed 2023년 12월 5일
Thank you so much!!!

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Contour Plots에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by