撼沧-优快云博客

原创使用ArcGIS提取HWSD中的土壤属性数据

基于世界土壤数据库的土壤数据集Siol Map Based Harmonized World Soil Database (HWSD)写在前面：网上关于该数据库的基本介绍以及如何加载到ArcGIS里方法教程较多，此处附上几个链接，读者自行参考。官方链接寒区旱区科学数据中心HWSD数据库的使用孟现勇-如何使用世界土壤数据库(HWSD)建立SWAT模型土壤数据库本文意在将HWSD数据库提取方法进行统一描述，方便自己复用和他人借鉴。HWSD拥有较多土壤属性数据，使用时往往需要将某条土壤属性

2022-01-17 16:18:59 12090 12

原创【Python】将微信收藏的文章批量导出为pdf

将微信收藏的文章批量导出为pdf

2022-01-17 15:48:05 5152

原创 R语言 | 儒略日与日期互相转换 the conversion between Julian day and year-month-day

基于package modiscloudlibrary(modiscloud)简介：R tools for processing Level 2 Cloud Mask products from MODISyearmonthday_to_julianday:Get the julian day for a year/month/day dateyear=2011; month=06; day=15yearmonthday_to_julianday(year, month, day)

2021-10-23 22:39:55 1037

原创批量修改文件名

使用Excel批量修改文件名提取文件夹内的文件名。参考批量提取文件夹内文件名；在Excel中打开文件名文件，第一列为旧文件名。第二列（B列）输入目标新文件名（带上文件名后缀）；输入改名公式在第三列（C列）输入公式：=“ren “&A2&” “&B2&””，需要注意，ren的后面有一个空格，"&A2&“和”&B2&“的中间也有一个空格，两端的引号”"需要闭合，所有字符输入的格式都采用英文格式；将改名公式导出为改名插件将C列中涉及到需要更

2021-08-02 17:17:37 4068

原创【python】无规律时间步长时序数据转为固定步长

