Obtaining the values within the boundary of my map
조회 수: 2 (최근 30일)
이전 댓글 표시
Hello everyone,
I have some plots that they're from a 180*360*12 matrix with a code like this:
for i = 1:12
load coast
figure(i)
hold on
imagesc(lambda-180,phi,EWT(:,:,i));
plot(long,lat,'k');
end
B_1 = shaperead('Aqueduct_river_basins_TIGRIS & EUPHRATES.shp'); %one of the shape files that ive loaded on my map
mapshow(B_1,'facecolor','none');
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/380651/image.jpeg)
you can see the boundaries on the map,the north africa and middle east region, now my question is: how could i derive EWT values within these boundaries, i want the average values on these areas at the end.
thanks for your considerations.
댓글 수: 3
채택된 답변
KSSV
2020년 10월 15일
Something like this:
[m,n,p] = size(EWT) ;
[X,Y] = meshgrid(lambda-180,phi) ;
B_1 = shaperead('Aqueduct_river_basins_TIGRIS & EUPHRATES.shp'); %one of the shape files that ive loaded on my map
N = length(B-1.X) ; % this should have different regions
iwant = cell(N,1) ;
for i = 1:N
% Get the indices insdie this polygon
idx = inpolygon(X,Y,B_1(i).X,B_i(i).Y) ;
% Get the mean of these points
thesum = 0 ;
for j = 1:p
T = EWT(:,:,i) ;
thesum = thesum+T(idx) ;
end
iwant{i} = thesum/p ;
end
The above code gives you inside the polygon points as scattered points. If you want grid, pick the chunk, use inpoltgon and make the outside indices data to NaN.
댓글 수: 0
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Elementary Polygons에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!