How to add map to an output

Question

Sangesh Pv 2023-9-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2027249-how-to-add-map-to-an-output

评论： Sangesh Pv 2023-9-30

采纳的回答： Dyuman Joshi

I am trying to add the output to a map but i can't figure out how to add it.

7 个评论
显示 5更早的评论隐藏 5更早的评论

Sangesh Pv 2023-9-29

在 MATLAB Online 中打开

subtowersutm.xlsx

data = readmatrix('subtowersutm.xlsx'); 
% Extract latitude, longitude, and RSRP values.
measurement_lat = data(:, 1);
measurement_lon = data(:, 2);
rsrp = data(:, 11); % Adjust the column index for RSRP data.
% Check if the dimensions of Lat and Lon match
if size(measurement_lat) ~= size(measurement_lon)
    error('Latitude and Longitude dimensions do not match.');
end
% Define the UTM zone for your area.
utm_zone = 32; 
% Step 2: Transform Coordinates using the built-in 'deg2utm' function
[utm_x, utm_y, utm_zone] = deg2utm(measurement_lat, measurement_lon);
Unrecognized function or variable 'deg2utm'.
% Define the pixel size and create the grid
pixel_size = 20; % Adjust as needed
x_grid = min(utm_x):pixel_size:max(utm_x);
y_grid = min(utm_y):pixel_size:max(utm_y);
% Step 3: Calculate average RSRP for each pixel
num_pixels_x = numel(x_grid) - 1;
num_pixels_y = numel(y_grid) - 1;
average_rsrp = zeros(num_pixels_y, num_pixels_x); % Initialize the grid.
for i = 1:num_pixels_x
    for j = 1:num_pixels_y
        % Define the current pixel polygon.
        polygon_x = [x_grid(i), x_grid(i + 1), x_grid(i + 1), x_grid(i)];
        polygon_y = [y_grid(j), y_grid(j), y_grid(j + 1), y_grid(j + 1)];
        
        % Check if measurements fall within the current pixel.
        in_polygon = inpolygon(utm_x, utm_y, polygon_x, polygon_y);
        
        % Calculate average RSRP for the measurements within the polygon.
        if any(in_polygon)
            rsrp_values_in_polygon = rsrp(in_polygon); % Replace with your RSRP data.
            % Convert dBm to watts, compute average, and convert back to dBm.
            rsrp_watts = 10 .^ (rsrp_values_in_polygon / 10);
            average_rsrp(j, i) = 10 * log10(mean(rsrp_watts));
        end
    end
end
% Step 4: Plot the heatmap of average RSRP
% You can use various plotting functions to visualize the heatmap.
heatmap(x_grid(1:end-1), y_grid(1:end-1), average_rsrp);
colorbar;
xlabel('UTM X');
ylabel('UTM Y');
title('Average RSRP Heatmap');

Dyuman Joshi 2023-9-29

在 MATLAB Online 中打开

subtowersutm.xlsx

What's the problem here?

What do you want to obtain? Can you provide an example as to what the expected output should be?

data = readmatrix('subtowersutm.xlsx');

% Extract latitude, longitude, and RSRP values.

measurement_lat = data(:, 1);

measurement_lon = data(:, 2);

rsrp = data(:, 11); % Adjust the column index for RSRP data.

% Check if the dimensions of Lat and Lon match

if size(measurement_lat) ~= size(measurement_lon)

error('Latitude and Longitude dimensions do not match.');

end

% Define the UTM zone for your area.

utm_zone = 32;

% Step 2: Transform Coordinates using the built-in 'deg2utm' function

[utm_x, utm_y, utm_zone] = deg2utm(measurement_lat, measurement_lon);

% Define the pixel size and create the grid

pixel_size = 20; % Adjust as needed

x_grid = min(utm_x):pixel_size:max(utm_x);

y_grid = min(utm_y):pixel_size:max(utm_y);

% Step 3: Calculate average RSRP for each pixel

num_pixels_x = numel(x_grid) - 1;

