GEE影像升尺度(10m->250m)

GEE影像升尺度(10m->250m)

代码

var ext = /* color: #d63000 *//* shown: false *//* displayProperties: [{"type": "rectangle"}] */ee.Geometry.Polygon([[[108.74625980473367, 28.562445155322063],[108.74625980473367, 26.566045171603736],[111.66862308598367, 26.566045171603736],[111.66862308598367, 28.562445155322063]]], null, false);var esa = ee.ImageCollection("ESA/WorldCover/v100").filterBounds(ext)
var esaProjection = ee.Image(esa.first()).projection()
//选择tree树的分类进行加载
var forest = ee.ImageCollection("ESA/WorldCover/v100").filterBounds(ext).mosaic().eq(10);
Map.addLayer(forest.selfMask(), {palette:'00ff00'}, 'forest');
print(forest)//设定默认的影像的分辨率
var forest = forest.setDefaultProjection(esaProjection)
print(forest)
var modis = ee.Image(ee.ImageCollection('MODIS/MYD11A2').first())
var modisProjection = modis.projection();
print('MODIS projection:', modisProjection);var modisProjectionAt250 = modisProjection.atScale(250)var forestMean = forest.reduceResolution({reducer: ee.Reducer.mode(),maxPixels: 1024}).reproject({crs: modisProjectionAt250});Map.addLayer(forestMean.selfMask(), {palette:['C5003D']}, 'forestmean');

注释

 ee.Geometry.Polygon([[[108.74625980473367, 28.562445155322063],[108.74625980473367, 26.566045171603736],[111.66862308598367, 26.566045171603736],[111.66862308598367, 28.562445155322063]]], null, false);

这是一个ee.Geometry.Polygon对象的实例,表示一个四边形区域。该对象由四个点的经纬度坐标组成,可以用于地理信息系统中的空间分析和处理。如果需要计算两个四边形之间的IOU,可以使用shapely.geometry.polygon库中的intersection()和area()方法来计算它们的交集和面积,然后根据IOU的定义计算它们的IOU值。具体实现可以参考引用中的代码示例。

另外,如果需要计算规则矩形框之间的IOU,可以使用引用中提到的两种方法之一,即计算两个矩形框的交集和并集的面积,然后根据IOU的定义计算它们的IOU值

var esa = ee.ImageCollection("ESA/WorldCover/v100").filterBounds(ext)

代码中的ee.ImageCollection(“ESA/WorldCover/v100”)表示获取ESA的WorldCover数据集,.first()表示获取该数据集的第一张图像。var visualization = { bands: [‘Map’], };表示可视化参数,其中bands表示使用哪个波段进行可视化。Map.centerObject(dataset);表示将地图中心设置为数据集的中心

var forest = ee.ImageCollection("ESA/WorldCover/v100").filterBounds(ext).mosaic().eq(10);

代码中的ee.ImageCollection(“ESA/WorldCover/v100”)表示获取ESA/WorldCover/v100图像集合,.filterBounds(ext)表示筛选出在ext范围内的图像,.mosaic()表示将图像集合拼接成一张图像,.eq(10)表示将图像中所有像素值为10的位置设置为1,其余位置设置为0,最终得到一个二值图像,表示森林区域。

Map.addLayer(forest.selfMask(), {palette:'00ff00'}, 'forest');

Map.addLayer(dataset, visualization, “Landcover”);表示将数据集添加到地图中进行可视化。

var modisProjectionAt250 = modisProjection.atScale(250)

建立一个250m分辨率的投影

var forestMean = forest.reduceResolution({reducer: ee.Reducer.mode(),maxPixels: 1024}).reproject({crs: modisProjectionAt250});

这段代码使用Google Earth Engine API计算了一个名为"forestMean"的变量。首先,它使用reduceResolution()方法将图像的分辨率降低到1024像素,并使用ee.Reducer.mode()方法计算每个像素的众数。然后,它使用reproject()方法将图像投影到名为"modisProjectionAt250"的投影上。最终,"forestMean"变量将是一个在新投影下的降低分辨率的图像,其中每个像素的值是原始图像中对应像素的众数。
在这里插入图片描述

