中国区域地图叠加绘制

已于 2024-02-21 14:29:55 修改
阅读量1.6k 收藏 11
点赞数 3
分类专栏： 地理信息数据处理 文章标签： matlab 开发语言
于 2023-05-05 17:50:29 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_38734327/article/details/130511032
版权

中国区域地图叠加绘制

0、声明

  未经声明的matlab版本是在2016版本上进行测试使用的，有小伙伴反映2022版本的matlab无法使用的情况，将内容更新在drawmap更新2中


1、数据准备

  地理信息地图在进行绘制的过程中，叠加上行政区域，更方便进行位置定位，更有利于针对问题进行具体分析。


  本文提供了两个Matlab函数，方便在画图时进行行政区域地图叠加操作。


  中国区域各种行政划分河流数据等shp文件，请自行下载，百度云链接链接：https://pan.baidu.com/s/1N6MIFYPHTz8Xpf3gCjvutw 提取码：082d


2、相关Matlab代码及使用方法

  叠加中国区域行政地图的函数如下:


  将下边函数复制，加载到matlab环境中，在需要进行叠加行政地图的区域使用drawmap(‘省’)即可完成行政区域的叠加操作，结果如下图所示。

function drawmap(Province, varargin)

% drawmap: Add administrative region information for a specified province on the canvas
% use whith default configuration    -->  drawmap('湖北省')
% use whithout global map            -->  drawmap('湖北省', 'global', 0)

% ------------------------------------------------------------ %
% Author:          Muyeqingfeng
% e-mail address:  muyeqingfeng@126.com
% Release date:    2023/05/05

ip = inputParser;
ip.addParamValue('global', 1);
parse(ip,varargin{:});  

% decoding the shp file wiht GBK
feature('DefaultCharacterSet', 'GBK');
path = 'shp_path';
China1=shaperead([path,'省.shp']);
China2=shaperead([path,'市.shp']);

idx1 = 0;
idx2 = [];
% Provincial  boundary infomation 
for i=1:length(China1)
    if strcmp(China1(i).x0xCA0xA1,Province), idx1=i;end
end
% Market boundary information
for j=1:length(China2)
    if strcmp(China2(j).x0xCA0xA1,Province), idx2=[idx2,j];end
end

if ip.Results.global == 1
    % the whole china info
    plot([China1(:).X],[China1(:).Y],'color',[0.2 0.2 0.2],'linewidth',1.0); hold on
end
% the special Province
plot([China1(idx1).X],[China1(idx1).Y],'color','r','linewidth',1.5); hold on
% all of the city of the special Province
plot([China2(idx2).X],[China2(idx2).Y],'color','k','linewidth',0.8); hold on
set(gca,'Layer','top');
end

函数调用结果展示

figure;drawmap('湖北省')

在这里插入图片描述

figure;drawmap('湖北省','global', 0)

在这里插入图片描述

3、图片固定位置信息标注

  以下提供了一个函数，在绘制图像时，在图像固定比例处进行信息标注，使得图像信息显得整齐一致

function text_info(varargin)
% text_info: Print information in a fixed position on the image
% use with the default configuration                               --->  text_info('text', 'the info need to print')
% change the position the display the info                         --->  text_into('xi', 0.1, 'yi', 0.2)
% change the color/fontsize of the display the info                --->  text_into('color', 'r', 'ftsiz', 20)

% ------------------------------------------------------------ %
% Author:          Muyeqingfeng
% e-mail address:  muyeqingfeng@126.com
% Release date:    2023/05/05

ip = inputParser;
ip.addParamValue('xi', 0.1);
ip.addParamValue('yi', 0.1);
ip.addParamValue('text', 'test');
ip.addParamValue('ftsiz', 30);
ip.addParamValue('color', 'k');
parse(ip,varargin{:});  

hh=get(gca);
X=hh.XLim;
Y=hh.YLim;
k3=[ip.Results.xi ip.Results.yi];
x_3=X(1)+k3(1)*(X(2)-X(1));
y_3=Y(1)+k3(2)*(Y(2)-Y(1));
text(double(x_3),double(y_3),ip.Results.text,'color',ip.Results.color,'fontweight','bold','fontsize',ip.Results.ftsiz-12);

end

函数调用结果展示