num_pixels_y = numel(y_grid) - 1;

average_rsrp = zeros(num_pixels_y, num_pixels_x); % Initialize the grid.

for i = 1:num_pixels_x

for j = 1:num_pixels_y

% Define the current pixel polygon.

polygon_x = [x_grid(i), x_grid(i + 1), x_grid(i + 1), x_grid(i)];

polygon_y = [y_grid(j), y_grid(j), y_grid(j + 1), y_grid(j + 1)];

% Check if measurements fall within the current pixel.

in_polygon = inpolygon(utm_x, utm_y, polygon_x, polygon_y);

% Calculate average RSRP for the measurements within the polygon.

if any(in_polygon)

rsrp_values_in_polygon = rsrp(in_polygon); % Replace with your RSRP data.

% Convert dBm to watts, compute average, and convert back to dBm.

rsrp_watts = 10 .^ (rsrp_values_in_polygon / 10);

average_rsrp(j, i) = 10 * log10(mean(rsrp_watts));

end

% Step 4: Plot the heatmap of average RSRP

% You can use various plotting functions to visualize the heatmap.

heatmap(x_grid(1:end-1), y_grid(1:end-1), average_rsrp);

colorbar;

xlabel('UTM X');

ylabel('UTM Y');

title('Average RSRP Heatmap');

Walter Roberson 2023-9-29

deg2utm appears to be from https://www.mathworks.com/matlabcentral/fileexchange/10915-deg2utm

Sangesh Pv 2023-9-29

在 MATLAB Online 中打开

% Load your data
data = readmatrix('subtowersutm.xlsx');
% Extract latitude, longitude, and RSRP values.
measurement_lat = data(:, 1);
measurement_lon = data(:, 2);
rsrp = data(:, 11); % Adjust the column index for RSRP data.
% Check if the dimensions of Lat and Lon match
if size(measurement_lat) ~= size(measurement_lon)
    error('Latitude and Longitude dimensions do not match.');
end
% Define the UTM zone for your area.
utm_zone = 32;
% Step 2: Transform Coordinates using the built-in 'deg2utm' function
[utm_x, utm_y, utm_zone] = deg2utm(measurement_lat, measurement_lon);
% Define the pixel size and create the grid
pixel_size = 20; % Adjust as needed
x_grid = min(utm_x):pixel_size:max(utm_x);
y_grid = min(utm_y):pixel_size:max(utm_y);
% Step 3: Calculate average RSRP for each pixel
num_pixels_x = numel(x_grid) - 1;
num_pixels_y = numel(y_grid) - 1;
average_rsrp = zeros(num_pixels_y, num_pixels_x); % Initialize the grid.
for i = 1:num_pixels_x
    for j = 1:num_pixels_y
        % Define the current pixel polygon.
        polygon_x = [x_grid(i), x_grid(i + 1), x_grid(i + 1), x_grid(i)];
        polygon_y = [y_grid(j), y_grid(j), y_grid(j + 1), y_grid(j + 1)];
        
        % Check if measurements fall within the current pixel.
        in_polygon = inpolygon(utm_x, utm_y, polygon_x, polygon_y);
        
        % Calculate average RSRP for the measurements within the polygon.
        if any(in_polygon)
            rsrp_values_in_polygon = rsrp(in_polygon); % Replace with your RSRP data.
            % Convert dBm to watts, compute average, and convert back to dBm.
            rsrp_watts = 10 .^ (rsrp_values_in_polygon / 10);
            average_rsrp(j, i) = 10 * log10(mean(rsrp_watts));
        end
    end
end
% Exclude cells with 0 value from the heatmap
average_rsrp(average_rsrp == 0) = NaN;
% Define the color for NaN values (white)
nanColor = [1, 1, 1]; % RGB color for white
% Plot the heatmap of average RSRP with NaN values represented as white
colormap('hot'); % Use the 'hot' colormap for heat colors
heatmap(x_grid(1:end-1), y_grid(1:end-1), average_rsrp, ...
    'MissingDataColor', nanColor); % Set the MissingDataColor property
