Month: October 2016

learn

How to download NDVI EVI form espa

弄清楚要研究区域的 Path/Row, 设置好研究时间段。在http://earthexplorer.usgs.gov/网站上把要下载的卫星影片筛选出来。 将筛选结果导出，单击“Click here to export your results”网站会给你发个邮件，在邮件中将链接打开，把要下载的文件名粘贴到“espa_request.txt”文件中。 打开并登陆https://espa.cr.usgs.gov/ordering/new/，在下图中单击“Choose File”并选择步骤3中的“espa_request.txt”文件，选择好其他参数。 单击submit订购，espa@usgs.gov会给你发两封邮件，一封是received，另外一封会慢点，因此需要点儿时间耐心等待。 用download_espa_order批量工具进行下载。

learn

Matlab中的类（Class）

对于类，我本没想搞明白，采用cftool工具拟合的结果返回值有两个fitresult和gof，其中fitresult是类，要保存它有些绕弯。只好读一些帮助理解的帖子，最后用coeffvalues(fitresult)把拟合系数变读取（访问）出来, 详见Matlab写入txt最后的例子，也许有更好的解决办法。将此好贴转载如下： 转载自：http://ycool.com/post/gxkf6q8 类的简介(class) 类是一种数据类型，与普通的数据类型不同的是类不仅包含数据，还包含对数据的操作，类把数据和数据操作方法封装在一起，作为一个整体参与程序的运行。类具有课继承性，创建一个新的类的时候，可以在一个基类中添加成员派生出新类。类的变量和类的实例是不同的，类的实例是动态分配的内存区域，通常称类的实例维“对象”，同一个类可以有不同的实例存在，他们个子有自己的数据，但是数据操作方法是相同的。类的变量可以看作是一个指针，指向类的实例。 在MATLAB中，为了更好地进行数据的封装，更加抽象地描述对象，也定义了类。图2.4.1给出了MATLAB中类的层次结构图，在建立类的结构的时候，可以根据这一层次结构进行设计和继承。 建立一个类 通常一个类应该包含四个基本的成员函数： 构造函数 ―― 与类名相同，可以在其中完成成员初始化的工作； 显示函数 ―― 名为display，用于显示成员的数据； 赋值函数 ―― 名为set，用于设置类成员的数值； 取值函数 ―― 名为get，用于读取类成员的函数。 […]

learn

Matlab写入txt

输入输出方式，即从数据文件读取数据或将结果写入数据文件。MATLAB提供了一系列低层输入输出函数，专门用于文件操作。 1）打开文件 在读写文件之前，必须先用fopen函数打开或创建文件，并指定对该文件进行的操作方式。fopen函数的调用格式为： fid=fopen（文件名，‘打开方式’） 说明：其中fid用于存储文件句柄值，如果返回的句柄值大于0，则说明文件打开成功。文件名用字符串形式，表示待打开的数据文件。常见的打开方式如下： ?      ‘r’：只读方式打开文件（默认的方式），该文件必须已存在。 ?      ‘r+’：读写方式打开文件，打开后先读后写。该文件必须已存在。 ?      ‘w’：打开后写入数据。该文件已存在则更新；不存在则创建。 ?      ‘w+’：读写方式打开文件。先读后写。该文件已存在则更新；不存在则创建。 ?      ‘a’：在打开的文件末端添加数据。文件不存在则创建。 ?      ‘a+’：打开文件后，先读入数据再添加数据。文件不存在则创建。 另外，在这些字符串后添加一个“t”，如‘rt’或‘wt+’，则将该文件以文本方式打开；如果添加的是“b”，则以二进制格式打开，这也是fopen函数默认的打开方式。 2）关闭文件 文件在进行完读、写等操作后，应及时关闭，以免数据丢失。关闭文件用fclose函数，调用格式为： […]

learn

去掉重复、寻找相同不再难

去掉重复、寻找相同这两件事在数据中会经常遇到，第一思路是通过比较循环来解决，要麻烦一些。但在Matlab中，把它们做为矩阵来处理，会很方便，并且效率也会高一些。   对于一个矩阵，去掉重复的行(no repetitions)，采用[C, ia, ic]=unique(A,’rows’); C——为A去除重复行后的矩阵； ia, ic 均为index，能够使A=C（ic），C=A（ia）成立。   对于A、B两个矩阵，若想找到A中和B相同的行，采用[Lia, Locb]=ismember(A,B,’rows’); Lia，Locb的维数 与矩阵A相同， Lia表示该行是否在B中存在，有了它就可以找到相同； Locb表示该行在B中的index； 用起来比较爽，具体看help。

travel

Big Falls Trail游记

计划要爬南加州最高的山，经准备，周日早起出发，早早到达san Bernardino Forest的Vivian停车场，因为没有permit，被管理人员劝阻。不过他给我们介绍了一个比较简单的Big Falls Trail，很是秀美。 停车场跨石河滩北行约500m，隔河岸，闻水声，如鸣佩环，心乐之。随行人，快走几步，下见小潭，水尤清冽。石砂以为底，近岸，怪石枯木嶙峋。为坻，为屿，为嵁，为岩。青树翠蔓，蒙络摇缀，参差披拂。寻源而上，斗折蛇行，若干瀑布或大或小，或上下阶梯状，或平行并列状，明灭可见。其岸势犬牙差互，不可知其源。 瀑布上下，有数名外国小孩磊石筑潭，不亦乐乎，有大人携小孩共游，有母抱婴携夫游者，虽语不甚通，但感觉相同，风景秀美，舒心惬意。 虽访山未至，却得玩水之乐。妻子共游，是以记之。 计划中的路线图石河滩瀑布1瀑布2枯木树根瀑布3瀑布4明灭可见小伙子也爬上来了

Uncategorized

关于潜热的理解Latent heat flux

以前只知道可以根据能量平衡的方法计算蒸发蒸腾量，“余项法”就是把潜热做为剩余项计算出。也知道，涡度相关计或者通量塔也可以测量潜热，由潜热就可以计算蒸发蒸腾量。那么潜热的单位是什么，和ET之间如何转换？ 潜热的单位一般为W/m2，表示单位面积上的功率，它与时间没有关系，一般每1小时或半小时测量一次，它每天的值就是每天测量n次的平均值。 Energy flux–Latent heat flux, Hl (W/m2) 蒸发就比较容易理解， Water flux–Evaporation rate, E (mm/day) 它们之间的关系为： Hl= ρlv E ρ= 1000 […]