How to create a chloropleth map using the mean of the shapefiles sub-fields from a raster file?

1 次查看（过去 30 天）
Abhishek Chakraborty
评论： Abhishek Chakraborty 2022-11-21
I am trying to create a chloropleth map which would show the means over the sub-regions of the shapefile using a raster file (the raster file is attached here and named as "Raster.mat". The raster file's latitude and longitude are attached here as "RasterLat.mat" and "RasterLon.mat", respectively. The shapefile is attached as "District_Boundaries.zip" with this question. Kindly help me with the same.
Thanks.

请先登录，再回答此问题。

回答（1 个）

KSSV
KSSV 2022-11-21
编辑：KSSV 2022-11-21
  1 个评论
Abhishek Chakraborty
Abhishek Chakraborty 2022-11-21
Hi. I tried running the code as per the suggestions but I am unable to get any results. Here is my code:
shapefile = 'C:\Users\District_Boundaries.shp' ;
S = shaperead(shapefile) ;
N = length(S) ;
for i = 1:N
plot(S(i).X,S(i).Y)
hold on
end
info=geotiffinfo("F:\India\Result.tif");
[x,y]=pixcenters(info); % x is Long and y is Lat
[X,Y]=meshgrid(x,y);
data = load("C:\Users\Abhishek\Downloads\Raster.mat") ;
[nx,ny,d] = size(data) ;
%%Extract data
iwant = cell(d,N) ;
for i =1:d
A = data(:,:,i) ;
for j = 1:N
idx = inpolygon(X(:),Y(:),S(i).X,S(i).Y) ;
iwant{i,j} = A(idx) ;
end
end

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Raster Data 的更多信息

标签

产品


版本

R2020b

Community Treasure Hunt

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

Start Hunting!

Translated by