無需底圖matlab實(shí)現(xiàn)畫任何矢量地圖

1、無需底圖matlab實(shí)現(xiàn)畫任何矢量地圖主要借助arcgis中的矢量shp線面和envi中的掩膜功能?？傮w思路是逐個要素的畫，讓某個圖層要素（比如海岸線）的所在矩陣位置的值先變成1其它變?yōu)?（envi中掩膜實(shí)現(xiàn)），然后1再去乘以不同的值（matlab中實(shí)現(xiàn)），以區(qū)分不同的圖層要素。在envi中將矢量轉(zhuǎn)成同一個圖層（圖一個波段的遙感影像圖）的roi就能保證每個要素所在的矩陣行列都一樣了。步驟：1、在arcgis中準(zhǔn)備好各種shp矢量面線，一般線要做個buffer轉(zhuǎn)化成面。analysis tools>>proximity>>buffer2、envi下打開矢量功能，file

2、>> open vector,選擇在memory 中間紅色圈是坐標(biāo)系統(tǒng) 默認(rèn)為arcgis矢量文件的坐標(biāo)，可以不用改，直接ok3、load矢量，然后選擇display，然后ok，到vector parameter界面4、在vector parameter界面，將矢量轉(zhuǎn)化成roi，file>>export layer to rois,然后選擇建立的roi給哪個圖層，這一點(diǎn)很重要，這一點(diǎn)也是保證你轉(zhuǎn)化成的掩膜文件，各種圖的大小都一樣，也就是在matlab中打開的矩陣行列都是一樣的。然后選擇第一個。5、將roi的主人（剛才給roi的那個圖層）打開，overly中打開roi，可

3、以看到剛才的矢量shp已經(jīng)變成roi了。6、接著用這個roi去建立掩膜。envi class中，basic tools>>masking>>build mask,選擇display，到達(dá)下面這個界面7、在masking definition這個界面中,用roi去建立掩膜，具體操作是，options>>import roi，選擇剛才的roi，然后可以momory，然后掩膜就建立好了。掩膜嘛 就是將值變成0和1，有要素的變成1 沒有的變成0嘛9、然后將這個掩膜保存，注意envi中的保存的區(qū)別，用envi class那個條的file菜單下保存功能，保存成tif嘛，

4、可以直接保存成灰度圖，在matlab中打開，矩陣的值就是在envi中波段的dn值。如果用顯示圖層的那個file菜單下的保存功能，那么保存出來的就是rgb圖，so，你懂得。10、接著在matlab中有要素的1去乘以不同的值，加以區(qū)別要素。例子：clc;clear;close all;% I=imread('H:modis2修改modis第二波段MOD09Q1.A2012241.sur_refl_b02GG.tif','tif'); %讀取圖像文件I=imread('I:混合像元分解后2012005.tif','tif'); %讀取圖

5、像文件% 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(I2>0 & 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; %注意是加起來哦