figure;drawmap('湖北省')
text_info('text', '湖北省')

在这里插入图片描述

figure;drawmap('湖北省', 'global', 0)
text_info('text', '湖北省', 'xi', 0.7,'yi', 0.9,'color', 'r')

在这里插入图片描述
在这里插入图片描述

drawmap更新

    drawmap更新，只画整个中国的图，不重点标出某一省份

function drawmap(varargin)

% drawmap: Add administrative region information for a specified province on the canvas
% use whith default configuration  -->  drawmap('湖北省')
% use whithout global map            -->  drawmap('湖北省', 'global', 0)

% ------------------------------------------------------------ %
% Author:             Muyeqingfeng
% e-mail address: muyeqingfeng@126.com
% Release date:    2023/05/05

ip = inputParser;
ip.addParamValue('global', 1)
ip.addParamValue('Province', 0);
parse(ip,varargin{:});  

feature('DefaultCharacterSet', 'GBK');
path = 'D:\My_Document\Code\Practice\data\全国shp\最新2021年全国行政区划\';
China1=shaperead([path,'省.shp']);
China2=shaperead([path,'市.shp']);

idx1 = 0;
idx2 = [];

if ip.Results.Province ~= 0
    % Provincial  boundary infomation 
    for i=1:length(China1)
        if strcmp(China1(i).x0xCA0xA1,ip.Results.Province), idx1=i;end
    end
    
    % Market boundary information
    for j=1:length(China2)
        if strcmp(China2(j).x0xCA0xA1,ip.Results.Province), idx2=[idx2,j];end
    end
end


if ip.Results.global == 1
    % the whole china info
    plot([China1(:).X],[China1(:).Y],'color',[0.2 0.2 0.2],'linewidth',1.0); hold on
end


if ip.Results.Province ~= 0
    % the special Province
    plot([China1(idx1).X],[China1(idx1).Y],'color','r','linewidth',1.5); hold on
    % all of the city of the special Province
    plot([China2(idx2).X],[China2(idx2).Y],'color','k','linewidth',0.8); hold on
end
set(gca,'Layer','top');
end

在这里插入图片描述

drawmap更新1

function drawmap(varargin)

% drawmap: Add administrative region information for a specified province on the canvas
% use whith default configuration  -->  drawmap('湖北省')
% use whithout global map            -->  drawmap('湖北省', 'global', 0)

% ------------------------------------------------------------ %
% Author:             Muyeqingfeng
% e-mail address: muyeqingfeng@126.com
% Release date:    2023/05/05

ip = inputParser;
ip.addParamValue('global', 1)
ip.addParamValue('Province', 0);
ip.addParamValue('City', 0);
ip.addParamValue('County', 0);
parse(ip,varargin{:});  

% feature('DefaultCharacterSet', 'UTF-8');
feature('DefaultCharacterSet', 'GBK');
path = './China/';
China1=shaperead([path,'省.shp']);
China2=shaperead([path,'市.shp']);
China3=shaperead([path,'县.shp']);

idx1 = 0;
idx2 = [];
idx3 = [];

if ip.Results.Province ~= 0
    % Provincial  boundary infomation 
    for i=1:length(China1)
        if strcmp(China1(i).x0xCA0xA1,ip.Results.Province), idx1=i;end
    end
    
    % Market boundary information
    for j=1:length(China2)
        if strcmp(China2(j).x0xCA0xA1,ip.Results.Province), idx2=[idx2,j];end
    end
    
    for k=1:length(China3)
        if strcmp(China3(k).x0xCA0xA1,ip.Results.Province), idx3=[idx3,k];end
    end
end


if ip.Results.global == 1
    % the whole china info
    plot([China1(:).X],[China1(:).Y],'color',[0.2 0.2 0.2],'linewidth',1.0); hold on
end


if  ip.Results.County ~= 0
% all of the city of the special Province
    plot([China3(idx3).X],[China3(idx3).Y],'color',[0.7,0.7,0.7],'linewidth',0.8); hold on
    set(gca,'Layer','top');
end

if  ip.Results.City ~= 0
% all of the city of the special Province
    plot([China2(idx2).X],[China2(idx2).Y],'color','k','linewidth',1); hold on
    set(gca,'Layer','top');
end

