当前位置: 首页 > news >正文

济南企业营销型网站建设价格商丘网络推广公司

济南企业营销型网站建设价格,商丘网络推广公司,问卷调查网站哪个好,wordpress+企业站模版最近发现好多人都在问一个问题,两张影像如何取其相交区域?其实这个问题简单来讲就是多张栅格影像进行叠加分析。在GEE中栅格影像不像矢量数据那样有直接的函数来做数据分析,需要我们自己手动写一些代码来实现这些操作。要实现这个功能有很多方…

最近发现好多人都在问一个问题,两张影像如何取其相交区域?其实这个问题简单来讲就是多张栅格影像进行叠加分析。在GEE中栅格影像不像矢量数据那样有直接的函数来做数据分析,需要我们自己手动写一些代码来实现这些操作。要实现这个功能有很多方法,这里简单说一下有几种方式:

第一种:先将栅格影像数据转换为矢量数据,然后在利用矢量数据计算。这种方式想想可以,但是千万别在实际中使用。

第二种:直接使用影像自带的方法mask()和updateMask()实现,这种方法简单快捷,下面我会通过具体的例子来说明一下如何实现这个功能。

    下面通过A:蓝色区域,B:粉色区域,来展示不同叠置分析具体是什么样子,相信大家对下面这个图非常清楚了,就在一一解释,下面通过具体例子来说明影像如何实现这几种情况。

公共代码:

var center = /* color: #0b4a8b */ee.Geometry.Point([115.78080896985244, 37.569622364096226]);
 
var img = ee.Image("LANDSAT/LC08/C01/T1_SR/LC08_123034_20160504")
 
            .multiply(0.0001)
 
            .normalizedDifference(["B5", "B4"])
 
            .rename("NDVI");
 
var vis = {
 
  min: 0,
 
  max: 1,
 
  palette: [
 
    'FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718', '74A901',
 
    '66A000', '529400', '3E8601', '207401', '056201', '004C00', '023B01',
 
    '012E01', '011D01', '011301'
 
  ],
 
};
 
Map.centerObject(center, 8);
 
Map.addLayer(img, vis, "NDVI", false);
 
var roi1 =
 
    ee.Geometry.Polygon(
 
        [[[115.58854822766494, 37.72402725521113],
 
          [115.58854822766494, 37.48467189570746],
 
          [115.92912439953994, 37.48467189570746],
 
          [115.92912439953994, 37.72402725521113]]], null, false),
 
    roi2 =
 
    ee.Geometry.Polygon(
 
        [[[115.78630213391494, 37.62185147647586],
 
          [115.78630213391494, 37.36688853978036],
 
          [116.1406112159462, 37.36688853978036],
 
          [116.1406112159462, 37.62185147647586]]], null, false);
 
var imgA = img.clip(roi1);
 
var imgB = img.clip(roi2);
 
Map.addLayer(roi1, {color: "00ffff"}, "roi1", false);
 
Map.addLayer(roi2, {color: "ff00ff"}, "roi2", false);
 
Map.addLayer(imgA, vis, "imgA");
 
Map.addLayer(imgB, vis, "imgB");
 
 

运行结果如下:

1、AB相交,取AB交集

下面的例子计算的就是A和B的交集,然后结果返回的是A的数据结果

  1. //A B交集

  2. var intersection = imgA.updateMask(imgB.mask());

  3. Map.addLayer(intersection, vis, "intersection");

运行结果:

2、AB相交,取AB并集

    取并集可以有很多中方式,下面介绍两种常用的方式,一种是利用mask()和updateMask()来做,另外一种方式则是利用构建ImageCollection方式来做。

第一种方式,需要说明的是下面这个代码返回相交区域结果是两张影像求和的结果

  1. //(1)A B并集

  2. var mask = imgA.mask().or(imgB.mask());

  3. var baseImg = ee.Image.constant(0)

  4.                 .updateMask(mask);

  5. var newImgA = baseImg.where(imgA.mask(), imgA);

  6. var newImgB = baseImg.where(imgB.mask(), imgB);

  7. var union = newImgA.add(newImgB);

  8. Map.addLayer(union, vis, "union");

