有效降水

水文学的有效降水:http://www.alanasmith.com/theory-Calculating-Effective-Rainfall.htm农田水利学的有效降水:http://www.fao.org/3/S2022E/s2022e03.htm laiEffective rainfall in irrigated agriculture;http://www.fao.org/3/X5560E/x5560e00.htm#Contents Rainfall, evaporation and effective rainfall:https://www.tocal.nsw.edu.au/tocals-new-e-farm/the-climate-of-tocal/rainfall,-evaporation-and-effective-rainfall 关于有效降水量的计算: 在CropWat中,有效降水的计算若选用USDA土壤保持局(USDA Soil Conservation Service)的方法,则可以分为月尺度计算和旬尺度计算,计算公式如下: Monthly step: Peff = […]

2019年春节——驾车返乡

昨日到达保定,今天(2019年2月1日)我们去了冉庄地道战纪念馆,因纪念馆闭馆修缮,我们没有参观,略显遗憾。但去参观了地道,地道很长,而且比较矮。若个子高的人进去的话容易碰到头,我家小孩还挺美的,走的很自如。我就走的难受了,哈腰驼背的。从保定去冉庄开车,大约50分钟。我们没有走高速,开了约1个小时。道路状况不是特别好,去的人也不多。他自己有一个专门儿的停车场。每辆校车收费五元。 下午去了清代省府第一衙————保定直隶总督署,每个人门票30元,觉得值得一看。直隶总督署是清代直隶总督的办公处所,占地30000余平米,分东中西三路,主体建筑集中在中路,是一座典型的北方衙署建筑群。该属历经清朝八帝182年历史,史料非常丰富,是目前我国保存最完整的清代省级衙署。在这里主要了解了做官的文化。里面保留了大量黄清代皇帝如道光、同治、乾隆、慈禧太后等的匾额。还有清代著名大臣、张之洞、李鸿章、曾国藩所题写的匾额。感觉挺有意思的。还包括一些做官的规则,对官员的要求,如何防止官员腐败,上任和离任的程序,官员的等级,各级官员及其家属穿衣的规定。历史上所产生的名人,永定河泛滥的历史,以及历代官员对永定河的治理。精正廉明 皇家御笔 战战兢兢即生时不忘地狱,坦坦荡荡虽逆境不忘天怀。曾国藩 曾国藩书法作品欣赏 http://www.sohu.com/a/148325686_224772

《通天大道宽又阔》改编

《通天大道宽又阔》改编 刚写完了几个字, 又做完了几道题。 模拟测试怎么就这么多! (白:俺老孙干它一场!) 整他个语也百分数也百。 你也开心,我也高兴, 考得那高分大学任我挑。 刚看过了几本书, 又上过了几门课, 作业上课怎么它就这么多。 (白:俺老孙去也!) 去你个眼要花来背要驼, 难也遇过,苦也吃过, 走出个人生大道宽又阔。 给二年级的儿子改编首歌曲,可惜自己不能唱,有人给唱了录音就好了。

徐总的插值工具

如何将点坐标的值插值成某一分辨率的asc文件? 假如是excel形式的文件,将excel的表转换成shape文件; 在Arcgis下,打开shape文件的属性表,添加X、Y字段; 右键X字段,选择calculate geometry(计算几何),Y字段同样,注意单位选择米。 导出成txt文件。 把X、Y和对应的值写入txt,用到了如下的matlab程序,以每年12个月为例: clear all; close all; clc; for mon=1:12; %表示月 %% Import the data % […]

郭总名言

1.不怕当最后一名,只要能够进入那一层次; 2.处处微妙,便无微妙; 3.争争抢抢不够吃,推推让让吃不够; 4.尚无立足之地,何谈栋梁之才; 5.原则有大小,事情无大小; 6.买一本30-50页的图集,只要有一页图纸能用,就已够本; 7.夫妻两个人要尿到一个壶里; 8.不要迷信他人,要有自己独立的思考; 9.不听老人言,吃亏在眼前,全听老人言,做不了神仙。

位山灌区初认识

