Drawing Tidal Ellipses on an m_map in MATLAB

21 次查看（过去 30 天）
은진
은진 2024-11-19，16:12
I want to draw tidal ellipses on an m_map, but the figure is not displaying correctly. Please help me.
% 사용자 지정 위도 및 경도 (덕적도)
latitude = 37 + 13/60 + 35/3600; % N 37° 13' 35"
longitude = 126 + 9/60 + 24/3600; % E 126° 09' 24"
% 데이터 파일 로드
fnm = 'cmems_mod_glo_phy_my_0.083deg_P1D-m_1731894708544.nc';
lon = ncread(fnm, 'longitude');
lat = ncread(fnm, 'latitude');
time = ncread(fnm, 'time');
dateT = datetime(1970, 01, 01, 00, 00, 00) + seconds(time);
u = ncread(fnm, 'uo');
v = ncread(fnm, 'vo');
% 지정된 위치에 가장 가까운 격자 점 찾기
[~, lat_idx] = min(abs(lat - latitude));
[~, lon_idx] = min(abs(lon - longitude));
% 해당 지점의 유속 데이터 추출
u_series = squeeze(u(lon_idx, lat_idx, :));
v_series = squeeze(v(lon_idx, lat_idx, :));
% 조석 성분 분석 (t_tide)
[NAME, FREQ, TIDECON, ~] = t_tide(u_series + 1i * v_series, ...
'interval', 1, 'start time', datenum(dateT(1)), 'latitude', latitude);
% 주요 조석 성분 설정
major_constituents = {'M2', 'S2', 'K1', 'O1'};
colors = {'r', 'g', 'b', 'k'}; % 각각 빨강, 초록, 파랑, 검정
NAME = strtrim(string(NAME)); % 조석 성분 이름 정리
% 지도 설정 (m_map)
figure;
m_proj('mercator', 'lon', [120 128], 'lat', [33 40]); % 범위 설정 (한반도)
m_gshhs_l('patch', [0.8 0.8 0.8]); % 육지 표시
m_grid('box', 'fancy', 'tickdir', 'out'); % 격자 추가
hold on;
% 주요 조석 성분별로 타원 그리기
for i = 1:length(major_constituents)
% 해당 조석 성분의 인덱스 찾기
idx = find(strcmp(NAME, major_constituents{i}));
if ~isempty(idx)
% 타원 매개변수 추출
major_axis = TIDECON(idx, 1); % 장축 반경 (확대)
minor_axis = TIDECON(idx, 3); % 단축 반경 (확대)
inclination = TIDECON(idx, 5) * pi / 180; % 경사각 (라디안)
x0 = longitude; % 타원의 중심 경도
y0 = latitude; % 타원의 중심 위도
% 타원 그리기
ellipsedraw(major_axis, minor_axis, x0, y0, inclination, colors{i});
end
end
% 제목 및 기타 설정
title('Tidal Ellipses for Major Constituents');
hold off;

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

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Oceanography and Hydrology 的更多信息

标签

产品


版本

R2024a

Community Treasure Hunt

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

Start Hunting!

Translated by