Random area of no data

1 次查看（过去 30 天）
Sophie
Sophie 2024-8-28
回答： Star Strider 2024-8-28
Hi,
I am creating a world plot to display data. However, there is a random line appearing in my plot and i dont know why. The figure didn't have this line and then i ran the code again and it randomly appeared and hasnt gone away since.
Is it something with my code?
clear;
clc;
data=load('date.mat', "temp_data_2030_245", "salinity_data_2030_245", "CO3_data_2030_245");
data1=load('245.mat',"co3_lat_245", "co3_lon_245");
Omegalimit = 2.82;
% Import data from files
temp_data = data.temp_data_2030_245;
salinity_data = data.salinity_data_2030_245;
CO3_data = data.CO3_data_2030_245;
lat=data1.co3_lat_245;
lon=data1.co3_lon_245;
CO3_data = CO3_data/1025; %convert to molkg^-1
%Ca
Ca_values = (salinity_data/35)*0.01028;
temp_data1=temp_data;
T = temp_data + 273.15;
S = salinity_data;
Ksp_values = -171.945 - (0.077993 .* T) + (2903.293 ./ T) + (71.595 .* log10(T)) ...
+ ((0.0017276 .* T - 0.068393 + (88.135 ./ T)) .* sqrt(S)) ...
- (0.10018 .* S) + (0.0059415 .* S.^1.5);
%Omega
omega = (CO3_data) .* Ca_values ./ 10.^Ksp_values;
omega = mean(omega, 3);
directory = '/Users/sophie/Documents/Masters/Research Project/01_Data';
fileName = 'WCMC008_CoralReef2021_Py_v4_1.shp';
filePath = fullfile(directory, fileName);
CurrentCoral = shaperead(filePath);
X = [CurrentCoral.X];
Y = [CurrentCoral.Y];
figure(3);
sgtitle('Aragonite Saturation State under SSP 245 in Summer');
subplot(3,1,1);
worldmap('World');
set(gcf, 'Color', 'w');
land = shaperead('landareas', 'UseGeoCoords', true);
geoshow(lat, lon, omega, 'DisplayType', 'texturemap');
mlabel off; plabel off; gridm off
hold on;
geoshow(land, 'DisplayType', 'polygon', 'FaceColor', 'k');
title('2030');
xlabel('Longitude');
ylabel('Latitude');
colormap('turbo');
set(gcf, 'Position', [100, 100, 800, 600]);
hold on;
plotm(Y, X, 'ro', 'MarkerSize', 1, 'MarkerFaceColor', 'r');

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

回答（1 个）

Star Strider
Star Strider 2024-8-28
I do not have your data, and I am not certain what ‘randomly’ implies.. Consider usinig the fillmissing function if the data creating the meridian line at about 50°E are actually missing data.
You did not list your MATLAB version/release. Be sure that you have downloaded and installed all the updates for your particluar version. This could be a bug that has been fixed.

类别

在 Help Center 和 File Exchange 中查找有关 Data Preprocessing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by