Data N' Method

Himawari-8静止卫星AHI云参数数据的下载

© 2023-2024 Jiheng HU, 原创内容，禁止转载。 要获取GOES系列静止卫星的L2+级云参数产品的FullDisk圆盘数据，需要清楚相关卫星的工作方式。 NOAA 第五代GOES静止卫星星座GOES-R系列静止卫星是目前NOAA在运行的第三代静止卫星任务，其编号分别为GOES-R/S/T三颗卫星。其后继的第五代静止卫星为U系列，预计于2024年以后发射。关于具体的GOES任务表，详见WIKIPEDIA.GOES ScheduleGOES静止卫星星座一般由两颗位于east和west点位上的业务卫星构成，赤经分别为75.2° W 和 137.2°W。目前GOES-R系列已经发射了三颗业务卫星，分别是GOES16、GOES17和GOES18，其发射前的名称分别为GOES-R/GOES-S/GOES-T。 GOES-16卫星是R系列的第一颗卫星，其与2017年12月18日机动到EAST点，取代原先的GOES-13星并开始业务化运行。自此，GOES-16一直是GOES-R系列唯一一颗EAST卫星。 GOES-17卫星于201 ...

MSG静止卫星云掩膜数据的下载

本文将包括以下解决方案： 如何自动化脚本下载？ 如何自定义ERA5再分析数据的水平分辨率？ ERA5再分析数据是大气科学领域常用的几种再分析数据之一，记得在3~5年以前，ERA5还不能普遍获取，只有新世纪年代的在分析数据集。目前，最新的ERA5数据集已经回溯到了1970年以前，无论是对于数值模拟还是年代际的气候变化研究来说，ERA5数据集都一个非常好用的数据集。相较于之前的ERA再分析数据集来说，ERA5有了更高效、更快速的数据分发系统，即Climate Data Store ，简称CDS。下载速度已经达到了MB级，会想起以前下载ERA-Interim/20C的痛苦回忆，不甚唏嘘。 准备工作：访问https://cds.climate.copernicus.eu/cdsapp#!/search?text=ERA5，注册账号，登录。 ERA5数据浏览器下载方式以下载ERA5小时分辨率的地表参数为例， 访问数据集页面ERA5-Land hourly data from 1950 to present ，会发现大量的复选框； 按需求勾选对应的地表参数，这里包含了个参数的分组， ...

使用gFortran 进行NetCDF文件的读取，需要netcdf库的支持。检查netcdf是否可用： $ ncdump 如果提示不存在命令,就需要自己安装了，可以参考Compilling WRF 安装步骤。Note: 此处介绍的安装方法可能不再适用，请参考最新的NETCDF/HDF5安装说明 — 2023-11-15 使用F90脚本读取NC文件脚本内容如下，读取Hamawari-8 L2CLP文件的经纬度参数： include 'netcdf.inc'!!main character*100 fnameout,tmpf,h8file,geoh8 logical alive!!!!H-8 参数 integer iopen,ioquire,iovar,latid,lonid,ncid,varid,len_file,ioclose integer, parameter :: dize=2401!定义维数 real*4 latitude(dize),longitude(dize)!!! open one h8-file for geolocat ...

Note: 本节介绍的安装方法可能不再适用，请参考最新的NETCDF/HDF5安装说明 — 2023-11-15 最近需要大批量读写一些FY卫星观测产品，决定使用fortran来读写，相对于NCL和Python会快不少，但是需要库依赖。网上的指导很少，也不全面，自己参照着摸索了半天，也算可以满足使用要求了，决定记录一下。 系统：CentOS release 6.7 (Final) 权限：普通用户 gfortran 版本： gcc version 4.4.7 20120313 (Red Hat 4.4.7-18) (GCC) HDF5版本: hdf5-1.8.22.tar.gz HDF5依赖库和依赖NetCDF库来读取NC文件一样，Fortran在读取HDF文件时需要HDF库的支持。HDFx库为诸多语言比如C,C++,Pyhton,Java,Fortran,IDL等提供HDF格式文件操作扩展。这里选择的是HDF5库，官网下载 https://portal.hdfgroup.org/display/support/Downloads.我选择的版本是hdf5-1 ...

背景Cumulative Distribution Function (CDF)累计分布函数 Probablity Distribution Function (PDF)概率分布函数GeoGebra演示 未完待续… 文献 VODCALIU water resources Soil Moisture 示例来源：https://www.researchgate.net/publication/271703152_Cumulative_Density_Function_CDF_matching_approach Source codefunction [SATCDF_ok,COEFF]=QMAPP(OBS,SAT)%--------------------------------------------------------------------------% CDF matching approach% The function allows to rescale the input data (SAT) to match the ...

以AE_L2A数据为例进行处理。下载的数据包括全球16~17条轨道。 初步裁剪： ;;;获得轨道的四至点行列下标ji = region_ind(Latitude_Low_Res_2d,Longitude_Low_Res_2d, latS, latN, lonW, lonE)jStrt = ji(0) ; lat startjLast = ji(1) ; lat lastiStrt = ji(2) ; lon startiLast = ji(3) ; lon lastLAT2D = Latitude_Low_Res_2d(jStrt:jLast,iStrt:iLast) grey dashed表示的是全部的轨道覆盖范围。阴影区是上述代码裁处的涉及到ROI(region of interest)的轨道。这样可以节省内存开销。进一步的，使用land> 50% 进行水体缺测

最近阅读被动微波天线信号积分时，亲自画了一张图像，展示立体角微分表达式的推导场景： 天线瞬时接受的信号来自于卫星瞬时视野IFOV内部地面上的不同位置的信号贡献效果： 天线温度为TA，地面上p点的亮温贡献为TB(p) 但是，由于卫星在积分时段内T内是移动的，有效视场EFOV必须是瞬时视场的时间积分，并且瞬时视场是时间的函数。那么序号为i的EFOV时刻积分区间[ti-T/2,ti+T/2],有效天线温度为： 虽然随着卫星移动和扫描，TB的采样位置会随时间变化，但是在极短的积分时间内(<8ms)，可以近视忽略亮温场的变化。那么在ti的TA可以表示为等效的天线增益图案和一个近视不变EFOV内亮温场的积分： 等效的天线增益图案可以表示为时间的积分： 按照Backus-Gilbert,Stogryn等人的处理思路，观测Resample过程中，粗分辨率的重建信号和内部多个i序号的高分辨率观测的关系应该是一个线性平均问题： 使用一组权重系数对信号进行重建，权重系数的约束条件为： NSIDC0630 ATB ...

文件说明mlse.land98.res1 为最新反演的地表微波比辐射率数据，第1个分辨率对应的卫星像斑直径为56km，格点距离10kmmlse.land98.res3 为最新反演的地表微波比辐射率数据，第3个分辨率对应的卫星像斑直径为21km，格点距离10kmMYD13C2.006 是月平均NDVI数据，存储格式是HDF格式，分辨率是5kmMYD13C1.005 是16天平均NDVI数据，存储格式是HDF格式，分辨率是5km 通道信息 res1：6.925v, 6.925h, 10.65v, 10.65h, 18.7v, 18.7h, 23.8v, 23.8h, 36.5v, 36.5h, 89.0v, 89.0h 单位GHz res3：18.7v, 18.7h, 23.8v, 23.8h, 36.5v, 36.5h, 89.0v, 89.0h 单位GHz 读取MLSE数据mlse是以TXT文件格式存储的，可以用NCL或者Fortran读取。 诸参数说明如下： ip integer型 每个AMSR-E数据点所拥有的唯一ID； ipp integer ...