attach a field to a shapefile and plot it

15 views (last 30 days)
I have a shapefile of India with state boundaries (34 polygons within the India shapefile).See the attached screenshot and the shapefile as .zip. I want to attach the mean rainfall data in a separate feild to each state in the shapefile and plot it in matlab, I am familiar with dealing with shapefiles in arcgis but would like to perform it in matlab as well.

Accepted Answer

Kelly Kearney
Kelly Kearney on 26 Sep 2018
Do you need to actually add a new field to the shapefiles themselves? Or just add it to the read-in structure in Matlab?
If the latter, you can easily assign new fields and values to your structure
India = shaperead('~/Downloads/shp/naturalearth_india_adm1.shp', 'usegeocoords', true);
for ii = 1:numel(India)
India(ii).rain = rand(1); % add your data here
end
You can plot shapefiles colored by a property by using the (extremely clunky, IMO) symbol spec parameter to geoshow:
spec = makesymbolspec('Polygon', ...
{'Shape_Area', [-11 31], 'Facecolor', parula})
geoshow(India, 'SymbolSpec', spec);
In this example, I've used the area of each polygon as the color-mapping data, but you can adjust it as needed. The down side to this plotting method is that the colors don't correspond to the colormap at all, so you need to manually sync your color limits to match the spec if you plan to use a colorbar or plot other data on the same color scale.
If you need more control over the output, you can instead plot without the symbol spec and set the CData property of the resulting patches:
h = geoshow(India);
set(h.Children, {'CData'}, {India.Shape_Area}');
set(h.Children, 'Facecolor', 'flat');
Using this method will link the patch color to the axis color limits and colorbar. Which behavior you prefer may depend on your application.
  2 Comments
Sourangsu Chowdhury
Sourangsu Chowdhury on 27 Sep 2018
Edited: Sourangsu Chowdhury on 27 Sep 2018
@Kelly Kearney Thank you for the help. I would ideally need to display the output as binned colors like for india.shp_area, 0-10 should have one uniform color, 11-20 should have a different color and so on.
Kelly Kearney
Kelly Kearney on 27 Sep 2018
Take a look at Example 3 (Using a Range of Attribute Values) in the makesymbolspec documenatation.

Sign in to comment.

More Answers (0)

Categories

Find more on Colormaps in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!