How to extract data through the shapefile (India) from global precipitation data.

Question

varunesh on 1 Jun 2018

0
Link

Direct link to this question

https://in.mathworks.com/matlabcentral/answers/403555-how-to-extract-data-through-the-shapefile-india-from-global-precipitation-data

Edited: varunesh on 6 Jun 2018

global_data.mat

Hello Guys,

I have shapefile for Indian region. Actually, I want to extract Indian region data through this shapefile from global data. Please help me for doing this. I am attaching here both (Shaplefile and Global data in .mat).

Thank you so much in advance.

Regards,

Varunesh

0 Comments
Show -2 older commentsHide -2 older comments

Sign in to comment.

Sign in to answer this question.

Answer 1

KSSV on 1 Jun 2018

2
Link

Direct link to this answer

https://in.mathworks.com/matlabcentral/answers/403555-how-to-extract-data-through-the-shapefile-india-from-global-precipitation-data#answer_322796

YOu need to read about inpolygon. This function helps you to extract what you want.

shapefile = 'India_Shapefile\ind_state_2000.shp' ;
S = load('global_data.mat') ;
C = shaperead(shapefile) ;
N = length(C) ;
lon = S.lon ; lat = S.lat ;
P = S.prc ;
[X,Y] = meshgrid(lon,lat) ;
iwant = cell(N,2) ;
for i =1:N
      Cx = C(i).X ; Cy = C(i).Y ;
      idx = inpolygon(X(:),Y(:),Cx',Cy') ;
      P1 = P(:,:,1)' ;
      iwant{i,1} = P1(idx) ;
      P2 = P(:,:,2)' ;
      iwant{i,2} = P2(idx) ;
  end