if ip.Results.Province ~= 0
    % the special Province
    plot([China1(idx1).X],[China1(idx1).Y],'color',[124, 0, 30]./255,'linewidth',1.5); hold on
end
  • 测试
set(gcf,'color','w','position',[0 0 800 1200])
axes('position', [0.2 0.2 0.6 0.4])
drawmap('Province','河南省','global',0,'city',1,'county',1)

在这里插入图片描述

drawmap更新2

  • Matlab version == 2022a
function drawmap(varargin)

% drawmap: Add administrative region information for a specified province on the canvas
% use whith default configuration  -->  drawmap('湖北省')
% use whithout global map            -->  drawmap('湖北省', 'global', 0)

% ------------------------------------------------------------ %
% Author:             Muyeqingfeng
% e-mail address: muyeqingfeng@126.com
% Release date:    2023/05/05

ip = inputParser;
ip.addParamValue('global', 1)
ip.addParamValue('Province', 0);
ip.addParamValue('City', 0);
parse(ip,varargin{:});  

% feature('DefaultCharacterSet', 'UTF-8');
% feature('DefaultCharacterSet', 'GBK');
path = '..\China\';
China1=shaperead([path,'省.shp']);
China2=shaperead([path,'市.shp']);

%{
matlab2016版本使用以下方式更改编码方式，就可以正常使用
feature('DefaultCharacterSet', 'UTF-8');
feature('DefaultCharacterSet', 'CBK');

但是2022版本好像不能解决上述问题，需要读取同路径下的dbf文件进行省份的匹配，然后作图
dbf也可以使用excel直接打开查看
%}
[num1,head1,text1] = xlsread("..\China\省.dbf");
[num2,head2,text2] = xlsread("..\China\市.dbf");

idx1 = 0;
idx2 = [];

if ip.Results.Province ~= 0
    % Provincial  boundary infomation 
    sheng = text1(2:end,1);
    shi = text2(2:end,1);
    for i=1:length(sheng)
        if strcmp(sheng(i),ip.Results.Province), idx1=i;end
    end
    
    % Market boundary information
    for j=1:length(shi)
        if strcmp(shi(j),ip.Results.Province), idx2=[idx2,j];end
    end
end


if ip.Results.global == 1
    % the whole china info
    plot([China1(:).X],[China1(:).Y],'color',[0.2 0.2 0.2],'linewidth',1.0); hold on
end


if ip.Results.Province ~= 0
    % the special Province
    plot([China1(idx1).X],[China1(idx1).Y],'color','r','linewidth',1.5); hold on
end

if  ip.Results.City ~= 0
% all of the city of the special Province
    plot([China2(idx2).X],[China2(idx2).Y],'color','k','linewidth',0.8); hold on
    set(gca,'Layer','top');
end
木叶清风666
关注
  • 3
    点赞
  • 踩
  • 11
    收藏
    觉得还不错? 一键收藏
  • 5
    评论
