difizuhvovs
这个作者很懒,什么都没留下…
展开
-
GEE|colab报错AttributeError: module ‘tornado.ioloop‘ has no attribute ‘_Selectable‘
JavaScript版的GEE可视化功能非常强大,可惜我太菜了实在用不惯js的接口,用map做循环对我来说太难了,在打印和导出上都遇到了很多困难。实验进行不下去了…日常想撞墙。还是决定试试python版的。最方便的还是使用colab编程,可以省去配环境的繁琐步骤。了解到geemap是一个python包,拓展了GEE python版本交互式可视化的功能。原创 2022-09-15 16:41:11 · 369 阅读 · 0 评论 -
GOOGLE/DYNAMICWORLD/V1
【代码】GOOGLE/DYNAMICWORLD/V1。原创 2022-08-24 18:49:10 · 425 阅读 · 1 评论 -
GEE|时间序列分析(五)
今天用哨兵数据做了一下,昨天的代码稍微改了一下(昨天的代码又忘记去云了,写了去云函数但是忘了调用。怎么感觉,不能说有点相似,只能说毫不相关呢。原创 2022-08-22 11:01:28 · 625 阅读 · 0 评论 -
GEE|时间序列分析(四)
和前面自协方差与自相关类似,不同的是互协方差和互相关是两个变量间的,而自协方差与自相关是一个变量的不同时间的。探究NDVI与降水的相关性(5天(CHIRPS每5天一幅影像)、30天(需要将整个降水列表求和并添加为一个波段)):第一步:把时间间隔设为34天,只保留同时有t-17、t-34的影像。第二步:取34天,这样的话,但凡不是前17天和前34天都有影像的时间就是空的(它的一个弊端):先解释下面这段代码的意思。first为影像集合中第一个影像(2021.1.8),执行迭代器(影像集合中每个影像都加进去原创 2022-08-21 21:23:40 · 628 阅读 · 0 评论 -
GEE|时间序列分析(三)
学习时间序列分析的第三天。前两天自作主张选了另一个数据集,不知道代码哪里改错了,结果有点离谱。还是换回教程上的数据集,先把流程跑通再改。原创 2022-08-21 10:15:03 · 605 阅读 · 1 评论 -
GEE|ee.ImageCollection.fromImages()用法示例(求MODIS上的NDVI年平均值)
ee.ImageCollection.fromImages()是把几个影像(image)组成影像集合(imageCollection)。代码来自GEE【ee.ImageCollection.fromImages()】的示例代码,做了改动。原创 2022-08-20 16:29:43 · 889 阅读 · 0 评论 -
GEE|ee.Filter.maxDifference()用法示例
分别打印了符合要求【外来入侵植物物种覆盖比例和非外来入侵植物物种覆盖比例相等】、【'InvAnnForbCover_AH’值和’NonInvAnnForbCover_AH’值不相等】、【'InvAnnForbCover_AH’值<’NonInvAnnForbCover_AH’值】、【'InvAnnForbCover_AH’值<='NonInvAnnForbCover_AH’值】、【'InvAnnForbCover_AH’值>'NonInvAnnForbCover_AH’值】、原创 2022-08-20 14:36:06 · 308 阅读 · 0 评论 -
GEE|时间序列分析(二)
还是参考了上篇文章的文档和视频(GEE官方的视频教程)。原创 2022-08-20 12:15:13 · 1105 阅读 · 0 评论 -
GEE|.updateMask()用法示例
哨兵二号的SWIR1波段可以区分干土和湿土、水体和陆地、扑灭森林火灾等。当B4、B3、B2全大于200时,即全为1时,才会成为掩膜。.updateMask()也可以掩膜多波段。的话,就是用landMask=1的部分做了掩膜。,即SWIR1大于100的值为1,其余为0。代码是GEE的官方文档里的示例代码。原创 2022-08-19 18:17:53 · 395 阅读 · 0 评论 -
GEE|时间序列分析(一)
qa波段的第三位和第四位分别表示云影和云,原创 2022-08-19 19:44:11 · 1022 阅读 · 3 评论 -
GEE|在地图上显示传入的shp数据
【代码】GEE|在地图上显示传入的shp数据。原创 2022-08-19 17:08:49 · 2119 阅读 · 0 评论 -
GEE|哨兵L1C和L2A求出来的NDVI到底差距多大
这个问题我之前问过我的导师,早些时候没有SR数据,大家都用TOA数据,后来有了SR数据,很多研究就开始采用SR数据了,理论上SR更准确,具体来说大气校正前后Blue波段反射率有较大差异,其他波段也有差异,但是这种差异相对比较小,所以计算NDVI之类的指数,影响不大,但如果计算有Blue波段的指数,那就有影响了。对于植被的话,差别还是有的,毕竟NDVI的范围只有[-1,1],差个十万八千里也不可能,但是经过这一对比实验,应该不用再纠结到底要不要选用大气校正后的数据了。还是用具体实验看看吧。原创 2022-08-19 12:09:14 · 1473 阅读 · 1 评论 -
GEE|typeof、ee.Algorithms.If、ee.Algorithms.IsEqual 语法
筛选2018-05-01至2018-06-01期间的云量小于10%的哨兵二号影像,s2_collection.size()为0,打印它的数据类型发现其为object(对象)。想知道s2_collection.size()是否为0,要用。ee.Algorithms.If——相当于if语句,ee.Algorithms.If(条件,如果正确则,如果错误则)ee.Algorithms.IsEqual——比较两个对象是否相等。,为0则打印*ture*,否则打印*false*。typeof——查看数据类型。......原创 2022-08-18 21:56:32 · 488 阅读 · 0 评论 -
GEE|求NDVI改进版
笔记原创 2022-06-01 19:29:50 · 194 阅读 · 0 评论 -
GEE|用随机森林进行农作物分类的一次尝试
前言:用别人已经分类好的结果来分类,小编也不知道自己在干什么呢!~_~前期准备:把2017-2019年的分类结果按掩膜提取到自己的研究区,并用重分类提取出玉米、大豆栅格。再取交集找到这三年都种玉米或都种大豆的地块。先上一下这部分结果:17年玉米种植区分布:18年玉米种植区分布:19年玉米种植区分布:取交集后:呃呃。。。怎么重叠部分这么少呢~_~相信大家懂得都懂,这里就看破不说破……………………………………………………既然如此那好像前面都有些多余,那就单选一年的试一下吧。把2019年的原创 2022-05-10 21:44:36 · 5723 阅读 · 3 评论 -
GEE|220408问题记录
把gee上简单的代码用python复杂实现(?)时,发现ndvi平均值求出来的结果不一样。GEE部分代码:var elev_img_mean=reprojected_20.reduceRegion({ reducer:ee.Reducer.mean(), geometry:roi.geometry()});print('The mean Elevation in Pengdun is(m):',elev_img_mean.get('NDVI_mean'))python部分代码:def原创 2022-04-08 09:37:42 · 1423 阅读 · 0 评论 -
GEE|下载研究区哨兵二号影像
任何一个遥感人不用GEE下载遥感影像我都会伤心的ok?!!众所周知哨兵二号影像可以在usgs上或者欧空局网站下载,usgs感觉还相对更方便些。但是都还是避免不了下载影像后先用SNAP预处理大气校正重采样再波段组合什么的才能得到L2A级的tif文件。现在有一个很方便的方法摆在我们面前,只需要十几行代码分分钟下载到预处理后的tif影像,还可以一步解决裁剪、投影、重采样的操作。GEE界面:总之就是先导入你研究区的shapefile和哨兵二号L2A,运行下面代码就好了(注释写得很详细了!!)。var b原创 2022-03-22 11:00:57 · 11645 阅读 · 65 评论 -
GEE|分级统计面积
因为只是存档给自己看的就不加那么多注释和说明了。总之,这段代码用语言描述就是,把[-1,1]分成[-1,0.2],[0.2,0.4],[0.4,0.6],[0.6,0.8],[0.8,1]五级,并统计这五类的面积占比。不知道自己这种在GEE中用“累加再相减”来划分等级的方法是不是最简单的,目前还是不知道怎么在GEE中做到像ArcGIS那样重分类。需要注意的是在GEE中将字符串转为数字要用:ee.Number.parse();在GEE中加减除用:.add() .subtract() .divide()。原创 2022-03-17 16:32:07 · 1225 阅读 · 0 评论 -
GEE|求NDWI
参考视频:https://www.youtube.com/watch?reload=9&v=Mptm1_TUdPY参考代码:https://code.earthengine.google.com/8c6c9d09d1f7e0a040fdce2ab0f729f6// 加载影像var image = ee.Image('LANDSAT/LC8_L1T_TOA/LC81270512015358LGN00');// 定义可视化参数//min、max、gamma好像是拉伸、调亮度饱和度的?var原创 2021-01-22 11:55:08 · 968 阅读 · 2 评论 -
GEE|获取某地区夜光数据并导出
VIIRS Stray Light Corrected Nighttime Day/Night Band Composites Version 1var image = ee.Image(VIIRS_Stray.filterDate("2019-10-01","2020-01-07").first());var viirscollection1 = image.clipToCollection(Victoria);var visParams = {bands:['avg_rad'],gain:100.原创 2021-02-06 11:11:54 · 2840 阅读 · 11 评论 -
GEE|导出到google drive的文件太大,每次下载都中断的解决方法
multcloud网址:https://www.multcloud.com/multcloud无需翻墙,账号注册很简单,一分钟即可成功。在邮箱点击激活邮件后登录multcloud,然后连接google drive(此时需要翻下墙)。连接好后就不用翻墙了。选中连接的google drive中要下载的文件,右键点击download,就可以下载了,不会像翻着墙在Google drive中下载总是下载到一半就失败了。...原创 2021-01-30 21:01:46 · 4769 阅读 · 3 评论 -
GEE|账号注册
第一次秒拒,不过第二次成功了!ㅠㅠ网址:https://earthengine.google.com/申请页面:第一次申请时,我的账户没有与学校邮箱绑定,Full name没有与注册时候的first name与last name一致,而且还勾选了最致命的“I am interested in commercial use of Earth Engine”!第二次痛定思痛,重新申请了一个账号,first name直接写的我的姓的英文,last name写的我的名的英文,注册到最后填生日的时候,那个原创 2021-01-20 16:32:03 · 8958 阅读 · 18 评论 -
GEE|求哨兵2号NDVI并保存到本地
第一步:导入研究区shapefile文件与哨兵2号L2A(经过了大气校正)数据集第二步:代码//哨兵2去云处理var s2_rmcloud = function(image) { var quality = image.select("QA60").unmask(); return image.updateMask(quality.eq(0));};//在原影像中加入一个新的波段,命名为NDVIvar s2_ndvi = function(image) { return image.a原创 2021-01-28 19:32:50 · 3419 阅读 · 28 评论 -
GEE|获取某地区DEM并导出
数据集:NASA SRTM Digital Elevation 30mvar Victoria_SRTM = SRTM.clipToCollection(Victoria);Map.centerObject(Victoria);Map.addLayer(Victoria,{color:'FFFF00'},'Victoria');var visParams = {min: 0,max: 3000,gamma: 1.4,};Map.addLayer(Victoria_SRTM, visParams,原创 2021-02-06 11:19:00 · 4002 阅读 · 5 评论 -
211124|GEE下载EVI
//返青期var month02 = ee.ImageCollection("MODIS/006/MOD13Q1") .filterBounds(roi) .filterDate("2009-02-10", "2009-03-10") .select('EVI');var month02_mean = month02.reduce(ee.Reducer.mean())原创 2021-11-24 10:21:02 · 2004 阅读 · 2 评论 -
GEE|获取某地区遥感影像
第一步:导入中国省级行政区划Assets->New这里可以看进度:第二步:导入Landsat 8 遥感影像第三步:导入你想获取地区的一个点使用GEE地图编辑栏的添加点线面的工具,在你想要获取的地区点一个点导入这三个后的界面:第四步:代码//筛选出点所在地区var cq = China.filterBounds(point).geometry();//筛选影像时间为2018-01-01至2019-12-31var Data_Collection = imageColle原创 2021-01-26 11:59:20 · 8321 阅读 · 44 评论 -
GEE|求某地最高/低点的高度
导入某地边界和dem数据集:var elev_img_max=dem.reduceRegion({ //改成min、mean也可以求最小值、平均值等 reducer:ee.Reducer.max(), geometry:pengdun.geometry()});print('The maxium Elevation in Pengdun is(m):',elev_img_max.get('elevation'))结果显示在右侧:...原创 2021-10-07 21:52:04 · 140 阅读 · 0 评论 -
GEE|直接输出某地区NDVI最大值
终于坚持到十点了,黄金周前夕还死撑到这个点的还有谁、?然而强扭的瓜不甜,实在搞不下去了,明天再看这个错怎么改吧,参考https://www.cnblogs.com/yhpan/p/14160221.htmlvar month02 = ee.ImageCollection("MODIS/006/MOD13Q1") .filterDate("2020-02-10", "2020-03-10") .select('NDV原创 2021-09-30 22:13:39 · 809 阅读 · 4 评论 -
GEE|利用裁剪获取某地区遥感影像的方法
第一步:导入一个矢量图形第二步:代码//去云处理 function maskL8sr(image) { var cloudShadowBitMask = (1 << 3); var cloudsBitMask = (1 << 5); var qa = image.select('pixel_qa'); var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0) .and(qa.bitw原创 2021-01-27 11:06:38 · 7138 阅读 · 27 评论