フィルターのクリア

how to generate correct logical mask of the size of shape file?

2 ビュー (過去 30 日間)
kuldeep
kuldeep 2024 年 9 月 25 日
コメント済み: kuldeep 2024 年 9 月 26 日
i am using following code to generate logical mask
function croppedimg = cropimage(data, R, shapefile)
% Load the shapefile and retrieve its info
S = shaperead(shapefile);
% Shapefile projections
info = shapeinfo(shapefile);
crs = info.CoordinateReferenceSystem;
[S.lon,S.lat] = projinv(crs,S.X,S.Y);
% read and reproject image file
p = R.ProjectedCRS;
[x,y] = worldGrid(R);
[lon,lat] = projinv(p,x,y);
% Create a logical mask
% Approach 1: Using 'inpolygon' function
xq = lon(:);
yq = lat(:);
in = inpolygon(xq,yq,S.lon,S.lat);
logical_mask = reshape(in,size(lon));
% Debugging: Display the logical mask summary
disp(['Total number of pixels inside polygon: ', num2str(numel(logical_mask))]);
disp(['Total number of pixels inside image: ', num2str(numel(data))]);
However it is showing the
Total number of pixels inside polygon: 16982229
Total number of pixels inside image: 16982229
please suggest me how to create correct logical mask of the size of shape file.
kuldeep

回答 (1 件)

DGM
DGM 2024 年 9 月 25 日
編集済み: DGM 2024 年 9 月 25 日
I don't have everything to run your example, but the last lines are a problem at the least.a
Use numel() to return the total number of elements regardless of their value.
Use nnz() if you want the number of true elements in a logical array.
A = magic(5);
mask = A>15
mask = 5x5 logical array
1 1 0 0 0 1 0 0 0 1 0 0 0 1 1 0 0 1 1 0 0 1 1 0 0
numel(mask)
ans = 25
nnz(mask)
ans = 10
  5 件のコメント
DGM
DGM 2024 年 9 月 26 日
If you're trying to crop out blobs from a logical mask, then you can use the 'Image' property from regionprops(), or you can use the 'BoundingBox' property and pass it to imcrop() instead.
That said, there's obviously something wrong because none of your query points are in the polygon. Without the inputs to this code, it's not apparent why that is.
kuldeep
kuldeep 2024 年 9 月 26 日
Thanks for your help. I am attaching the shape file and data file.

サインインしてコメントする。

製品


リリース

R2024a

Community Treasure Hunt

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

Start Hunting!

Translated by