小葱拌豆腐

原创 word页眉问题

设置页眉为章标题时不显示原因：回车键导致该页首个大纲级别为空名，删掉即可

原创 arcgis取反+掩膜提取

目的去掉全域中的A区域，获取非A区域0.获取A区域shparcgis 画图问题1.利用全区域shp-A区域shp=非A区域shp(交集取反)GIS如何获取两个图层叠加后去掉公共区域的部分2.利用非A区域裁剪(按掩膜提取)...

原创 python读取txt矩阵数据（转换为excel格式利用pandas读取）

目标：读取txt格式的矩阵数据数据描述：m*n矩阵，以空格分隔（开头存在空格，末尾存在\n）方法：复制到excel中，分列，利用pandas读取如何去除开头空格：利用excel以空格分列后，参考 ：excel快速删除空值单元格，数据上移，数据左移，之后利用pandas读取即可。也可以直接读，在python中处理开头空格问题，但逐行操作太慢了，参考Python3 中把txt数据文件读入到矩阵中的方法，其中报错：ValueError: could not convert string to f

原创 ENVI报错：Unable to recognize this file as a standard format.

ENVI 报错：无法识别文件Unable to recognize this file as a standard format.lf it is an external format with a known type,please specify the type.原始文件TIFF，经bandmath后保存为‘xx.tif’，bandmath生成文件默认为envi格式，因此无法识别，保存注意文件格式。...

原创 pycharm: connot find declaration to go to

背景：github 上下载项目后没有保存到Source root，出现connot find declaration to go to解决：将文件夹添加到interpreter path参考pycharm ctrl+点击无法跳转问题（connot find declaration to go to）

原创 arcgis鹰眼图问题

1.插入比例尺的图要有坐标系，不然比例尺有问题2.对原始图/鹰眼图插入比例尺时要激活对应的图层：图层-右键-激活鹰眼图制作流程：参考：【ArcGIS教程】如何用ArcGIS做一张专题地图

原创 arcgis 画图问题

问题：arcgis ：利用图片格式底图画线，出图1.自定义-工具条-绘图 调出绘图工具2.选择 ’线‘，在底图上画线3.绘制-将图形转化为要素（shp文件）（勾选转换后删除图形，删除没用的图形）4.出图...

原创 ENVI:There are no available ROls or EVFs associated with this input file.

利用ENVI subset Data from ROIs 工具时报错：There are no available ROls or EVFs associated with this input file.原因：原影像没有空间参考（基准）

原创 ENVI5.3 No Data像元值

在ENVI中被设置为data ignore value 的值将在影像中显示为No Data:不会参与之后的处理（波段运算、大气校正等等）可以通过编辑头文件的方式来更改No Data 值或者去掉No Data 值。多用于去掉黑色背景，如：...

原创 运用ENVI bandmath对多波段进行计算

原创 envi 打开影像报错：‘HISTOGRAM:illegal binsize or max/min‘.The result maybe invalid

利用ENVI打开影像时报错：'HISTOGRAM:illegal binsize or max/min'.The result maybe invalid原因在于影像中存在非正常值（极大值/极小值），使得影像不能正常显示。ARCGIS可以正常打开，存在一些非正常值。可以通过arcgis掩模裁剪去掉这些值。参考：喵家小抓的博客...

原创 python 生成可执行文件.exe

库PyInstaller只生成可执行文件：1.终端进入py文件目录2.pyinstaller -F xxx.py之后exe文件就生成在此目录下的dist文件中参考：Python 程序打包成 exe 可执行文件

原创 K-fold vs. Monte Carlo cross-validation（K折交叉验证与蒙特卡洛交叉验证（MCCV））

1.K折交叉验证单次k折交叉验证每一个数据顶多作为一次测试集（因为数据划分已经固定），对k-fold交叉验证运行的结果进行平均可以得到算法性能的(几乎)无偏估计，但方差很高；想要获得更多的切分方式：N次K折交叉验证（每次进行k折交叉验证时有不同的数据划分,共K*N个模型，取平均。实现：单次k折：sklearn.model_selection -kfold多次k折：sklearn.model_selection-repeatedkfold2.MCCV蒙特卡罗交叉验证(MCCV) =重复随机子

原创 R^2 score is not well-defined with less than two samples

使用sklearn.metrics中的‘r2’（决定系数R^2）时，报错：UndefinedMetricWarning: R^2 score is not well-defined with less than two samples.说明Y值只有一个（所以R^2不适用于留一法交叉验证）

原创 ARCGIS中生成KML文件

1.google地球批量导入坐标点2.解决：arcgis指定的表没有object-id字段3.图层输出比例：1

原创 Building Autoencoders in Keras

目的：利用autoencoder进行降维，并尽量保留主要信息（仅针对特征），mnist数据集数据处理读取数据、归一化（仅针对特征）import matplotlib.pyplot as pltimport numpy as npfrom keras import regularizersfrom keras.datasets import mnistfrom keras.layers import Dense, Inputfrom keras.models import Modelfrom

