- 博客(31)
- 资源 (2)
- 收藏
- 关注
RSS订阅原创 关于C指针的初始化理解
指针（变量的储存地址）在C语言的函数中有重要的应用，有必要对指针有更清晰的认知。下面给一个将两个数升序排列的小例子来说明指针的基本理解。（思路均有注释，可参考）假设要排序的变量名字为a，b。 很多人理解 int *p 认为*p是一个整体，实际上这里应该分开读：int * p，p是一个int类型的指针变量。而*p 是一个int变量。对此项的理解关系到对指针的认识。如下例子中，可关注两点:1.为什么在sort函数中运行swap函数没有传入带*的变量。2. a，b传入局部函数是如何可以被彻底互换.
2021-01-12 11:02:41
278
原创 如何自动修正VIC流向
大概800行代码，供各位参考。#include <stdio.h>#include <stdlib.h>#include <string.h>#include "flowgen.h"BOOLEAN Verbose = FALSE; /* flag for optional verbose */char FlowChar[] = {' ', '|', '/', '-', '\\', '|', '/', '-', '\\'};/*char F
2020-11-27 08:59:01
309
原创 C语言动态内存分配需注意的问题
// 数组动态内存分配时，注意查看自己机器的位数。若是32位则4字节，64位8字节。// parr= (int * )malloc(4 * len); // 4 或者8#include <stdio.h>#include <math.h>#include <malloc.h>int * parrvoid f(int *q){ *q = 200;}int main(void){ int len = 20； parr...
2020-11-26 13:21:04
212
原创 GMT安装指南-补充版本
VSCODE 配合linux简直不能太爽！！！！！强烈推荐remote 插件。 在一定程度上弥补了gmt计算能力弱的缺点。安装需要的包如上。gmt，visul studio，gs, sumatra的安装网上均有教程，基本一路点确认即可。本文重点： VS code的配置，以及如何在 cmd 使用 awk 等 linux 小工具。首先安装pdf 和 png viwer。在VS COD...
2019-12-21 20:46:05
2477
2
原创 matplotlib 堆积柱状图-正负叠加问题
没看到很好的解决方案，希望有好的方案的在下方提供链接。目前相对较好的方案：import numpy as npimport matplotlib.pyplot as plt# Juwairia's data: a = [4,-6,9]b = [2,7,1]c = [3,3,1]d = [4,0,-3]data = np.array([a, b, c, d])d...
2019-11-08 16:42:27
3668
1
原创 GRDC径流数据处理
GRDC是国际径流共享网站。下载下来是txt格式如下图：提供一个python小脚本，输入文件，处理成excel年-月-径流的格式。import pandas as pdimport numpy as nppth = 'F:/STA-RUNOFF/AMU/NUSUR.txt' # orignal filesave_pth = 'F:/STA-RUNOFF/AMU/' ...
2019-09-27 11:08:43
4593
24
原创 如何利用 GDAL-Python裁剪栅格数据（DEM等）2019
ARCGIS裁剪数据量比较大的栅格很容易卡死，故用gdal。可批量裁剪多个shp的组合。翻阅了一些博客，这么一个简单的功能写得很复杂。提供代码供有需要的各位参考。from osgeo import gdalinput_raster = r"F:\JUNK\dem90.tif"# or as an alternative if the input is already a gda...
2019-06-02 10:50:03
7251
12
原创 全球降水数据处理专题之--APHRODITE
最近出一批降水数据处理专题，主要提供代码及处理步骤，数据说明请自己查阅。另外，最近qq问问题的太多，故暂时停止qq免费答疑。有需求，请付费。APHRODITE 官网：http://aphrodite.st.hirosaki-u.ac.jp/download/http://aphrodite.st.hirosaki-u.ac.jp/faq_en.html#Q01 FAQ 这些问题先读一...
2019-05-09 15:44:49
6342
2
原创 MODIS10A2积雪数据处理
首先从https://nsidc.org/data/mod10a2 下载数据。 现在modis数据的下载已经非常方便了，以前的教程可忽略。 下载下来为hdf根式，正弦投影。 第一步转换成wgs84坐标系统。转换投影可参见： https://blog.csdn.net/hengcall/article/details/81508700补充下HEG投影工具如何使用：https...
2018-12-10 20:31:53
6858
2
原创 Modis AMSR雪当量处理（待更新）
处理modis数据最麻烦的就是投影问题，翻阅nasa网站后总算有些眉目。记录下来，以供大家参考。例子采用的AMSR 雪当量月数据，这个数据采用的EASE-GRID投影方式，分辨率是25km*25km。整个北半球是721*721的pixels。（具体怎么投影成这样的请见官网）处理思路：把原文件换到wgs84坐标，接下来在matlab或其他就很好做矩阵运算了。首先hdf文件里雪当量是按照m...
2018-11-21 10:36:00
1892
1
原创 NDVI 数据处理，及估算植被覆盖度（一）
若是研究区在中国，则在地理空间数据云有中国区合成好的数据，直接下载用即可。数据格式为TIF。坐标为WGS84，不用再转换投影。下载网址为：http://www.gscloud.cn/sources/?cdataid=265&pdataid=10若研究区在国外，需从nasa下载modis标准数据产品，进行投影转换等工作。后续有时间会补充如何转换投影等内容，可参见《modis AMSR 雪...
2018-11-09 16:57:32
81113
23
原创 根据经纬度自动寻找最近的点，并返回相应属性
import pandas as pdimport numpy as npfrom scipy.spatial.distance import cdistimport scipy.stats as statsfrom pylab import *STA_tmin = 'F:/AAA/STA-T-TMAX-TMIN-CSV/STA-TMEAN.csv'PG_tmin = 'F:/AAA...
2018-09-21 22:05:29
2643
1
原创 NDVI 数据处理，及估算植被覆盖度（二）arcgis 栅格转ascii
公告：最近忙其他事情，暂时停止更新。经过上一步的处理，得到了11年*12月 共132各tif图像。接下来把tif图像转成我们习惯使用的文本格式，即经度，纬度，ndvi。共三列。当然也可以根据个人偏好制作为nc文件等（制作nc文件可参见我的博客）。第一步：把tif文件转为ascii格式的文本文件。所谓ascii格式的文本文件可参见博客：https://blog.csdn.net/kiki...
2018-09-13 20:01:35
32088
14
原创 ARCGIS arcpy批量裁剪命令
全部自动化批量处理，当然右键batch也可以处理，但是一旦有上百个文件，那就运动手指了。代码如下，注释都有，根据需要改动参数即可。 import arcpy# Set the current workspacearcpy.env.workspace = ("F:/ZPATR-TIME-JOB/ndvi-china-/") # "D:\\Test2\\Grids.gdb"...
2018-09-09 18:40:36
7497
3
原创 python统计列表中个元素出现各次数
from collections import CounterCounter（a）a = [9, 9, 3, 1, 1, 2]result = Counter(a)print(result)或者for x in a: if a.count(x)>1: print(x)#
2018-09-04 17:45:35
15576
1
原创 狗十三
对于一个年轻人来说，如果他很早就洞察人事，谙于世故，如果他很快就懂得如何与人交接、周旋，胸有成竹地步入社会，那么不论从理智还是道德的角度来考虑，这都是一个不好的迹象。这预示着他的本性平庸。相反，一个年轻人对世人的行为方式感到诧异和惊讶，并且与他们交往过程中显得笨拙，乖僻，则显示出他有着好贵的品质。 ...
2018-08-25 18:07:22
1324
原创 CMORPH数据处理3
之前是先把二进制文件转成矩阵，再转成nc。写了很多代码，实在复杂。现在用CDO 一句命令就可以直接转为nc。关于CDO的用法可参见官网。很强大的语句。我是在windows的cmd窗口运行CDO的，写了一个bat文件来批量执行任务。安装CDO等请参见气象家园。直接放出bat文件，从1998-2016，有需要的改成自己的目录即可使用。bat文件的具体制作过程不在这里说了。下载好，修改成...
2018-08-17 22:23:20
2764
9
原创 如何制作nc文件（小插曲）
借处理CMORPH降水数据，写一篇关于python制作nc文件的文章。制作nc的模块Matlab， R ，Python 都有，任意选择一个使用即可。之前尝试看过CDO软件，是一个很强大的软件，命令简洁高效，有兴趣可自行google。https://wenku.baidu.com/view/867650d7647d27284a7351ff.html原始矩阵是一个1440*720的二维矩阵，以这个...
2018-08-17 12:18:21
4589
5
原创 CMORPH数据处理2 （请见3 此篇繁琐）
采用数据是v1.0 BLD（融合站点的）daily数据，分辨率是0.25*0.25.数据说明和下载网址见前一篇博客。CMORPH的官方文档暂时没有找到很详细的说明。发现很多人对这个数据的处理都存在疑问，参考了论坛和官网，记录下来供大家参考。欢迎提出更好的思路。首先数据是以二进制字符计录的，4byte float32，这个很重要，找了很久才找到。二进制文件需要转换成能看的数字。文档的re...
2018-08-16 23:05:37
5133
5
原创 CMORPH 数据处理1（下载指南）
共三种版本1：raw 2：CRT 3:BLD 多采用后两者（后面两个是经过站点校正的），下载网址如下：CRT – bias-corrected (over-estimates in the tropics?) 1998/01/01 – 2013/12/31 ftp://ftp.cpc.ncep.noaa.gov/precip/CMORPH_V1.0/CRT/0.25deg-DLY...
2018-08-16 20:06:23
4523
4
原创 GSMaP 降水数据下载指南
1打开官网：http://sharaku.eorc.jaxa.jp/GSMaP/index.htm 点击注册，填好邮箱等。2：收到确认邮件，会给一个ftp地址，进入即可下载。3：可提供数据一览注意;00Z-23Z代表从0点到23点的降水。根据需要，自行选择。 ...
2018-08-16 16:39:42
6009
7
原创 modis AMSR 雪当量数据处理
查看更多文章请点击头像。Modis 处理最麻烦的就是投影转换了，翻阅了下nasa的官网，记录下来，供大家参考。http://nsidc.org/data/AE_MoSno 数据页面。文件格式为25*25km的月雪当量数据，投影为EASE-GRID nasa的投影方式。北半球一共721*721个pixel。注意雪当量以mm计，为了节约储存空间，hdf文件里的雪当量数据被除了一个2，后续处理注...
2018-08-08 22:41:06
2877
7
原创 UNIX paste命令怎么避免文本中间出现制表符
paste默认连接符号是tab^w, 若需要更换为空格则 paste -d ‘ ’ file1 file2 即可。paste的两个或多个文件需提前进行dos2unix的处理。
2018-08-07 11:13:53
484
原创 ERA 处理成ascii
ERA的数据格式是我最喜欢的，读取速度够快。私人记录备忘，非教程。clear alllat_max = 46;lat_min = 32;lon_max = 80;lon_min = 60;save_path = 'F:\AAA\interim\vic\';Vname = 'tp';files = dir('*.nc');days = length(files);%---...
2018-08-04 17:39:43
374
原创 把vic ASCII格式文本转为csv
代码记录：Change vic format file to standar file. 语言：python IDE;pycharm """""==============================================================get annual time series form vic format filenotice: format...
2018-08-04 14:17:56
1896
2
原创 era 数据下载需要注意的问题
首先附上官网的解释，有兴趣可看：https://confluence.ecmwf.int/pages/viewpage.action?pageId=56658233这张图说的非常清楚：对steps选择的疑问： 从 00:00 12:00 出发, 分别有：3, 6, 9, 12 这几个步长可选，例如选0，3，代表从0点开始，到三点的值。如果下载降水就是累积降水，下载温...
2018-08-03 16:44:53
5690
1
原创 ERA 数据下载指南 （python）
说明：数据是用python下载的，安装python请百度。第一步：注册账号，尽量用gmail邮箱。第二步：在c盘 用户/用户名 这个目录下创建一个 .ecmwfapirc的文件（以点开头,以点结尾）第三步：打开https://confluence.ecmwf.int/display/WEBAPI/Access+ECMWF+Public+Datasets/点击会出现以下内容，把内...
2018-08-03 15:06:16
13189
5
原创 用pandas读取大文件
读取上GB大文件的方式有很多，但是并不是都有像pandas这样强大的数据清洗功能。数据太大是，可以分块处理，及时手动删除内存，最后在合并数据。 今天在读取一个超大csv文件的时候，遇到困难：首先使用office打不开 然后在python中使用基本的pandas.read_csv打开文件时：MemoryError 最后查阅read_csv文档发现可以分块读取。read_csv...
2018-05-20 15:48:30
11337
原创 对比两组经纬度，找出最相近的点，直接看代码
两组数据作为例子，函数可直接使用。 import pandas as pdfrom scipy.spatial.distance import cdistdata1 = {'Lat': pd.Series([50.6373473,50.63740441,50.63744285,50.63737839,50.6376054,50.6375896,50.6374239,50.637...
2018-05-15 22:09:41
4338
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人