遥感
文章平均质量分 92
陈修一
加里敦大学地理系毕业生,编程练习时长三个月的个人练习生。
主要研究方向:人口地理、城市地理、大数据、遥感。
展开
-
【Earth Engine】基于GEE协同Sentinel-1/2使用随机森林回归实现高分辨率相对财富(贫困)制图
最近在做一些课题,需要使用Sentinel-1/2进行机器学习制图。然后想着总结一下相关数据和方法,就花半小时写了个代码。然后再花半小时写下这篇博客记录一下。这篇博客主要参考数据(相对贫困指数,RWI)来自这个GEE社区网站,有缺数据的可以直接在这上面找,要用的时候调用一下就行了。这是这个数据的一个交互地图预览。工作完全是在GEE平台上写的,如上面所说,这个工作跟我的课题内容关系不大,纯粹是拍脑袋的需求然后拍脑袋写的代码。原创 2023-12-21 21:54:20 · 1507 阅读 · 0 评论 -
【Python】使用Python调用Fragstats批量进行万级及以上数据的景观指数运算
本文针对的场景是,已经通过arcpy制作了渔网(fishnet)并且把一个tif根据渔网格网分裂成成千上万个带编号的小方格tif,接下来就要计算这些带编号的小放个的各种景观指数。本文假设同学们已经会使用界面化的fragstats,并且会制作批量计算的fbt脚本文件(也就是import batch需要导入的文件)。本文包括两部分:1,拆分批量运算的脚本文件(默认为fbt文件,但本文中为了方便使用csv格式的脚本);2,根据脚本文件集批量计算景观指数。原创 2022-09-23 20:00:36 · 4357 阅读 · 13 评论 -
【Earth Engine】基于GEE合成Landsat4/5/7/8/9影像并进行NDVI、NDWI和MNDWI等指数计算
调用符合要求的所有Landsat系列影像合成,并计算各种基础指数。原创 2022-06-13 17:59:36 · 10021 阅读 · 23 评论 -
【Python】批量爬取OSM建筑瓦片数据并整合为shp格式数据
目录1 简介2 效果展示3 思路及代码3 完整代码1 简介最近在做一项课题,涉及到建筑足迹(Building footprints)数据。所以想看看现在比较常用的都有什么产品。目前查到的比较权威的有微软的和谷歌的,虽然以上两个数据的覆盖率都挺全的,但是无奈几乎都只有图形信息,没有其他字段。所以我还找到了OSM建筑数据作为补充,OSM数据虽然遗漏的建筑比较多,但有比较丰富的字段信息,如下。比较麻烦的是,OSM建筑数据不是一键下载的,而是geojson的瓦片数据。但还好它提供了url。所以,这篇文原创 2022-05-16 11:09:02 · 3047 阅读 · 4 评论 -
【Python】一个简陋的基于混淆矩阵计算遥感分类精度(总体精度、Kappa系数、用户精度、生产者精度等)的代码
目录1 简介2 案例和代码说明3 完整代码1 简介有一篇文章返修了,由于文章的重点不在分类所以我之前就随便写个Kappa系数上去,没想到审稿人居然要求我提供其他参数ಥ_ಥ可是我只大概存了个各类型的分类。。。后来整理整理也只能得到一个混淆矩阵。。。说实话分类我做得也不多,作为纯ArcGIS党以前只知道有软件可以跑各种精度或GEE用现成代码跑。。本着靠自己不靠别人的想法,这次我就想写一个基于混淆矩阵计算总体精度(Overall Accuracy,OA)、Kappa系数(Kappa Coefficient)、原创 2022-02-28 13:55:39 · 8477 阅读 · 3 评论 -
【Python】基于Python统计文件目录内每景栅格图像均值并输出成表
目录简介效果展示代码简介乔博士最近在做一个研究,涉及到过去十几年某些地区某些指标逐日的提取。一个指标5000+副影像需要获取均值,一共四五个指标。。。然而ArcGIS上手动做还是ArcPy做都不太现实,所以就问问我有没有相关的代码程序:我是有的,然而因为我之前涉及的没有那么多影像,最多也就几十景求均值,所以现成的代码是输入在终端里的。所以乔博士提出了这么个要求:所以,这次的代码目标是提取某个文件夹里所有影像的均值,并且输出成excel可以读的表。效果展示首先,得把tif扔进一个文件夹里:原创 2021-07-14 16:04:50 · 2354 阅读 · 5 评论 -
【Earth Engine】基于GEE进行非监督学习
目录简介与摘要采样点的选取合成多波段影像分类器参数设置影像的显示与结果后记简介与摘要最近导师给了个新方向,其中要用到长时间和大空间的高分辨率遥感影像分类数据,其中有一类还是季节性池塘。。这。。这么多要求。。实在没有免费的并且符合要求的数据。。于是我寻思还是自己做吧,所以马上学习了下GEE。。。在第一阶段,我使用landsat合成全年的影像,比如landsat5,只用到1-7波段,输入分类器的波段实在太少了,分类精度很低,而且根本提不出来季节性的地物。第二阶段,我终于意识到用全年合成影像是不行的,原创 2021-04-14 21:26:37 · 1844 阅读 · 3 评论 -
【Earth Engine】基于GEE对季节性地物进行分类(多源数据叠图+监督分类)
目录简介与技术流程数据及其预处理要调用的包代码函数:计算斜率和变异系数代码:计算栅格图像的斜率代码:计算栅格图像的变异系数代码:函数的调用与结果后记简介与技术流程最近导师给了个新方向,其中要用到长时间和大空间的高分辨率遥感影像分类数据。这。。这么多要求。。实在没有免费的并且符合要求的数据。。于是我寻思还是自己做吧,所以马上学习了下GEE。。。关于GEE的申请什么的我这里就不多说了,这篇博客直接进入应用。总的来说,就是把多幅影像叠加在一起,逐像元构建一组时间序列,然后计算它的变异系数。比方说,原创 2021-03-16 20:47:02 · 11831 阅读 · 7 评论 -
【Envi】基于单窗算法的地表温度反演实验操作记录
前天刚下的Envi,高强度学习了两天摸索出了点门道,不一定准确。本文谨作为个人记录,如有问题欢迎指正。写这篇博客的时候,我是边做边写的,想到啥写啥,可能会比较混乱一点。后面有空的话我会整理一下。参考博客在实践操作这方面,相关博客给的启发更大一点;而文献大部分只参考了公式或者是经验数值,太多了太杂了就不一一写出来了。记 / ENVI5.3温度反演-单窗算法-操作要点地表温度热红外遥感反演-以单窗算法为例地表温度反演操作总结技术流程1.数据预处理1.1数据获取首先在地理空间数据云下载lan原创 2020-08-10 15:14:37 · 35738 阅读 · 80 评论 -
【Python】基于Python计算长时间遥感栅格图像的像元值变化度(斜率)和变异系数
目录简介与技术流程数据及其预处理要调用的包代码函数:计算斜率和变异系数代码:计算栅格图像的斜率代码:计算栅格图像的变异系数代码函数的调用后记简介与技术流程之前看到一篇文章,用变异系数(CV)计算年际初级净生产力(NPP)的变异系数,研究一个地区的产量稳定性。下面是这篇文章(Hong C, Jin X, Ren J, et al. 2019. Satellite data indicates multidimensional variation of agricultural production in原创 2021-01-09 16:05:38 · 9595 阅读 · 37 评论