I am creating a polygon using coordinate points with geoplot, and I would like to fill the resultant shape in a partially transparent blue color. I am able to create the filled shape, and I can overlay the outline on the map with geoplot, but I cannot do both.

%% Plot

% Create a figure
figure;

% Create a geographic axes
gx = geoaxes;
hold(gx, ‘on’);

% Define the indices of visible points
visibleIndices = maxRadiusPerAngle > 0;

% Close the loop by adding the first point at the end
latTargetsClosed = [latTargets(visibleIndices) latTargets(visibleIndices(1))];
lonTargetsClosed = [lonTargets(visibleIndices) lonTargets(visibleIndices(1))];

% Plot the outline of the polygon using geoplot
geoplot(gx, latTargetsClosed, lonTargetsClosed, ‘b-‘, ‘LineWidth’, 2);

% Plot the central point
geoplot(gx, latCenter, lonCenter, ‘ro’, ‘MarkerSize’, 8, ‘MarkerFaceColor’, ‘r’);

% Plot the visible points
geoplot(gx, latTargets(visibleIndices), lonTargets(visibleIndices), ‘go’, ‘MarkerSize’, 6, ‘MarkerFaceColor’, ‘g’);

% Set the basemap
geobasemap(gx, ‘streets’);
title(‘Furthest Outward Points with Positive LOS’);

% Define the Web Mercator projection
projCRS = projcrs(3857); % EPSG:3857 Web Mercator

% Convert geographic coordinates to projected coordinates
[x, y] = projfwd(projCRS, latTargetsClosed, lonTargetsClosed);

% Create a new Cartesian plot for the filled polygon, ensuring it overlays the geographic axes
figure;
hold on;

% Create a filled polygon in the Cartesian plot
fill(x, y, ‘b’, ‘FaceAlpha’, 0.3, ‘EdgeColor’, ‘none’);

% Set Cartesian axes limits to match the geographic plot
axis equal;
xlabel(‘X (meters)’);
ylabel(‘Y (meters)’);
title(‘Filled Polygon in Projected Coordinates’);

% Adjust the limits to fit the geographic plot
latlim = gx.LatitudeLimits;
lonlim = gx.LongitudeLimits;
[xlim, ylim] = projfwd(projCRS, [latlim(1), latlim(2)], [lonlim(1), lonlim(2)]);
set(gca, ‘XLim’, xlim, ‘YLim’, ylim);

% Ensure the Cartesian plot is visible
set(gcf, ‘Visible’, ‘on’);I am creating a polygon using coordinate points with geoplot, and I would like to fill the resultant shape in a partially transparent blue color. I am able to create the filled shape, and I can overlay the outline on the map with geoplot, but I cannot do both.

%% Plot

% Create a figure
figure;

% Create a geographic axes
gx = geoaxes;
hold(gx, ‘on’);

% Define the indices of visible points
visibleIndices = maxRadiusPerAngle > 0;

% Close the loop by adding the first point at the end
latTargetsClosed = [latTargets(visibleIndices) latTargets(visibleIndices(1))];
lonTargetsClosed = [lonTargets(visibleIndices) lonTargets(visibleIndices(1))];

% Plot the outline of the polygon using geoplot
geoplot(gx, latTargetsClosed, lonTargetsClosed, ‘b-‘, ‘LineWidth’, 2);

% Plot the central point
geoplot(gx, latCenter, lonCenter, ‘ro’, ‘MarkerSize’, 8, ‘MarkerFaceColor’, ‘r’);

% Plot the visible points
geoplot(gx, latTargets(visibleIndices), lonTargets(visibleIndices), ‘go’, ‘MarkerSize’, 6, ‘MarkerFaceColor’, ‘g’);

% Set the basemap
geobasemap(gx, ‘streets’);
title(‘Furthest Outward Points with Positive LOS’);

% Define the Web Mercator projection
projCRS = projcrs(3857); % EPSG:3857 Web Mercator

% Convert geographic coordinates to projected coordinates
[x, y] = projfwd(projCRS, latTargetsClosed, lonTargetsClosed);

% Create a new Cartesian plot for the filled polygon, ensuring it overlays the geographic axes
figure;
hold on;

% Create a filled polygon in the Cartesian plot
fill(x, y, ‘b’, ‘FaceAlpha’, 0.3, ‘EdgeColor’, ‘none’);

% Set Cartesian axes limits to match the geographic plot
axis equal;
xlabel(‘X (meters)’);
ylabel(‘Y (meters)’);
title(‘Filled Polygon in Projected Coordinates’);

