- 博客(157)
- 资源 (3)
- 收藏
- 关注
RSS订阅原创 用sentinel-2光学影像通过GAN生成Sentinel-1影像
参考的代码是：https://github.com/ZirvedaAytimur/GAN_Electro-Optic_to_SAR_Images。记录在学习用sentinel-2光学影像通过GAN生成Sentinel-1影像时遇到的问题。
2023-03-07 22:23:11
402
原创 ArcMap|以表格显示分区统计 ERROR:999999 解决方法
解决方法是：把区域字段换成FID居然就可以了。如果还不行的话 ，要记得输出表的位置也改成系统默认的，路径都是英文。问题描述：当我按照下图操作时，报了999999的错误。
2023-02-10 10:14:11
3343
2
原创 RS|下载RapidEye(5m分辨率)/Planetscope (3m)的网站申请流程（planet.com）
参考：高分辨率遥感影像在哪免费下载？- Endless的回答 - 知乎 https://www.zhihu.com/question/373272521/answer/2786776699。第三步：提交，出现此弹窗，等待邮件回复。（收到邮件回来更新）第三步一开始我没有弹出该窗口，在浏览器做以下修改后弹窗成功。
2023-02-07 21:26:49
1280
1
原创 ArcMap|多值提取至点 ERROR 001143 后台服务器抛出异常
报这个错以为是我的点太多了，通过查阅资料把后台处理中的【启用】取消勾选后arcmap就开始闪退。然后发现我的点和栅格一个是地理坐标系一个是投影坐标系。结果把点缩小到100个也还是报错。
2022-10-24 11:12:08
3018
3
原创 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
444
原创 GEE|时间序列分析（五）
今天用哨兵数据做了一下，昨天的代码稍微改了一下（昨天的代码又忘记去云了，写了去云函数但是忘了调用。怎么感觉，不能说有点相似，只能说毫不相关呢。
2022-08-22 11:01:28
729
原创 GEE|时间序列分析（四）
和前面自协方差与自相关类似，不同的是互协方差和互相关是两个变量间的，而自协方差与自相关是一个变量的不同时间的。探究NDVI与降水的相关性（5天(CHIRPS每5天一幅影像)、30天（需要将整个降水列表求和并添加为一个波段））：第一步：把时间间隔设为34天，只保留同时有t-17、t-34的影像。第二步：取34天，这样的话，但凡不是前17天和前34天都有影像的时间就是空的（它的一个弊端）：先解释下面这段代码的意思。first为影像集合中第一个影像（2021.1.8），执行迭代器（影像集合中每个影像都加进去
2022-08-21 21:23:40
734
原创 GEE|时间序列分析（三）
学习时间序列分析的第三天。前两天自作主张选了另一个数据集，不知道代码哪里改错了，结果有点离谱。还是换回教程上的数据集，先把流程跑通再改。
2022-08-21 10:15:03
714
1
原创 GEE|ee.ImageCollection.fromImages()用法示例（求MODIS上的NDVI年平均值）
ee.ImageCollection.fromImages()是把几个影像（image）组成影像集合（imageCollection）。代码来自GEE【ee.ImageCollection.fromImages()】的示例代码，做了改动。
2022-08-20 16:29:43
1019
原创 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
365
原创 GEE|.updateMask()用法示例
哨兵二号的SWIR1波段可以区分干土和湿土、水体和陆地、扑灭森林火灾等。当B4、B3、B2全大于200时，即全为1时，才会成为掩膜。.updateMask()也可以掩膜多波段。的话，就是用landMask=1的部分做了掩膜。，即SWIR1大于100的值为1，其余为0。代码是GEE的官方文档里的示例代码。
2022-08-19 18:17:53
506
原创 GEE|哨兵L1C和L2A求出来的NDVI到底差距多大
这个问题我之前问过我的导师，早些时候没有SR数据，大家都用TOA数据，后来有了SR数据，很多研究就开始采用SR数据了，理论上SR更准确，具体来说大气校正前后Blue波段反射率有较大差异，其他波段也有差异，但是这种差异相对比较小，所以计算NDVI之类的指数，影响不大，但如果计算有Blue波段的指数，那就有影响了。对于植被的话，差别还是有的，毕竟NDVI的范围只有[-1,1]，差个十万八千里也不可能，但是经过这一对比实验，应该不用再纠结到底要不要选用大气校正后的数据了。还是用具体实验看看吧。
2022-08-19 12:09:14
1888
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
594
原创 Time-weighted dynamic time warping (TWDTW) method
文章地址：High-resolution map of sugarcane cultivation in Brazil using a phenology-based method验证样本不一定必须要实地采集：
2022-08-11 20:06:23
3856
原创 WorldCereal下玉米和小麦全球作物日历
文章地址：Global crop calendars of maize and wheat in the framework of the WorldCereal project获取作物日历的网站：
2022-08-09 21:51:16
606
原创 ArcMap|用字段计算器对不同类别的id赋值
一开始我的id有四类，分别是1、2、3、4，类型是空的。然后我想分别将1、2、3、4的类型赋为水稻、其它、玉米、大豆。具体操作步骤是：①右键【类型】，打开字段计算器。②这是操作界面。③在【预逻辑脚本代码】输入以下代码。这个代码就是一个简单的python函数。④在【类型】输入以下代码。⑤点击确定即可。...
2022-06-23 21:10:24
1887
原创 DL|循环神经网络部分
一、序列模型时间序列预测方式有自回归模型（与过去所有时间有关）、马尔可夫模型（只与x个过去数据有关）、潜变量模型（h+x->新的h，再循环）。多层感知机实现马尔可夫模型：%matplotlib inlineimport torchfrom torch import nnfrom d2l import torch as d2lT = 1000 # 总共产生1000个点time = torch.arange(1, T + 1, dtype=torch.float32)x = to
2022-05-26 16:47:58
174
原创 LeetCode|3. 无重复字符的最长子串
最近没什么事干，又开始做leetcode了，对一个专业与写代码沾边的非科班出身边缘人来说，每次下定决心刷leetcode都是一次从入门到放弃的过程。。思路：用dfs暴力搜索先找到所有子串，然后用集合判断是否有重复字符，再比较长度。语法：获取字符串中某一字符：s.charAt(i)字符串无重复字符个数：import java.util.HashSet;import java.util.Set;public class text { public static void main(Str
2022-05-12 12:12:58
244
原创 GEE|用随机森林进行农作物分类的一次尝试
前言：用别人已经分类好的结果来分类，小编也不知道自己在干什么呢！~_~前期准备：把2017-2019年的分类结果按掩膜提取到自己的研究区，并用重分类提取出玉米、大豆栅格。再取交集找到这三年都种玉米或都种大豆的地块。先上一下这部分结果：17年玉米种植区分布：18年玉米种植区分布：19年玉米种植区分布：取交集后：呃呃。。。怎么重叠部分这么少呢~_~相信大家懂得都懂，这里就看破不说破……………………………………………………既然如此那好像前面都有些多余，那就单选一年的试一下吧。把2019年的
2022-05-10 21:44:36
6280
3
原创 python|批量修改文件名
参考：https://www.runoob.com/note/27030import ospath = "E:/网页下载/change/"fileList = os.listdir(path)for i in range(len(fileList)): oldname = path + os.sep + fileList[i] # os.sep添加系统分隔符 newname = path + os.sep + fileList[i][:8] + fileList[i][-4:]
2022-04-29 14:41:58
736
原创 RS|numpy.ndarray把-3.40282347e+38设为nodata
参考：https://www.pythonheidong.com/blog/article/486038/253098ead152c6bec4a7/问题描述：我的研究区只有左图中带颜色的部分，但是最开始用matplotlib着色的时候，那些nodata的全被当0着成了红色。把ndvi.tif（代码中的gray）的ndarray打印出来，发现nodata值全显示为-3.40282347e+38，但是gray[gray==-3.40282347e+38]=np.nan无效。解决方法：gray[(gr
2022-04-29 10:54:13
1893
原创 arcpy|批量投影
改下目录名称就可# -*- coding: utf-8 -*-# @Time: 2022/4/26 18:40import arcpy# 解决中文乱码import sysreload(sys)sys.setdefaultencoding('utf8')# 要投影的shp所在的目录arcpy.env.workspace = "E:/▲▲▲/2022-项目/220422-肃州玉米/肃州/"arcpy.env.overwriteOutput = True# 投影后的文件要存放的目录ou
2022-04-28 20:08:18
767
原创 Error|assert (os.path.isfile(mxd) or (mxd.lower() == “current“)), gp.getIDMessage(89004, “Invalid M
参考：https://www.cnblogs.com/youzi-xuchongyou/p/9241891.html错误：因为路径中有中文mymxd = "E:/myPy/in_data/china400/简称.mxd"解决方法：路径前加中文mymxd = u"E:/myPy/in_data/china400/简称.mxd"
2022-04-27 14:46:22
848
原创 python|格式化输出
%s 字符串%d 十进制整数%i 十进制整数%f 浮点数（默认小数点后6位）默认右对齐，加“-"变左对齐例如：x,y,z="john",20,1.85print("5Name:%8s Age:%4d Height:%5.2f"%(x,y,z))# 5Name: john Age: 20 Height: 1.85print("6Name:{:8s} Age:{:4d} Height:{:5.2f}".format(x, y, z))# 6Name:john Age: 20
2022-04-26 18:05:09
682
原创 python|读取txt文件的三种饭饭
import fileinputfile = "./in.txt"for line in fileinput.input(file): print(line,end="")fileinput.close()path = "./in.txt"file = open(path, "r", encoding="utf-8")file_text = file.readlines()for line in file_text: print(line, end="")file.clos
2022-04-23 18:36:33
913
原创 python|读取文件夹下所有文件名
import ospath = "D:/anaconda3/envs/lianxi/study/220415"file = os.listdir(path)for i in file: print(i)
2022-04-16 22:11:03
491
原创 DL|深度学习记录（计算机视觉部分）
一、数据增广%matplotlib inlineimport torchimport torchvisionfrom torch import nnfrom d2l import torch as d2ld2l.set_figsize()img = d2l.Image.open('/content/drive/MyDrive/DL/data/cat.jpg')d2l.plt.imshow(img)def apply(img, aug, num_rows=2, num_cols=4, s
2022-04-14 16:17:23
2935
原创 python|根据坐标生成shapefile
参考：利用百度地图坐标拾取系统获取经纬度，并将其坐标转换保存成SHP文件txt文件长这样（最后一行不是第一行也可以）：代码：import osfrom osgeo import osr, ogrdef coordinateTrans(txtPath): folder = str(os.path.split(txtPath)[:-1][0]) driver = ogr.GetDriverByName("ESRI Shapefile") out_shp_path = os.
2022-04-13 14:55:13
825
原创 python|生成游戏机超级玛丽风格音乐（上世纪的味道）
小姜说等她生日的时候给她写个生日快乐歌MIDI To RTTL:http://midi.mathewvp.com/midi2RTTL.php单音轨下载网址参考：https://www.bilibili.com/video/BV1LK4y1L7eR?share_source=copy_web%matplotlib inlinefrom pylab import *from IPython.display import Audiofrom scipy.signal import square..
2022-04-11 21:30:56
679
原创 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
1463
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人