MATLAB Answers

How to Crop India shape file from Global shapefile

15 views (last 30 days)
Poulomi Ganguli
Poulomi Ganguli on 30 Aug 2020
Commented: Poulomi Ganguli on 31 Aug 2020
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.
  4 Comments
Poulomi Ganguli
Poulomi Ganguli on 31 Aug 2020
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.

Sign in to comment.

Accepted Answer

KSSV
KSSV on 31 Aug 2020
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')
  3 Comments
Poulomi Ganguli
Poulomi Ganguli on 31 Aug 2020
could you please show how it could be done here?

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!