NetCDF文件合并

MeteoInfo

哞哞三年半 发表于 2013-3-26 16:41

                               
登录/注册后可看大图

经纬度是北半球，共1981-2000共二十年的夏天的hgt,现在想把这20年合并在一起

是Surfer的二进制格点数据格式吗？

哞哞三年半

MeteoInfo 发表于 2013-3-26 16:51

                               
登录/注册后可看大图

是Surfer的二进制格点数据格式吗？

亲！忽略surfer，就是普通的二进制格点数据，配个ctl能在grads中画图的那个~~~~~~~~~~~~

哞哞三年半

MeteoInfo 发表于 2013-3-26 16:51

                               
登录/注册后可看大图

是Surfer的二进制格点数据格式吗？

所以应该怎么合并呢 觉得用fortran好麻烦~

MeteoInfo

哞哞三年半 发表于 2013-3-26 17:31

                               
登录/注册后可看大图

所以应该怎么合并呢 觉得用fortran好麻烦~

如果每个文件里只有一个格点数据（只有一个变量、一个时次、一个高度），可以简单的用DOS的copy命令合并，比如：

copy /b *.dat newdata.dat

当然ctl文件需要根据合并后的数据情况重新编写。

如果每个文件里的格点数据比较复杂（比如多变量、多时次、多高度），则需要写程序来合并数据，以保证合并后的数据符合GrADS二进制数据格式。

royalosyin

合并NC文件，如果是二次开发，可以使用MeteoInfoC中的函数查看要合并的nc文件信息，然后使用nco中的命令做实际工作（http://nco.sourceforge.net/）。
上个函数，贴点代码。 这里代码的功能是将新产生的小文件不停地append已有的文件后面。        /// fileout = fileout + filein;
        private void UseNCOToConcatenateFiles(string filein, string fileout)
        {
            // Prepare the process to run
            ProcessStartInfo start = new ProcessStartInfo();

            // Enter in the command line arguments, everything you would enter after the executable name itself
            StringBuilder sb = new StringBuilder();
            sb.Append(" -O --no_tmp_fl ");
            sb.Append(fileout);
            sb.Append(" ");
            sb.Append(filein);
            sb.Append(" -o ");
            sb.Append("test.nc");

            start.Arguments = sb.ToString();
            // Enter the executable to run, including the complete path
            start.FileName = @"nco\ncrcat.exe ";

            // Do you want to show a console window?
            start.WindowStyle = ProcessWindowStyle.Hidden;
            start.CreateNoWindow = true;

            // Run the external process & wait for it to finish
            using (Process proc = Process.Start(start))
            {
                this.Cursor = Cursors.WaitCursor;

                proc.WaitForExit();

                // Retrieve the app's exit code
                int exitCode = proc.ExitCode;

                if (exitCode != 0)
                {
                    string errmessage = "There occured an error.\n Please try again.";
                    MessageBox.Show(errmessage, "Alert");
                }
                else
                {
                    #region Using ncrcat, have to use the following extra file operations.
                    //Sometimes, they are very dangerous. So it'd better backing up firstly.
                    File.Copy(fileout, fileout + ".backup.nc", true);
                    File.Copy("test.nc", fileout, true);
                    File.Delete("test.nc");
                    #endregion End of extra file operations.

                    string message = Path.GetFileName(filein) + " has been concatenated to " + Path.GetFileName(fileout) + " successfully!";
                    MessageBox.Show(message, "Info");
                }

                this.Cursor = Cursors.Default;
            }//end of using
        }



合并前，要用到MeteoInfoC中的nc函数查看是否已经合并过等等。
MeteoInfoC中的函数我就不说了，自己看文档吧

royalosyin

我用nco+MeteoInfoC合并了两个700多M的文件，没有崩溃

MeteoInfo

royalosyin 发表于 2013-3-28 04:55

                               
登录/注册后可看大图

我用nco+MeteoInfoC合并了两个700多M的文件，没有崩溃

用MeteoInfo脚本也很方便：http://bbs.06climate.com/forum.p ... 298&fromuid=106

曾记否

每天学习一点，感谢分享！

a1s1

“文件合并分为两种情况，一种是合并不同时次的数据文件，要求各文件的维、变量等设置都一样，只是时间不同；”请问这种情况用MeteoInfo如何操作合并啊？
刚才用文中方法合并后发现只有第一个文件的数据保留着，其它的文件数据都没有了。

a1s1

a1s1 发表于 2014-1-2 08:26

                               
登录/注册后可看大图

“文件合并分为两种情况，一种是合并不同时次的数据文件，要求各文件的维、变量等设置都一样，只是时间不同 ...

哈，更新了文件，合并没问题了