2018年11月23日,第一次来到山东聊城的位山灌区。感触颇深,写小文以记之。 取水口位于长江水和黄河水的交叉点,引黄第二大灌区,引水闸坐落在泰山岩石基础上,三大干渠自流灌溉——位山灌区。引黄河水,水高沙高,水来沙来,水走沙留,调水调沙,沙走人走,相生相克,既恨又爱。60载工程,河道变迁,或淤或冲,曾经断流,曾经贫穷,该输沙的输沙,该沉沙的沉沙,该灌溉的灌溉,它依然在那里,默默地镌刻着这片土地。如今,新时代,新理念,重生态,爱环境,堆沙为山,化沙为宝,沉沙池变成湿地公园,传统灌区成为风景区,高效管理智慧决策,黄河逐渐褪去浑浊,河床不再升高,震旦雅雀悉数来归,人民安居乐业,经济迅速发展,一切都在向好之中。 引黄灌区具有不同的特征,上游的引黄灌区,如甘肃、宁夏,一般情况下是大引大排;中游的引黄灌区,如河道灌区,一般为中引小排;而下游的引黄灌区则是小引不排。这样也许与它们的降水量和地下水位有关,甘肃的降水量少,只有排水才能把盐分带走,河套灌区的降水量,所以也需要一定的排水量,而位山灌区则因为本身有一定的降水量,大约为500-600mm,这样本身起到了洗盐排盐的作用。一般认为,地下水埋深大于2m则不容易会产生盐碱化。 位山灌区引水闸位置闸前黄河泥沙西东输沙渠原来的沉砂池,如今的湿地公园聊城的早晨,雾霭缭绕二干渠穿过市区 位山灌区降水比较集中,7-9月降水500m。引黄河水量5.6亿m3,灌溉面积515万亩。地下水用的少。灌区分级管理。274km干渠。地下水埋深在3-4m,井渠配套。地势高差十多米。深县多蔬菜大棚,也有个引黄灌区,引水口在河南濮阳。主要作物为小麦,玉米,有果树。 截止2018年,秸秆还田已经实施了5-6年。但秸还田以后,土壤有变酸的趋势,小前期扎根好,后期须根多,不容易向下扎根。因此在秸秆还田之后,需要冬灌,密实土壤,让小麦向土壤深处扎根,以防止后期倒伏,冬灌定额50-60方/亩,然而在秸秆还田实施之前是不需要冬灌的。灌浆水100多方/亩。产量为小麦1100-1200斤,玉米1400斤。田间灌水方式为小白龙。灌区没有退水,所以存在潜在盐碱化的风险,目前地下水位深,盐碱化明显。 在干渠只衬砌边坡而渠底没有衬砌的地方,渠道的渗漏对底线水的补给作用明显。比如,在干渠周围的同一块儿地里,如果采用井水灌溉,如果干渠有水是灌溉每亩地需要40分钟,若干渠没水时则需要1.5小时。在穿过城市的干渠段,因为景观用水的需求,渠道底部也进行了衬砌,但是如果不再放水更新,大约1周左右,渠道的水也就见底了。位山灌区灌溉试验站36.385085, 116.132421,位于一干上,面积越40亩租用耕地,当地老百姓比较淳朴。 位山二干渠36.647907, 116.049960高营节制闸,清华大学在此地修建了通量塔,已经进行了10多年的水汽通量观测。 总之,一切都在向好之中,感谢调研组和陪同讲解人员,祝愿灌区,祝愿我国的水利事业。    

用矩阵散点做等高线图和梯度场图

matlab提供了可以根据矩阵做等高线的命令contour,和求梯度的函数gradient,可以根据梯度用quiver命令将梯度做成箭头图。但要注意一些参数的设置。 一般情况下,例如: x=-2:0.2:2; y=-2:0.2:2; [X,Y]=meshgrid(x,y); Z=X.^3+Y.^2; contour(X,Y,Z); [DX,DY]=gradient(Z,0.2,0.2); hold on; quiver(X,Y,DX,DY); 作图如下: 设置一下不同的坐标,纵坐标的范围为 [-20 20],间距为2。 x=-2:0.2:2; y=-20:2:20; [X,Y] = meshgrid(x,y); […]

位图格式图片如何转化为Arcgis下的shape文件

本文以jpg格式的图片为例,jpg格式是位图,它没有坐标系,但有时候又非常有必要把它做成shape文件。本文介绍两种方法,一种是通过以下6个步骤来进行的,另外一种简单一些,只有2个步骤。 1.把jpg图片放到有坐标的Arcgis里面,采用Georeferencing工具进行配准。 2.借助矢量化工具比如R2V等进行矢量化,将矢量化后的结果在CAD中对所有的线进行分层。 3.在步骤1里面选取两个特征点A,B并查询它们的坐标。 4.在CAD文件中输入A,B两个参照点的坐标,用AL对齐命令将整个目标图形旋转并缩放到正确的位置。 5.将目标图形加载到Arcgis里面,用export把相应的面或者线层导出并另存为shape文件。 注意:CAD目标文件要干净,没有多余的信息,最好分层做。采用Arcgis10以上的版本,9.3版本的有可能出问题。 6.将shape文件在catalog里面添加坐标和投影系统。就可大功告成。 另外一种简单的方法是,将jpg图片在Arcgis下配准以后,直接新建面对象,采用编辑面对象的工具Editor直接绘制边界。  

散点图矩阵(Scatter plot matrix)和散布矩阵(Scatter matrix)

1)散点图矩阵对于matlab来说很容易实现,采用plotmatrix实现就行。参看帮助: x = randn(50,3); y = x*[-1 2 1;2 0 1;1 -2 3;]’; plotmatrix(y,’*r’) 2) 散布矩阵有用python的Matplotlib库来实现的,matlab实现起来有点儿小麻烦,不过可以借助scatter命令来进行一些设置也能完美实现,现在以一个17*17的实对称矩阵来说,程序如下: clear all figure1=figure; % […]