ToB企服应用市场:ToB评测及商务社交产业平台

标题: 将NC栅格表示时间维度的数据提取出来的方法 [打印本页]
作者: 徐锦洪    时间: 2024-10-22 15:51
标题: 将NC栅格表示时间维度的数据提取出来的方法
  本文先容基于Python语言，逐一读取大量.nc格式的多时相栅格文件，导出此中所具有的全部时间信息的方法。
  .nc是NetCDF（Network Common Data Form）文件的扩展名，表示一种常用的科学数据存储格式。NetCDF是一种自形貌的、可移植的二进制文件格式，用于存储科学和工程范畴的大型数据集；由于其自身的特性，.nc数据被广泛应用于景象学、海洋学、地球科学、气候研究、大气科学、地理信息系统等范畴。
  首先，明确一下本文的需求。现在有一个文件夹，此中具有大量的.nc格式的栅格文件，如下图所示。

  此中，每一个.nc格式的文件都具有多个时相（或者说是多个维度），而不仅仅只是一个时相。我们希望，读取这个文件夹中的全部.nc格式文件，并获取其所表示的每一个时相。
  明确了需求后，我们就可以开始具体的操作。首先，本文所需用到的代码如下。
  1. # -*- coding: utf-8 -*-
  2. """
  3. Created on Sun Dec 31 20:28:03 2023
  5. @author: fkxxgis
  6. """
  8. import os
  9. import netCDF4
  10. from netCDF4 import Dataset
  12. def list_nc_dates(folder_path):
  13.     nc_dates = []
  15.     for file_name in os.listdir(folder_path):
  16.         if file_name.endswith(".nc"):
  17.             file_path = os.path.join(folder_path, file_name)
  18.             try:
  19.                 dataset = Dataset(file_path)
  20.                 time_var = dataset.variables["time"]
  21.                 time_values = time_var[:]
  22.                 time_units = time_var.units
  23.                 time_calendar = time_var.calendar
  25.                 dates = []
  26.                 for value in time_values:
  27.                     date = netCDF4.num2date(value, units=time_units, calendar=time_calendar)
  28.                     dates.append(date.strftime("%Y-%m-%d %H:%M:%S"))
  30.                 nc_dates.append((file_name, dates))
  31.             except Exception as e:
  32.                 print(f"Error reading file {file_name}: {str(e)}")
  34.     return nc_dates
  36. folder_path = "F:/Data_Reflectance_Rec/soil_1"
  37. nc_dates = list_nc_dates(folder_path)
  39. for nc_file, dates in nc_dates:
  40.     for date in dates:
  41.         print(date)
复制代码
  这段代码团体思路也很明确。
  首先，我们导入所需的模块。在这里，需要导入Python的os模块，用于处理文件和文件夹路径操作；同时导入netCDF4库，并接着从netCDF4库中导入Dataset类，用于打开和读取.nc文件。在这里，如果需要设置netCDF4库，大家可以参考文章设置h5py、netCDF4库的方法：Anaconda环境。
  接下来，我们定义了一个名为list_nc_dates的函数，继承一个文件夹路径作为参数。在函数中，首先创建一个空列表nc_dates，用于存储每个.nc文件及其对应的日期列表；随后，使用os.listdir()函数遍历文件夹中的所有文件，通过查抄文件名是否以.nc末了来筛选出.nc文件。紧接着，对于筛选出来的.nc文件，使用os.path.join()函数构建其完整路径。
  其次，使用Dataset类打开.nc文件，并将打开的文件对象赋值给dataset变量；随后，获取.nc文件的时间，在本文的.nc数据中，也就是名为time的变量，并将时间变量的值读取到time_values变量中。接下来，分别获取时间变量的单元与时间范例。
  随后，我们创建一个空列表dates，用于存储日期字符串。遍历时间变量的每个值，使用netCDF4.num2date()函数将时间值转换为日期对象。紧接着，将日期对象转换为指定格式的字符串，并将其添加到dates列表中。此外，这里还将.nc文件名和对应的日期列表作为元组添加到nc_dates列表中，方便我们后期对日期的核对。函数的最后，返回包含每个.nc文件及其对应日期的列表。
  在函数外部，我们设置文件夹路径，随后即可调用list_nc_dates函数，将文件夹路径传递给它，并将返回的效果赋值给nc_dates变量。最后，通过循环，打印每个日期即可。
  执行上述代码，即可出现如下图所示的效果（效果很长，就截取一部分）。由于在本文中，每一个.nc格式文件的每一个维度（即每一个时相）都是准确到天的，所以下图天数后的时、分、秒都是00。当然，如果大家的.nc格式文件维度很多，时相打印出来的话也欠好完全表现，所以可以考虑将时间信息导出为表格文件等；例如，可以将每一个date都放在DataFrame中，随后导出为.csv文件。

  至此，大功告成。

免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！更多信息从访问主页：qidao123.com:ToB企服之家，中国第一个企服评测及商务社交产业平台。



欢迎光临 ToB企服应用市场:ToB评测及商务社交产业平台 (https://dis.qidao123.com/) Powered by Discuz! X3.4