Author: admin

travel

美西游记——第一天（20160720 Joshua Tree NP and LAS）

Riverside 雨季中，突生懒惰心。宅居，品茶，看完特朗普就职演说后，想整理下历时20天的第一次大规模美国旅游照片。从旅游第一天2016年7月20日到今天2017年1月20日，恰半年已过。遂以天为序，回忆记录这带夫人孩子一起在美西旅游的20天旅程。是以成《美西游记》，以为念。 第一天 Joshua tree NP 7月20日，暑假自驾游正式如期进行。我没有开过这么长时间的车，自己的腰也有毛病，只能带着一颗忐忑的心边开边体验了。出发当天临时决定对Joshua tree NP 南进北出。这是当天的路线图。 去公园的路上遇到大面积风力发电机，下了10号高速，不远就是公园名字的牌子。稍作休息，夫人准备的面包夹草莓酱，小朋友吃的很开心。这里到公园门口（公园南侧的visitor center）还有一段距离，由于人少车少，夫人替我开了一段车，让她感受下在美国开车的感觉，permit已经拿到。 在南Visitor center80刀买了一张国家公园年卡，服务人员态度蛮好，给小朋友送了一本书和铅笔，里面有各种问题和图画，小朋友在游览公园时，看到什么动物或者植物，可以在书上做上标记。等从公园出去时，可以在visitor center换一枚有纪念意义的勋章。进入公园，就一路大太阳，无处庇荫，隔一段路，就会有一个view point，有关于景点的介绍，来来往往的车辆在此停泊，下车拍照，但空气温度很高，没有会在外面停留太长时间。 进公园就是这种景色 这张照片很震撼，如果你有耐心放大来看。在太阳的炙烤下，还有蜥蜴类的动物从沙子轻巧地掠过。 不过，晚上在这个公园是欣赏星空拍星星的绝佳去处。路上还碰到了花园，供游人去走一走，近距离观察这种植物。这里会有因为缺水而疯狂的蜜蜂，漫山遍野的Joshua tree。蜜蜂警告Joshua […]

learn

调用拟合好的函数

Matlab 2016要调用一个拟合好的函数非常简单： f = fittype(‘a*x^2+b*exp(n*x)’)  %给出函数 f =      General model:        f(a,b,n,x) = a*x^2+b*exp(n*x) c = cfit(f,1,10.3,-1e2)         %给出相应参数，返回函数值 c […]

learn

[转载]matlab化简符号表达式

化简符号表达式计算机毕竟还是挺笨的, 经过一系列的符号计算后, 得到的结果可能只有它自己才能看懂， Matlab提供大量函数以用于符号表达式的化简. collect(f): 函数用途是合并多项式中相同的项, 如:syms x t f=(1+x)*t+x*t; collect(f) expand(f):展开多项式,syms x f=x*(x*(x-1)+3)+2; expand(f); horner(f) 对转换多项式为Horner形式, 这种形式的特点是乘法嵌套, 其有着不错的数值计算性质.syms x; […]

learn

长短矩阵的对齐

短矩阵向长矩阵对齐 A=[]; % 长系列序号 B=[]; % 短系列序号 Z=cell(1,1); % 和短系列同长度的矩阵，可包含汉字 [Lia, Locb]=ismember(A,B,’rows’); % 把短系列的矩阵向长系列看齐，缺的值以[]填写 for i=1:749;     if Locb(i)==0         […]

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。