FY-4B LST原始数据预处理及质量控制

通过Python代码实现

import os
import time
import numpy as np
import netCDF4 as nc
from scipy.interpolate import griddata
import rasterio
from rasterio.crs import CRS
from tqdm import tqdm
from concurrent.futures import ThreadPoolExecutor, as_completed
from datetime import datetime

# =============================================================================
# 配置区
# =============================================================================
BASE_INPUT_FOLDER = r\”H:\\FY-4B\\Raw_Data\” # 原始数据根目录
BASE_OUTPUT_FOLDER = r\”H:\\FY-4B\\Data_逐时_tif_clip\” # 输出根目录
NUM_THREADS = 4 # 线程数（可改6/8加速）

# 固定输出参数（研究区范围裁剪）
OUT_RES = 0.04
TARGET_ROWS = 901
TARGET_COLS = 1601
BASE_MINX = 72.0
BASE_MAXY = 54.0

# 经纬度表路径
RAW_1330E = r\”H:\\FY-4B\\Raw_Data\\FY4B-_DISK_1330E_GEO_NOM_LUT_20220323000000_4000M_V0001.raw\”
RAW_105E = r\”H:\\FY-4B\\Raw_Data\\FY4B-_DISK_1050E_GEO_NOM_LUT_20240227000000_4000M_V0001.raw\”

# 需要处理的月份列表