[转][修]利用matlab绘制地图上的点、线、面

一.绘制点

%生成背景地图地图  

h = worldmap('France');

%读取和显示大陆架  

landareas = shaperead('landareas.shp','UseGeoCoords', true);  

geoshow (landareas, 'FaceColor', [1 1 .5]);  

%埃菲尔铁塔的坐标  

TowerLon = 2.28;%经度坐标  

TowerLat = 48.85;%纬度坐标  

%用红色标记绘制%埃菲尔铁塔  

geoshow(TowerLat, TowerLon, 'Marker','.','MarkerEdgeColor','red')  

                                  %点显示             红色点

textm(TowerLat,TowerLon + 0.5, 'Eiffel Tower');%放置文本  

绘制的图像如下:

二.绘制线段、河流

%生成背景地图地图  

worldmap world;  

framem on;%set the map axes property Frame to 'on' or 'off'. The default setting for map axes is 'off'.

%绘制线条  

linem([15; 0; -45; -25],[-100; 0; 70; 110],'r-') 

linem([15; -30; -60; -65],[-100; -20; 100; 150],'b-')  

linem([15; 20; 40; 20],[-100; -20; 40; 50], 'g-')  

 

绘制的线条结果如下:

绘制线条在地图中是很常用的,geoshow的内部就是调用linem来绘制河流的:

clc;clear;

clf;%Clear current figure window

%生成背景地图地图  

h = worldmap('France');

landareas = shaperead('landareas.shp','UseGeoCoords', true);  

geoshow (landareas, 'FaceColor', [1 1 .5]);  

%河流数据

rivers = shaperead('worldrivers', 'UseGeoCoords', true);  

%rivers =   128x1 struct array(128条河流) with fields:  

     %   Geometry  几何结构：线型

     %   BoundingBox  河流端点的经纬度

     %   Lon  纬度（点集）

     %   Lat  经度（与经度对应）

     %   Name  名字

%显示河流  

geoshow(rivers, 'Color', 'blue') ·

rivers(8)  %查看rivers的第八个元素的信息 

绘制的图形如下:

三.绘制多边形、州

%多边形的顶点坐标  

%注意:数据中多边形的初始顶点和结束顶点必须重合 

x = [40 55 33 10 0  5  10 40 NaN 10 25 30 25 10 10 NaN 90 80 65 80 90 NaN];  

y = [50 20 0  0  15 25 55 50 NaN 20 10 10 20 30 20 NaN 10 0  20 25 10 NaN];  

mapshow(x,y,'DisplayType','polygon') %可以更改'DisplayType'的参数，具体查看帮助

%mapshow绘制'polygon'，注意多边形顶点的环绕方向：

%1/顺时针为正,里面被填充->"面"  

%2/逆时针为负,里面不被填充->"洞"  

 

使用geoshow来绘制美国的地图,每个州都是一个多边形结构:view plain

 

figure; ax = usamap('conus');  

axis off                            %不显示坐标轴

states = shaperead('usastatelo', 'UseGeoCoords', true);  

    %states =   51x1 struct array（51个州） with fields:

    %每个州都有如下属性，例如取出states(1)

    %        Geometry: 'Polygon'      几何属性是：多边形

    %     BoundingBox: [2x2 double]   最小、最大边界点的坐标

    %             Lon: [1x105 double] 点的纬度

    %             Lat: [1x105 double] 点的经度、这些离散点闭合组成一个州

    %            Name: 'Alabama'      州的名称

    %        LabelLat: 32.282672      

    %        LabelLon: -86.9206  

    %     PopDens2000: 87.6  

indexConus = 1:numel(states);%州的数目  

stateColor = [0.5 1 0.5];    %设置颜色  

%画出所有的州,方法一：

geoshow(ax, states(indexConus), 'FaceColor', stateColor) %更精确

%画出所有的州,方法二：

%geoshow(ax, states, 'FaceColor', stateColor)  

%只画出一个州

%geoshow(ax, states(1), 'FaceColor', stateColor)   

setm(ax, 'Frame', 'off', 'Grid', 'off',...    %Set properties of map axes and graphics objects 

     'ParallelLabel', 'off', 'MeridianLabel', 'off')  

         %平行圈                   %子午圈

绘制的图形如下:

原文链接：https://blog.csdn.net/rumswell/article/details/16927565