- 博客(34)
- 资源 (1)
- 收藏
- 关注
RSS订阅原创 GEE之S-G滤波
//影像去云等预处理/** * Function to mask clouds using the Sentinel-2 QA band * @param {ee.Image} image Sentinel-2 image * @return {ee.Image} cloud masked Sentinel-2 image *///ROIvar pt=ee.Geometry.Point([119.85872739185788,33.06337427687241]);var fc=ee..
2023-06-30 16:54:54
893
3
原创 python 带置信区间的折线图
今天来尝试一下学以致用，资料参考于https://www.jianshu.com/p/4d7d7e4ff4f8jupyter notebook对新手非常友好import seaborn as snstips = sns.load_dataset("tips")tips两个非常有用的快捷键（shortcut）Tabshift+Tab可以逐步查看sns.regplotsns.regplot(x="total_bill",y="tip",data=tips)尝试打开自己的工作表时
2022-05-07 20:25:45
2708
原创 python之package
推荐一个网址https://pypi.org/search for “openpyxl” on this websitecopy the commander, and paste it on the terminal of PyCharm不知道为什么我的terminal打开后有红色字体的提示，先不管了，如果报错，再来解决这个问题）结果，我应该在anaconda安装过了果然。先看看能不能用吧1.右键当前工程文件“HelloWorld"2.open in explorer(系统不同，命
2022-05-06 20:38:25
707
原创 python之class
当构建一个类时，相当于用语言去描述真实世界的一个对象，如Person，这个Person有名字name，年龄age，会执行动作saying，我最开始的代码是这样的有一个Maryclass Person(): def __init__(self, name,age) : Person.name = name Person.age= age def saying(self): print("hi, I am " + self.name+".
2022-05-05 19:46:10
196
原创 python之function
将具有特定功能的多行语句组合成一个function，以后就可以重复利用这个function啦，具体代码如下def greeting(): print("start") name= input ("what is your name?") print (name +", you are such a lucky dog!") print("see you")start= input ("do you wanna play a game? yes or no: ")if s
2022-05-05 18:17:57
301
原创 python的第一个尝试
今天安装了PyCharm，下载网址为https://www.jetbrains.com/pycharm/download/#section=windows安装后，界面如下：点击 “new project”，修改保存路径（可选）click on “new project”, and then change the saving location (optional)在新建前，查看一下base interpreter是否正确，应该选择自己安装python的版本（应该是3.0以上版本了），然后crea
2022-05-04 20:37:18
814
原创 1.简单介绍R语言
下载网址https://cran.r-project.org/运行的界面如下：查看帮助文档，可以键入help.start()并按下enter键执行先随便尝试一个例子> x<-rnorm(100)> y<-rnorm(x)> plot(x,y)用#来注释，用<-来赋值而不是用=...
2021-07-17 14:43:07
88
原创 GEE按属性选择矢量
//在ROI矢量文件中选出小麦print(table);var wheat=ee.FeatureCollection(table) .filterMetadata('CLASS_NAME',"equals",'wheat');print(wheat);
2021-05-15 15:47:29
1406
原创 GEE时序NDVI
//define ROIvar ROI= ee.FeatureCollection(table);//imagevar startDate='2020-10-01', endDate='2021-05-10';var rapeseed = ee.ImageCollection('COPERNICUS/S2_SR') .filterDate(startDate, endDate) .filterBounds(ROI);//cal
2021-05-15 15:45:51
2598
原创 GEE之HSV变换
官方的hsv变换实例是对单幅影像进行操作的，对影像集进行HSV变换，需要先定义function，再map，否则会出现整幅影像的值都是（0，0，1）的异常情况。核心功能如下：先定义一个function，再遍历function HSV_Trans(image){var rgb = image.select('B4', 'B3', 'B2');var huesat = rgb.rgbToHsv();return huesat;}var HSVset=dataset.map(HSV_Trans);
2021-05-15 15:44:42
170
原创 GEE显示行政边界
参考无形的风整个思路就是首先，上传矢量边界，我一般都是传到asset，然后一键导入。然后，加载需要的影像，进行去云、日期和地区筛选等操作。接着，设置影像的显示参数。然后，通过构造一个空的image来加载矢量边界。具体代码如下：//前面部分是加载影像，可以略过，记得要上传ROI/** * Function to mask clouds based on the pixel_qa band of Landsat 8 SR data. * @param {ee.Image} image in
2021-05-15 10:22:02
2851
4
原创 GEE上传矢量文件失败
错误提示：Task details: Ingest table: “projects/earthengine-legacy/assets/users/20152600056/8”State: FailedStarted: 2h ago (2021-04-29 19:09:31 +0800)Runtime: 17sAttempts: 1Id: RITWCSG5XCJJ7RYUFOZU66DXError: Projection exception. Ensure the projection is
2021-04-29 20:44:56
5693
8
原创 GEE 使用矢量裁剪影像
先将需要的矢量文件上传到GEE//矢量文件var fc=ee.FeatureCollection(table);//sentinel-2影像，筛选时间var dataset=ee.ImageCollection('COPERNICUS/S2_SR') .filterDate('2020-03-01','2020-04-30');//对这段时间的影像集作中值处理var median=dataset.median();//裁剪var clip=median.clipToCo
2021-04-22 09:00:05
7906
13
原创 关于ArcMap中打不开dat文件的解决办法
我记得ArcMap是可以打开dat文件的，但最近显示不出来，也没有办法拖进去。解决办法：存储路径中不要有中文，命名方式最好也都改成英文。
2021-04-13 09:23:12
4604
原创 excel提取文字、字母、符号
我收到了一份非常乱的数据，需要把地块名（字母+数字+符号“-”）和植被（中文）分开，不能用excel的分列功能中的固定宽度和字符分隔来实现参考知乎用户https://zhuanlan.zhihu.com/p/89321644首先要在excel中打开VBA窗口，点击“插入-模块”将下列代码复制粘贴，代码如下：Function TQ(rng As Range, Optional i As String = "数字")With CreateObject("vbscript.regexp")
2020-11-26 11:17:49
3192
原创 监督分类
打开文件，以5，4，3波段（近红，红，绿）组合增强显示。右键图层，单击new region of interest可以用各种性状工具除了plant，如果还想再增加ROI类型，有两种方法最后的分类样本选取如下options->compute ROI separability, select all items, ok参数的值为0~2.0，大于1.9说明样本之间可分离性好，属于合格样本，小于1.8，需要重新选择样本，小于1，考虑将两类样本合成一类样本（在ROI tool对话框中，
2020-11-02 19:47:05
10076
原创 空间增强和图像融合
一、空间增强卷积滤波（convolutions）Kernel Size（卷积核大小）以奇数来表示，如 3×3、5×5 等，有些卷积核不能 改变大小，包括 Sobel 和 Roberts。默认卷积核是正方形，如果需要使用非正 方形，选择 Option￫Square kernel。Image Add Back输入一个加回值（add back）。将原始图像中的一部分“加回” 到卷积滤波结果图像上，有助于保持图像的空间连续性。该方法经常用于图像 锐化。 “加回”值是原始图像在结果输出图像中所占
2020-10-29 22:17:16
2078
1
原创 光谱增强处理
光谱增强是基于多光谱数据度波段进行变换达到图像增强处理，如可以直接利用多光谱数据的不同波段组合方式改善显示效果；或采用主成分变换、独立成分变换、色彩空间变换和色彩拉伸等光谱处理方法进行光谱增强。...
2020-10-28 22:41:26
3873
6
原创 辐射增强
辐射增强处理是通过对单个像元的灰度值进行变换处理来增强处理，如直方图拉伸、直方图匹配、去除条带噪声处理。交互式直方图拉伸display->custom stretch单波段影像：拉伸前后对比多波段影像切换拉伸方法坏道填补（replace bad lines)去条带处理（destripe data)要求数据是原始格式（平行条带），没有被旋转或地理坐标定位。...
2020-10-27 09:40:03
1858
原创 arcmap中无法导入excel的问题
当把excel拖入arcmap时，出现下列提示框。方法一：添加驱动见博客方法二：用wps打开excel，另存为dbf文件，然后拖入arcmap中。我记得还有第三种方法，暂时没想起来，后续想起来了再补回来
2020-10-26 11:15:28
4125
原创 几何校正（image to image）
这次用envi classic ，打开文件：在主菜单中，map->registration->image to image, 选择窗口2的spot影像作为基准图像，TM影像作为待校正图像。然后就可以开始选择控制点了，一般都会选择道路交叉点、桥梁、建筑物等变化不大的地物类型。在控制点选择的窗口上，点击show list,可以看到已选控制点的误差，如果觉得精度可以（看具体需求，本科作业是不能超过1），且控制点分布均匀，就可以校正了。如果精度不行，如#4、#5、#7点，误差较大，需要删除
2020-10-25 15:44:43
4654
1
原创 ArcGIS新建要素丢失
最近在新建的点（或面图层），添加了几个点后，点击save，就会在属性表中丢失该点，但在图层中还能看到该点，就是在属性表中找不到该点的属性。只有选中了该点，才能在属性表的“已选”视图当中看到该点属性。在“全部”视图中仍然看不到该点属性。前几周，我还笨笨地删掉该点重新新建，再改属性。今天偶然发现了一个解决办法。保存编辑以后，关闭属性表，重新打开属性表，即可。...
2020-10-24 20:12:59
5132
原创 裁剪与拼接
1.图像裁剪file->save as,这里以spatial subset（空间裁剪）为主，spectral subset略过不讲。裁剪有多种方式，一种是规则的矩形裁剪，一种是通过其他文件（ROI或矢量文件）裁剪。（1）通过矩形裁剪设置好裁剪范围后，就可以点击OK了。点击ok后，在窗口可以设置保存的文件格式，设置好输出路径及文件名即可。（2）通过其他文件裁剪设置输出文件名及路径即可拼接打开文件文件列表如图，有很多个波段，如果想打开多波段影像，可以打开带有“_MTL"的文
2020-10-21 21:25:59
1902
原创 ASD光谱仪数据的查看
使用ASD地物光谱仪测得光谱数据后，进行室内处理。地物光谱仪获得的数据如下图：所需软件：View SpecPro Graph打开文件安装好软件后，打开软件，file->open，选择想要打开的数据，在本例中，全选，打开设置默认文件夹setup->input directoriy,设置打开数据默认的文件夹，软件会询问默认的输出文件夹是否设为一样，可以选择yes，也可以另选路径。查看光谱选中数据（可选中一条数据或多个数据），点击菜单栏view->Graph Data，可以查
2020-10-21 15:20:46
8419
8
原创 文件打开等基本操作（一）
本科学的envi是经典版的4.8，现在需要用到5.0以上版本，然而两种版本的界面差异太大，不太习惯，因此重温一遍，熟悉一下envi5.3的操作界面。envi5.3界面打开影像点击“file”->open,选择打开的文件，可以打开无后缀的文件，也可以打开后缀为hdr的头文件。选好文件后，点击打开，会弹出“data manager”窗口，在此窗口，可以看到影像的波段，比如，示例文件中就有六个波段，分别是band1, band 2, band 3, band 4, band 5, band 7.
2020-10-16 17:35:17
3034
原创 第二期 GEE授权遇到的问题及解决办法
前几天有人问我授权代码的事情，我配置好开发环境后，就没继续学，也没敲代码测试一下，还不知道授权代码是什么，后来看官方文档，才知道必须授权才能运行。错误原因如果不授权，直接在jupyter notebook上直接敲代码时，会要求先授权，授权会出错。当我按照官方文档，对Conda进行GEE授权时，同样出错，如图：请教了大佬，大佬甩给我一个链接https://groups.google.com/g/google-earth-engine-developers/c/GeO3d6Utyto问题大概在于，
2020-08-30 23:51:28
9144
原创 高分屏软件字体过小
部分软件字体过小的原因自从买了13.3寸轻薄本后，软件的字体过小，尤其是工具栏，如图：在网上搜了，原因大概是屏幕分辨率过高，而屏幕尺寸小，导致像素点小而密集。设置方法1.更改字体大小设置->显示->更改文本、应用等项目的大小，选择合适大小即可。2.更改软件的DPI在Windows找到相应的软件，右键，打开文件所在位置，此时的文件是快捷方式所在的文件，选中对应的软件，再次右键，打开文件所在位置，此时打开的是软件安装的位置，右键->属性->兼容性->更改高DPI设
2020-08-25 19:14:14
3603
原创 GEE学习资源
1.官方文档https://developers.google.com/earth-engine/guides/getstarted2.B站可关注以下用户：geemap https://space.bilibili.com/527404442/底数e https://space.bilibili.com/173893517/3.GitHub:https://github.com/giswqs/geemap4.CSDN论坛https://www.csdn.n...
2020-08-22 10:23:35
700
原创 第一期 GEE开发环境配置
参考B站用户geemap的视频和网页GitHub、CSDN博文并结合自己的实际操作补充内容。安装geemap之前的准备工作：1. 安装anaconda有两种方式可以安装geemap，一种是pip，另一种是conda， 我选用了anaconda。在win10上安装anaconda的时候，要以管理员身份运行，并且要添加到path environment variable(第一个框默认为空，记得勾选），没有规定一定要安装在C盘，也可以安装在D盘。安装好anaconda后，运行Anaconda Promp
2020-08-21 14:49:31
3222
7
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人