Matlab大津法提取水体，并输出带有参考坐标系的tif影像

368 阅读 0 评论 243 点赞

我是靠谱客的博主无心羊，这篇文章主要介绍Matlab大津法提取水体，并输出带有参考坐标系的tif影像，现在分享给大家，希望可以做个参考。

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
close all;
clear all;
clc;

F=imread('C:UsersMyDesktopNDSI.tif');  
f=F(:,:,1);
histogram(f);
%thresh = multithresh(f,1);
thresh = graythresh(f);

g1=f;
g1(find(f>thresh))=1;
g1(find(f<=thresh))=0
imshow(g1)

%面积计算
S=numel(g1);%像素点总数
s=sum(sum(g1));%白色点总数
ratio=(S-s)/S;%白色面积比

%小连通域去除
BW2 = bwareaopen(g1,10,8); %P为自定义面积，删除小于该面积的连通区域
imshow(BW2);

%[A,R] = geotiffread(filename),从filename指定的GeoTIFF文件中读取一个地理参考的灰度、RGB或多光谱图像或数据网格到A中，并创建一个空间参考对象，R。
[image,R]=geotiffread('C:UsersMyDesktopNDSI.tif');
%Info=imfinfo(filepath):获取图片信息并判断是否为tif
info=geotiffinfo('C:UsersMyDesktopNDSI.tif');
%输出 第一个参数是输出路径 第二个参数是图像 第三个参数是图像的地理信息，第四、五个参数是投影信息
%geotiffwrite(outFile, treeTop, geo,'CoordRefSysCode',CoordRefSysCode)
geotiffwrite('C:UsersMyDesktoppig',BW2,R,'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);