% Adjust the limits to fit the geographic plot
latlim = gx.LatitudeLimits;
lonlim = gx.LongitudeLimits;
[xlim, ylim] = projfwd(projCRS, [latlim(1), latlim(2)], [lonlim(1), lonlim(2)]);
set(gca, ‘XLim’, xlim, ‘YLim’, ylim);

% Ensure the Cartesian plot is visible
set(gcf, ‘Visible’, ‘on’); I am creating a polygon using coordinate points with geoplot, and I would like to fill the resultant shape in a partially transparent blue color. I am able to create the filled shape, and I can overlay the outline on the map with geoplot, but I cannot do both.

%% Plot

% Create a figure
figure;

% Create a geographic axes
gx = geoaxes;
hold(gx, ‘on’);

% Define the indices of visible points
visibleIndices = maxRadiusPerAngle > 0;

% Close the loop by adding the first point at the end
latTargetsClosed = [latTargets(visibleIndices) latTargets(visibleIndices(1))];
lonTargetsClosed = [lonTargets(visibleIndices) lonTargets(visibleIndices(1))];

% Plot the outline of the polygon using geoplot
geoplot(gx, latTargetsClosed, lonTargetsClosed, ‘b-‘, ‘LineWidth’, 2);

% Plot the central point
geoplot(gx, latCenter, lonCenter, ‘ro’, ‘MarkerSize’, 8, ‘MarkerFaceColor’, ‘r’);

% Plot the visible points
geoplot(gx, latTargets(visibleIndices), lonTargets(visibleIndices), ‘go’, ‘MarkerSize’, 6, ‘MarkerFaceColor’, ‘g’);

% Set the basemap
geobasemap(gx, ‘streets’);
title(‘Furthest Outward Points with Positive LOS’);

% Define the Web Mercator projection
projCRS = projcrs(3857); % EPSG:3857 Web Mercator

% Convert geographic coordinates to projected coordinates
[x, y] = projfwd(projCRS, latTargetsClosed, lonTargetsClosed);

% Create a new Cartesian plot for the filled polygon, ensuring it overlays the geographic axes
figure;
hold on;

% Create a filled polygon in the Cartesian plot
fill(x, y, ‘b’, ‘FaceAlpha’, 0.3, ‘EdgeColor’, ‘none’);

% Set Cartesian axes limits to match the geographic plot
axis equal;
xlabel(‘X (meters)’);
ylabel(‘Y (meters)’);
title(‘Filled Polygon in Projected Coordinates’);

% Adjust the limits to fit the geographic plot
latlim = gx.LatitudeLimits;
lonlim = gx.LongitudeLimits;
[xlim, ylim] = projfwd(projCRS, [latlim(1), latlim(2)], [lonlim(1), lonlim(2)]);
set(gca, ‘XLim’, xlim, ‘YLim’, ylim);

% Ensure the Cartesian plot is visible
set(gcf, ‘Visible’, ‘on’); geoplot, mapping toolbox MATLAB Answers — New Questions

​

I wrote a Python script and got a path selection at a certain point as [‘Satellite_185’, ‘Satellite_113’, ‘Satellite_162’, ‘Satellite_108’, ‘Satellite_87’, ‘Satellite_68’, ‘Satellite_305’, ‘Satellite_335’, ‘Satellite_384’]. Then I wrote the following MATLAB code, but it gets stuck at the simulation scene step. I suspect it’s because of the time issue, but I can’t fix it.
% Create Satellite Scenario
startTime = datetime(2021,12,10,18,27,57); % 10 December 2021, 6:27:57 PM UTC
stopTime = startTime + hours(3); % 10 December 2021, 9:27:57 PM UTC
sampleTime = 60; % Seconds
sc = satelliteScenario(startTime,stopTime,sampleTime,"AutoSimulate",false);

% Load the entire constellation, but we’ll only use the specified relay satellites
satellites = satellite(sc,"largeConstellation2.tle");

% Known relay satellite names or identifiers (assuming you have them)
% If you only have indices, you’ll need to map them to satellite objects somehow
% Here we assume you have the satellite objects directly

% relaySats = {…}; % This should be filled with the actual satellite objects
% Example: relaySats = {satellites(10), satellites(20), satellites(30)};
% But since you likely have names or other identifiers, you’ll need to find them

% Or, if you have indices into the ‘satellites’ array
relaySatsIndices = [185, 113, 162, 108, 87, 68, 305, 335, 384]; % Example indices
relaySats = satellites(relaySatsIndices);

