Matlab写入txt

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

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

去掉重复、寻找相同这两件事在数据中会经常遇到,第一思路是通过比较循环来解决,要麻烦一些。但在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。

批量读取文件名read file name in batch

目前总结出3种方法对卫片文件名进行整理: 1.用excel合成; 2.在EarthExplorer网站等搜索引擎中 export your results,将搜索后的文件名导出成csv文件; 3.用批处理命令读取文件名并写进zxd.txt文件中,保存成bat文件运行,或者直接用cmd运行: dir /b J:\Anderson_lab\for_Xudong\Landsat8_surface_reflectance\*band1*.tif >> E:\zxd.txt

在USGS批量下载卫片download_espa_order

1.在USGS网上预订(http://earthexplorer.usgs.gov/),这个需要一段时间大概1-2天,预订好后会通过邮件给你发一个链接,不过要打开链接,还是需要登录用户名。 2.使用基于python的批下载工具 Bulk Download Client,下载 Bulk Download Client的链接为: https://github.com/USGS-EROS/espa-bulk-downloader/blob/master/espa_bulk_downloader_v2.0.1.zip 3.解压缩文件夹 4.将解压缩后的文件都copy到python27下 5.仔细阅读README,里面有download_espa_order.py的调用方式和详细的使用说明,举例实例如下,在CMD下输入: C:\ python27\ python download_espa_order.py -e zxxddd58@yahoo.com -o zxxddd58@yahoo.com-0101609202926 […]

几个比较好的网站nice website

1.一个关于下水道历史的网站:http://www.sewerhistory.org 介绍各国下水道的发展历程和历史文化。 2.北方水网站:http://www.northernwater.org/Default.aspx 介绍科罗拉多河流域的各种水工程项目,节水措施等。 3.加州灌溉管理信息系统:http://wwwcimis.water.ca.gov/SpatialData.aspx 有各种数据及报告。 4.美国土壤类型统计:http://casoilresource.lawr.ucdavis.edu/gmap/ 美国的土壤类型分布数据。 5.摄影入门:http://www.360doc.com/content/14/0405/22/358457_366661735.shtml 6.查询高程:http://elevationmap.net/ 7.作物种类:https://nassgeodata.gmu.edu/CropScape/ 8.ETcb,Kcb,Fc数据,:http://ec2-54-197-48-121.compute-1.amazonaws.com/dgw/sims/ 9.METRIC ET data: http://eeflux-level1.appspot.com/ 10.中华珍宝馆:http://ltfc.net/ 11.灾害统计网:http://www.gdacs.org/   webinars […]

Matlab read txt in batch processing(Matlab批量读取txt文件)

假如有636个txt文件,需要读取其中的特定的行或者列。文件名如下: 可以先读取一个正确的,让matlab自动生产读取的位置和m文件;然后将所需要读取的文件名字和循环次数(所需要读取的文件个数)加到m文件中并修改之。 最终,Matlab 批量读取txt数据的程序为: clear; clc; dir =’C:\In\’; fyj=[196101,196201,196301,196401,196501,196601,196701,196801,196901,197001,197101,……]; zhr=reshape(fyj,636,1); for i=1:636 ii=zhr(i); FileNames1=’SURF_CLI_CHN_PRE_MON_GRID_0.5-‘; FileNames2=num2str(ii); FileNames3=’.txt’; %% Initialize variables. filename […]