无需底图matlab实现画任何矢量地图.doc

无需底图matlab实现画任何矢量地图主要借助arcgis中的矢量shp线面和envi中的掩膜功能。总体思路是逐个要素的画，让某个图层要素（比如海岸线）的所在矩阵位置的值先变成1其它变为0（envi中掩膜实现），然后1再去乘以不同的值（matlab中实现），以区分不同的图层要素。在envi中将矢量转成同一个图层（图一个波段的遥感影像图）的roi就能保证每个要素所在的矩阵行列都一样了。步骤：1、在arcgis中准备好各种shp矢量面线，一般线要做个buffer转化成面。analysis toolsproximitybuffer2、envi下打开矢量功能，file open vector,选择在memory 中间红色圈是坐标系统 默认为arcgis矢量文件的坐标，可以不用改，直接ok3、load矢量，然后选择display，然后ok，到vector parameter界面4、在vector parameter界面，将矢量转化成roi，fileexport layer to rois,然后选择建立的roi给哪个图层，这一点很重要，这一点也是保证你转化成的掩膜文件，各种图的大小都一样，也就是在matlab中打开的矩阵行列都是一样的。然后选择第一个。5、将roi的主人（刚才给roi的那个图层）打开，overly中打开roi，可以看到刚才的矢量shp已经变成roi了。6、接着用这个roi去建立掩膜。envi class中，basic toolsmaskingbuild mask,选择display，到达下面这个界面7、在masking definition这个界面中,用roi去建立掩膜，具体操作是，optionsimport roi，选择刚才的roi，然后可以momory，然后掩膜就建立好了。掩膜嘛 就是将值变成0和1，有要素的变成1 没有的变成0嘛9、然后将这个掩膜保存，注意envi中的保存的区别，用envi class那个条的file菜单下保存功能，保存成tif嘛，可以直接保存成灰度图，在matlab中打开，矩阵的值就是在envi中波段的dn值。如果用显示图层的那个file菜单下的保存功能，那么保存出来的就是rgb图，so，你懂得。10、接着在matlab中有要素的1去乘以不同的值，加以区别要素。例子：clc;clear;close all;% I=imread(H:modis2修改modis第二波段MOD09Q1.A2012241.sur_refl_b02GG.tif,tif); %读取图像文件I=imread(I:混合像元分解后2012005.tif,tif); %读取图像文件% I22=rgb2gray(I);I2=double(I);figure(1);imshow(I2,);sizeI2= size(I2);a1=zeros(sizeI2);a2=a1;a3=a1;a1(I2=0)=2;a2(I2=0.83)=1;a3(I20 & I2=0.83)=0;hyym=imread(I:混合像元分解后海洋掩膜.tif,tif); %读取图像文件hyymd=double(hyym);hyym2=hyymd*(-1);bjym=imread(I:混合像元分解后行政边界研磨.tif,tif); %读取图像文件bjymd=double(bjym);bjym2=bjymd*(4);axym=imread(I:混合像元分解后海岸线掩膜.tif,tif); %读取图像文件axymd=double(axym);axym2=axymd*(16);B=a1+a2+a3+hyym2+bjym2+axym2; %注意是加起来哦 下面的怎么显示