% Add Ground Stations
gsSource = groundStation(sc,42.3001,-71.3504, "Name","Source Ground Station");
gsTarget = groundStation(sc,17.4351,78.3824, "Name","Target Ground Station");

% Build the path using the known relay satellite sequence
pathNodes = {gsSource, relaySats, gsTarget};

% Visualize Path
sc.AutoSimulate = true;
ac = access(pathNodes{:});
ac.LineColor = "red";

% Determine Intervals When Calculated Path Can Be Used (if needed)
% intvls = accessIntervals(ac); % Uncomment this line if you need the access intervals

% Create satellite scenario viewer
v = satelliteScenarioViewer(sc,"ShowDetails",false);
% Note: Setting MarkerSize directly on the satellites object may not work as expected.
% Instead, adjust the visualization properties in the viewer or use a loop to set markers.

% Play the scenario
play(sc);
The code above references part of the code from this website:https://www.mathworks.com/help/satcom/ug/multihop-path-select-through-sat-constellation.html

This picture is my stuck interface after running in matlab.I wrote a Python script and got a path selection at a certain point as [‘Satellite_185’, ‘Satellite_113’, ‘Satellite_162’, ‘Satellite_108’, ‘Satellite_87’, ‘Satellite_68’, ‘Satellite_305’, ‘Satellite_335’, ‘Satellite_384’]. Then I wrote the following MATLAB code, but it gets stuck at the simulation scene step. I suspect it’s because of the time issue, but I can’t fix it.
% Create Satellite Scenario
startTime = datetime(2021,12,10,18,27,57); % 10 December 2021, 6:27:57 PM UTC
stopTime = startTime + hours(3); % 10 December 2021, 9:27:57 PM UTC
sampleTime = 60; % Seconds
sc = satelliteScenario(startTime,stopTime,sampleTime,"AutoSimulate",false);

% Load the entire constellation, but we’ll only use the specified relay satellites
satellites = satellite(sc,"largeConstellation2.tle");

% Known relay satellite names or identifiers (assuming you have them)
% If you only have indices, you’ll need to map them to satellite objects somehow
% Here we assume you have the satellite objects directly

% relaySats = {…}; % This should be filled with the actual satellite objects
% Example: relaySats = {satellites(10), satellites(20), satellites(30)};
% But since you likely have names or other identifiers, you’ll need to find them

% Or, if you have indices into the ‘satellites’ array
relaySatsIndices = [185, 113, 162, 108, 87, 68, 305, 335, 384]; % Example indices
relaySats = satellites(relaySatsIndices);

% Add Ground Stations
gsSource = groundStation(sc,42.3001,-71.3504, "Name","Source Ground Station");
gsTarget = groundStation(sc,17.4351,78.3824, "Name","Target Ground Station");

% Build the path using the known relay satellite sequence
pathNodes = {gsSource, relaySats, gsTarget};

% Visualize Path
sc.AutoSimulate = true;
ac = access(pathNodes{:});
ac.LineColor = "red";

% Determine Intervals When Calculated Path Can Be Used (if needed)
% intvls = accessIntervals(ac); % Uncomment this line if you need the access intervals

% Create satellite scenario viewer
v = satelliteScenarioViewer(sc,"ShowDetails",false);
% Note: Setting MarkerSize directly on the satellites object may not work as expected.
% Instead, adjust the visualization properties in the viewer or use a loop to set markers.

% Play the scenario
play(sc);
The code above references part of the code from this website:https://www.mathworks.com/help/satcom/ug/multihop-path-select-through-sat-constellation.html

This picture is my stuck interface after running in matlab. I wrote a Python script and got a path selection at a certain point as [‘Satellite_185’, ‘Satellite_113’, ‘Satellite_162’, ‘Satellite_108’, ‘Satellite_87’, ‘Satellite_68’, ‘Satellite_305’, ‘Satellite_335’, ‘Satellite_384’]. Then I wrote the following MATLAB code, but it gets stuck at the simulation scene step. I suspect it’s because of the time issue, but I can’t fix it.
% Create Satellite Scenario
startTime = datetime(2021,12,10,18,27,57); % 10 December 2021, 6:27:57 PM UTC
stopTime = startTime + hours(3); % 10 December 2021, 9:27:57 PM UTC
sampleTime = 60; % Seconds
sc = satelliteScenario(startTime,stopTime,sampleTime,"AutoSimulate",false);

% Load the entire constellation, but we’ll only use the specified relay satellites
satellites = satellite(sc,"largeConstellation2.tle");

