Google Earth Engine(GEE)实例代码学习十二——霍夫变换(Hough Transform)提取几何图形

GEE了利用霍夫变换提取几何图形

本文利用GEE提供的Hough Transform函数探测线或圆等几何图形。先来看看霍夫变换是什么？

1、霍夫变换检测直线原理
霍夫变换，英文名称Hough Transform，作用是用来检测图像中的直线或者圆等几何图形的。
一条直线的表示方法有好多种，最常见的是 y=mx+b 的形式。 假设有一幅图像，经过滤波，边缘检测等操作，怎么把这张图片中的直线提取出来。基本的思考流程是：如果直线 y=mx+b 在图片中，那么图片中，必需有N多点在直线上（像素点代入表达式成立），只要有这条直线上的两个点，就能确定这条直线。该问题可以转换为：求解所有的(m,b)组合

再来看看实例代码

// 利用Landsat8影像，计算NDVI
var image = ee.Image('LANDSAT/LC08/C01/T1_TOA/LC08_033032_20170719');
var ndvi = image.normalizedDifference(['B5', 'B4']);
// 先利用边缘检测出边界
var canny = ee.Algorithms.CannyEdgeDetector({
  image: ndvi,
  threshold: 0.4
}).multiply(255);

// 对上述边缘检测利用霍夫变换提取直线
//inputThreshold表示大于输入值像元被考虑；
//lineThreshold我理解是线上大于100个点被探测为直线
var h = ee.Algorithms.HoughTransform({
  image: canny,
  gridSize: 256,
  inputThreshold: 50,
  lineThreshold: 100
});

// 来进行图像显示
Map.setCenter(-103.80140, 40.21729, 13);
Map.addLayer(image, {bands: ['B4', 'B3', 'B2'], max: 0.3}, 'source_image');
Map.addLayer(canny.updateMask(canny), {min: 0, max: 1, palette: 'blue'}, 'canny');
Map.addLayer(h.updateMask(h), {min: 0, max: 1, palette: 'red'}, 'hough');

来看看运行结果
蓝色为边缘检测得到结果
红色为霍夫变换检测结果，明显比边缘检测得到少