我是菜鸡，我不敢睡-优快云博客

原创 Google Earth engine 机器学习算法分类回归算法实现

回归算法，这里以随机森林算法为例：以土地利用分类为例从var tem=ee.ImageCollection("OREGONSTATE/PRISM/AN81d") .filterDate('2016-08-01','2016-09-01') .first().select('tmax')var elev=ee.Image("CGIAR/SRTM90_V4")var lucc=ee.ImageCollection("MODIS/006/MCD12Q1").fil...

2021-06-05 15:26:39 2405 1

原创 GEE Expected a homogeneous image collection, but an image with an incompatible band was encountere

这个问题时主要出现在对ImageCollection进行运算时，其中的一个imaged的数据类型或值域范围不同导致。我们首先看错误是否为数据类型不同，即，int,float,double.如果是。只需要image.toInt(),toDouble(),toFloat()即可。而如果是值域范围不同，这里我们以NDVI 为例，其正常范围是-1到1，但是当存在异常值时，数据的值域范围就不是这样的了，这时候我们采用最简单的方法，将大于1和小于-1的分别赋值为1，-1.并且设置数据的类型和值域范围。代码如下...

2021-05-27 09:12:40 2125 2

原创 Google Earth Engine(GEE) 多元线性回归拟合 ee.Reducer.linearRegression

在GEE中想要进行多元线性拟合可以使用ee.Reducer.linearRegression方法。这里举个例子，建立温度，降雨，NDVI 与EVI之间的关系。数据时间范围2000-2020。这里我们需要将数据处理为以一个包含21的image的imagecollection，每个image包含四个波段，如上面顺序。ee.Reducer.linearRegression采用最小二乘方法拟合参数，以最小化残差为目标结果返回两个波段，一个是系数，一个是残差（均方根误差）具体实现如下：首先将数据处理为im

2021-05-08 09:12:03 6226 8

原创 Python ESA 土壤水分数据nc转为tif格式

ESA 土壤水分下载网址：http://www.esa-landcover-cci.org/?q=node/164数据处理结果是每天一个数据，结果如下在arcgis中的信息如下：处理代码如下：import numpy as npimport netCDF4 as ncfrom osgeo import gdal,osr,ogrimport osimport globdef nc2tif(path,name,save_path): os.chdir(pa...

2021-05-03 19:51:22 1686 3

原创 Google earth engine 获取栅格数据的唯一值，unique values

caiyong reduceRegion方法，使用reducer:ee.Reducer.frequencyHistogram().来统计不同值的分布情况，然后将结果转化为字典形式，获取其keys().运行结果如下：var da1=LUCCS.filterDate('2016-01-01','2016-12-31').select('LC_Type1').first()var re=get_unique_value(da1,geometry,1000)print(re)var lucc_v

2021-04-26 16:47:02 1239

原创 arcpy 批量采样点数据转为点shp文件

在实际科研生活中，通常我们需要吧我们的采样数据导入到arcgis，当采样点只有一批时，我们可以通过使用arcgis软件，添加采样数据（csv格式）--display xy--导出为点shp文件。但是当我们想要批量将采样数据转化为点shp时，上述方法过于笨重。下面我们介绍通过批处理的方法，批量将点数据转为shp点文件。主要思想，我们首先建立一个空的点shp文件--添加字段--使用游标逐行插入的方法实现单个采样点到点shp文件的转化。表格主要内容如下：表格如以下形式。代码如下:我

2021-04-22 21:02:20 1576 3

原创基于Python 利用桑基图（Sankeydiagram）表示转移矩阵（土地利用）

桑基图桑，基能量分流图，也叫桑基能量平衡图。它是一种特定类型的流程图，图中延伸的分支的宽度对应数据流量的大小，通常应用于能源、材料成分、金融等数据的可视化分析。因1898年Matthew Henry Phineas Riall Sankey绘制的“蒸汽机的能源效率图”而闻名，此后便以其名字命名为“桑基图”。看到这么炫酷的图，是不是心动了呢对于土地利用文献中大多使用转移矩阵的方式。本博客实现2000年，2005年，2010年，2015年，2019年土地利用之间的转移情况使用桑基图可视化。首选.

2020-11-17 19:59:21 7550 10

原创 arcgis DEM 提取流域范围（详细教程）(空间分析--Hydrology)

arcgis DEM 提取流域范围（详细教程）(空间分析–Hydrology)1：洼地填充 Fill（生成无洼地的DEM）2：计算流向 Flow Direction(基于无洼地的DEM 计算流向)3：计算流量 Flow Accumulation (基于无洼地的DEM 计算流量)4：基于Con函数根据径流量大小提取河网结构（栅格），确定径流量的阈值，将大于该阈值的区域赋值为1，这里径流量的大小是指的相元个数。4：基于提取的河网结构（栅格），生成河网链接，5：对河网进行分级。6：将生成的带