原创 利用BayesianOptimization库对模型进行贝叶斯调参(XGBOOST)

构建目标函数目标函数(参数优化最大化的函数)：-1*mse，使用10折交叉验证，取10次的平均MSE作为最终目标函数.def xgb_cv(max_depth, learning_rate, n_estimators, min_child_weight, subsample, colsample_bytree, reg_alpha, gamma): val = cross_val_score(estimator=XGBRegressor(max_depth=int(max_depth),

原创 MNIST手写数字辨识

多层感知器相关库import matplotlib.pyplot as pltimport numpy as npfrom keras.datasets import mnistfrom keras.layers import Densefrom keras.models import Sequentialfrom keras.utils import np_utils加载数据file = mnist.load_data()# tuple类型# 包含训练集：特征(60000,28,2

原创 ENVI去除NDVI/EVI异常值、NAN(not a number,无效值)、无穷大(inf)

NDVI/EVI范围：（-1，1）利用ENVI 去除异常值公式：b1为NDVI/EVI，小于-1赋予-1，大于1赋予1，其他正常(b1 lt -1)*(-1)+(b1 gt 1)*1+(b1 ge -1 and b1 le 1) *b1

原创 利用集成学习（Adaboost（DTs））对高光谱影像数据集（mat）分类

1.数据下载地址：Hyperspectral Remote Sensing Scenes包含：Indiana pines，Pavia University，KSC，Salinas的mat文件以Indiana pines为例：Indian_pines.mat、Indian_pines_gt.matBand:220pixel:145*145class:16类别：2.方法adaboost（基学习器：决策树）评价指标：OA(总体精度、混淆矩阵)3.主要函数3.1.相关的库：import

原创 matlab利用geotiffread读取tif文件报错：‘错误使用 tifflib， 无法打开 TIFF 文件’

在使用matlab读取TIF格式文件时：[a,R]=geotiffread('C:\2009.tif')报错：错误使用 tifflib，无法打开 TIFF 文件 'C:\2009.tif'原因：geotiff函数不能读取的情形：手动修改后缀，包括修改后缀为.tif，envi中band math等保存时修改为.tifgeotiff函数可以读取以下途径获得的TIF类型文件：1.envi中另存为tif的文件，即File-Save as-Save as (ENVI,TIFF......)，选择

原创 使用MODIS REPROJECTION TOOL(MRT)进行MODIS影像批处理（拼接+转投影）

1.准备数据MRT进行批处理时，所有路径必须为英文，要求数据格式为HDF，且命名规则符合官方命名规则（非常重要），如：‘MYD09A1.A2018001.h26v05.006.2018011145322.hdf’，如果将文件命名为其他，则无法识别出哪些影像可以拼接在一起，哪些影像不可以拼接在一起我的数据是某地区的6天的MYD09A1影像(反射率影像，分辨率500m，8天合成)，需要两幅拼在一起...

原创 模型评估与选择——参数、超参数、训练集、验证集、测试集

1.模型选择与评估1.1 模型参数、超参数了解什么是模型参数，什么是超参数是很必要的，不同类型的参数需要在不同的阶段进行调节（参数——训练集；超参数——验证集）模型参数：模型内部的配置变量，可以通过数据估计（学习）模型参数的值超参数：模型的外部配置变量，需要手动设置，好的超参数可以进一步提高模型的性能模型参数的特点：1.进行模型预测时需要模型参数；2.可定义模型功能；3.可使用数...

原创 XGBOOST模型简单调参（分类）

1.模型选择与评估1.1 模型参数、超参数了解什么是模型参数，什么是超参数是很必要的，不同类型的参数需要在不同的阶段进行调节（参数——训练集；超参数——验证集）模型参数：模型内部的配置变量，可以通过数据估计（学习）模型参数的值超参数：模型的外部配置变量，需要手动设置，好的超参数可以进一步提高模型的性能模型参数的特点：1.进行模型预测时需要模型参数；2.可定义模型功能；3.可使用数...

原创 OBS( Open Broadcaster Software )录制视频黑屏问题&录制方法

黑屏问题：设置-系统-显示-图形设置-浏览-obs.exe-选项-节能简单录制步骤0.添加场景1.添加来源方法一.窗口捕获（选择特定窗口录制）方法二.显示器捕获（全屏录制）2.音量调节3.开始/停止录制4.输出路径/格式设置5.查看已录制视频：文件-显示录像...

原创 pypi 镜像使用帮助——临时使用镜像网站

https://mirrors.tuna.tsinghua.edu.cn/help/pypi/

原创 IDL where函数与array_indices函数

WHERE简介按照指定条件查找数组元素，并返回数组中满足条件的元素的下标。语法：result=where(array_expression [,count] [,complement=variable],[,ncomplement=variable] )array_expression ：数组表达式count:符合条件的元素个数complement:不满足条件的元素下标ncomp...

原创 IDL size函数

简介数组基本操作函数：size，返回数据结构的基本信息（维数，大小，类型）语法：result=size(expression[,/n_dimensions|,/dimensions|,/type|,/tname|,/n_elements])其中：expression为标量、数组或者表达式；n_dimensions设定函数只返回维数；dimensions设定函数只返回每维大小；t...

原创 Pycharm下安装Tensorflow问题（利用anaconda）

目的：安装了Pycharm、Anaconda，并且已经配置好环境，想安装tensorflow包。以前安装tensorflow的时候直接pip install ，利用python自带IDE可以使用。这里的future warning 是由于numpy版本过高，两种解决方法，参考：链接Pycharm里也可以运行想在pycharm下的anaconda中安装tensorflow，做了以下...

原创 为numpy数组增加一个维度的方法

skearn中要求输入的特征数组至少二维，而有时特征数组不符合要求，需要增维如a=np.array([1,2,3])方法0.利用数组长度扩维b=a.reshape(((len(a),1)))#(3,1)1.直接扩维c=a.reshape(-1,1)#(3,1)-1的意义:根据另一个参数的维度计算出数组的另外一个shape属性值。[[1] [2] [3]]c=a.re...

原创 报错’TypeError: only integer scalar arrays can be converted to a scalar index‘

使用numpy.concatenate（）时报错：TypeError: only integer scalar arrays can be converted to a scalar index想要实现numpy数组按列拼接，使用numpy.concatenate函数函数说明我写成了：x=np.linspace(1,100,100)ones=np.ones(len(x))X=np.co...

原创 赫尔默特方差分量估计Python

基本思想先对各类观测值定初权，进行预平差，利用预平差得到的相应各类观测值的残差平方和，依据一定的原则对各类观测值的验前方差作出估计，重新定权。通过这样不断计算迭代，使得不同类观测值的单位权中误差趋向一致，从而达到最佳平差效果。赫尔默特方差分量估计的计算步骤（二元为例）实际应用有一组模拟的GPS和GLONASS观测数据，采用伪距单历元单点定位，根据经验，ＧＰＳ的伪距中误差为０.３ｍ，GLO...

原创 numpy.cumsum()函数

语法：numpy.cumsum(a, axis=None, dtype=None, out=None)按照所给定的轴参数（axis）返回元素的梯形累计和，axis=0，按照列累加。axis=1，按照行累加。axis不给定具体值，就把numpy数组当成一个一维数组举个简单例子（二维）import numpy as npa=np.arange(0,12).reshape(3,4)print...

原创 python中 return self的作用

举个例子class A: def __init__(self): self.c=0 def count(self): self.c+=1 print(self.c)a=A()a.count()#1当我们想对实例对象的方法进行连续调用时a.count().count()AttributeError: 'NoneType'...

转载 Adaboost python实现

AdaBoost算法详解与python实现

原创 利用sklearn实现adaboost，以单一分类树为例

0普通决策树0.1生成数据sklearn提供了随机数据生成的功能:make_regression 生成回归模型的数据make_multilabel_classification生成分类模型数据用make_blobs生成聚类模型数据用make_gaussian_quantiles生成分组多维正态分布的数据import matplotlib.pyplot as pltimport nu...

原创 KNN算法python实现（鸢尾花，iris）

读取数据集import operatorimport randomimport numpy as npimport sklearn.datasets as sdfrom sklearn.metrics import classification_report, confusion_matrixfrom sklearn.model_selection import train_test...

原创 稳健估计,P范数最小法

前言数据探测和稳健估计均可以探测出粗差，但是数据探测仅仅在数据中仅一个观测量含有粗差时好用，如果想要检验多个粗差，需要逐步剔除已检验粗差，再继续检验。而稳健估计在处理包含有多个粗差时有一定优势。稳健估计稳健估计分为三大类：M估计，R估计，L估计，M估计最常用，它是一种广义的极大似然估计，有分为选权迭代法和P范数最小法。M估计的准则其中ρ为M估计的函数，p为权M估计的稳健性与ρ的选择有...

原创 决策树的简单实现与可视化

0.原理略1.数据http://www.ai-xlab.com/my/course/135课时1-课程资料是CSV文件，依据已有数据推断是否买电脑数据属性：age，income，student，credit_rating标签属性：Yes/No共14笔数据2.需要的库csv，sklearn，graphviz其中graphviz需要独立安装并自行设置环境变量，否则报错:make...

原创 外文文献的检索、与下载

0检索文献网站：谷歌学术（用不了）谷粉学术（推荐）谷粉学术支持关键词、作者搜索：如deep learning可以显示时间、引用次数、类型及名称（book or journal）1下载文献选择一篇文献，点击进入期刊网站，下载方法有两种：DOI和URL1.DOI找到它的DOI复制到sci-hub然后进入下载界面：点击save有的时候需要输入验证码下载完毕2.URL...