Obtaining the values within the boundary of my map

조회 수: 6 (최근 30일)
Benyamin Nemati
Benyamin Nemati 2020년 10월 14일
편집: Benyamin Nemati 2020년 10월 27일
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');
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
KSSV 2020년 10월 14일
inpolygon is your function..go for it.
Benyamin Nemati
Benyamin Nemati 2020년 10월 15일
how can i use this function on my Z values (EWT(:,:,i)) as xq and yq ?
by considering my code, what is your suggestion?

댓글을 달려면 로그인하십시오.

채택된 답변

KSSV
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개)

카테고리

Help CenterFile Exchange에서 Elementary Polygons에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by