函数和映射
函数是一组执行特定任务的指令:
function functionName(Arguments) { statements; };
调用函数
var result = functionName(input);
在集合上映射函数
//map在GEE中非常常用 var result = input.map(functionName);
在集合上映射函数将函数应用于集合中的每个元素。
几何上的常见操作
寻找几何的面积
var geoArea = geometry.area(maxError);
默认情况下,地球引擎中的所有单位都以米为单位。
求一条线的长度
var linLen = lineString.length(maxError);
求几何的周长
var geoPeri = geometry.perimeter(maxError);
减少几何中的顶点数
//这个非常有用,因为当你上传的矢量边界过于复杂或者节点数太多的化会报错,这是一个不错的选择 var simpGeo = geometry.simplify(maxError);
寻找几何的质心
var centrGeo = geometry.centroid(maxError);
在几何体周围创建缓冲区
//这里面的第二个参就是你允许存在的最大误差 var buffGeo = geometry.buffer(radius, maxError);
查找几何的边界矩形
//这个一般会用到边界的最大框选,filterbounds一般会用来进行边界的筛选也是通过矩形的方式 var bounGeo = geometry.bounds(maxError);
寻找可以包围几何的最小多边形
//这个一般用的比较少,根据个人情况而定 //这个就是矩形的强化班,多一些线性包裹着你的研究区 var convexGeo = geometry.convexHull(maxError);
查找两个或多个几何图形之间的公共区域
//就是arcgis中的相交 var interGeo = geometry1.intersection(geometry2, maxError);
查找包含两个或更多几何图形的区域
<span style="color:var(--devsite-code-color)"><code><span style="color:var(--devsite-code-keywords-color)">var</span> unGeo = geometry1.<span style="color:var(--devsite-code-keywords-color)">union</span>(geometry2, maxError); </code></span>
示例:几何运算
让我们使用 Earth Engine 上可用的美国公共县要素集合的几何图形在美国康涅狄格州运行一些这些操作:
1.我们首先放大感兴趣的区域并通过从相应的特征中提取它们来加载/创建感兴趣的几何图形。
// 设置影像的中心点 Map.setCenter(-72.6978, 41.6798, 8); // 加载美国地图 var countyData = ee.FeatureCollection('TIGER/2018/Counties'); //过滤位于康涅狄格州的县。餐厨就是选择这个波段,这个州编号是09 var countyConnect = countyData.filter(ee.Filter.eq('STATEFP', '09')); // 获取康涅狄格州所有县几何的联合。 var countyConnectDiss = countyConnect.union(100); // 使用康涅狄格州的第一个县创建一个圆形区域 // 矢量集合建立圆形缓冲区半径5km,误差100 var circle = ee.Feature(countyConnect.first()) .geometry().centroid(100).buffer(50000, 100); // 使用指定的颜色和图层名称将图层添加到地图。 Map.addLayer(countyConnectDiss, {color: 'red'}, 'CT dissolved'); Map.addLayer(circle, {color: 'orange'}, 'Circle');
2.使用该bounds()
函数,我们可以找到包含几何最南、最西、最东和最北点的矩形。
var bound = countyConnectDiss.geometry().bounds(100); // 使用指定的颜色和图层名称将图层添加到地图。 Map.addLayer(bound, {color: 'yellow'}, 'Bounds');
3.同样,但不限于矩形,凸包 ( convexHull()
) 是一个覆盖几何体末端的多边形。
var convex = countyConnectDiss.geometry().convexHull(100); // 使用指定的颜色和图层名称将图层添加到地图。 Map.addLayer(convex, {color: 'blue'}, 'Convex Hull');
4.继续进行一些基本操作以组合多个几何图形,交集 ( intersection()
) 是两个或多个几何图形共有的区域。
var intersect = convex.intersection(circle, 100); // 使用指定的颜色和图层名称将图层添加到地图。 Map.addLayer(intersect, {color: 'green'}, 'Circle and convex intersection');
5.联合 ( union()
) 是包含两个或多个特征的区域。
// 设置最大的误差 var union = convex.union(circle, 100); // 使用指定的颜色和图层名称将图层添加到地图。 Map.addLayer(union, {color: 'purple'}, 'Circle and convex union');
6.我们还可以找到difference()
两个几何图形之间的空间差异 ( )。
var diff = convex.difference(circle, 100); Map.addLayer(diff, {color: 'brown'}, 'Circle and convex difference');
最后的图形结果基本上是这样的也就是我们所说的子集、并集、交集
7.最后,我们可以计算并显示我们几何图形的面积、长度、周长等。
// 找到计算的矢量集合并计算面积 var ar = countyConnectDiss.geometry().area(100); print(ar); // 查找线几何的长度(因为这是一个多边形,所以你得到零) var length = countyConnectDiss.geometry().length(100); print(length); // 找到特征的周长。 var peri = countyConnectDiss.geometry().perimeter(100); print(peri);