备注

在这个例子中,以ESA的2020年土地覆盖数据集按阈值法处理得到森林的10m分辨率分布,,然后用modis的250M分辨率产品制作投影,将10m的每个栅格用众数计算聚合为250m分辨率的一个大栅格。即在这个250m上,表示森林的像元(1)越多,认为250M也表示森林。
这么大一个海南岛,大部分都是森林?所以地区大并不代表经济很发达,土地覆盖类型也很重要。平原地带能够发展城镇化才能提高经济。但是也不能为了钱过度破坏环境去砍树。
在这里插入图片描述

参考博文

Google Earth Engine(GEE)——将原始影像进行升尺度计算(以海南省为例)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/212944.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【MySQL】之死锁问题及其解决方案

前言 数据库死锁问题是我们老生常谈的问题了,在我们实际开发过程中经常会遇到,为了尽量避免出现死锁,我们需要了解出现死锁的场景。同时,如果线上出现了死锁之后怎么去分析、排查和解决,下面我就这两点介绍一下。 一、…

ubuntu22.04 怎么开启SSH服务

在 Ubuntu 22.04 LTS 中,默认情况下不会自动启动 SSH 服务。如果你想通过 SSH 访问你的 Ubuntu 系统,你需要手动安装 SSH 服务器,并确保 22 端口(SSH 的默认端口)是开放的。以下是必要的步骤: 安装 SSH 服…

Java 多线程之同步(锁)相关类总结

文章目录 一、概述二、volatile 可见性/有序性三、synchronized 互拆锁/排他锁/非观锁四、DCL(Double-Checked Locking)五、CAS(Compare and Set)六、ReentrantLock 可重入锁/公平/非公平锁七、ReentrantReadWriteLock 读写锁/共享…

Day56力扣打卡

打卡记录 数对统计&#xff08;DP状态压缩&#xff09; 参考文献 #include <bits/stdc.h>using namespace std;void solve(){int n;cin >> n;map<int, int> mapp;vector<int> a(n);for (auto& x : a){cin >> x;mapp[x] ;}vector<array&…

使用WebyogSQLyog使用数据库

数据库 实现数据持久化到本地&#xff1a; 使用完整的管理系统统一管理&#xff0c; 数据库&#xff08;DateBase&#xff09;&#xff1a; 为了方便数据存储和管理&#xff08;增删改查&#xff09;&#xff0c;将数据按照特定的规则存储起来 安装WebyogSQLyog -- 创建数…

101基于matlab的极限学习机ELM算法进行遥感图像分类

基于matlab的极限学习机ELM算法进行遥感图像分类&#xff0c;对所获取的遥感图片进行初步分类和最终分类。数据可更换自己的&#xff0c;程序已调通&#xff0c;可直接运行。

如何使用 Explain 分析 SQL 语句?

如何使用 Explain 分析 SQL 语句&#xff1f; MySQL中EXPLAIN命令是我们分析和优化SQL语句的利器。 如何使用EXPLAIN来分析SQL语句&#xff0c;接下来有15个例子&#xff0c;一起学习呗 1. EXPLAIN的基本使用 EXPLAIN可以用于分析MySQL如何执行一个SQL查询&#xff0c;包括如…

ElasticSearch之cat repositories API

命令样例如下&#xff1a; curl -X GET "https://localhost:9200/_cat/repositories?vtrue&pretty" --cacert $ES_HOME/config/certs/http_ca.crt -u "elastic:ohCxPHQBEs5*lo7F9"执行结果输出如下&#xff1a; id type repo1 fs repo2 s3查…

python+gdal地理坐标转投影坐标

1 前言 地理坐标系&#xff0c;是使用三维球面来定义地球表面位置&#xff0c;以实现通过经纬度对地球表面点位引用的坐标系。 地理坐标系经过地图投影操作后就变成了投影坐标系。而地图投影是按照一定的数学法则将地球椭球面上点的经维度坐标转换到平面上的直角坐标。 2 流程…

基于STM32的四位数码管计数器设计与实现

