- 博客(27)
- 资源 (1)
- 收藏
- 关注
RSS订阅
原创 基于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
5053
4
原创 基于GEE与哨兵1号影像数据提取水体
在进行光学数据提取水体时，常会发现部分时间、部分区域内由于云的存在而出现大面积的空白区域，从而使得我们提取水体面积过程中存在不精确的问题。因此，基于微波数据进行云的提取就成为我们较好的一个选择，通过阅读各种文献以及实践，我在原有基于GEE平台和哨兵1号影像提取水体代码的基础上，进行了部分补充，从而使得提取水体的精度更加准确。代码如下所示：// An highlighted blockvar slope_lib = require('users/wushupu18/Nanhu:slope_correcti
2021-05-27 09:45:57
11143
18
原创 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
623
1
原创 对于由面数据转换点数据的实验记录
1.对于只有一个面属性的一条属性表，我们要用multipart to singlepart命令把整个面打散。2.打散之后，需要用另一个完备数据的进行spatial join操作，使得单个属性表具有名称属性。3.合并具有相同名称的数据，这个命令为dissolve4.使用公式统计shp面积，并将其转换为点数据（Feature to point）...
2021-09-29 20:25:32
447
原创 批量对文件夹中的文件进行重命名（筛选特定条件命名）
#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
722
转载 批量对文件进行重命名
#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
201
原创 按照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
389
原创 关于栅格反向擦除的公式
Con(IsNull(“2000_dissolve.tif(擦除栅格)”),“2000.tif.tif(被擦除栅格)”,0)公式来源于：https://blog.csdn.net/u010885059/article/details/82348514
2021-09-24 14:57:25
285
原创 如何批量创建空的文件夹
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
428
原创 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
909
原创 对于ArcMap中的不同尺度的影像进行处理
对于不同空间范围影像进行交互处理时，我们需要在环境设置里面对空间范围进行处理，设置最大范围为处理空间数据的范围。
2021-08-28 12:06:58
208
原创 利用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
2330
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
3450
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
591
原创 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
109
原创 对矢量数据属性表进行合并
当我们获得的研究区中含有多个小的行政区时，而我们需要的是小的行政区的上一级行政区域时，可以使用Dissolve工具进行合并。
2021-08-13 21:30:36
1532
原创 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
1177
原创 关于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
380
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
2183
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
2939
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
72
原创 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
206
原创 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
5979
3
原创 基于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
10027
21
原创 关于气候序列的显示
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
152
原创 哨兵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
4228
1
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人