how can i use utm2degfunchtion

Question

Sangesh Pv 2023-10-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2032614-how-can-i-use-utm2degfunchtion

评论： Sangesh Pv 2023-10-17

i am trying to convert the utm coordinates back to degree but i cant figure out. i get multiple errors when i use the function. i need to add this function

[Lat,Lon] = utm2deg(utm_x,utm_y,utm_zone) in the end so that the utm coordinates can be converted back to deg before plotting in the map

% 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 = 10; % 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
% Set the MissingDataColor property and the colormap
heatmap(x_grid(1:end-1), y_grid(1:end-1), average_rsrp, ...
    'MissingDataColor', nanColor, Colormap=parula);
colorbar;
xlabel('UTM X');
ylabel('UTM Y');
title('Average RSRP Heatmap (Excluding 0 Values)');

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

Steven Lord 2023-10-13

Do you receive warning and/or error messages? If so the full and exact text of those messages (all the text displayed in orange and/or red in the Command Window) may be useful in determining what's going on and how to avoid the warning and/or error.

请先登录，再进行评论。

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

Answer 1

Gyan Vaibhav 2023-10-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2032614-how-can-i-use-utm2degfunchtion#answer_1334949

在 MATLAB Online 中打开

Hi Sangesh, 

To convert UTM coordinates back to degrees, you can add the “utm2deg” function to your code. Here's how you can modify your code to include the “utm2deg” function:

% 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); 
% Step 3: Convert UTM coordinates back to degrees using the 'utm2deg' function 
[Lat, Lon] = utm2deg(utm_x, utm_y, utm_zone); 
% Define the pixel size and create the grid 
pixel_size = 10; % Adjust as needed 
x_grid = min(utm_x):pixel_size:max(utm_x); 
y_grid = min(utm_y):pixel_size:max(utm_y); 
% Step 4: 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 
heatmap(x_grid(1:end-1), y_grid(1:end-1), average_rsrp, ... 
    'MissingDataColor', nanColor, 'Colormap', parula); 
colorbar; 
xlabel('UTM X'); 
ylabel('UTM Y'); 
title('Average RSRP Heatmap (Excluding 0 Values)'); 
% Function to convert UTM coordinates to degrees 
function [Lat, Lon] = utm2deg(utm_x, utm_y, utm_zone) 
% Implement the conversion logic here 
% ... 
end 