如何利用 GDAL-Python裁剪栅格数据(DEM等)2019

ARCGIS裁剪数据量比较大的栅格很容易卡死,故用gdal。 可批量裁剪多个shp的组合。 翻阅了一些博客,这么一个简单的功能写得很复杂。提供代码供有需要的各位参考。 from osgeo import gdal input_raster = r"F:\JUNK\dem90.t...

2019-06-02 10:50:03

阅读数 148

评论数 1

全球降水数据处理专题之--APHRODITE

最近出一批降水数据处理专题,主要提供代码及处理步骤,数据说明请自己查阅。 另外,最近qq问问题的太多,故暂时停止qq免费答疑。有需求,请付费。 APHRODITE 官网:http://aphrodite.st.hirosaki-u.ac.jp/download/ http://aphrodi...

2019-05-09 15:44:49

阅读数 222

评论数 0

新年开工

今年将持续不定期更新博客,欢迎交流。 3周不来办公室,甚好!

2019-02-12 16:15:41

阅读数 106

评论数 1

MODIS10A2积雪数据处理

  首先从https://nsidc.org/data/mod10a2 下载数据。 现在modis数据的下载已经非常方便了,以前的教程可忽略。   下载下来为hdf根式,正弦投影。 第一步转换成wgs84坐标系统。 转换投影可参见: https://blog.csdn.net/hengca...

2018-12-10 20:31:53

阅读数 646

评论数 1

Modis AMSR雪当量处理(待更新)

处理modis数据最麻烦的就是投影问题,翻阅nasa网站后总算有些眉目。记录下来,以供大家参考。 例子采用的AMSR 雪当量月数据,这个数据采用的EASE-GRID投影方式,分辨率是25km*25km。整个北半球是721*721的pixels。(具体怎么投影成这样的请见官网) 处理思路:把原文...

2018-11-21 10:36:00

阅读数 272

评论数 0

NDVI 数据处理,及估算植被覆盖度(一)

若是研究区在中国,则在地理空间数据云有中国区合成好的数据,直接下载用即可。数据格式为TIF。坐标为WGS84,不用再转换投影。下载网址为:http://www.gscloud.cn/sources/?cdataid=265&pdataid=10 若研究区在国外,需从nas...

2018-11-09 16:57:32

阅读数 17003

评论数 7

考虑海拔的IDW的插值

import pandas as pd from geopy.distance import vincenty import numpy as np loworhigh = 'high' # 输入需要哪段 ,low 代表第一段,high代表第二段。 year = 1951 # 是否要选...

2018-10-02 18:22:23

阅读数 419

评论数 0

根据经纬度自动寻找最近的点,并返回相应属性

import pandas as pd import numpy as np from scipy.spatial.distance import cdist import scipy.stats as stats from pylab import * STA_tmin = 'F:/AAA/S...

2018-09-21 22:05:29

阅读数 655

评论数 0

NDVI 数据处理,及估算植被覆盖度(二)arcgis 栅格转ascii

公告:最近忙其他事情,暂时停止更新。 经过上一步的处理,得到了11年*12月 共132各tif图像。 接下来把tif图像转成我们习惯使用的文本格式,即经度,纬度,ndvi。共三列。当然也可以根据个人偏好制作为nc文件等(制作nc文件可参见我的博客)。 第一步:把tif文件转为ascii格式的...

2018-09-13 20:01:35

阅读数 5952

评论数 5

ARCGIS arcpy批量裁剪命令

  全部自动化批量处理,当然右键batch也可以处理,但是一旦有上百个文件,那就运动手指了。 代码如下,注释都有,根据需要改动参数即可。   import arcpy # Set the current workspace arcpy.env.workspace = (&quo...

2018-09-09 18:40:36

阅读数 847

评论数 0

python统计列表中个元素出现各次数

from collections import Counter Counter(a) a = [9, 9, 3, 1, 1, 2] result = Counter(a) print(result) 或者 for x in a: if a.count(x)&g...

