How to Crop India shape file from Global shapefile

Hello:

I have two shapefiles: global shapefiles of landuse map and India administrative boundary shapefile. Is it possible to clip/crop the India region using shapefile of India from global shapefile?

The shapefile for India is available in this link: http://otlet.sims.berkeley.edu/imls/world/IND/IND_ADM0.shp

I have attached the global shape file (in Zip) here.

[KSSV]

Yes it should be very much possible.....attach your shape file here...

[Poulomi Ganguli]

I have attached it.

[KSSV]

There is no .dbfread file in the zip folder. And the link given downloads on .shp file. To read shape file it needs all supporting files.

[Poulomi Ganguli]

There is a .dbf file in the zip folder. For India shapefile with all supporting files are available in this link: https://otlet.berkeley.edu/imls/world/IND/

Kindly consider files associated with IND_ADM0.shp

For global shape file, please download from this link, if the zipped file does not load associated files: https://datacatalog.worldbank.org/dataset/world-maps-k%C3%B6ppen-geiger-climate-classification/resource/274ce998-385d-403b-8196

Please let me know if you have any further issues.

Answer 1

[KSSV]

shpFile1 = "c1976_2000.shp" ; 
shpFile2 = "IND_ADM1.SHP" ; 

S1 = shaperead(shpFile1) ; 
S2 = shaperead(shpFile2) ; 

P1 = [[S1(:).X]' [S1(:).Y]'] ; 
P2 = [[S2(:).X]' [S2(:).Y]'] ; 
N = length(S2) ; 
%% Get the points from gloabl file from given Indian file 
% Remove NaN's from P1
idx = isnan(P1(:,1)) ; 
P = cell(N,1) ;   % required matrix 
% Get the points lying inside shapefile2 from shape shapefile1
for i = 1:N
  i
  P2i = [[S2(i).X]' [S2(i).Y]'] ; 
  idx = isnan(P2i(:,1)) ; 
  P2i(idx,:) = [] ; 
  idx = inpolygon(P1(:,1), P1(:,2),P2i(:,1),P2i(:,2)) ; 
  P{i} = P1(idx,:) ; 
end
P = cell2mat(P) ; 
plot(P2(:,1),P2(:,2),'.b')
hold on
plot(P(:,1),P(:,2),'.r')

[Poulomi Ganguli]

This is good but can we crop the global file that lies within India map? It would be good if we can retain the polygon that lies within the India shapefile instead of points.

[KSSV]

With the points use boundary. Read about the function boundary.

[Poulomi Ganguli]

could you please show how it could be done here?