写在前面日常可能会遇到时间步长无规律的数据，需要转化为固定时间步长，此时需要进行重采样或插值。例子在选择好时间间隔后，可以用pandas的resample来操作。import pandas as pdfrom pandas.tseries.offsets import Hour,Secondimport numpy as npimport datetimedata = pd.read_csv('fuel.csv')data['date'] = pd.to_datetime(data['d

2021-07-19 14:46:07 1780 1

转载 C语言中的数据类型

写在前面：该文章内容转自wiki，介绍的非常清楚。在C语言中，数据类型可以分为两类：基础数据类型和复合数据类型。目录1 基础数据类型2 复合数据类型• 2.1结构• 2.2联合• 2.3枚举1 基础数据类型注意：以下是典型的数据位长和范围。编译器可能使用不同的数据位长和范围。请参考具体的参考手册。在标准头文件limits.h 和 float.h中说明了基础数据的长度。float，double和long double的范围就是在IEEE_754标准中提及的典型数据。另外，C99添加了新.

2021-05-25 16:07:11 1585

原创使用C语言创建文件夹

使用C语言创建文件夹的代码可以如下：#include <io.h>#include <direct.h>// 创建文件夹void CreateFolder(){ //文件夹名称 char folderName[] = "RunData"; // 文件夹不存在则创建文件夹 if (_access(folderName, 0) == -1) { _mkdir(folderName); }} //然后在main函数中调用CreateF

2021-05-25 15:52:14 11316 1

原创基于R语言进行栅格数据统计及Raster包简介

基于R语言进行栅格数据统计及Raster包简介写在前面：日常数据处理中常需要对单个栅格图层进行统计，包括计算其均值、最值和总和等等。使用R语言可方面快速地实现该需求，主要借助的package为Raster包。library(raster)x = raster('data/data.tif')cellStats(x,stat='mean')x表示待计算的栅格；stat表示需要统计的内容，包括sum, mean, min, max, sd（标准差）, skew（偏度）和 rms（均方根）；

2021-05-25 15:41:58 12042

原创 ArcGIS许可管理器License server administrator无法启动的解决方法

ArcGIS许可管理器License server administrator无法启动的解决方法在安装ArcGIS软件，甚至日常使用启动软件时，ArcGIS许可管理器License Manager经常会无法启动，点多少次都没反应，更别谈往后的破解了。这里谈一下自己遇到问题后用过的能用的方法。

2021-01-07 22:50:16 120713 26

原创移除word文档的域代码 remove all field codes of word

去除word中因使用mendeley或endnote等文献管理软件而产生的域代码。网上很多说使用 CTRL+SHIFT+F9，但亲测不行，可能是软件版本问题。使用以下方法解决：CTRL+A 全选文章；CTRL+6 去除全文域代码。参考：https://community.endnote.com/t5/EndNote-Styles-Filters-and/Removing-field-codes-adjusts-spacing-in-Word-document/td-p/152892http

2021-01-07 22:26:09 17234

原创 python绘制1:1对角线（1:1 line）

本文基于python绘制1:1对角线并x、y数据之间的R方。效果图：蓝色线表示拟合线，黑色线表示1:1线。详细代码：import matplotlib.pyplot as pltimport numpy as npimport pandas as pdfrom pylab import mplfrom matplotlib.pyplot import MultipleLocatormpl.rcParams['font.sans-serif'] = ['SimHei']mpl.rcPa

2021-01-06 09:42:52 9688 2

原创基于python进行小波分析，频率谱分析

该方法基于python进行时间序列的小波分析并出图（包括功率谱图和小波分解后的图）。默认的小波为morlet小波。该代码由 Evgeniya Predybaylo 博士提供：https://github.com/chris-torrence/waveletshttps://atoc.colorado.edu/research/wavelets/Copyright © 1995-2004, Christopher Torrence and Gilbert P.CompoPython versio

2020-12-28 14:50:03 24073 26

原创基于R语言进行热浪指数的计算_HWMid

写在前面当前对于高温热浪事件指标的定义缺乏统一标准，不同学者有不同的算法。Russo（2015）提出了一种每日热浪幅度指数（HWMId，the heat wave magnitude index daily）计算方法，适合对全球的热浪强度进行计算和对比分析。该方法也被wiki百科对热浪指数的介绍中采用（wiki）。该方法作者提供了相关的R语言包和函数，本文对此进行原理和代码的介绍。放上wiki上的一张图（Zampieri, M., 2016）1. 原理介绍HWMId定义为每年最大热浪事件的幅

2020-12-26 22:57:58 4431 11

原创基于R语言读取NC数据并提取单个波段（netcdf转为tiff或其他格式）

写在前面netcdf是气象、生态、地理领域非常常用的一种数据格式，读写的方法很多，arcgis、arcpy、python、matlab等等均有相应的方法。本文基于R语言实现netcdf的数据读取，提取数据中的单个波段（某一时间的数据），展示和写出。注：本文仅展示基础的代码，多个波段或多个数据的读写操作加上相应的循环即可。具体代码加载包library(raster)library(rasterVis)library(ncdf4)library(lattice)打开nc数据，并展示其中

2020-12-26 00:33:43 10942 2

原创基于R语言或Matlab直接读取压缩包中的数据

写在前面以gosif网站上的tiff数据为例。该数据是地学上常用的一种数据，希望读取该网站上的gosif数据：http://data.globalecology.unh.edu/data/GOSIF_v2/8day/但数据全是压缩包格式的，解压后内存占用又会变大。所以需要直接从压缩包中读取数据。基于R语言library(R.utils)library(stringr)library(raster)library(rasterVis)# download and read the dat

2020-12-23 12:47:01 3434

原创 Voronoi图 | 泰森多边形

概念Voronoi图，又叫泰森多边形、沃罗诺伊图或Dirichlet图，是一个关于空间划分的基础数据结构。它是由一组由连接两邻点直线的垂直平分线组成的连续多边形组成。N个在平面上有区别的点，按照最邻近原则划分平面；每个点与它的最近邻区域相关联。Delaunay三角形是由与相邻Voronoi多边形共享一条边的相关点连接而成的三角形。Delaunay三角形的外接圆圆心是与三角形相关的Voronoi多边形的一个顶点。应用Voronoi图的应用非常广泛。在计算几何中，重心Voronoi图(CVT)方法被用

2020-12-09 15:15:30 10036

原创数字城市、智能城市、智慧城市、数字中国和智慧社会简析

1. 百度百科上各术语的概念数字城市： “数字城市”（英文：digital city） [1] 以计算机技术、多媒体技术和大规模存储技术为基础，以宽带网络为纽带，运用遥感、全球定位系统、地理信息系统、遥测、仿真-虚拟等技术，对城市进行多分辨率、多尺度、多时空和多种类的三维描述，即利用信息技术手段把城市的过去、现状和未来的全部内容在网络上进行数字化虚拟实现。城市：智能城市是把基于感应器的物联网和现有互联网整合起来，通过快速计算分析处理，对网内人员、设备和基础设施实施，特别是交通、能源、商业、安全、医疗等

2020-12-09 15:02:22 3263

原创 c语言中define和typedef的区别

c语言中define和typedef的区别define只是简单的文本替换，typedef是声明一种新的类型，等同自带的基本类型。它们完全是两个东西，是由不同的“程序”处理的：#define是宏，处理的时候位于编译前阶段，宏处理器基本上对你的C/C++程序不会有任何的感知。它只处理宏的语法。而编译阶段的“程序”得到的是宏处理完的结果。typedef是编译阶段的一部分。它的意义是单一的。用宏来定义类型，只是它恰好替换完字符串，可以达到那样的效果。但没有任何措施能够阻止你用它干别的事情。

2020-12-09 14:57:40 842

原创 c语言结构体的定义和使用

一、什么是结构体在实际问题中，一组数据往往具有不同的数据类型；例如在学生信息登记表中，姓名为字符型，学号为整型或字符型，年龄为整型，性别为字符型，成绩为整型或实型。因为数据类型不同，显然不能用一个数组来存放。在C语言中，可以使用结构体（Struct）来存放一组不同类型的数据。结构体(struct)指的是一种数据结构，是C语言中聚合数据类型(aggregate datatype)的一类。结构体可以被声明为变量、指针或数组等，用以实现较复杂的数据结构。结构体同时也是一些元素的集合，这些元素称为结构体的成员(

2020-12-09 14:55:12 74476 12

原创 c语言产生随机数的方法

随机数种子什么是随机数种子电脑不是人，不能随便就凭空想象一个数字出来。所谓的随机数,不过是通过一个函数经过计算出来的数字而已。因此,随机数不过是函数计算出来的值，而一个函数要得到这个值Y，肯定需要一个输入X，所以X就是那个种子。简单点说 Y = function(x)；随机数 = Y，随机数种子 = X。第二，rand()%Y+1是什么意思。首先你要弄懂%是什么意思,是取余运算。举个例子 10 % 3 = 1 ，10 % 5 = 0;就是得到余数。余数的范围永远是从 0到被除数-1 不可能等

2020-12-09 14:48:20 3088 1

原创 DIR命令帮助文档

帮助文档：help DIR:显示目录中的文件和子目录列表。DIR [drive:][path][filename] [/A[[:]attributes]] [/B] [/C] [/D] [/L] [/N] [/O[[:]sortorder]] [/P] [/Q] [/R] [/S] [/T[[:]timefield]] [/W] [/X] [/4] [drive:][path][filename]:指定要列出的驱动器、目录和/或文件。 /A:显示具有指定属性的文件。属性:D:目录;

2020-12-09 14:46:23 630 1

原创【python】使用python绘制地图时手动添加比例尺

写在前面使用basemap、cartopy等包在python中绘制地图添加比例尺有投影限制，不太方便。本博文使用最基础的matplotlib画线实现比例尺的添加，为大家提供一种可选思路。思路：提供了三种可选风格，根据自己需要进行选择。原理为根据1°=111千米的原理计算需要绘制的长度，因此本代码主要适用于低纬度地区的wgs84投影的绘图。高纬度地区投影变形太大，使用这种方法会不准确。大家也可以根据这种思路自己修改代码。详细代码：#-----------函数：添加比例尺------------

2020-12-08 23:15:14 6752

原创【python】使用python绘制地图时添加指北针

本博文基于matplotlib，手动实现在python中绘制地图时添加指北针。代码如下：from mpl_toolkits.basemap import Basemapimport matplotlib.pyplot as pltimport matplotlib.patches as mpatches#-----------函数：添加指北针--------------def add_north(ax, labelsize=18, loc_x=0.88, loc_y=0.85, width=0

2020-12-08 22:54:21 5039

转载【python】【Gif制作】使用多张图片合成gif动图

使用多张图片合成gif动图

2020-12-07 10:18:06 5577

原创【R语言】使用r语言进行单个或多个栅格数据间的计算

本文介绍如何使用r语言进行单个或多个栅格数据间的计算，类似于ArcGIS中的栅格计算器。单个栅格数据：library(raster)calc(x, fun, filename='', na.rm, forcefun=FALSE, forceapply=FALSE, ...)元素：x：Raster* objectfun：functionfilename：character. Output filename (optional)Examples# NOT RUN {r <- rast

2020-12-07 10:08:31 8170

原创【R语言】栅格数据快速制图-levelplot函数（2）

rasterVis包中的levelplot函数除了出单张图，也可以快速随多张地图进行制作。单张地图制图参考：【R语言】栅格数据快速制图-levelplot函数此处主要使用raster包中的stack函数将多幅图合为影像集。#栅格数据包括.tiff、.nc、.flt等常见格式library(raster)library(rasterVis)library(lattice)#设定工作空间setwd('E:\\hancang')#读取文件夹中的全部文件files = list.files(

2020-12-07 09:53:21 5762 5

翻译【R语言】栅格数据快速制图-levelplot函数

使用R语言中的levelplot函数可以快速方便地对栅格数据进行制图，非常实用，下面做简单地讲解实例。栅格制图# 加载packagelibrary(raster)library(rasterVis)library(lattice)基础展示:根据坐标轴展示的栅格图像；两个边缘分布展示为行列的平均值；颜色图例；levelplot(rprob)对图像的栅格值进行对数变换(zscaleLog=TRUE)，并添加水平等值线(contour=TRUE)。levelplot(rprob,

2020-12-06 23:05:47 6918 2

原创自然资源部标准地图底图转矢量Shapefile并配准

当前许多期刊投稿要求凡涉国界图件（国内部分地区、全国、世界部分地区、全球）必须使用自然资源部标准地图底图（下载网址：http://bzdt.ch.mnr.gov.cn），所用底图边界要完全无修改。但该网站并未提供可编辑的.shp矢量文件，除p图外，无法直接应用于自己的数据，在arcgis或其他软件进行作图。因此，本文介绍了如何将标准地图底图转为Shapefile矢量图层，及如何将自己的数据与标准底图进行配准，使shp边界对齐，达到在arcgis里出图的目的。

2020-12-04 10:18:02 38061 24

原创 C Primer Plus第八章课后答案

其他章节答案/*Project 1*/#include<stdio.h>int main(void){ int ch; int n = 0; while ((ch = getchar()) != EOF) n++; printf("%d",n); return 0; } /*Project 2*/#include&l...

2020-10-15 20:44:47 506 2

原创 C Primer Plus第七章课后答案

其他章节答案第七章课后答案/*Project 1*/#include <stdio.h>int main(void){ int space = 0; int low = 0; int chars = 0; char ch; printf("Please enter:"); while((ch = getchar()) != '#...

2020-10-14 18:30:52 708

原创 C Primer Plus第六章课后答案

其他章节答案//Project 1#include <stdio.h>#define SIZE 26int main(void){ char letter[SIZE]; char l = 'a'; int n; for (n = 0;n < SIZE;n++,l++) { letter[n] = l; }...

2020-10-14 18:30:37 930 1

原创 C Primer Plus第五章课后练习答案

其他章节答案//Project 1#include <stdio.h>#define S_M_PER 60int main(void){ int min_1,hour,min_2; printf("Please enter the minute: "); while(scanf("%d",&min_1)>0) { ...

2020-10-12 13:39:19 591

原创 C Primer Plus第四章课后练习答案

其他章节答案/*The answer for Unit 4*//*//Project 1#include <stdio.h>int main(void){ char name[40],surname[40]; printf("Please enter your name and surname.\n"); scanf("%s%s",name,surn...

2020-10-12 13:38:48 594

原创 C Primer Plus第三章课后练习答案

其他章节答案/*The third Unit programming answer*///Project 1#include <stdio.h>int main(void){ int i=2147483647; float f=3.4E38,F=0.1234E-10; printf("Integer overflow:\n\ti=%d,\n\ti+1=...

2020-10-11 09:52:55 857

原创 C Primer Plus第二章课后练习答案

其他章节答案/*Answer of the fourth question */#include <stdio.h>void jolly(void); //Note that there are semicolons here!void deny(void);int main(void) //Note that there are't semicolons here!{ ...

2020-10-11 09:52:15 1347

原创 C Primer Plus第九章课后答案

其他章节答案/*Project 1*/#include<stdio.h>double min(double,double);int main(void){ double x,y; printf("请输入一个double类型的值："); scanf("%lf",&x); /*double类型要用%lf*/ printf("请再输入一个do...

2020-10-10 22:49:19 597

原创 C Primer Plus第十章课后答案

其他章节答案/*Project 1*/#include<stdio.h>#define MONTHS 12#define YEARS 5int main(void){const float rain[YEARS][MONTHS] = { {4.3,4.3,4.3,3.0,2.0,1.2,0.2,0.2,0.4,2.4,3.5,6.6}, {8.5,8.2...

2020-10-10 22:45:36 785

原创文献图复现-方块图-python

来源Guo, W., et al., Vegetation can strongly regulate permafrost degradation at its southern edge through changing surface freeze-thaw processes. Agricultural and Forest Meteorology, 2018. 252: p. 10-17.解释用数字表示矩阵内包含的点数，颜色表示另外一个数据。复现图复现代码（python）impor

2020-08-03 13:51:59 1362 1

原创文献图复现-散点图-python

文献图来源Guo, W., et al., Vegetation can strongly regulate permafrost degradation at its southern edge through changing surface freeze-thaw processes. Agricultural and Forest Meteorology, 2018. 252: p. 10-17.解释颜色和形状表示类别，灰色线表示标准差，适合用来展示不同类型对象的分布。复现图复现

2020-08-03 00:55:51 734

原创【R语言】根据经纬度从tif图像中提取兴趣点

library(stringr)library(raster)library(sp)#set work spacegetwd()setwd("G:/data")#初始设置，提取点lai <- matrix(nr=50,nc=2)point<-cbind(128.5,45.9) #lon,lati <- 1n <- 1repeat { #读取栅...

2020-04-03 10:45:28 4061 2