2 hydromonkey

尚未进行身份认证

分享让世界变得更美好

等级
TA的排名 11w+

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

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

2019-06-02 10:50:03

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

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

2019-05-09 15:44:49

新年开工

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

2019-02-12 16:15:41

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

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

处理modis数据最麻烦的就是投影问题,翻阅nasa网站后总算有些眉目。记录下来,以供大家参考。例子采用的AMSR雪当量月数据,这个数据采用的EASE-GRID投影方式,分辨率是25km*25km。整个北半球是721*721的pixels。(具体怎么投影成这样的请见官网)处理思路:把原文件换到wgs84坐标,接下来在matlab或其他就很好做矩阵运算了。首先hdf文件里雪当量是按照m...

2018-11-21 10:36:00

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

若是研究区在中国,则在地理空间数据云有中国区合成好的数据,直接下载用即可。数据格式为TIF。坐标为WGS84,不用再转换投影。下载网址为:http://www.gscloud.cn/sources/?cdataid=265&pdataid=10若研究区在国外,需从nasa下载modis标准数据产品,进行投影转换等工作。后续有时间会补充如何转换投影等内容,可参见《modisAMSR雪...

2018-11-09 16:57:32

考虑海拔的IDW的插值

importpandasaspdfromgeopy.distanceimportvincentyimportnumpyasnploworhigh='high'#输入需要哪段,low代表第一段,high代表第二段。year=1951#是否要选择年份dem=pd.read_table('F:/AAA/PG-T-LAPSE-COR/llde...

2018-10-02 18:22:23

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

importpandasaspdimportnumpyasnpfromscipy.spatial.distanceimportcdistimportscipy.statsasstatsfrompylabimport*STA_tmin='F:/AAA/STA-T-TMAX-TMIN-CSV/STA-TMEAN.csv'PG_tmin='F:/AAA...

2018-09-21 22:05:29

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

ARCGIS arcpy批量裁剪命令

 全部自动化批量处理,当然右键batch也可以处理,但是一旦有上百个文件,那就运动手指了。代码如下,注释都有,根据需要改动参数即可。 importarcpy#Setthecurrentworkspacearcpy.env.workspace=("F:/ZPATR-TIME-JOB/ndvi-china-/")#"D:\\Test2\\Grids.gdb"...

2018-09-09 18:40:36

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

fromcollectionsimportCounterCounter(a)a=[9,9,3,1,1,2]result=Counter(a)print(result)或者forxina:ifa.count(x)>1:print(x)# 

2018-09-04 17:45:35

狗十三

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

2018-08-25 18:07:22

CMORPH数据处理3

之前是先把二进制文件转成矩阵,再转成nc。写了很多代码,实在复杂。现在用CDO一句命令就可以直接转为nc。关于CDO的用法可参见官网。很强大的语句。我是在windows的cmd窗口运行CDO的,写了一个bat文件来批量执行任务。安装CDO等请参见气象家园。直接放出bat文件,从1998-2016,有需要的改成自己的目录即可使用。bat文件的具体制作过程不在这里说了。下载好,修改成...

2018-08-17 22:23:20

如何制作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

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

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

2018-08-16 23:05:37

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

共三种版本1:raw 2:CRT 3:BLD  多采用后两者(后面两个是经过站点校正的),下载网址如下:CRT–bias-corrected(over-estimatesinthetropics?) 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

GSMaP 降水数据下载指南

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

2018-08-16 16:39:42

modis AMSR 雪当量数据处理

查看更多文章请点击头像。Modis处理最麻烦的就是投影转换了,翻阅了下nasa的官网,记录下来,供大家参考。http://nsidc.org/data/AE_MoSno 数据页面。文件格式为25*25km的月雪当量数据,投影为EASE-GRIDnasa的投影方式。北半球一共721*721个pixel。注意雪当量以mm计,为了节约储存空间,hdf文件里的雪当量数据被除了一个2,后续处理注...

2018-08-08 22:41:06

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

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

2018-08-07 11:13:53

ERA 处理成ascii

ERA的数据格式是我最喜欢的,读取速度够快。私人记录备忘,非教程。clearalllat_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

查看更多

勋章 我的勋章
  • 持之以恒
    持之以恒
    授予每个自然月内发布4篇或4篇以上原创或翻译IT博文的用户。不积跬步无以至千里,不积小流无以成江海,程序人生的精彩需要坚持不懈地积累!
  • 勤写标兵Lv1
    勤写标兵Lv1
    授予每个自然周发布1篇到3篇原创IT博文的用户。本勋章将于次周周三上午根据用户上周的博文发布情况由系统自动颁发。