- 博客(29)
- 收藏
- 关注
RSS订阅原创 2月科研——arcgis计算植被差异
ArcGIS中，设置高于或低于某个值的像元为 -9999，然后将这些地方设为空——目的：去除异常值和黑色背景值。2、得到格网后，统计0.5°网格内500m像元的平均值，分区统计（zonal statistics）3、zonal statistics：分别计算均值、sum，再用合成波段。以此类推，再累加得到7-10月的总差异，分析其与SPEI的一致性。2、栅格计算器，Con，有值的部分全部赋值为1，以便在网格中计数。1、栅格计算器，将500m有值的像元统一赋值为1。2、统计格网中500m像元的数量。
2024-07-19 17:13:21
928
原创 GEE计算百分位数、报错、重命名波段、SIF基础
GEE计算数据的第XX百分位数的数值：第一段代码为具体某个百分位数，第二段是多个等差的百分位数。这个问题时主要出现在对ImageCollection进行运算时，其中的一个image的。把 Export 中的region参数修改为：roi。
2024-05-03 10:17:49
364
4
原创 arcgis操作：copy raster、高清导出图、渲染色带保存
打开新数据属性-classified，右上角import刚保存的lyr文件，就会让新数据的分级及颜色与原有数据的一致。以及下面这个分级save和load，也是很好用的功能。：设为tiff、LZW格式，doi在300以上。这种方法会比同样doi、png格式导出的大小低一半。新方法：copy raster，这是个好工具，可以把nodata值设为空或其他值。2）保存色带：调整好后点击确定，选择save as layer file。
2024-05-03 10:16:30
631
原创 GEE learning：图像波段合成、LC、地类面积统计、内存超限
填补的时候，遇到 Composites_2022 的分辨率是111319.49079327357m，而Collection2022和 lswi_interp_2022的分辨率都是500m的情况，实在想不明白为什么会这样。monthly_2022每张image的 system:time_start属性都被设置为一样的了，怎么改？MYD09A1和MOD09A1的区别：MOD和MYD的区别是卫星不同，前者搭载Terra（上午过境），后者搭载Aqua（下午过境）。A1和Q1的主要区别是波段和空间分辨率不同。
2024-02-20 20:58:00
1011
原创 GEE绘制直方图、提取像元值和经纬度及个数
SIF 的地类聚合位置放错，导致 .gt(0.5) 无效；而且2022年的cropland并没有 .gt(0.5)
2024-01-29 22:25:47
661
原创 中国1千米分辨率逐日全天候地表土壤水分数据集（2003-2022）HDF5数据处理
0-10cm处的土壤水分，比例为0.001，单位是cm3/cm3。和m3/m3单位的换算应该是1：1【？不确定 全网居然没资料……如果有小伙伴能回答就好了】；逐日的SM加和得到每月的【同样不确定是加和还是求平均。个人感觉和降水应该差不多 是加和】3、栅格计算器设-9999为 nodata，然后逐日加和得到对应月份的土壤水分。arcgis- 模型构建器。arcgis- 模型构建器。记录一下这个数据集的处理过程~
2024-01-29 22:25:20
1120
14
原创 GEE learning: 重投影reproject、重采样resample、聚合aggregate的区别
例如，如果有30m和500m的数据，需要统一尺度，就需要对30m的影像进行下采样。遥感影像产品使用的一种全球正弦投影，属于投影坐标系。需要转换为墨卡托投影（了默认的最邻近法重采样（可能出现网格，插值效果不好）和三次卷积方法的效果。GEE重投影时默认的是最近邻法，可以用。：如果设置scale 默认最邻近法重采样。还可以用其他重采样方法。当目标不是在重投影时进行。），从而消除上述像元倾斜的状态。，而是将小尺度像元采用不同的投影。重采样后需要重投影。为一个更大尺度的像元，可结合。
2023-11-20 17:18:03
1939
1
原创 R——《数据科学中的R语言》学习【基础篇】
向量是R中最基础的数据结构，其衍生出了矩阵、列表（小火车）、数据框（excel表，是特殊的列表，其中向量等长。广泛应用于数据分析）。3、c() 函数构造向量：向量是有序结构，四种类型（数值型、字符串型、逻辑型、因子型）。返回FALSE FALSE TRUE FALSE TRUE FALSE。2 == 0，判断是否为偶数（返回TRUE / FALSE）如 is.na(c(1, 2, NA, 4, NA, 6))2、查看对象的类型和长度的命令。判断元素是否属于向量。判断元素是否为缺失值。
2023-10-12 21:15:58
111
原创 GEE learning - gif动图、ImageCollection操作、SG滤波
multiyear的SPEI值：https://code.earthengine.google.com/700d5e7bdd87b9ab934be1312c34d892。2022年的SPEI值：https://code.earthengine.google.com/85ac123db3eb8448f389908469746a4d。）：窗口大小是用于局部拟合的数据点的数量。窗口大小越大，平滑效果越明显，但过大的窗口可能会导致信号失真。：从很常用的SPEI网站（时间分辨率1月、3月……
2023-10-11 22:16:36
425
原创 GEE Training：qualityMosaic和max的区别 | iterate迭代（3）
qualityMosaic和max的区别 | iterate迭代的基本框架和代码
2023-08-30 14:51:16
340
1
原创 GEE Training——Guideline篇（1）
GEE指南解读&需要注意的点（image、imagecollection overview）
2023-08-17 17:18:47
68
1
原创 坐标系与常见投影
例如：“WGS_1984_UTM_Zone_50N”这个坐标为投影坐标系，其中“WGS_1984_”为地理坐标系，“UTM”为投影类型。常用于我国各省市投影）、Lambers兰伯特投影——“实际处理中，只有地理坐标一致的数据，投影成一个坐标系后才能较好匹配。：北京54、西安80、WGS1984、CGCS2000——所以不管三七二十一，首先设定地理坐标成一致 然后再投影。平面坐标系，用xy表示位置，单位通常为米。球形坐标系，由经纬度表示位置。（我国遥感影像常用）、Albert投影坐标系。（常对应WGS84）
2023-07-31 10:25:42
757
1
原创 Zotero操作技巧
—导出txt并复制，在zotero新建一个文件夹（方便区分），点击菜单栏文件。2、条目面板：查看-布局试图。1、批量从WOS导入文献↑。
2023-07-31 10:24:54
271
1
原创 毕设思路和过程梳理【4月27日及以后】
对高温热害造成结实率下降的中稻，要抓住有限的灌浆时间，科学增施穗粒肥，喷施叶面肥，增强光合作用和灌浆速率，努力增加粒重，弥补结实率下降损失；安徽、湖北、江西、湖南、福建等地因旱情持续或发展，导致无水源保障地区的晚稻长势弱、灌浆进度缓慢，影响棉花裂铃吐絮、柑橘果实膨大及果实着色。然后想要试着剔除异常值，代码写好了发现分辨率500m的数据有问题：1000m分辨率数据导出是没问题的（像元的值呈正态分布），但是500m分辨率的数据有80%的像元的值都分布在同一个狭小的区间内，如（-0.1， 0.1）。
2023-05-16 23:36:14
253
原创 MCD43A4_deprecated / 作物底图的古早代码
MCD43A4 天顶BRDF调整的反射率(NBAR) 产品 (MODIS/Terra Nadir BRDF-Adjusted Reflectance Daily L3 Global 500m SIN Grid) 在当地太阳正午时为每个 MODIS 光谱带 (1-7) 计算. 由于视角效应将从定向反射率中移除，这将产生更稳定和一致的产品。NBAR 值直接用作生产全球 MODIS 土地覆盖产品 MCD12Q1的先进技术分类器的主要输入。网友有的说可用，有的说不能用了。
2023-04-20 16:57:53
312
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人