clear; close; clc

latlim = [25 50]; lonlim = [-105 -67];  % Breiten- und Laengengrad-Bereich definieren
                                        % Teilbereich klappt

% latlim = [-90 90]; lonlim = [-180 180];  gesamte Darstellung haut nicht
% hin, 

axe = usamap(latlim, lonlim);

landareas = shaperead('landareas.shp','UseGeoCoords', true, 'BoundingBox', [lonlim', latlim']);
geoshow(landareas);