- 博客(27)
- 资源 (1)
- 收藏
- 关注

原创 基于GEE平台土地类型分类
GEE可以进行大尺度的地物类型分类,主要原理就是根据不同地物类型的光谱反射率的不同,目前地物分类的算法有很多种,比如随机森林算法等等,本文章主要基于随机森林算法进行大尺度地物分类,代码如下所示:var samples=ee.FeatureCollection(table);function maskL8sr(image){ var timeStart = image.get('system:time_start'); var srImageList = ee.ImageCollection('
2021-06-08 13:14:45
6208
4

原创 基于GEE与哨兵1号影像数据提取水体
在进行光学数据提取水体时,常会发现部分时间、部分区域内由于云的存在而出现大面积的空白区域,从而使得我们提取水体面积过程中存在不精确的问题。因此,基于微波数据进行云的提取就成为我们较好的一个选择,通过阅读各种文献以及实践,我在原有基于GEE平台和哨兵1号影像提取水体代码的基础上,进行了部分补充,从而使得提取水体的精度更加准确。代码如下所示:// An highlighted blockvar slope_lib = require('users/wushupu18/Nanhu:slope_correcti
2021-05-27 09:45:57
12579
19

原创 GEE学习资料汇总
相关学习资源基本学习材料:官方地址:https://earthengine.google.com/平台地址:https://code.earthengine.google.com/数据资源介绍:https://earthengine.google.com/datasets/查看数据资源:https://explorer.earthengine.google.com/#workspace入门介绍文档:https://developers.google.com/earth-engine/完整API文档
2021-05-18 20:58:12
708
1
原创 对于由面数据转换点数据的实验记录
1.对于只有一个面属性的一条属性表,我们要用multipart to singlepart命令把整个面打散。2.打散之后,需要用另一个完备数据的进行spatial join操作,使得单个属性表具有名称属性。3.合并具有相同名称的数据,这个命令为dissolve4.使用公式统计shp面积,并将其转换为点数据(Feature to point)...
2021-09-29 20:25:32
512
原创 批量对文件夹中的文件进行重命名(筛选特定条件命名)
#coding=gbkimport osimport sysdef rename(): path = input("请输入路径(例如D:\\\\picture):") name=input("请输入开头名:") startNumber = input("请输入开始数:") fileType = input("请输入后缀名(如 .jpg、.txt等等):") filelsit=os.listdir(path) print(filelsit) cou
2021-09-29 17:40:27
864
转载 批量对文件进行重命名
#coding=gbkimport osimport sysdef rename(): path=input("请输入路径(例如D:\\\\picture):") name=input("请输入开头名:") startNumber=input("请输入开始数:") fileType=input("请输入后缀名(如 .jpg、.txt等等):") print("正在生成以"+name+startNumber+fileType+"迭代的文件名") count
2021-09-27 10:42:54
249
原创 按照mask提取tif区域
# -*- coding: utf-8 -*-import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = "F:\\Third_article\\final_data\\country_pro"#shp格式数据储存的路径inraster = "F:\\Third_article\\final_data\\worldpro\\reclass\\grass_fence.tif"#tif格式数据储存的路径shps = a
2021-09-25 17:41:40
461
原创 关于栅格反向擦除的公式
Con(IsNull(“2000_dissolve.tif(擦除栅格)”),“2000.tif.tif(被擦除栅格)”,0)公式来源于:https://blog.csdn.net/u010885059/article/details/82348514
2021-09-24 14:57:25
339
原创 如何批量创建空的文件夹
import ospath="F:\\creat\\"name="9月"end="日"for i in range(1,21): print(i) filename=path+name+str(i)+end print(filename) os.makedirs(filename) print(filename+"has finished")print("mission end")
2021-09-24 10:08:44
480
原创 ArcPy按照属性提取
# -*- coding: utf-8 -*-import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = "F:\\Third_article"#tif格式数据储存的路径outputpath = "F:\\Third_article\\NDVI\\extract_ndvi"#提取属性数据输出路径rasters = arcpy.ListRasters("*.tif")#遍历tif文件for raster in ra
2021-09-23 14:28:53
1035
原创 对于ArcMap中的不同尺度的影像进行处理
对于不同空间范围影像进行交互处理时,我们需要在环境设置里面对空间范围进行处理,设置最大范围为处理空间数据的范围。
2021-08-28 12:06:58
273
原创 利用arcgis中的cell_statistics统计某个文件夹中的影像均值
import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace="C:\\Users\\DELL\Desktop\\test\\tif\\output"outputpath="C:\\Users\\DELL\\Desktop\\test\\tif\\output\\mean"rasters=arcpy.ListRasters("*.tif")listraster=[]for raster in rasters:
2021-08-26 22:36:28
2550
1
原创 关于GEE中merge命令和mosaic命令的说明
merge简单说就是合并,即数据集合的合并,mosaic即拼接,这两个命令常用作拼接不完整的影像时进行使用。可以简单的说:b4=b1.merge(b2).merge(b3).mosaic()b4返回的是b3的值,如果把merge命令中b2和b3位置进行调换时,mosaic()b4返回的就是b2的值了。所以在写代码时,弄懂每个简单的命令很重要。参考代码如下所示:https://code.earthengine.google.com/2a4fed8cb7c1b5fac9cf1ebb04f7669a
2021-08-25 11:36:40
3715
1
原创 运用python进行批量裁剪
# -*- coding: utf-8 -*-import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = "J:\\python_test\\tif"outputpath = "J:\\python_test\\out"rasters = arcpy.ListRasters("*.tif")inMaskData = "J:\\python_test\\clip_shp\\mask.shp"for raster
2021-08-22 11:10:49
633
原创 python函数块联系
分别是两种解法,第一种为按照步骤进行分块写,第二种定义了def函数体m=int(input('m='))n=int(input('n='))fm=1for pm in range(1,m+1): fm *=pmprint fmfn=1for pn in range(1,n+1): fn *=pnprint fnfmn=1for pmn in range(1,m-n+1): fmn*=pmnprint fmnprint (fm//fn//fmn)def f
2021-08-15 12:42:50
144
原创 对矢量数据属性表进行合并
当我们获得的研究区中含有多个小的行政区时,而我们需要的是小的行政区的上一级行政区域时,可以使用Dissolve工具进行合并。
2021-08-13 21:30:36
1677
原创 Arcpy遍历文件夹的代码
下面展示一些 遍历与投影操作。import arcpyimport osworkspace = r"C:\Users\DELL\Desktop\test\tif"feature_classes=[]walk = arcpy.da.Walk(workspace, type="TIF")prj="C:\\Users\\DELL\\Desktop\\test\\tif\\2002.tif"for dirpath, dirnames, filenames in walk: for filen
2021-07-28 15:23:04
1324
原创 关于ArcModel模型的构建
我们在做遥感栅格数据处理时,特别是长时间尺度数据,常常用到ArcPy进行批处理。这里运用ArcPy进行简单的栅格数据演示:import arcpyfrom arcpy import envimport os#out_folder = r"D:\Date\NewStart\Lake\Supplement\Lake30_80\reclass\out"data = r'D:\Date\NewStart\Lake\Supplement\Lake30_80\reclass\'env.workspace
2021-07-26 21:50:20
468
1
原创 关于栅格计算空值处理的思路
假如我们有两个tif文件,如果其中一个tif文件中某像元A为空值,而另一个tif文件中的对应空间位置的像元B有值,我们在进行栅格计算是,很容易得到A+B为nodata的情况。下面这个方法可以解决这种问题。Con(IsNull(“A”),0,“A”)+Con(IsNull(“B”),0,“B”)虽然这个方法能够解决空栅格值相加的问题,但是我们在计算多个栅格值均值的时候,其遇到的问题则更为复杂(即空栅格值0也参与到了运算),欢迎交流此类情况的处理。...
2021-07-25 15:52:44
2433
1
原创 使用Arcpy进行批量栅格数据相加求取均值
import arcpyfrom arcpy.sa import *arcpy.CheckOutExtension("spatial")arcpy.gp.overwriteOutput = 1# customarcpy.env.workspace = "D:\\Date\\NDVI_new"# customoutpath = "D:\\Date\\NDVI_new\\mean\\"# customoutfilename = "mean";n = 21;Sum = 0files = a
2021-07-06 19:35:48
3109
6
原创 Python函数练习(基本函数以及点的创建)
def myAdd(a,b): c=a+b return cdef myClass(num): if num>0 and num<10: return "A" elif num>=10 and num<20: return "B" else:return "C"d=myAdd(3,5)e=myAdd(14,15)print dprint ef=myClass(100)print fimport
2021-07-05 10:13:25
100
原创 ArcPy练习代码块
#if的用法#-- coding:utf8 --x=10if x>1 and x<5: print "x not here"elif x>5 and x<7: print "x not here too"elif x>7 and x<9: print "x not here tooo"else: print"you find it"#coding=utf-8strStudy=["xinjiangsuo","nanhusuo
2021-07-04 18:13:59
245
原创 S-G滤波算法应用于NDVI数据
M-G滤波算法主要是用来进行数据平滑使用,一种基于平滑时间序列数据和最小二乘原理的卷积算法,采用移动窗口的加权平均算法,但其加权系数不是简单的常数窗口,而是通过在滑动窗口内对给定高阶多项式的最小二乘拟合得出,关于M-G滤波算法网上有很多的文献,这里主要介绍在GEE平台上进行应用,代码如下:var oeel=require('users/OEEL/lib:loadAll');modis=modis.filterDate('2020-01-01','2020-12-31').filterBounds(roi
2021-06-29 15:28:53
6596
4
原创 基于GEE平台提取水体
本文主要介绍如何利用GEE平台与哨兵2号影像提取水体。水体的提取主要是基于NDWI指数进行,当然,通过改变波段运算,也可以根据MNDWI进行提取。主要代码如下:下面展示一些 内联代码片。//首先对数据进行去云处理 var s2_rmcloud = function(image) { var quality = image.select("QA60").unmask(); return image.updateMask(quality.eq(0));};//NDWIvar s2_ndwi
2021-05-21 22:17:38
10962
22
原创 关于气候序列的显示
GEE平台提供了大量的数据集,其中包括大量的气候数据集合,本文主要介绍如何显示某地区的气候时间序列。其中涉及的内容主要有两方面,第一为气候序列的获取(简单),第二为气候序列的显示(较简单),下篇文章可能会涉及Chart表的显示,欢迎各位批评指正。话不多说,上代码:var imageCollection=ee.ImageCollection("ECMWF/ERA5/MONTHLY")var pr = imageCollection.filterDate('2010-01-01','2020-12-31'
2021-05-16 00:04:16
208
原创 哨兵2号GEE去云处理
本篇文章的主要处理主要为利用GEE平台进行哨兵2号影像数据的去云处理,哨兵2号数据空间分辨率较高,常常用于小尺度的遥感分析工作,这里利用GEE平台优势,对哨兵2号数据进行去云处理,去云方法有很多种,此篇文章主要是用QA波段进行去云处理,后期会逐步对其他去云方法进行上传,欢迎交流。 下面为详细代码。function rmCloudByQA(image) { var qa = image.select('QA60'); var cloudBitMask = 1 << 10;...
2021-05-13 23:04:57
4478
1
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人