[已解决]matlab代码

聆風祈

求助大佬，想知道这个降水grd文件中数据是怎样的，比如是站点日降水还是流域年降水之类的。如果是日降水怎样汇聚成年降水并输出到excel中方便和径流做相关分析。

最佳答案

月排行榜 / 总排行榜

isdkz

2023-4-17 18:59:54

聆風祈 发表于 2023-4-17 17:31
大佬这个处理完出现了9216行数据，是72个纬度一组经度还是128个经度一组纬度啊，然后他那个描述是从多少 ...

抱歉之前的代码输出格式可能会让您感到困惑。在这个情况下，输出的数据是按照128个经度一组纬度排列的。也就是说，前128行对应于第一个纬度（18度），接下来的128行对应于第二个纬度（18.5度），依此类推。按照ctl文件中的描述，数据从纬度18度、经度72度开始，每个纬度间隔0.5度，每个经度间隔0.5度。

要使输出的数据更易于阅读和理解，我们可以在导出到Excel之前调整其结构。下面是一个修改后的代码示例：

1. % 输出到Excel
   
2. xlswrite('annual_precipitation.xlsx', permute(annual_precipitation, [2, 1, 3]));

复制代码

这个代码示例使用permute函数交换了经度和纬度的顺序，这样输出的数据就会按照72个纬度一组经度排列。也就是说，前72行对应于第一个经度（72度），接下来的72行对应于第二个经度（72.5度），依此类推。这种方式可能更符合直观的理解。

跳转到最佳答案楼层

聆風祈

*dset D:\tempt\11.grd
dset d:/wind/precn-1961-2020.grd
title Sea Surface Temperature Analysis
undef -9999
ydef 72 linear 18 0.5
xdef 128 linear 72 0.5
zdef 1 linear 0 1
tdef 720 linear 00Z01JAN1961 1mo
vars 1
pre=>pre  0  t,y,x rainfall
endvars
ctl文件里面是这样的

歌者文明清理员

根据提供的信息，降水数据是以格点形式存在的，并且时间分辨率是1个月。因此，该数据表示的是地表降水场，覆盖整个研究区域，不是局部站点数据。如果您需要计算流域年降水量，可以先将每个格点的日降水加总成月降水数据，然后再将月降水加总成年降水。最后，您可以将结果输出到Excel中进行相关分析。

isdkz

根据您提供的ctl文件，这是一个降水数据的网格文件，包括1961年到2020年每月的降水数据。文件包含了72行（纬度）和128列（经度）的网格数据，每个网格单元的大小分别为0.5度。现在，为了计算日降水到年降水的累积并将结果导出到Excel，我们需要首先解析这个grd文件，然后对数据进行处理。下面是一个简单的MATLAB代码示例：

1. % 读取grd文件
   
2. filename = 'D:\tempt\11.grd';
   
3. fileID = fopen(filename, 'r');
   
4. data = fread(fileID, [128*72, 720], 'float32');
   
5. fclose(fileID);
   
6. 
   
7. % 转换数据维度
   
8. data = reshape(data, [128, 72, 720]);
   
9. 
   
10. % 计算每年降水总量
    
11. year_num = 2020 - 1961 + 1;
    
12. annual_precipitation = zeros(128, 72, year_num);
    
13. 
    
14. for i = 1:year_num
    
15.     start_month = (i-1) * 12 + 1;
    
16.     end_month = i * 12;
    
17.     annual_precipitation(:, :, i) = sum(data(:, :, start_month:end_month), 3);
    
18. end
    
19. 
    
20. % 输出到Excel
    
21. xlswrite('annual_precipitation.xlsx', reshape(annual_precipitation, [128*72, year_num]));
    

复制代码

注意：这个代码示例假设了每个月的降水量是一个月的累积降水量。根据您的数据源，您可能需要进行适当的调整以确保结果的正确性。这个示例仅用于展示如何进行计算和导出。

首先，代码读取了grd文件，然后将数据重新整形为128x72x720的矩阵。接着，代码计算每年的降水总量并将结果存储在一个新的矩阵中。最后，使用xlswrite函数将结果写入到一个Excel文件中，以便进行后续的相关分析。

聆風祈

isdkz 发表于 2023-4-17 09:38
根据您提供的ctl文件，这是一个降水数据的网格文件，包括1961年到2020年每月的降水数据。文件包含了72行（ ...

大佬这个处理完出现了9216行数据，是72个纬度一组经度还是128个经度一组纬度啊，然后他那个描述是从多少经纬度开始啊

isdkz

聆風祈 发表于 2023-4-17 17:31
大佬这个处理完出现了9216行数据，是72个纬度一组经度还是128个经度一组纬度啊，然后他那个描述是从多少 ...

抱歉之前的代码输出格式可能会让您感到困惑。在这个情况下，输出的数据是按照128个经度一组纬度排列的。也就是说，前128行对应于第一个纬度（18度），接下来的128行对应于第二个纬度（18.5度），依此类推。按照ctl文件中的描述，数据从纬度18度、经度72度开始，每个纬度间隔0.5度，每个经度间隔0.5度。

要使输出的数据更易于阅读和理解，我们可以在导出到Excel之前调整其结构。下面是一个修改后的代码示例：

1. % 输出到Excel
   
2. xlswrite('annual_precipitation.xlsx', permute(annual_precipitation, [2, 1, 3]));

复制代码

这个代码示例使用permute函数交换了经度和纬度的顺序，这样输出的数据就会按照72个纬度一组经度排列。也就是说，前72行对应于第一个经度（72度），接下来的72行对应于第二个经度（72.5度），依此类推。这种方式可能更符合直观的理解。

聆風祈

isdkz 发表于 2023-4-17 18:59
抱歉之前的代码输出格式可能会让您感到困惑。在这个情况下，输出的数据是按照128个经度一组纬度排列的。 ...

感谢大佬