2020-11-09 16:50:20 30978 4

原创迁移学习案例基于VGG16的迁移学习

由于深度学习网络的训练复杂，耗时，并且对硬件要求较高，因此对于目标相似的任务，并且自己的数据集数目较小的时候，比如分类任务，不同的分类任务之间的对象具有较为相似的特征，我们就可以使用已经预训练的模型来应用到自己的项目任务中，当然对于已经预训练的模型，需要经过微调才能应用到自己的分类任务中来。微调主要包括训练参数的变动，全连接层神经元个数，当我们自己的数据集很少时，我们可以仅仅改变最后一层全连接层的神经元的个数即可，当我们的数据量大的时候，我们可以选择训练模型的后几层参数。这里我们采用VGG16来实现简单

2020-10-18 23:12:46 2331

原创利用卷积神经网络对CIFAR-10数据集分类

CIFAR-10　　该数据集共有60000张彩色图像，这些图像是32*32，分为10个类，每类6000张图。这里面有50000张用于训练，构成了5个训练批，每一批10000张图；另外10000用于测试，单独构成一批。测试批的数据里，取自10类中的每一类，每一类随机取1000张。抽剩下的就随机排列组成了训练批。注意一个训练批中的各类图像并不一定数量相同，总的来看训练批，每一类都有5000张图。python 版本数据集下载网址：）http://www.cs.toronto.edu/~kriz/cifar-

2020-10-17 17:50:53 2442 2

原创 deep dream 案例

什么是DeepDream？DeepDream是谷歌发布的对卷积神经网络（CNN）进行可视化的方法，旨在更好的理解深度学习的过程。其通过某一网络（网络结构固定）根据提供的图片，来调整输入的图片，以生成目标的图片，其原理图如下：其中在生成图片的过程中，网络的结构是固定的，是预训练好的，其改变的是输入数据（下图的随机噪声）以生成目标图片（下图的海星）其中具体原理大家可自行搜索学习。本案例是慕课上的深度学习案例，使用GoogleNet的Inception预训练模型来生成目标图片，预训练模型的下载地址： h

2020-10-17 00:58:23 541 1

原创 module ‘scipy.misc‘ has no attribute ‘toimage‘ module ‘scipy.misc‘ has no attribute ‘imresize

1:module ‘scipy.misc’ has no attribute ‘toimage’2:module ‘scipy.misc’ has no attribute ‘imresize第一种解决方法：卸载当前的scipy, 重装较低版本scipy==1.2.1第二种解决方法（使用PIL，cv2等库进行替换）rezise函数：cv2.resize(img,(h,w)))toimage函数：cv2.imwrite(’./inception/’+img_name,img_array)...

2020-10-17 00:39:08 228

转载 arcgis提示ImageServer can not be started

在运行arcgis水文模块fill等工具时，出现错误提示ImageServer can not be started。解决方法参考https://blog.youkuaiyun.com/wxqqq666/article/details/91359769进入环境设置将parallel processing 设置为0.就可以运行了，但是感觉运行速度变慢了...

2020-09-02 13:09:12 1100

转载 Geopandas 安装

Geopandas 依赖库从 http://www.lfd.uci.edu/~gohlke/pythonlibs 下载 Fiona , GDAl , pyproj , Shapely其他包可以通过 pip install geopandas 安装参考：https://blog.youkuaiyun.com/SmallTankPy/article/details/71194862

2020-08-24 11:17:46 542

原创 arcgis 根据条件批量给另一个字段赋值

arcgis 根据条件批量给以字段赋不同的值

2020-04-23 22:39:22 23676 3

原创使用Pyqt5制作GUI界面，并使用pyinstaller打包成exe文件

基于pyqt5制作GUI界面，实现简单的图片的打开，缩小和放大功能，使用pyinstaller打包成exe文件

2020-04-13 10:06:37 5175 5

原创 python 爬虫小案例８６８４网站爬取北京公交路线站点信息。

基于Ｐｙｔｈｏｎ　在８６８４网站爬取北京公交路线的信息，包括站点，站点的经纬度，路线的描述，票价等信息。

2020-04-08 10:25:57 5433 4

原创基于百度API使用经纬度反查位置信息，百度API密钥申请

基于百度API使用经纬度反查位置信息，包括百度API密钥的申请

2020-04-08 04:41:04 2071 1

原创 PYTHON API获取某地点的经纬度

查询网址http://api.map.baidu.com/lbsapi/getpoint/index.htmlimport jsonimport requestsfrom urllib.request import urlopen, quoteaddress = '焦作市云台山'url= 'http://api.map.baidu.com/geocoder?output=json&a...

2020-04-07 08:30:23 891 1

原创 python 翻页式爬取百度图片（两种实现方法）

博客是基于两种方法实现对百度图片的翻页式爬取

2020-04-07 07:42:29 2527

原创遗传算法原理，交叉、变异、适应度函数的设置