% Known relay satellite names or identifiers (assuming you have them)
% If you only have indices, you’ll need to map them to satellite objects somehow
% Here we assume you have the satellite objects directly

% relaySats = {…}; % This should be filled with the actual satellite objects
% Example: relaySats = {satellites(10), satellites(20), satellites(30)};
% But since you likely have names or other identifiers, you’ll need to find them

% Or, if you have indices into the ‘satellites’ array
relaySatsIndices = [185, 113, 162, 108, 87, 68, 305, 335, 384]; % Example indices
relaySats = satellites(relaySatsIndices);

% Add Ground Stations
gsSource = groundStation(sc,42.3001,-71.3504, "Name","Source Ground Station");
gsTarget = groundStation(sc,17.4351,78.3824, "Name","Target Ground Station");

% Build the path using the known relay satellite sequence
pathNodes = {gsSource, relaySats, gsTarget};

% Visualize Path
sc.AutoSimulate = true;
ac = access(pathNodes{:});
ac.LineColor = "red";

% Determine Intervals When Calculated Path Can Be Used (if needed)
% intvls = accessIntervals(ac); % Uncomment this line if you need the access intervals

% Create satellite scenario viewer
v = satelliteScenarioViewer(sc,"ShowDetails",false);
% Note: Setting MarkerSize directly on the satellites object may not work as expected.
% Instead, adjust the visualization properties in the viewer or use a loop to set markers.

% Play the scenario
play(sc);
The code above references part of the code from this website:https://www.mathworks.com/help/satcom/ug/multihop-path-select-through-sat-constellation.html

This picture is my stuck interface after running in matlab. large satellite constellation MATLAB Answers — New Questions

​

Simplied my problem so it’s easier to solve. Let’s pretend I have MatrixA that has 2 columns, and I want to sort the rows from highest to lowest depending on the value of the cell in the second column. How would I do this?

I also have a VectorB that starts at 10 and would like to add the value of each row in column 2 of the sorted MatrixA. How would I do this? Note this needs to come after MatrixA is sorted from highest to lowest.

Any suggestions matlab wizards out there? Really want to avoid using loops to make the code inefficient.

%Just creating example of the matrix I want to sort you can ignore this
MatrixA = cell(3,2);
MatrixA{1,1} = "D1";
MatrixA{2,1} = "D2";
MatrixA{3,1} = "D3";
MatrixA{1,2} = 5;
MatrixA{2,2} = 15;
MatrixA{3,2} = 10;

% Creating Vector B based on the already defined values for MatrixA which aren’t in order
VectorA = [10 10+MatrixA{1,2} 10+MatrixA{2,2} 10+MatrixA{3,2}];

% How would I sort MatrixA from highest to lowest here and the values of VectorB to be based on the sorted MatrixA?Simplied my problem so it’s easier to solve. Let’s pretend I have MatrixA that has 2 columns, and I want to sort the rows from highest to lowest depending on the value of the cell in the second column. How would I do this?

I also have a VectorB that starts at 10 and would like to add the value of each row in column 2 of the sorted MatrixA. How would I do this? Note this needs to come after MatrixA is sorted from highest to lowest.

Any suggestions matlab wizards out there? Really want to avoid using loops to make the code inefficient.

%Just creating example of the matrix I want to sort you can ignore this
MatrixA = cell(3,2);
MatrixA{1,1} = "D1";
MatrixA{2,1} = "D2";
MatrixA{3,1} = "D3";
MatrixA{1,2} = 5;
MatrixA{2,2} = 15;
MatrixA{3,2} = 10;

% Creating Vector B based on the already defined values for MatrixA which aren’t in order
VectorA = [10 10+MatrixA{1,2} 10+MatrixA{2,2} 10+MatrixA{3,2}];

% How would I sort MatrixA from highest to lowest here and the values of VectorB to be based on the sorted MatrixA? Simplied my problem so it’s easier to solve. Let’s pretend I have MatrixA that has 2 columns, and I want to sort the rows from highest to lowest depending on the value of the cell in the second column. How would I do this?

I also have a VectorB that starts at 10 and would like to add the value of each row in column 2 of the sorted MatrixA. How would I do this? Note this needs to come after MatrixA is sorted from highest to lowest.

Any suggestions matlab wizards out there? Really want to avoid using loops to make the code inefficient.

%Just creating example of the matrix I want to sort you can ignore this
MatrixA = cell(3,2);
MatrixA{1,1} = "D1";
MatrixA{2,1} = "D2";
MatrixA{3,1} = "D3";
MatrixA{1,2} = 5;
MatrixA{2,2} = 15;
MatrixA{3,2} = 10;