colorbar;
xlabel('UTM X');
ylabel('UTM Y');
title('Average RSRP Heatmap (Excluding 0 Values)');
% updated code

请先登录，再进行评论。

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

Answer 1

Dyuman Joshi 2023-9-30

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2027249-how-to-add-map-to-an-output#answer_1322629

在 MATLAB Online 中打开

When using heatmap, colormap has to be specified in the call. See below -

% Load your data
data = readmatrix('subtowersutm.xlsx');
% Extract latitude, longitude, and RSRP values.
measurement_lat = data(:, 1);
measurement_lon = data(:, 2);
rsrp = data(:, 11); % Adjust the column index for RSRP data.
% Check if the dimensions of Lat and Lon match
if size(measurement_lat) ~= size(measurement_lon)
    error('Latitude and Longitude dimensions do not match.');
end
% Define the UTM zone for your area.
utm_zone = 32;
% Step 2: Transform Coordinates using the built-in 'deg2utm' function
[utm_x, utm_y, utm_zone] = deg2utm(measurement_lat, measurement_lon);
% Define the pixel size and create the grid
pixel_size = 20; % Adjust as needed
x_grid = min(utm_x):pixel_size:max(utm_x);
y_grid = min(utm_y):pixel_size:max(utm_y);
% Step 3: Calculate average RSRP for each pixel
num_pixels_x = numel(x_grid) - 1;
num_pixels_y = numel(y_grid) - 1;
average_rsrp = zeros(num_pixels_y, num_pixels_x); % Initialize the grid.
for i = 1:num_pixels_x
    for j = 1:num_pixels_y
        % Define the current pixel polygon.
        polygon_x = [x_grid(i), x_grid(i + 1), x_grid(i + 1), x_grid(i)];
        polygon_y = [y_grid(j), y_grid(j), y_grid(j + 1), y_grid(j + 1)];
        
        % Check if measurements fall within the current pixel.
        in_polygon = inpolygon(utm_x, utm_y, polygon_x, polygon_y);
        
        % Calculate average RSRP for the measurements within the polygon.
        if any(in_polygon)
            rsrp_values_in_polygon = rsrp(in_polygon); % Replace with your RSRP data.
            % Convert dBm to watts, compute average, and convert back to dBm.
            rsrp_watts = 10 .^ (rsrp_values_in_polygon / 10);
            average_rsrp(j, i) = 10 * log10(mean(rsrp_watts));
        end
    end
end
% Exclude cells with 0 value from the heatmap
average_rsrp(average_rsrp == 0) = NaN;
% Define the color for NaN values (white)
nanColor = [1, 1, 1]; % RGB color for white

Colormap specified inside the heatmap() call.

% Plot the heatmap of average RSRP with NaN values represented as white

% Set the MissingDataColor property and the colormap

heatmap(x_grid(1:end-1), y_grid(1:end-1), average_rsrp, ...

'MissingDataColor', nanColor, Colormap=hot);

%% ^^^^^^^^^^^^

colorbar;

xlabel('UTM X');

ylabel('UTM Y');

title('Average RSRP Heatmap (Excluding 0 Values)');

3 个评论
显示 1更早的评论隐藏 1更早的评论

Dyuman Joshi 2023-9-30

Which map are you trying to plot? Could you attach the code for it?

Sangesh Pv 2023-9-30

在 MATLAB Online 中打开

% Create a web map using the Mapping Toolbox
webmap;
% Display the map using a web map service (you may need an internet connection)
basemap = wmsfind('nasa', 'SearchField', 'serverurl');
wmslayer = wmsinfo(basemap(1).ServerURL);
layers = wmslayer(1).LayerName;
addWebMap(webmap, layers);

% I tried with this code but it gives as 2 seperate output i dont think so it works, do you think it will be better to convert this utm coordinates back to lat lon so that geoscatter or geoshow command can be used ?

% I just need it to be on top of the streets which the data has been taken with and not a white background in the grided view

请先登录，再进行评论。