遗传算法（Genetic Algorithm，GA）由霍兰德教授在20世纪70年代提出，是以自然选择和遗传变异为理论依据的全局性概率搜索优化算法模型。采用遗传算法寻优时需要将问题的候选解进行编码，即一个候选解对应一个编码，编码通常采用二进制，用“0”、“1”表示，为“1”则表示该候选解被选中，所有候选解组合在一起定义为染色体；在迭代进化的过程中，通过构造适应度函数，计算每个个体的适应度，适应度值越...

2020-04-06 06:11:48 22308

原创论文排版之插入图表目录，交叉引用, 出现图0.1的解决方法

笔者论文的排版需要插入图表目录，且图表的命名需要根据带有章节号比如图1.1，图1.2，图2.1.下面是具体的操作的步骤。第一部分插入图表目录由于大部分人在排版时，各级标题的第1章，1.1，1.2是自己手动敲上去的，不是自动的编号，所以首先需要进行给论文编号，如下图：在（开始-编号，可以随便插入一种，这里是给一级标题编号，然后在点击改样式，在定义新编号方式将编号改为第1章或需要的格式...

2020-04-02 17:07:44 13105

原创基于DEM提取河谷网络、流域范围和坡长

1：下载DEM（地理空间数据云、USGS等），使用fill进行洼地填充。2：基于得到的经过洼地填充的DEM，采用flow direction计算流向3：使用flow accumulation 工具计算得到汇流量。河谷网络的提取，基于得到的汇流面积使用使用栅格计算器通过设定阈值将汇流量数据分为0，1，0可代表流量小于阈值的汇流量，1可表示汇流量大于阈值的汇流量，即河谷网络。采用stre...

2020-03-05 09:36:52 4276

原创 arcgis WGS1984地理坐标系转krasovsky_1940_Albers

首先采样arcgis define projection（该工具为定义参考系，不是投影坐标系），将地理参考定义为GCS_Krasovsky_1940，而后才有project raster（该工具定义投影坐标系）设置输出的坐标系统为krasovsky_1940_Albers。如果地理参考系和投影坐标系均为WGS1984系列，则只需要project raster工具即可，当上述两者不一致时，需要适...

2019-12-21 18:26:16 23116 7

原创哨兵2数据预处理大气校正分为2017之前与2017年之后

哨兵2数据下载可以从USGS和欧空局下载，下载网站自行搜索。由于2016年年底，欧空局改变了对哨兵2 数据的命名。哨兵2数据其他博客说2017年之前的数据，可以采用envi5.3打开，同时 envi app store 中也提供了哨兵LIC数据的辐射定标插件，之后采用flassh进行大气校正，具体插件在app store 中搜索 sentinel 2即可，笔者采用这个方法一致有问题，感兴趣...

2019-07-21 17:59:42 6839 1

原创 endnote 只剩.data文件，怎么找回以前的.enl

假如是a.DATA.启动记事本，将其另存为 a.enl。路径和.data文件在一个根目录下（这个很重要）。只需要重新创建一个与.data同名的.enl文件就可以了。使用endnote重新打开 a.enl即可...

2019-07-18 19:02:11 7379

原创利用ENVI实现landsat 5地表温度反演

1：将landsat 5数据进行大气校正2：计算NDVI,3:利用混合相元法计算植被覆盖度fv=代表裸地的NDVI值，代表全部植被覆盖的NDVI值根据自己的实际情况确定上述两个值，比如为-0.2，为0.7利用envi的波段运算（band math）工具输入（b1 gt 0.7）*1+(b1 lt -0.2)*0+(b1 gt -0.2 and b1 lt 0.7)*(b1+...

2019-04-18 21:46:23 26817 38

原创 arcgis制作超链接

Arcgis批量制作超链接1：在属性表中新建属性字段（text），键入超链接地址，2：更改超链接基础基础路径File→Map Document properties→Haperlink base 键入存放图片的路径3：设置属性表超链接字段。layer properties-----display---support hylinks using fields 设置超链接字...

2019-04-04 17:31:46 6245

原创基于arcpy处理 modis hdf格式转geotif

1：modis转tif.有一下几种方法MRT(MODIS REPROJECTION TOOL)HEG(The HDF-EOS To GeoTIFF Conversion Tool)IDL调用MRT2：利用arcgis中集成的arcpy工具包，实现hdf转tif格式ExtractSubDataset_management (in_raster, out_raster, {sub...

2018-10-24 16:05:30 3603

原创 matconvnet 配置。

本人matlab2014a+vs2010,matconvnet-1.0-beta20.win101：在网上下载matconvnet，网址为：http://www.vlfeat.org/matconvnet/2:在matlab下将工作空间切换到 \matconvnet-1.0-beta20\matconvnet-1.0-beta20目录下。3：在matlab命令行输入mex -setup;4:在输...

2018-06-01 10:25:56 272