MATLAB Answers

How to Crop India shape file from Global shapefile

7 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

Show 1 older comment
KSSV
KSSV on 31 Aug 2020
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
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
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
KSSV on 31 Aug 2020
With the points use boundary. Read about the function boundary.
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!