% Creating Vector B based on the already defined values for MatrixA which aren’t in order
VectorA = [10 10+MatrixA{1,2} 10+MatrixA{2,2} 10+MatrixA{3,2}];

% How would I sort MatrixA from highest to lowest here and the values of VectorB to be based on the sorted MatrixA? matlab, matlab code, mathematics, matrix array, matrix, matrix manipulation, matrices, array, arrays, cell array, cell arrays, vector, vectors, vectorization, cell, excel, importing excel data, matlab function, matlab coder, struct, script, speed, solve, sort, filter, fsolve MATLAB Answers — New Questions

​

Hi experts,
I want to calculate the power of the spatial frequency range 20-30 cycles per picture in an picture. I am not sure how I can accomplish this using MATLAB code. FYI. Cycles per picture is a spatial frequency unit, also called cycles per image.
I downloaded following code online and made some modifications. Could you check the code for me and tell me whether I am correct?

% 1. Read the image
img = imread(‘your_image.jpg’);
img_gray = rgb2gray(img); % Convert to grayscale if the image is colored

% 2. Compute the Fourier transform of the image
img_fft = fft2(double(img_gray));
img_fft_shifted = fftshift(img_fft);

% Get the size of the image
[rows, cols] = size(img_gray);

% Compute the frequency coordinates
u = (-rows/2:(rows/2-1)) / rows;
v = (-cols/2:(cols/2-1)) / cols;
[U, V] = meshgrid(u, v);
D = sqrt(U.^2 + V.^2);

% 3. Calculate the power in the 20-30 cycles per image range
freq_low = 20 / rows;
freq_high = 30 / rows;
mask = (D >= freq_low) & (D <= freq_high);
energy = sum(sum(abs(img_fft_shifted .* mask).^2));

% Output the power
fprintf(‘Power in the range 20-30 cycles per image: %fn’, energy);Hi experts,
I want to calculate the power of the spatial frequency range 20-30 cycles per picture in an picture. I am not sure how I can accomplish this using MATLAB code. FYI. Cycles per picture is a spatial frequency unit, also called cycles per image.
I downloaded following code online and made some modifications. Could you check the code for me and tell me whether I am correct?

% 1. Read the image
img = imread(‘your_image.jpg’);
img_gray = rgb2gray(img); % Convert to grayscale if the image is colored

% 2. Compute the Fourier transform of the image
img_fft = fft2(double(img_gray));
img_fft_shifted = fftshift(img_fft);

% Get the size of the image
[rows, cols] = size(img_gray);

% Compute the frequency coordinates
u = (-rows/2:(rows/2-1)) / rows;
v = (-cols/2:(cols/2-1)) / cols;
[U, V] = meshgrid(u, v);
D = sqrt(U.^2 + V.^2);

% 3. Calculate the power in the 20-30 cycles per image range
freq_low = 20 / rows;
freq_high = 30 / rows;
mask = (D >= freq_low) & (D <= freq_high);
energy = sum(sum(abs(img_fft_shifted .* mask).^2));

% Output the power
fprintf(‘Power in the range 20-30 cycles per image: %fn’, energy); Hi experts,
I want to calculate the power of the spatial frequency range 20-30 cycles per picture in an picture. I am not sure how I can accomplish this using MATLAB code. FYI. Cycles per picture is a spatial frequency unit, also called cycles per image.
I downloaded following code online and made some modifications. Could you check the code for me and tell me whether I am correct?

% 1. Read the image
img = imread(‘your_image.jpg’);
img_gray = rgb2gray(img); % Convert to grayscale if the image is colored

% 2. Compute the Fourier transform of the image
img_fft = fft2(double(img_gray));
img_fft_shifted = fftshift(img_fft);

% Get the size of the image
[rows, cols] = size(img_gray);

% Compute the frequency coordinates
u = (-rows/2:(rows/2-1)) / rows;
v = (-cols/2:(cols/2-1)) / cols;
[U, V] = meshgrid(u, v);
D = sqrt(U.^2 + V.^2);

% 3. Calculate the power in the 20-30 cycles per image range
freq_low = 20 / rows;
freq_high = 30 / rows;
mask = (D >= freq_low) & (D <= freq_high);
energy = sum(sum(abs(img_fft_shifted .* mask).^2));

% Output the power
fprintf(‘Power in the range 20-30 cycles per image: %fn’, energy); digital image processing, fft2, fftshift MATLAB Answers — New Questions

​