[转]利用matlab绘制地图上的点线面
Posted 司徒鲜生
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[转]利用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
以上是关于[转]利用matlab绘制地图上的点线面的主要内容,如果未能解决你的问题,请参考以下文章
用C#代码写,ArcEngine开发,在加载的地图上绘制多边形后,导出绘制区域内的所有矢量要素(点线面)