GEE|假彩色目视解译山东省玉米、水稻、小麦等样本集制作代码
相关背景:
山东种植的玉米主要是夏玉米,采用与冬小麦轮作的种植制度。一般在6月冬小麦收获后播种;9 月底至 10 月收获,收获夏玉米后便播种冬小麦。与之对应的,播种期在春季、收获在秋天的春玉米则占比很小。如果是春玉米,山东地区则一般于4月初播种,8 月末收获。玉米全生育期分为播种、出苗、三叶、七叶、拔节、抽雄、开花、灌浆、 乳熟、成熟等主要发育时期。玉米生育期是指玉米从播种到新种子成熟所经历的天数。生育期的长短因品种、播种期、光照、温度等环境条件差异而有所不同,一般品种早熟、播种晚、温度高等情况下,玉米生育期相对短一些。
主要思想:
样本部分通过图像预处理即按一定的云量覆盖进行筛选后进行去云,再按照最小值、中值、均值等方式进行影像融合,同时按照能较为清晰地反映各地物差别的波段组合函数,选择时间尺度为各个作物生育期内遥感影像,生成 real、 Corn_NIRSW1SW2、Rice_SW1NIRR、wheat_NIRRG 共四个图层。再根据这些图层经过目视解译标注玉米cron_samples和其他类地物 other_samples 两类样本点,其中,其他类包括但不限于水稻、小麦、花生、大豆、山、河流水域、城镇、交通用地、建筑用地等。每类样本点标注不少于 300 个后,生成数据点相 同的两个数据集 corn 和 other 用于后续进行随机森林分类。并将 corn 的 value 值设置为 1,other 设置为 2,方便后续分类训练。
成果展示:
代码:
//===========================加载山东省行政边界数据===========================
var shpfile= ee.FeatureCollection('users/xianda19/China_shp').filter(ee.Filter.eq('ID','Shandong'))
//===========================确定地图中心与放缩比例===========================
Map.centerObject(shpfile)
//===========================建立Landsat5去云函数============================
function maskL5sr(image) {var cloudShadowBitMask = 1 << 3;var cloudsBitMask = 1 << 5;var qa = image.select('pixel_qa');var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0).and(qa.bitwiseAnd(cloudsBitMask).eq(0));return image.updateMask(mask).divide(10000); //
}
//=================引入Landsat5地表反射率数据,并进行云量筛选====================
var composition = function(startdate,enddate,cloudiness,compositiontype,bands){
var image = ee.ImageCollection('LANDSAT/LT05/C01/T1_SR').filterDate(startdate,enddate)
.filterBounds(shpfile).filter(ee.Filter.lt('CLOUD_COVER',cloudiness))
//===================================去云===================================
image = image.map( maskL5sr)
//===============================多时相影像融合==============================
//0—使用最小值融合、1—使用均值融合、2—使用中值融合
if(compositiontype === 0){image = image.min()}
else if(compositiontype == 1){image = image.mean()}
else if(compositiontype == 2){image = image.median()}
else image = image.max()
//================================裁剪======================================
image = image.clip(shpfile)
//================================显示图层==================================
if(bands == 543) {Map.addLayer(image,{min:0, max:0.3, bands:["B5", "B4", "B3"]},'Rice_SW1NIRR'+startdate)}
else if(bands == 457) {Map.addLayer(image,{min:0, max:0.3, bands:["B4", "B5", "B7"]},'Corn_NIRSW1SW2'+startdate)}
else if(bands == 432) {Map.addLayer(image,{min:0, max:0.3, bands:["B4", "B3", "B2"]},'wheat_NIRRG_'+startdate)}
else if(bands == 321) {Map.addLayer(image,{min:0, max:0.3, bands:["B3", "B2", "B1"]},'real'+startdate)}
}
// // 目视解译识别小麦
composition('1990-03-01', '1990-04-30',100,2,432)
composition('1990-08-12', '1990-11-30',100,2,543)
// // 目视解译识别水稻
composition('1990-08-12', '1990-11-30',100,2,543)
// 目视解译识别玉米
composition('1990-07-12', '1990-09-30',100,2,457)
//真彩色合成
composition('1990-06-12', '1990-09-30',100,2,321)
//==============================下载样本数据集 =============================
Export.table.toDrive({
collection:corn_samples,
description: "corn",
fileNamePrefix: "corn_samples",
fileFormat: "KML"
});Export.table.toDrive({
collection:other_samples,
description: "other",
fileNamePrefix: "other_samples",
fileFormat: "KML"
});
GEE|假彩色目视解译山东省玉米、水稻、小麦等样本集制作代码相关推荐
- python爬取玉米、小麦、水稻信息数据到本地为网页形式和mysql数据库中
1.创建Scrapy项目 scrapy startproject ExGrain 2.进入项目目录,使用命令genspider创建Spider scrapy genspider exgrain ex- ...
- WorldCereal下玉米和小麦全球作物日历
文章地址:Global crop calendars of maize and wheat in the framework of the WorldCereal project 获取作物日历的网站: ...
- GEE哨兵二号去云不成功的原因(代码修改)
基于GEE的哨兵二号QA波段去云,但是去不成功 原因:因为筛选出来的影像云量几乎都超过50% 官方的官方的解答:文档里是提取云量20以下的进行去云的 源码: // 去云 function maskS2 ...
- 第九届(2018)蓝桥杯 山东省赛解题报告(题目+分析+代码)
1标题:第几天 2000年的1月1日,是那一年的第1天. 那么,2000年的5月4日,是那一年的第几天? 注意:需要提交的是一个整数,不要填写任何多余内容. [答案]:125 2标题:明码 汉字的字形 ...
- 应用层加密方_加密应用层数据之前要问的6个问题
应用层加密方 In the era of Facebook data scandals and Edward Snowden whistleblowers, data has never been a ...
- PHP全栈学习笔记7
图形图像处理技术,gd库的强大支持,PHP的图像可以是PHP的强项,PHP图形化类库,jpgraph是一款非常好用的强大的图形处理工具. 在PHP中加载GD库 gd官方网址下载: http://www ...
- 育种数据分析之放飞自我----关注数1万+的感言
前几天,我的公众号关注人数突破了10000,这么小众的公众号,还有这么多人信任我,值的记录一下. 公众号2016年开通,佛系更新,大都是自己的学习笔记,还有放飞自我的段子. 育种数据分析 很多人问我公 ...
- md5算出来不一样_天天吃味精,您知道是用什么做的吗?说出来可能和您想的不一样...
天天吃味精,您知道是用什么做的吗?说出来可能和您想的不一样.亲爱的好朋友们,大家好,我是大厨江一舟,今天又到了,给大家分享美食的时刻了,你们准备好了吗? 味精是咱们天天都吃的一道调味料.可以说.很多的 ...
- 山东田野稻花香 国稻种芯·中国水稻节:威海荣成引种的旱稻
山东田野稻花香 国稻种芯·中国水稻节:威海荣成引种的旱稻 新闻中国采编网 中国新闻采编网 谋定研究中国智库网 中国农民丰收节国际贸易促进会 国稻种芯·中国水稻节 中国三农智库网-功能性农业·农业大健康 ...
最新文章
- BootLoader引导程序制作及移植(一)
- c语言中return的作用_C语言简介
- php 枚举cookie内容,php如何修改cookie内容?
- 基础篇--Java IO--概览
- XML—— 验证XML文档
- leetcode 368. 最大整除子集(dp)
- 帆软决策报表JS实现点击超链切换TAB页
- 学生信息的顺序表存储
- BUMO 3.0原理以及设计.
- ShaderForge插件介绍详解
- 2008年上半年程序员考试上午真题自我汇总
- web前端登录页面设计技巧
- 一套完整的后台管理系统
- 下载腾讯视频为mp4格式
- 改造摄像头模块成开发板
- Arch Linux 下安装 st (终端模拟器)
- [wirteup] unctf-web-can_you_hacked_me
- 我的世界Bukkit插件开发(4)——命令
- OSGi模块化框架详解
- 我眼中的匈牙利命名法
热门文章
- 开发一个基于 Android系统车载智能APP
- 【git】You have not concluded your merge (MERGE_HEAD exists).
- xcode13 swift语言 ios开发 快捷代码优化方式(代码重构)例子
- 因子模型:套利定价理论APT
- 10个实用的Excel技巧
- php jwt使用案例,PHP JWT基础知识及其简单示例
- java h d,Java HijrahChronology zonedDateTime(TemporalAccessor)用法及代码示例
- python dataframe 写入到doc文件_将Python Pandas DataFrame写入Word文档
- NLP入门(八)使用CRF++实现命名实体识别(NER)
- 【sdx62】通过dtsi控制gpio高低电平实例