2018-09-04 17:45:35

阅读数 2685

评论数 0

狗十三

对于一个年轻人来说,如果他很早就洞察人事,谙于世故,如果他很快就懂得如何与人交接、周旋,胸有成竹地步入社会,那么不论从理智还是道德的角度来考虑,这都是一个不好的迹象。这预示着他的本性平庸。相反,一个年轻人对世人的行为方式感到诧异和惊讶,并且与他们交往过程中显得笨拙,乖僻,则显示出他有着好贵的品质。...

2018-08-25 18:07:22

阅读数 849

评论数 0

CMORPH数据处理3

之前是先把二进制文件转成矩阵,再转成nc。写了很多代码,实在复杂。现在用CDO 一句命令就可以直接转为nc。 关于CDO的用法可参见官网。很强大的语句。 我是在windows的cmd窗口运行CDO的,写了一个bat文件来批量执行任务。 安装CDO等请参见气象家园。 直接放出bat文件,从1...

2018-08-17 22:23:20

阅读数 608

评论数 6

如何制作nc文件(小插曲)

借处理CMORPH降水数据,写一篇关于python制作nc文件的文章。制作nc的模块Matlab, R ,Python 都有,任意选择一个使用即可。之前尝试看过CDO软件,是一个很强大的软件,命令简洁高效,有兴趣可自行google。https://wenku.baidu.com/view/8676...

2018-08-17 12:18:21

阅读数 809

评论数 3

CMORPH数据处理2 (请见3 此篇繁琐)

采用数据是v1.0 BLD(融合站点的)daily数据,分辨率是0.25*0.25.数据说明和下载网址见前一篇博客。 CMORPH的官方文档暂时没有找到很详细的说明。 发现很多人对这个数据的处理都存在疑问,参考了论坛和官网,记录下来供大家参考。欢迎提出更好的思路。 首先数据是以二进制字符计录...

2018-08-16 23:05:37

阅读数 1188

评论数 4

CMORPH 数据处理1(下载指南)

共三种版本1:raw  2:CRT  3:BLD   多采用后两者(后面两个是经过站点校正的),下载网址如下: CRT – bias-corrected (over-estimates in the tropics?) 1998/01/01 – 2013/12/31 ftp://ftp.c...

2018-08-16 20:06:23

阅读数 815

评论数 3

GSMaP 降水数据下载指南

1打开官网:http://sharaku.eorc.jaxa.jp/GSMaP/index.htm  点击注册,填好邮箱等。 2:收到确认邮件,会给一个ftp地址,进入即可下载。 3:可提供数据一览 注意; 00Z-23Z代表从0点到23点的降水。根据需要,自行选择。  ...

2018-08-16 16:39:42

阅读数 876

评论数 4

modis AMSR 雪当量数据处理

查看更多文章请点击头像。 Modis 处理最麻烦的就是投影转换了,翻阅了下nasa的官网,记录下来,供大家参考。http://nsidc.org/data/AE_MoSno 数据页面。文件格式为25*25km的月雪当量数据,投影为EASE-GRID nasa的投影方式。北半球一共721*721个...

2018-08-08 22:41:06

阅读数 642

评论数 3

UNIX paste命令怎么避免文本中间出现制表符

paste默认连接符号是tab^w,  若需要更换为空格则 paste -d ‘  ’  file1 file2  即可。 paste的两个或多个文件需提前进行dos2unix的处理。  

2018-08-07 11:13:53

阅读数 140

评论数 0

ERA 处理成ascii

ERA的数据格式是我最喜欢的,读取速度够快。私人记录备忘,非教程。 clear all lat_max = 46; lat_min = 32; lon_max = 80; lon_min = 60; save_path = 'F:\AAA\interim\vic\'; Vname = 'tp'...

2018-08-04 17:39:43

阅读数 97

评论数 0

提示
确定要删除当前文章?
取消 删除
关闭
关闭