返回结果如下,点击的是箭头指示的地方,可以看到union的影像是两张影像的和,非相交的地方则是对应自己的结果。

第二种方式,构造ImageCollection,然后计算并集

  1. //(2)A B并集

  2. var imgAB = ee.ImageCollection.fromImages([imgA, imgB]);

  3. var union1 = imgAB.mosaic();

  4. Map.addLayer(union1, vis, "union1");

  5. var union2 = imgAB.sum();

  6. Map.addLayer(union2, vis, "union2");

其中:union1是普通的拼接而得到的结果,每一个像素结果都是唯一的;union2是求和计算拼接的结果,运算结果和第一种方式类似。

3、AB相交,AB的交集取反

    计算这个关键就是如何删除交集,主要逻辑就是考虑清楚下面代码mask2就可以,这个是交集的掩膜,那么只要对这个掩膜取反就可以了。

  1. // AB的交集取反

  2. var mask = imgA.mask().or(imgB.mask());

  3. var baseImg = ee.Image.constant(0)

  4.                 .updateMask(mask);

  5. var newImgA = baseImg.where(imgA.mask(), imgA);

  6. var newImgB = baseImg.where(imgB.mask(), imgB);

  7. var mask2 = imgA.mask().and(imgB.mask());

  8. var symmetricDifference = newImgA.add(newImgB)

  9.                                  .updateMask(mask2.not());

  10. Map.addLayer(symmetricDifference, vis, "symmetricDifference");

运行结果如下:

4、AB相交,取属于A但非B

具体代码如下:

  1. //取属于A但非B

  2. var difference = imgA.updateMask(imgB.mask().not());

  3. Map.addLayer(difference, vis, "difference");

运行结果:

    总结一下就是影像要实现叠置分析,最关键就是明白影像中的mask()、updateMask()、unmask()等方法使用,通过这些方法结合and、or、not等逻辑判断就可以实现自己想要的各种结果。

来源请引用:地理遥感生态网科学数据注册与出版系统

http://www.15wanjia.com/news/47307.html

相关文章:

  • 东莞做汽车有没有买票的网站软文推广怎么做
  • 网站开发建设推荐十大经典事件营销案例
  • 浪起科技做的网站怎么样税收大数据
  • 重庆网站建设中心吸引人气的营销方案
  • 视频网站开发视频中国市场营销网
  • 建立网站的风险天津百度网络推广
  • 免费搭建商城网站东莞百度快速优化排名
  • O2O网站制作需要多少钱怎么在百度上推广自己的店铺
  • 夸克搜索sem与seo的区别
  • 飓风算法受影响的网站有哪些cpa广告联盟
  • php网站调试环境搭建湖南seo优化服务
  • 2022国内外重大新闻事件10条seo优化顾问服务阿亮
  • 建设银行网站进不去关键词优化排名软件案例
  • iis7网站建设上海百度公司地址
  • 精英学校老师给学生做的网站武汉seo学徒
  • 网站建设公司推荐金石下拉g宁德市中医院
  • php商城建站系统幽默软文经典案例300
  • 网站优化说明深圳百度公司地址在哪里
  • 苏州网站建设logo软文发稿网站
  • 计算机网络工程师中级职称宁波seo外包优化
  • 政府网站开发文档百度一下百度一下
  • 建设网站需要机房吗网络销售是什么工作内容
  • 花店网站建设百度关键词推广条件
  • 网络营销软件哪个好杭州seo排名收费
  • 中石化建设工程有限公司网站最新seo新手教程
  • 如何让广域网访问利用公网ip和本地服务器建设的网站如何制作付费视频网站
  • 汕头 网站设计营销型企业网站的功能
  • 北京做网站的公司汽车网络营销推广方案
  • 缩短网址做钓鱼网站常德网站设计
  • 英文网站建设公司百度扫一扫识别图片在线