在matlab中实现遥感影像和shp文件的结合显示

clc;close all;clear;
road=shaperead(‘boston_roads.shp‘);  %读取shape文件
figure,
mapshow(‘boston.tif‘);  %读取tif遥感影像并显示
axis image manual off; %关闭本地坐标系统
surveyFeetPerMeter = unitsratio(‘sf‘,‘meter‘); %将米级单位转换成surveyFeet单位
colorTab=‘rgbcymkw‘; %颜色代码表，不同等级的公路显示不同颜色
for i=1:numel(road)
    x = surveyFeetPerMeter * [road(i).X]; %将shp文件坐转换为tif文件坐标
    y = surveyFeetPerMeter * [road(i).Y];
    mapshow(x,y,‘Color‘,colorTab(road(i).CLASS)); %根据公路类别来进行不同颜色显示
end
print(‘-dpng‘,‘-r600‘,‘d:/boston.png‘);  %以600的分辨率保存结果图片

时间： 2024-10-05 05:02:06

在matlab中实现遥感影像和shp文件的结合显示的相关文章

在matlab中进行遥感影像地理坐标的相互转换

在matlab中进行图像处理,一般使用的都是图像本地坐标,以左上角(1,1)开始.处理完成后,如果要将结果在带地理坐标的遥感影像中显示,或者需要输出成shp文件,就需要涉及到本地坐标和地理坐标的转换,我们可以使用pix2map函数来实现这个功能. clc;clear;close all; [pic, R] = geotiffread('boston.tif'); %读取tif图像 figure(1),imshow(pic);% pic即为本地图像,下面对本地图像进行操作 gray=rgb2gra

在matlab中将处理结果输出为shp文件

在matlab中读入shp文件很简单,一个函数shaperead就可以了,但输出为shp文件就稍微麻烦一些了.shp文件实际上就是一个struct,因此得到处理结果后,要先将数据变成struct结构,然后调用shapewrite函数就可以了. clc;clear;close all; [pic, R] = geotiffread('boston.tif'); %读取tif图像 gray=rgb2gray(pic); bw=im2bw(gray,0.1); %阈值分割 bw=1-bw; bw1 =

Matlab中怎样将.dat格式的文件进行加载转化为一个矩阵??

在网上下载了一个数据集rating.dat,格式为.dat的,要在matlab中使用,必须将其转化为一个.mat格式的矩阵,下面写一下具体步骤: 第一步:打开matlab,界面如图: 然后点击“打开”选项:界面如下 例如选择rating.dat 文件 点击上图右下角“打开”,得到如下界面: 然后在“列分隔符”里进行选择,由于原文本文件格式如下图: 所以选择“列分隔符”中的::,然后选中所有数据,再点击"导入",选择其中的“导入数据”,则得到结果为m*4的矩阵:

Matlab中的数组保存为csv文件，在VC里面使用

由于图像处理用到Matlab的一些算法,得到一些数组数据(三维重构的深度数据). 在VC中使用OpenGL画图,需要使用这些数据(深度数据)做显示用. 在Matlab下使用语句:csvwrite('depthMap',depthMap); //后面的depthMap是数组变量,前面的depthMap是csv文件名. 函数ReadCsv是读取csv文件(二维数组)到指定变量depthMap. 变量的定义: double depthMap[240][320] = {0}; //深度图 调用函数: /

Matlab中legend位置

转自:http://blog.sina.com.cn/s/blog_7db803c10102weyk.html Matlab中legend默认的位置在NorthEast,如图所示: %Matlab中legend的位置设置 clc clear close all Npoint = 100; x = linspace(0,4*pi,Npoint); y1 = sin(x); y2 = cos(x); H = plot(x,y1,x,y2); legend('sin(x)','cos(x)'); 然而

Matlab计算多幅遥感影像交集面积占比

一个文件夹中有同一时间不同的遥感影像,有LST NDVI.反射率等,而想要知道这些像元数量相同,填充值和背景值均为nan的影像的交集部分占该影像研究区(已经经过研究区矢量裁剪)的比例,以及不同类别遥感影像的覆盖面积比例,主要目的是想以多种数据交集覆盖面积比例来判断改天是否适合数据反演 1 函数 函数基本结构 function[输出形参表: output1, ...,outptn] = 函数名(输入形参表: input1, ... , inputn) 注释说明部分 函数体代码部分 end 函数名:

基于单分类器的高分辨率遥感影像道路提取

本人硕士阶段做了很久的高分辨率遥感影像道路提取,颇有心得,在此,本人将最新的研究成果进行开源... 大家都知道,传统的基于机器学习的分类方法通常需要正负样本的同时参与,才能得到目标类,但是负样本的勾选,通常很困难,也非常难获得,根据文献- <Elkan, Charles, and Keith Noto. "Learning classifiers from only positive and unlabeled data." Proceedings of the 14th ACM

遥感影像滤波处理软件 — timesat3.2

最近因为要做遥感影像的滤波处理,经过女神推荐,决定用Timesat,可是该软件3.1版本只适合xp系统以及2011的matlab,后来在官网上找到了最新的3.2版本.支持64位操作系统以及2014的matlab.大家可以直接上官网(http://www.nateko.lu.se/TIMESAT/timesat.asp )注册下载,但是有一个问题是,官网注册的时候,会用到Google公司的rechaptcha验证码系统,因为大天朝的墙太高了,所以是刷不出来的,可以选择挂vpnFQ下载.如果不想FQ

提取遥感影像有效区域

利用某些算法处理遥感影像时,需要排除无效区域,因此找出有效区域范围变的至关重要,下面根据我的经验谈谈有效区域提取的具体方法及实现. 创建金字塔图层 由于遥感影像数据量比较大,因此为了提高处理效率,提取有效区域需要在金字塔图层进行.金字塔的创建可以利用GDAL进行,也可以自己进行缩放. 提取有效区域shp文件 为了提取有效区域,我研究了两种方法可以实现,分别是栅格矢量化和轮廓跟踪,下面具体解释两种方法的区别: 栅格矢量化 该方法原理比较简单,首先按照指定的无效值,将图像进行二值化,得到掩模图,然后