Rachel0417的博客

原创 matlab逐像元计算栅格数据10年间的变化趋势代码

python的linregress和matlab的polyfit，都是用的最小二乘法，直接简单粗暴地将[1,2,3,4,…N]的时间序列与像元值进行的回归拟合。和那种高级的滑动时间序列预测啥的还是不一样的。

原创 GEE多行注释疑似失灵Ctrl+Shift+/

GEE多行注释疑似失灵Ctrl+Shift+/

原创 PermissionError: [Errno 13] Permission denied:报错解决

PermissionError: [Errno 13] Permission denied:解决报错

原创 栅格数据有属性表，但表是空的

如图，DEM数据可以正常显示，但由于属性表没有内容，所以没办法进行转矢量等一系列涉及到值的操作直接构建属性表就🆗

原创 arcgis去除遥感影像黑边

采用“复制栅格”工具，在下图开始之前，还尝试过1 不填忽略背景值 参数，不修改nodata值（默认为3开头的一串数字）2 设置nodata值为0，不填忽略背景值 参数以上两步均失败，黑边依然存在3 继续尝试，单个波段单个地进行，设置如下结束后进行波段合成。去除失败4 先是一整张图像进行下图的设置，试了两次，软件都崩了。然后单个波段单个地进行，参数改为去除成功。单个合成完后，波段合成即可因为我是一步一步做的，就是每一次尝试都在之前的结果上进行，所以我不知道是因为我在前面的步骤中修改了

原创 ENVI ERROR:Radiometriccalibration failed: calibration processing failed.

不是中文路径的问题，作为一个ENVI小白也不知道怎么查这种问题重启无济于事，于是怀疑是软件本身的补丁出了问题，就开始卸载软件结果卸不掉？？出现下图情况准备把QQ登上进行OCR然后百度时，QQ先提醒我清理个人空间，我一看E盘好家伙，只有8M了，程序运行得了个啥？？？太笨了太笨了，以后一定多多关注可怜E盘的情况因为其本身就没有多大，然后之前又已经操作了两轮，生成了新的大文件，给本就存储困难的空间雪上加霜这个助手让我感觉自己受到了猫饭...

原创 arcpy中波段合成与波段提取的坑

波段合成这是个大坑，目前我没有找到好的解决办法它是什么问题呢？就是arcgis中的波段合成，生成的是临时图层，需要手动导出虽然帮助中给了波段合成的代码，但至少对于我来说，是根本无法实现！！！生成的临时图层根本无法转成栅格文件进行存储总之就是不要使用arcpy进行批量合成波段import arcpyfrom arcpy import envfrom arcpy.sa import *import osarcpy.CheckOutExtension("spatial")arcpy.gp.ov

原创 生态遥感指数的mk分析的准备工作

在获得逐年逐月的裁剪后的rsei图层后，首先需要将需要分析的波段进行波段合成（这里不得不吐槽我的正版arcgis10.8，无论做几个波段的波段合成，结果出来后在arcgis文件夹里的图层都是band1，而且是一个东西，值全是一样的，但使用信息查看工具看到的数值又是正确的，根本就无法再提取原始波段；但在另一台电脑上就是正常显示，这就是正版的妖力吗？？？？）安装mk分析python包，网址为https://www.cnpython.com/pypi/pymannkendall，pip即可接下来开始做分析的

原创 多波段图像进行主成分分析时，遇到栅格数量错误情况的解决方法

在老师的帮助下做遥感生态指数（RSEI）的主成分分析，有一个月的图层死活有问题，报如下错误：因为开始PCA之前需要抠水体，裁剪的过程中，就没有掩膜栅格，不知道是哪一步出了问题，重来了无数次，用了各种各样的掩膜去裁剪或提取都无济于事，只好又求助老师老师说，是因为合成前某个原始图层的nodata值有问题，所以合成后的合成图像没有Nodata值就是说每个栅格图像的nodata值并不是就是真的没有值，而是系统自己生成的一个很小很小的值但这个出问题的图层的nodata值有问题，所以转换图层的时候就会出错

原创 关于掩膜与捕捉栅格的重要性

当你的栅格不论是用矢量裁剪、还是栅格裁剪、或者按掩膜提取处理后的数据，都与其他数据产生偏差的时候一定记得看看环境中是否设置处理范围、掩膜与捕捉栅格这三个选项我已经吃了不下20次亏了…时间全耗在这些上面...

原创 arcgis网络分析——设施服务水平分析笔记

原创 arcgis10.6 模型构建器笔记

原创 arcgis pro导入三维模型踩过的坑

如何通过arcgis pro导入各类格式的三维模型，并发布至arcgisonline平台的方法，官方在说明文档中已经写得比较清楚了，这里记录一下我当时总是导入失败的原因我自己的情况中秋的时候企图发布三个模型，从晚上九点折腾到了凌晨五点终于发上去了/(ㄒAㄒ)/~~现在想起还是应该来记录一下不会使用sketch up软件，对软件中的坐标系等基础设置一窍不通，因此可能导致了后续一系列问题第一个模型导入失败是因为操作不当，在su中使用obj/3dmax导出后，再用“导入3D文件”导入pro，此时视图中仍然

原创 执行arcgis缓冲区时可能遇到的问题

1 最开始的要素类没有定义投影，然后报错：ERROR 999999: 执行函数时出错。输入数据的范围超出了经度 ±360 或纬度 ±90执行(Buffer)失败。就是没有投影的结果，最好不要用地理坐标系做2 ArcGIS做多环缓冲区报错：超出输出几何属性范围出现 坐标或测量值超出范围 的报错时，当要素的坐标超出要素类的属性域范围时会发生此错误。这种情况会在创建新要素或编辑现有要素的坐标时发生。坐标必须位于要素类 x,y 属性域以及 z 或 m 属性域（如果要素类可存储这些值）的范围内。属性域在

