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

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