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

最新推荐文章于 2024-03-31 17:40:32 发布

weixin_42408281

最新推荐文章于 2024-03-31 17:40:32 发布

阅读量1.4k

点赞数

本文链接：https://blog.csdn.net/weixin_42408281/article/details/125359630

版权

图像处理二值化连通域去除 GeoTIFF 地理信息系统

关键词由CSDN通过智能技术生成

close all;
clear all;
clc;

F=imread('C:\Users\My\Desktop\NDSI.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:\Users\My\Desktop\NDSI.tif');
%Info=imfinfo(filepath):获取图片信息并判断是否为tif
info=geotiffinfo('C:\Users\My\Desktop\NDSI.tif');
%输出 第一个参数是输出路径 第二个参数是图像 第三个参数是图像的地理信息，第四、五个参数是投影信息
%geotiffwrite(outFile, treeTop, geo,'CoordRefSysCode',CoordRefSysCode)
geotiffwrite('C:\Users\My\Desktop\pig',BW2,R,'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);