MATLAB: Problem in reading .shp file

mappingMapping Toolboxshaperead

Hi,

I have a problem in reading .shp file.

My .shp file contains all the major cities in Australia.

I want to read all the major cities located in the state of Queensland only.

I've tried to write code

shapedir = fullfile (pwd, 'Shapefile');
main_cities = shaperead (fullfile (shapedir, 'MainCities.shp'), 'usegeocoords', true);

But after this, I do not know how to do it.

When I opened the variable in the variable editor, I find a field named 'state', containing the name of different states in each struct.

Can someone help me? thank you

Best Answer

Hi, Raymond

You can perform that by adjusting the selector parameter.

I'll give you sample code :

shapedir = fullfile(pwd,'Shapefile');
main_cities = shaperead(fullfile(shapedir,'MainCities.shp'), ...
                       'usegeocoords',true,'selector', ...
                       {@(name) strcmp(name,'Queensland'),'state'});

Related Solutions

MATLAB: Issue creating a map of the USA

You can see an example here.

figure
ax = usamap('all');
set(ax, 'Visible', 'off')
states = shaperead('usastatelo', 'UseGeoCoords', true);
names = {states.Name};
indexHawaii = strcmp('Hawaii',names);
indexAlaska = strcmp('Alaska',names);
indexConus = 1:numel(states);
indexConus(indexHawaii|indexAlaska) = []; 
stateColor = [0.5 1 0.5];

Display the three regions.

geoshow(ax(1), states(indexConus),  'FaceColor', stateColor)
geoshow(ax(2), states(indexAlaska), 'FaceColor', stateColor)
geoshow(ax(3), states(indexHawaii), 'FaceColor', stateColor)

Hide the frame.

for k = 1:3
    setm(ax(k), 'Frame', 'off', 'Grid', 'off',...
      'ParallelLabel', 'off', 'MeridianLabel', 'off')
end

Add labels to each state.

lat = [states(indexConus).LabelLat];
lon = [states(indexConus).LabelLon];
textm(lat, lon, {states(indexConus).Name},'HorizontalAlignment', 'center','Parent',ax(1));
latA = [states(indexAlaska).LabelLat];
lonA = [states(indexAlaska).LabelLon];
textm(latA, lonA, {states(indexAlaska).Name},'HorizontalAlignment', 'center','Parent',ax(2));
latH = [states(indexHawaii).LabelLat];
lonH = [states(indexHawaii).LabelLon];
textm(latH, lonH, {states(indexHawaii).Name},'HorizontalAlignment', 'center','Parent',ax(3));

MATLAB: Read a shapefile with name that was provided from sprintf

    getshp = shaperead(shp2open,'UseGeoCoords',true)

Best Answer

Related Solutions

MATLAB: Issue creating a map of the USA

MATLAB: Read a shapefile with name that was provided from sprintf

Related Question