function data = readdata(directory)
file = dir(directory); % 获取文件夹directory中的所有文件名
n = length(file);
j = 1;
for i = 1:n
point = strfind(file(i).name, '.'); % 寻找文件扩展名起始标志点'.'
if ~strcmp(file(i).name(point:end), '.txt') % 如果当前文件不是文本文件,则跳过该文件
continue;
end
dirname = [directory, '\', file(i).name]; % 将文件夹目录和文件名组集成当前文件的完整目录
mydata = importdata(dirname);
station =mydata(:,1);
weidu=mydata(:,2);
jingdu=mydata(:,3);
gaodu=mydata(:,4);
nian=mydata(:,5);
yue=mydata(:,6);
ri=mydata(:,7);
jiangshui=mydata(:,10);
%% 1961
jiangshui1961=jiangshui(nian == 1961 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1962
jiangshui1962=jiangshui(nian == 1962 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1963
jiangshui1963=jiangshui(nian == 1963 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1964
jiangshui1964=jiangshui(nian == 1964 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1965
jiangshui1965=jiangshui(nian == 1965 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1966
jiangshui1966=jiangshui(nian == 1966 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1967
jiangshui1967=jiangshui(nian == 1967 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1968
jiangshui1968=jiangshui(nian == 1968 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1969
jiangshui1969=jiangshui(nian == 1969 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1970
jiangshui1970=jiangshui(nian == 1970 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1971
jiangshui1971=jiangshui(nian == 1971 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1972
jiangshui1972=jiangshui(nian == 1972 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1973
jiangshui1973=jiangshui(nian == 1973 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1974
jiangshui1974=jiangshui(nian == 1974 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1975
jiangshui1975=jiangshui(nian == 1975 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1976
jiangshui1976=jiangshui(nian == 1976 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1977
jiangshui1977=jiangshui(nian == 1977 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1978
jiangshui1978=jiangshui(nian == 1978 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1979
jiangshui1979=jiangshui(nian == 1979 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1980
jiangshui1980=jiangshui(nian == 1980 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1981
jiangshui1981=jiangshui(nian == 1981 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1982
jiangshui1982=jiangshui(nian == 1982 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1983
jiangshui1983=jiangshui(nian == 1983 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1984
jiangshui1984=jiangshui(nian == 1984 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1985
jiangshui1985=jiangshui(nian == 1985 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1986;
jiangshui1986=jiangshui(nian == 1986 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1987
jiangshui1987=jiangshui(nian == 1987 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1988
jiangshui1988=jiangshui(nian == 1988 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1989
jiangshui1989=jiangshui(nian == 1989 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1990
jiangshui1990=jiangshui(nian == 1990 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1991
jiangshui1991=jiangshui(nian == 1991 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1992
jiangshui1992=jiangshui(nian == 1992 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1993
jiangshui1993=jiangshui(nian == 1993 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1994
jiangshui1994=jiangshui(nian == 1994 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1995
jiangshui1995=jiangshui(nian == 1995 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1996
jiangshui1996=jiangshui(nian == 1996 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1997
jiangshui1997=jiangshui(nian == 1997 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1998
jiangshui1998=jiangshui(nian == 1998 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 1999
jiangshui1999=jiangshui(nian == 1999 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2000
jiangshui2000=jiangshui(nian == 2000 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2001
jiangshui2001=jiangshui(nian == 2001 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2002
jiangshui2002=jiangshui(nian == 2002 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2003
jiangshui2003=jiangshui(nian == 2003 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2004
jiangshui2004=jiangshui(nian == 2004 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2005
jiangshui2005=jiangshui(nian == 2005 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2006
jiangshui2006=jiangshui(nian == 2006 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2007
jiangshui2007=jiangshui(nian == 2007 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2008
jiangshui2008=jiangshui(nian == 2008 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2009
jiangshui2009=jiangshui(nian == 2009 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2010
jiangshui2010=jiangshui(nian == 2010 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2011
jiangshui2011=jiangshui(nian == 2011 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2012
jiangshui2012=jiangshui(nian == 2012 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2013
jiangshui2013=jiangshui(nian == 2013 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 2014
jiangshui2014=jiangshui(nian == 2014 & yue == 1 & jiangshui >= 50 & jiangshui < 1000000000000);
%% 存储数据
data{j, 1} = file(i).name % 将文件名保存在元胞数组data的第1列
data{j, 2} = mean(station); % 将该的数据的第一列的平均值保存在元胞数组data的第2列
data{j, 3} = mean(weidu);
data{j, 4} = mean(jingdu);
data{j, 5} = mean(gaodu);
data{j, 6} = sum(jiangshui1961);
data{j, 7} = sum(jiangshui1962);
data{j, 8} = sum(jiangshui1963);
data{j, 9} = sum(jiangshui1964);
data{j, 10} = sum(jiangshui1965);
data{j, 11} = sum(jiangshui1966);
data{j, 12} = sum(jiangshui1967);
data{j, 13} = sum(jiangshui1968);
data{j, 14} = sum(jiangshui1969);
data{j, 15} = sum(jiangshui1970);
data{j, 16} = sum(jiangshui1971);
data{j, 17} = sum(jiangshui1972);
data{j, 18} = sum(jiangshui1973);
data{j, 19} = sum(jiangshui1974);
data{j, 20} = sum(jiangshui1975);
data{j, 21} = sum(jiangshui1976);
data{j, 22} = sum(jiangshui1977);
data{j, 23} = sum(jiangshui1978);
data{j, 24} = sum(jiangshui1979);
data{j, 25} = sum(jiangshui1980);
data{j, 26} = sum(jiangshui1981);
data{j, 27} = sum(jiangshui1982);
data{j, 28} = sum(jiangshui1983);
data{j, 29} = sum(jiangshui1984);
data{j, 30} = sum(jiangshui1985);
data{j, 31} = sum(jiangshui1986);
data{j, 32} = sum(jiangshui1987);
data{j, 33} = sum(jiangshui1988);
data{j, 34} = sum(jiangshui1989);
data{j, 35} = sum(jiangshui1990);
data{j, 36} = sum(jiangshui1991);
data{j, 37} = sum(jiangshui1992);
data{j, 38} = sum(jiangshui1993);
data{j, 39} = sum(jiangshui1994);
data{j, 40} = sum(jiangshui1995);
data{j, 41} = sum(jiangshui1996);
data{j, 42} = sum(jiangshui1997);
data{j, 43} = sum(jiangshui1998);
data{j, 44} = sum(jiangshui1999);
data{j, 45} = sum(jiangshui2000);
data{j, 46} = sum(jiangshui2001);
data{j, 47} = sum(jiangshui2002);
data{j, 48} = sum(jiangshui2003);
data{j, 49} = sum(jiangshui2004);
data{j, 50} = sum(jiangshui2005);
data{j, 51} = sum(jiangshui2006);
data{j, 52} = sum(jiangshui2007);
data{j, 53} = sum(jiangshui2008);
data{j, 54} = sum(jiangshui2009);
data{j, 55} = sum(jiangshui2010);
data{j, 56} = sum(jiangshui2011);
data{j, 57} = sum(jiangshui2012);
data{j, 58} = sum(jiangshui2013);
data{j, 59} = sum(jiangshui2014);
data{j, 60} = numel(jiangshui1961);
data{j, 61} = numel(jiangshui1962);
data{j, 62} = numel(jiangshui1963);
data{j, 63} = numel(jiangshui1964);
data{j, 64} = numel(jiangshui1965);
data{j, 65} = numel(jiangshui1966);
data{j, 66} = numel(jiangshui1967);
data{j, 67} = numel(jiangshui1968);
data{j, 68} = numel(jiangshui1969);
data{j, 69} = numel(jiangshui1970);
data{j, 70} = numel(jiangshui1971);
data{j, 71} = numel(jiangshui1972);
data{j, 72} = numel(jiangshui1973);
data{j, 73} = numel(jiangshui1974);
data{j, 74} = numel(jiangshui1975);
data{j, 75} = numel(jiangshui1976);
data{j, 76} = numel(jiangshui1977);
data{j, 77} = numel(jiangshui1978);
data{j, 78} = numel(jiangshui1979);
data{j, 79} = numel(jiangshui1980);
data{j, 80} = numel(jiangshui1981);
data{j, 81} = numel(jiangshui1982);
data{j, 82} = numel(jiangshui1983);
data{j, 83} = numel(jiangshui1984);
data{j, 84} = numel(jiangshui1985);
data{j, 85} = numel(jiangshui1986);
data{j, 86} = numel(jiangshui1987);
data{j, 87} = numel(jiangshui1988);
data{j, 88} = numel(jiangshui1989);
data{j, 89} = numel(jiangshui1990);
data{j, 90} = numel(jiangshui1991);
data{j, 91} = numel(jiangshui1992);
data{j, 92} = numel(jiangshui1993);
data{j, 93} = numel(jiangshui1994);
data{j, 94} = numel(jiangshui1995);
data{j, 95} = numel(jiangshui1996);
data{j, 96} = numel(jiangshui1997);
data{j, 97} = numel(jiangshui1998);
data{j, 98} = numel(jiangshui1999);
data{j, 99} = numel(jiangshui2000);
data{j, 100} = numel(jiangshui2001);
data{j, 101} = numel(jiangshui2002);
data{j, 102} = numel(jiangshui2003);
data{j, 103} = numel(jiangshui2004);
data{j, 104} = numel(jiangshui2005);
data{j, 105} = numel(jiangshui2006);
data{j, 106} = numel(jiangshui2007);
data{j, 107} = numel(jiangshui2008);
data{j, 108} = numel(jiangshui2009);
data{j, 109} = numel(jiangshui2010);
data{j, 110} = numel(jiangshui2011);
data{j, 111} = numel(jiangshui2012);
data{j, 112} = numel(jiangshui2013);
data{j, 113} = numel(jiangshui2014);
j = j+1;
end