How to Crop India shape file from Global shapefile

Question

Poulomi Ganguli 2020-8-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/586508-how-to-crop-india-shape-file-from-global-shapefile

评论： Poulomi Ganguli 2020-8-31

1976-2000_Global.zip

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 个评论
显示 2更早的评论隐藏 2更早的评论

KSSV 2020-8-31

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 2020-8-31

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 2020-8-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/586508-how-to-crop-india-shape-file-from-global-shapefile#answer_487256

在 MATLAB Online 中打开

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')