How to find number of pixels and pixel ID based on Latitude and Longitude Information?
5 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Amjad Iqbal
am 8 Apr. 2024
Kommentiert: Amjad Iqbal
am 18 Apr. 2024
Dear MATLAB Research Community,
I have to find number of pixels and pixel IDs based on Latitude and Longitude Information given below.
Here is the link where I have uploaded the data for reference and code.
I have sorted the way to read the data which is for one day of scan and the way to read latitude and longitude information as well.
This file loop_read_1_day_data.m is containing code.
I have to find the pixel numbers and ID of that pixel from complete data for this range [Lat, long] = [ -72.345785°, 170.044070°] in decimal.
Then I have to find that ID in brightness image to find brightness temperature of that pixels.
2 Kommentare
Pratyush Swain
am 18 Apr. 2024
Hi Amjad,
Correct me if I am wrong but this is what i have understood from your query. You have to find the set of indices in form (row,col) of the values from latitude and longitude information which are the closest to this range [Lat, long] = [ -72.345785°, 170.044070°] and use the same set of indexes to find the brightness values. Is this your requirement ?
Akzeptierte Antwort
Pratyush Swain
am 18 Apr. 2024
Hi Amjad,
Please refer to this workflow which can help you in finding common rows and cols for pixel values satisfying the custom range conditions:
% Specifying latitude and longitude range as per your requirement %
lat_range = [-72.3250, -72.2000];
lt_range = [169.95, 170.05];
% Creating global arrays which will store the lat and long values
% satisfying range criteria and also the brightness values at their
% intersection points
all_lat_values = [];
all_long_values = [];
all_brightness_values = [];
% Iterate through all files
for i = 1:14
% Implementation for reading brightness image and lat & long
% ------------------------------------------------ %
% Implementation for resizing brightness img,lat and long matrix to
% common dimensions
% ------------------------------------------------ %
% Finding latitudes and longitudes rows & cols satisfying range
% criteria %
[row_lat, col_lat] = find(lati_resized > lat_range(1) & lati_resized < lat_range(2));
[row_long, col_long] = find(lt_resized > lt_range(1) & lt_resized < lt_range(2));
% Forming latitude and longitude indices array consisting of rows and
% cols values
lat_indices = [row_lat, col_lat];
long_indices = [row_long, col_long];
% Finding matching indices between lat and long indices array by using the intersect function %
% "Stable" parameter insures the matching indices come in ordered way %
[common_indices, ~] = intersect(lat_indices, long_indices, 'rows', 'stable');
% Segregating matched rows and matched columns %
matched_rows = common_indices(:, 1);
matched_cols = common_indices(:, 2);
% Forming linear indices using the sub2ind function %
linear_indices = sub2ind(size(B), matched_rows, matched_cols);
% These linear indices are applicable to all matrices since all were
% resized to common dimensions
lat_values = lati_resized(linear_indices);
long_values = lt_resized(linear_indices);
% Finally Using the indices to find brightness values %
matched_brightness_values = B(linear_indices);
% Append the values to the global arrays %
all_lat_values = [all_lat_values; lat_values];
all_long_values = [all_long_values; long_values];
all_brightness_values = [all_brightness_values; matched_brightness_values];
end
You can also visualize the findings using a scatter plot between latitude,longitude and brightness values:
figure;
scatter3(all_lat_values, all_long_values, all_brightness_values, 36, all_brightness_values, 'filled');
colorbar;
colormap(jet);
title('3D Plot of Brightness Temperatures vs. Latitude and Longitude');
xlabel('Longitude');
ylabel('Latitude');
zlabel('Brightness Temperature');
view(3);
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1672636/image.png)
Some of the functions I used above and their uses are as follow :
1- find() : Find all the indices satisfying a criteria in a matrix
2- intersect(): Finds data common to two arrays.
3- sub2ind(): Returns the linear indices corresponding to the row and column subscripts.
For more information on functions and their parameters, please refer to following resources:
I hope this example implementation helps you in your use case. I have also attached the code file for your reference. I also have tested it on the data you had provided.
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Geographic Plots finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!