专栏目录
使用Matlab绘制山东省各个市地理分区示意图
2301_77925375的博客
07-13 963
我们可以在Matlab中使用自带的地图数据集，也可以从网上下载相应的地理数据集。这里，我们使用了`parula`函数生成一个包含5种颜色的颜色映射，然后遍历每个市，根据所属地理分区用不同颜色绘制市。这里，`cn_province.shp`是自带的中国行政区划数据集文件，`Shandong`是我们要获取的省份名称。这里，`FaceColor`参数设置地图颜色为灰色，`LineWidth`参数设置线宽为1。% 根据地理分区用不同颜色绘制市。% 根据地理分区用不同颜色绘制市。% 获取市的名称和所属地理分区。
【GIS】基于MATLAB的超详细中国地图绘制（附代码）.zip
03-03
【GIS】基于MATLAB的超详细中国地图绘制（附代码）.zip
调用百度地图API：绘制多个地区的轮廓图
最新发布
knighthood2001
05-22 830
以下代码中，通过更改其中的地名，打开就会高亮具体的地区，并框选轮廓。// 创建多边形覆盖物 fillOpacity设置填充物的不透明度。如果你需要多个地区，你也可以改源代码。'未能获取到该地区的边界数据''未能获取到该地区的边界数据'// 要绘制的地区及对应颜色。// 绘制某个地区的轮廓图。// 绘制多个地区的轮廓图。// 获取行政区域边界。// 创建多边形覆盖物。// 添加覆盖物到地图。// 获取行政区域边界。// 添加覆盖物到地图。
省级分布图matlab实现（有九段线）
m0_50309522的博客
07-20 720
clc; clear; close all; Chprovince_0 = shaperead('省界111.shp', 'UseGeoCoords', true);%读入31省界数据，建立结构体 Chprovince_1 = shaperead('台湾和香港.shp', 'UseGeoCoords', true); Chprovince_2 = shaperead('九段线.shp', 'UseGeoCoords', true) Chprovince_3 = shaperead('南海诸岛及其...
MATLAB绘制中国地图超全教程
slandarer的博客
02-25 9995
依旧需要用到Mapping Toolbox不会安装的可以看我上一篇 虽然我们只读取shp文件，但需要保证文件夹里还有shx文件及dbf文件 各省边界线绘图 provinces=shaperead('bou2_4l.shp','UseGeoCoords',true); % 绘图 worldmap('China'); geoshow(provinces) 省份填色图 provinces=shaperead('bou2_4p.shp','UseGeoCoords',true); % 绘图 worldmap
MATLAB—地图
hellogril前端库
11-09 3108
如果数据已经被投影了就选择false，如果还未被投影，可以选择true以生成投影。3） geoshow(_,Name,Value,)输入相关属性并设定参数，生成所需要的地图。1） geoshow(lat,lon)直接输入经纬度，将以该经纬度为中心，生成一小片地图。1） facecolor：定义了地图表面的颜色，需要输入三个参数且均在0至1之间。geoshow是用来显示地图数据的函数，非常重要，大部分的地图都使用该函数进行显示。
map1_绘制中国地图_
09-29
本知识点将深入探讨如何使用编程语言，如MATLAB，来绘制中国地图，特别是针对"map1_绘制中国地图_"这个主题。在这个过程中，我们将涉及到`map1.m`这个脚本文件，它很可能是实现地图绘制功能的核心代码。 首先，我们...
中国分省行政区矢量图（shp）
03-23
【中国分省行政区矢量图（shp）】是一个专门针对中国行政区域划分的数据集，主要服务于地理信息系统（GIS）的应用。这个数据集包含了中国所有省份、直辖市、自治区以及特别行政区的边界信息，以矢量格式存储，可以...
中国行政区矢量图.rar
08-20
综合以上信息，这个压缩文件提供了一套完整的中国省级行政区域的地理矢量数据，可以用于各种用途，比如绘制高精度的地图、进行人口分布分析、规划交通网络、评估自然灾害的影响范围等。在GIS软件中，如ArcGIS或QGIS...
中国城市地名及区域划分数据集.zip
10-09
《中国城市地名及区域划分数据集》是一个包含中国各个城市及其行政区域划分的重要资源，对于研究中国地理、社会经济、人口分布以及进行数据分析的人员来说具有极高的价值。这个数据集以压缩文件形式提供，其中的核心...
matlab奇技淫巧——绘制三维地图
qq_43557686的博客
10-18 3200
matlab绘制三维地图（地球仪）
MATLAB轻松绘制地图路线——已知及未知坐标下的处理方法(2)
Wayne的CSDN博客
05-02 4565
11B交通巡警平台地图绘制_简化版
MATLAB绘制地图超详细教程
slandarer的博客
02-23 2万+
worldmap和usamap是axesm的子类，worldmap是用于生成世界地图坐标区域，usamap用于生成美国地图坐标区域，本文先在worldmap函数基础上讲解如何导入各种数据绘制各种类型图片，略提一下如何使用usamap，再讲解axesm的各种属性以满足更多元化的地图绘制需求，并在最后总结地图绘制其他可用函数。 Mapping Toolbox工具箱安装 请通过如下步骤安装工具箱： 点击附加功能： 搜索Mapping Toolbox，并点击进入介绍界面： 点安装一路点确定： wor
matlab怎么合并两个shp地图_没想到使用R画地图数据这么简单
weixin_42512474的博客
12-18 703
用R画地图数据首先，先下载中国地图的GIS数据，这是一个压缩包，完全解压后包含三个文件(bou2_4p.dbf、bou2_4p.shp和bou2_4p.shx)，将这三个文件解压到同一个目录下。用R绘制地图比较简单。比如画一下全国范围的区域，可以用如下代码：library(maptools)mydat = readShapePoly("china-province-border-data.tar/...
在matlab中对中国地图中的不同省份按照高度进行渲染
weixin_34344677的博客
01-04 2302
直接上优化后的代码和效果图 colour_totall=128; % 载入地图数据 --各省的多边形数据 shp_data=shaperead('maps/bou2_4p.shp', 'UseGeoCoords', true);% 省 %{ // matlab 2016b已经解决这个问题了 % 检查数据发现中午名称为乱码，可能是由于字符编码引起的，需要纠正文件导入时出现的...
MATLAB绘制一幅中国地图
say_book的博客
07-30 9949
     今天博主跟大家讲一下如何用MATLAB制作一幅中国地图，那废话不多说，我们先看一下最终效果吧。 mercator墨卡托圆柱投影地图    lambert兰伯特圆锥投影地图  一张中国地图大概包括以下要素： 中国矢量数据（面、省界线等shp数据）、地图四要素（图名或标题、指北针、比例尺、图例），最好加一个海南岛及南海诸岛的鹰眼图。    大概思路： 清理内存 地图投影 加载面文件 加载省界线文件（包含国界（含九段线）
Matlab下地形图绘图包m_map安装与使用
热门推荐
yangshun_cug
01-27 4万+
m_map是Matlab下用于绘制地形图的工具箱，和GMT有些相似，被广泛应用于地学研究成图。那么，应该如何下载和安装m_map呢？
利用MATLAB创建栅格地图(代码可复制)
algorithmzzy的博客
10-25 6996
分享两段MATLAB创建栅格地图的代码~
Matlab：使用纬度和经度数据创建地图
希望我的博客，能帮上你解决学习中工作中所遇到的问题
02-27 1430
Matlab：使用纬度和经度数据创建地图
splitContainer1.Panel2叠加绘制100个图像
05-27
要在 `splitContainer1.Panel2` 控件中叠加绘制100个图像，您可以在 `splitContainer1_Panel2_Paint` 函数中使用 `Graphics` 对象的 `DrawImage` 方法，并将 `ImageAttributes` 对象的 `ColorMatrix` 属性设置为透明度小于 1 的值，从而实现图像叠加效果。 以下是示例代码： ```csharp private void splitContainer1_Panel2_Paint(object sender, PaintEventArgs e) { // 加载图像文件 Image image = Image.FromFile("image.jpg"); // 创建 ImageAttributes 对象 ImageAttributes imageAttributes = new ImageAttributes(); // 创建 ColorMatrix 对象，将透明度设置为0.5 ColorMatrix colorMatrix = new ColorMatrix(); colorMatrix.Matrix33 = 0.5f; imageAttributes.SetColorMatrix(colorMatrix); // 计算每个图像在控件中的位置和大小 int imageWidth = splitContainer1.Panel2.Width / 10; int imageHeight = splitContainer1.Panel2.Height / 10; Rectangle imageRect = new Rectangle(0, 0, imageWidth, imageHeight); // 在控件中循环绘制每个图像 for (int i = 0; i < 100; i++) { // 计算当前图像在控件中的位置 int x = (i % 10) * imageWidth; int y = (i / 10) * imageHeight; imageRect.Location = new Point(x, y); // 绘制当前图像 e.Graphics.DrawImage(image, imageRect, 0, 0, image.Width, image.Height, GraphicsUnit.Pixel, imageAttributes); } } ``` 在上面的代码中，我们首先加载了一个图像文件，并创建了一个 `ImageAttributes` 对象。然后，我们创建了一个 `ColorMatrix` 对象，并将其中的 `Matrix33` 属性设置为 0.5，即将透明度设置为 0.5。接下来，我们计算了每个图像在控件中的位置和大小，并使用一个 for 循环在控件中叠加绘制每个图像。在绘制每个图像时，我们使用 `e.Graphics.DrawImage` 方法，并将 `ImageAttributes` 对象作为最后一个参数传递，从而实现图像叠加效果。 请注意，叠加绘制多个图像可能会对性能产生影响，具体取决于图像的大小和计算机的性能。如果您需要叠加绘制更多的图像，可以考虑使用双缓冲技术来提高性能。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 5
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值