weixin_37659245的博客

原创 菜鸡的自救之路（四）

今天是2021年9月1日。按照国内的计算方法的话，同时也可以看作是我博士二年级的第一天。值此特殊时刻，书写“菜鸡自救之路”第四部。距离上一篇总结已经过去了一年零7个月，这段时间我都做了什么呢？完成毕业论文，顺利毕业；正式开始我的PhD生涯；发表一篇SCI一区论文，投稿一篇一区，待投稿一篇。同时，手头有三个研究工作在并行，预期今年可以完成。...

原创 菜鸟的自救之路（三）

距离上一篇总结已经是5个月了，也到了2020年。小结：科研进展：8月在迷茫，9月、10月修改文章并投出去，11月被接受。虽然是一个开源期刊，但自己还算满意。总算是给自己的“科研之路”开了一个头。但是这么长的时间里，其实没有任何新的实验进展，这是我最为难过的地方。留学申请：拿到了欧洲两个大学的Offer，老师都是领域内最尖端的大佬。同时也在和英国几个TOP学校的老师保持联系、修改RP，因为还...

原创 菜鸡的自救之路（二）

今天是7月28日，相比上一篇文章[《菜鸡的自我拯救》]已经过去4个月了，在此做个小结，重做计划。(https://blog.csdn.net/weixin_37659245/article/details/88928659)

原创 菜鸡的自救之路

博主研二在读，浑浑噩噩，一无所学。一朝面临毕业求职，惶惶恐恐，虽无大志，但也不甘堕落。于是今日开启菜鸡的自我拯救之路，争取通过数月乃至一年的努力，不求成为大神，起码脱离菜鸡称号，等来年春招，在来年春招面临残酷的竞争时不至于任人宰割。首先明确需要点亮的技能树都有哪些：基础算法数据结构数据挖掘机器学习算法深度学习算法其他（爬虫、可视化等）毕业要紧，最理想的情况当然是在满足科研项目的...

原创 GEE批量下载：控制台代码

function runTaskList() { var tasklist = document.getElementsByClassName('awaiting-user-config'); for (var i = 0; i < tasklist.length; i++) tasklist[i].children[2].click(); } function confirmAll() { var ok = document.getElemen..

原创 第五篇个人总结

今天是2022年2月11日。截至今天，博士期间发表3篇一区论文，投稿两篇。手头有3篇论文在进行中，等待投稿。以上8个工作，算作博士研究的第一阶段吧。下一个阶段，希望能够做一些有挑战性的工作。例如对地理定律的探索。这时就深刻地意识到自己其实过往有些“温水煮青蛙”，当然我是被煮的那一个。不知不觉中，已经太久没有输入新技术了。记得《菜鸡的自救之路》这个系列创建的初衷就是研二的某一天，下定决心要学习编程，一改往日代码小白的身份。最后的结果差强人意，学了点编程、机器学习的皮毛，但也总算是完成了不少工作，尤其是毕业

原创 用R语言实现spatial join

今天实现了一个任务：有一个library(rgdal)library(gstat) # Use gstat's idw routinelibrary(sp) # Used for the spsample functionlibrary(raster) # Used to clip out thiessen polygonsimos=read.csv('D://research//GHM//imos.csv')coordinates(imos) <- c("lon","

原创 Python：多个tiff转csv，并按照像元合并成一个csv

1）定义tiff转csv的函数from osgeo import gdalimport pandas as pdimport numpy as npgdal.AllRegister()def tifftocsv(path): filePath = path dataset = gdal.Open(filePath) adfGeoTransform = dataset.GetGeoTransform() nXSize = dataset.RasterXSize

原创 Python: 计算时间序列影像的斜率图和均值图

1）定义斜率计算的函数from osgeo import gdalimport numpy as npfrom sklearn import linear_modelimport copy# 计算斜率def calculate_slope(data): reg = linear_model.LinearRegression() reg.fit(np.array(range(len(data))).reshape(-1, 1), np.array(data).reshape(-1,

原创 R语言学习：批量tiff转PDF并且组图

library(raster)library(ggplot2)library(rgdal)library(sf)library(viridis) library(ggthemes) library(ggpubr)setwd('D:/research/5_AOD/results/Analysis/Slope/test/')list <- list.files()#list=list[12:23]j=0for(i in list){j=j+1 # print(i)

原创 修改别人R包，并编译

当前有一个任务，需要使用他人的R包，并且进行一定程度上的修改，重新编译使用。经过探索，终于…成功了。1）下载R包源文件可能你电脑上已经有安装好的该R包，这时候最好是提前删除掉，防止后面有了两个同名R包出现冲突。用remove.packages(‘XXX’)。之后直接从Google搜索R包名称下载。然后在本地解压缩。解压缩后的文件结果如下图：2) 修改源文件打开Rstudio，打开需要修改的文件。3）建立新project，编译Rstudio- File- New project, 选择：E

原创 R语言ggplpt绘制地理图

近期在用ggplot做论文的图，深感这个绘图包的博大精深。加载包和数据：## plot raw data of soil moisturelibrary(ggplot2)library(sf)library("rnaturalearth")library("rnaturalearthdata")# 加载数据data_fr<-read.csv('D:\\research\\2_STH\\data_fr_0430_class.csv')data_th<-read.csv('D:

原创 问题解决：R语言将num格式列读取成chr

原因是该列存在NULL值，但是NULL值又不是用"NULL"字符表示的，因此R语言傻乎乎的就以为这一列是字符型变量，而不是数字型。（str()函数查看数据的格式）解决办法：我对R十分不熟悉，因此就用python先读取了文件, fillna(0)，导出。然后从R中重新读取，就没问题了。...

原创 GEE统计区域气溶胶差值

var shp=tablevar shp1 = shp.filterMetadata('Result_GWR','less_than','30000');var shp2=shp.filterMetadata('Result_GWR','not_less_than','30000');//Map.addLayer(shp1)// var images = imageCollection4.filterBounds(shp1)// .filter(ee.Filter.date('20.

原创 基于R语言GD包的Risk Map制作（批量生成）

短短的代码…踩了多少坑过去完全没有接触过R语言，很多本该知道的基础知识都不清楚，遇到了盲点就要花很长的时间解决。先放总代码：```{r}factors2011<-data[,c(4,7,9:39)]name2011<-names(factors2011)name2011<-name2011[-1]for(i in 1:length(name2011)){ b <- ndvigdm$Risk.mean[[name2011[i]]]$meanrisk

原创 使用R包GD实现地理探测器算法

首先载入数据data<-read.table(file="D:/GWRR/统计总表-终版.csv",header=TRUE,sep=',')然后由于我的研究需要两个年份的，factors2011<-data[,c(4,7,9:39)]factors2016<-data[,c(5,7,9:39)]定义分组discmethod <- c("equal","natural","quantile")discitv <- c(4:6)运行ndvigdm <

原创 使用Arcpy批量进行可达性分析

# Name: MakeODCostMatrixAnalysisLayer_Workflow2.py# Description: Find the travel time to the closest hospital from each census# tract and join the travel time and hospital name to the input# tracts.# Requirements: Network An

原创 对科研实验的流程管理思考

近期在做一个研究，流程图与方法都早早确定了，唯一的工作就是尽快做出结果。但是我预期认为很简单的操作，却遭遇了很多问题，也走了一些弯路，主要是做了一些重复的过程，浪费了宝贵的时间。因此，在这里总结下我得到的教训，以及在其中的收获。

原创 使用R语言的spgwr包进行地理加权回归（GWR）

今天需要做要给GWR回归，使用了R语言中的spgwr包。做的工作主要有；1）安装R语言直接去R官网就好。2）安装spgwr包> install.packages("spgwr")发现报错Warning in install.packages : 'lib = "C:/Program Files/R/R-4.0.3/library"' is not writableWarning in install.packages : cannot create dir 'D:\onedri

原创 使用arcpy批量计算路网可达性（网路分析方法）

# Name: MakeODCostMatrixAnalysisLayer_Workflow2.py# Description: Find the travel time to the closest hospital from each census# tract and join the travel time and hospital name to the input# tracts.# Requirements: Network An

原创 计算各行政区划内的人口加权质心

本研究的目的是使用人口栅格数据计算每个行政区划内的人口加权质心。

原创 初学ArcGIS Pro:遇到并解决的一些问题

SHP数据修改显示属性（style）的时候，没有可选择的样式原因是没有连接上System styleInsertAddAdd system style有时候发现显示是选上了

原创 问题解决：ArcGis的shp连接表格时，显示不出对应字段名称

如图，想要用CSV匹配shp（含同一字段名称），但是在join时却显示不出字段名称，实际上也匹配不上数据解决方法：将改CSV文件导出为xls文件（老版excel格式），即可匹配

原创 时间序列聚类

参考这篇文章

原创 GIS空间分析方法（一）：学习计划与目录

近期愈发感觉自己对于一些经典的空间分析方法都不甚了解,极大地限制了科研思路。因此打算抽空系统地补一下这方面的知识，目录暂定如下：地理加权回归克里金插值地理加权的克里金回归插值Moran I index混合度、区位熵...

原创 与师长谈科研

因为对未来要面临的博士生涯充满未知，故而有些焦虑。在我是生命中，这样的事情已经发生了太多次。在这样的处境中，我通常是寻找我将要前往的处境中的长辈去咨询、探讨。近期和一些优秀的青椒聊了一些，把个人的感悟总结如下：

转载 EXCEL复制选取筛选后的单元格（而不是全部）

1）对筛选后的数据进行选定。选定需要的数据区，如图所示需要选定B2:C11,但其中的B4:C5并不是我们想要的，所需要使用“选定可见单元格的操作”选定数据区域后，使用CTRL+G打开定位选择，如图所示：2）使用定位条件选择可见单元格。点击进入后，可看如图所示，选择可见单元格，后点击确定。3）选择复制和粘贴。点击确定后可发现源数据区域变成了一层灰色~代表数据的选定状态，此时使用CTRL...

原创 EXCEL多列vlookup匹配

关键有二：$字符该符合保证了后面的字符不会随着拉动而自动变化Column(A1)函数该函数返回输入变量所在的列名

原创 GIS学习（二）：插值分析与密度分析的区别

密度分析与插值分析的区别1 定义1 .1插值分析在实际工作中，由于成本的限制、测量工作实施困难大等因素，我们不能对研究区域的每一位置都进行测量（如高程、降雨、化学物质浓度和噪声等级）。这时，我们可以考虑合理选取采样点，然后通过采样点的测量值，使用适当的数学模型，对区域所有位置进行预测，形成测量值表面。插值之所以可称为一种可行的方案，是因为我们假设，空间分布对象都是空间相关的，也就是说，彼此接...

原创 GIS学习（一）：密度分析、可达性分析、空间句法

1)密度分析2) 路网可达性分析所需数据：成本栅格、目标点成本栅格一般由路网矢量转换栅格、重分类（nodata为高成本值）权重的设置：人为设置权重、空间句法设置权重3)空间句法人与环境相互影响可以对路网进行空间句法分析，用作可达性分析的权重...

原创 我对科研的一点看法

在读文献、做实验的过程中，会逐渐认识到，真正决定一个人科研产出上限的，并不是接触到了多么前沿的技术，而是对本领域基础知识的理解程度。本科时专业课上走过的神、偷过的懒，早晚有一天会成为限制自己idea的桎梏。对于绝大多数的博士生来说，那种开天辟地的创新是不现实的，大多数让人眼前一亮的研究可能也仅仅是在本领域基础理论上的一点小小的改进，甚至只是小小的结合。就我本人对科研肤浅的理解以及薄弱的产出，...

原创 问题解决：Python作图时坐标轴label显示不全

问题描述在用matplotlib或sns作图时，横坐标轴是每隔两个显示一个，如图：代码如下：import seaborn as sns​cmap = sns.cubehelix_palette(start = 1.5, rot = 3, gamma=0.8, as_cmap = True)​sns.heatmap(outdata02,linewidths = 0.05, cmap=...

原创 如何利用好手机信号数据

现有新疆乌鲁木齐市中心的手机信号数据，是一个真的特别棒的数据。一定要利用好这个数据，做一些研究，并且完成毕业设计。数据描述手机信号数据分为两种：居住人口定义：1）居住地观测时段：21:00到次日8:00；2）用户每日在观测时间段内被观测到的秒数，进行月度累加，并进行排名，取排名最高的为用户的居住地；3）一个月内出现天数超过10天。包含字段：归属地（本地/外地）、富裕度、年龄...

原创 POI高德坐标转WGS84

因为下载的高德POI数据，相比WGS84是有一定的坐标偏移，因此需要进行校正一下。首先，文件夹是这样的：CSV文件格式为：代码如下：#!/usr/bin/python# -*- coding: UTF-8 -*- import pandas as pdimport jsonimport mathimport osimport csv## 转换函数：x_pi = 3...

原创 sklearn的feature_importances_含义是什么？

参考这篇回答：https://stackoverflow.com/questions/15810339/how-are-feature-importances-in-randomforestclassifier-determined/15821880#15821880其中有Sk-learn作者的回复：There are indeed several ways to get feature “...

原创 Pandas显示全部行、列

尤其是在jupyter中，显示不全有时会很不方便import pandas as pd#显示所有列(参数设置为None代表显示所有行，也可以自行设置数字)pd.set_option('display.max_columns',None)#显示所有行pd.set_option('display.max_rows',None)#设置数据的显示长度，默认为50pd.set_option('...

原创 问题解决：Word多了一页空白页

原创 ArcGIS/ENVI学习（一）：近期一些操作总结

1）计算城市比例目的：使用清华10m的全球土地覆盖数据，得到城市像元所占的比例影像Arcgis操作：重分类：将本来属于不透水层的8，重分类为1，其余类别（像元值）重分类为0块统计：对一个33的矩形，计算mean值。即可得到30m30m的城市区块比例影像（注意，此时影像的分辨率仍旧为10m，只是每一个30m*30m的区块都是同样的数值）重采样：重采样至所需的分辨率2）计算某影像的标...

原创 问题解决：pdf打不开，出现Acrobat failed to connect to a DDE server问题

1、开始菜单——输入services.msc回车，启动服务页面2、找到Network List Service Properties，对其右键选择属性3、设置启动类型为自动，应用即可。

原创 Python学习：pandas读取大文件的一些经验

1、设置CSV的上限2、加载进度条（读取太慢时，可以判断程序是否还在运行）3、循环读进pandas比如，想要从全国的POI数据中，读取四川省的：import pandas as pdimport numpy as npimport osfrom tqdm._tqdm_notebook import tnrangefilepath='D:\\高德POI数据全国\\2018-POIC...