✅作者简介&#xff1a;热爱科研的嵌入式开发者&#xff0c;修心和技术同步精进&#xff0c; 代码获取、问题探讨及文章转载可私信。 ☁ 愿你的生命中有够多的云翳,来造就一个美丽的黄昏。 &#x1f34e;获取更多嵌入式资料可点击链接进群领取&#xff0c;谢谢支持&#xff01;…

Docker Compose(容器编排)——9

目录 什么是 Docker Compose生活案例为什么要 Docker ComposeDocker Compose 的安装Docker Compose 的功能Docker Compose 使用场景Docker Compose 文件&#xff08;docker-compose.yml&#xff09; 文件语法版本文件基本结构及常见指令Docker Compose 命令清单 命令清单如下命…

垃圾回收器CMS和G1的区别

CMS和G1的区别 区别一&#xff1a; 使用范围不一样 CMS收集器是老年代的收集器&#xff0c;可以配合新生代的Serial和ParNew收集器一起使用 G1收集器收集范围是老年代和新生代。不需要结合其他收集器使用 区别二&#xff1a; STW的时间 CMS收集器以最小的停顿时间为目标的收…

C++11(下)

可变参数模板 C11的新特性可变参数模板能够创建可以接受可变参数的函数模板和类模板. 相比C98/03, 类模版和函数模版中只能含固定数量的模版参数, 可变模版参数无疑是一个巨大的改进, 然而由于可变模版参数比较抽象, 使用起来需要一定的技巧, 所以这块还是比较晦涩的.掌握一些基…

Vue 3项目的运行过程

概述&#xff1a; 使用Vite构建Vue 3项目后&#xff0c;当执行yarn dev命令启动服务时&#xff0c;项目就会运行起来&#xff0c;该项目会通过src\main.js文件将src\App.vue组件渲染到index.html文件的指定区域。 文件介绍&#xff1a; src\App.vue文件 Vue 3项目是由各种组件…

递归实现指数型枚举

title: 递归实现指数型枚举 date: 2023-12-10 19:29:20 tags: 递归 catgories: 算法进阶指南 —> 传送门 题目大意 从 1 ~ n n n 这 n n n 个整数随机选取任意多个&#xff0c;输出所有可能的选择方案 思路 这等价于每个整数可以选或者不选&#xff0c;所有的方案总数共有…

Spring Boot的日志

打印日志 打印日志的步骤: • 在程序中得到日志对象. • 使用日志对象输出要打印的内容 在程序中得到日志对象 在程序中获取日志对象需要使用日志工厂LoggerFactory,代码如下: package com.example.demo;import org.slf4j.Logger; import org.slf4j.LoggerFactory;public c…

STM32——继电器

继电器工作原理 单片机供电 VCC GND 接单片机&#xff0c; VCC 需要接 3.3V &#xff0c; 5V 不行&#xff01; 最大负载电路交流 250V/10A &#xff0c;直流 30V/10A 引脚 IN 接收到 低电平 时&#xff0c;开关闭合。

Go Fyne 入门

Fyne是一个用于创建原生应用程序的UI工具包&#xff0c;它简单易用&#xff0c;并且支持跨平台。以下是一个简单的Fyne教程&#xff0c;帮助你入门&#xff1a; 1. 安装Fyne 首先&#xff0c;确保你已经安装了Go语言。然后&#xff0c;在终端中运行以下命令来安装Fyne&#x…

android-xml语法

xml解析器 Android的XML文件语法是由Android系统中的解析器解析的。具体来说&#xff0c;Android使用了一个名为"Android Asset Packaging Tool (AAPT)"的工具来解析和处理XML文件。AAPT负责将XML文件编译为二进制格式&#xff0c;并在构建过程中将其打包到Android应…

第2节:Vue3 模板语法

Vue3 的模板语法主要包括以下几个部分&#xff1a; 插值表达式&#xff1a;使用双大括号 {{ }} 包裹变量&#xff0c;可以直接在模板中显示变量的值。 <div>{{ message }}</div>指令&#xff1a;以 v- 开头&#xff0c;后面跟一个自定义的名字&#xff0c;用于操…