原创 C语言防止非法字符输入——ret的用法

PPT上猜数字的防止非法字符输入的写法：ret = scanf("%d", &guess);while (ret != 1){ while (getchar() != '\n'); printf("please guess:"); ret = scanf("%d", &guess);}counter++;1 ret的用法百度某回答说：“scanf函数成功输入的数值个数 我们可以通过判断if(scanf("%d",&man)==1)判断是否

原创 cmd与vscode终端的中文不一致的问题

今天运行C语言的作业，vscode的终端里的中文始终是乱码，但是cmd里的中文又是正常的简体中文。到网上查了资料明白两者应该保持统一才行，但无论是把两者都改为utf-8还是gbk都不能改变两个框框显示不一样的事实。尝试了直接使用旧版命令提示符，总算是一样了。修改方法如下图：（现在是旧版cmd）在这里也能看到cmd的编码方式。网页上的chcp 56001只成功了一次，不清楚原因TT勾选过后重启vscode就可以啦...

原创 打印后加换行符的网页参数报错：TypeError: can‘t concat str to bytes

背景：获得网页post请求结果的各项参数一开始不知道这些参数是干什么的。。想打印的时候换行，就挨着加了换行符试了一下post_response=requests.post(url,headers=headers,data=None,json=None)#打印出来是显示响应状态码，通过.text 方法可以返回是unicode 型的数据，一般是在网页的header中定义的编码形式print(post_response)print(get_response.text+'\n')print

原创 python模块中urllib和requests发送请求的区别

此次实践中获取网页使用的代码urllib:importurllib.requesturl=['http://wz.lanzh.95306.cn/mainPageNoticeList.do?method=init&id=2000001&cur=1']#通常我们爬取网页，在构造http请求的时候，都需要加上一些额外信息，什么Useragent，cookie等之类的信息，或者添加代理服务器。往往这些都是一些必要的反爬机制headers={ 'User-Agent.

原创 python监控网站学习笔记2——网页中文乱码

一 查询网页编码方式右键——检查——console——输入：document.charset——回车二 requests的编码方式爬取的所有网页无论何种编码格式，都转化为utf-8格式进行存储，与源代码编码格式不同所以出现乱码.unicode是一种二进制编码。1.编码方式为“utf-8”的网页也可能会乱码，那么可以也试试下面的一些函数（对不住，没学过计算机网络，只能先试出来再说了TAT2.一些乱码解决方式（1）教务网:编码方式为GBKa.xxx.encoding = "utf-8

原创 python监控网站学习笔记1——web基础

一 HTML基础知识二 一些名词抓包分析：对截取到的数据进行分析解析json数据：json数据模拟登录IP代理池

原创 ArcScene三维地图导出失败

在ArcScene里用图层的属性进行一顿拉伸，全都成了钟乳石，看着还挺稀奇，见多识广的老师就只会说丑/(ㄒoㄒ)/~~地形图还好，我其实也没搞懂自己怎么就去拉斜率的图了，被说丑是活该然后就疯狂找能把.wrl转成.obj、.stl之类的能导进PPT里的格式，全都失败了！！！这些都不得行！！！！！找到了一个介绍文件类型的网站，...

原创 QGIS初体验

一、设置中文就找了好久！！！也太蠢了吧！！！二、数据怎么导入。。看到了一个灯光图教程，不错，就用它做一个NDVI植被之光的图好了。（瞎78乱编的名字只能导入栅格数据集？但是我不认识一个栅格数据集是由哪些文件组成的（这只是暂时的！！理论知识量与期末倒计时成反比）好的，看不懂，打开栅格数据就需要代码来实现了吗？对不起，是我太天真了TAT这就滚去看使用手册...

原创 18年NDVI斜率图生成

想自己用python实现多幅栅格图像多年变化率图的制作，但网上始终没有找到合适的代码。基础太差了55555贴老师写的多波段图像实现多年降水与气温的月NDVI相关系数图代码：import numpy as npfrom pyrsgis import rasterfrom scipy.stats.stats import pearsonrRasterFile1 = 'D:/2ndcmpttions/遥感小分队/相关系数图/NDVI波段/NDVIAPR.tif'RasterFile2 =

原创 2020-12-17

C还有两个小时就看完好老师的C了。从7.16到现在，一个计算机入门课程上了半年也是让人很惭愧。目录好老师的课主要包含了7个方面：基础知识如何学习运算符数据类型的介绍基本的输入输出函数强制类型转化进制与补码一点ASC||码流程控制选择循环数组数组的意义数组的分类函数函数的定义定义函数的方法注意的问题常用的系统函数递归指针结构体链表前期听的还是比较仔细，但是后期重点的时候反而因为时间不够变得马虎起来，好在百分之七十还是听懂了；接下来希望在阅读大话数据结构

转载 GIt初体验（假的）

写在前面第一次写博客。从没想到自己会想来干这样一个曾经以为多此一举到后来觉得我不配的事情，但对越来越杂乱的生活来说，每天对自己进行一定的总结变得必不可少。常常会在某个瞬间反思起自己在过去的一个小时、一天、一周做了什么？我得到或失去了什么？转眼间就要迎来年底——十二月，在这特殊的一年，我有什么收获？虚度了多少时间？不能再等了，现在就要开始思考自己为每一件事情耗费的时间是否值得，是否可以花费更少；现在是2020年的十一月底，这也许就是我的进步吧。没有怎么用过博客，不知在这个大佬云集的平台是否接受一个计算机课