ui外包网站,长沙竞价网站建设价格,深圳推广公司哪家好,旅行社网站开发Sentinel-5P是欧空局#xff08;Europe Space Agency#xff0c;ESA#xff09;于2017年10月13日发射的一颗全球大气污染监测卫星。卫星搭载了对流层观测仪#xff08;Tropospheric Monitoring Instrument#xff0c;TROPOMI#xff09;#xff0c;可以有效的观测全球各…Sentinel-5P是欧空局Europe Space AgencyESA于2017年10月13日发射的一颗全球大气污染监测卫星。卫星搭载了对流层观测仪Tropospheric Monitoring InstrumentTROPOMI可以有效的观测全球各地大气中痕量气体组分包括NO2、O3、SO2、 HCHO、CH4和CO等重要的与人类活动密切相关的指标加强了对气溶胶和云的观测。
⛄Sentinel-5P数据介绍
数据产品
TROPOMI是目前世界上技术最先进、空间分辨率最高的大气监测光谱仪成像幅宽达2600km每日覆盖全球各地成像分辨率达7km×3.5kmESA提供了L1B和L2两种级别的数据下载。L2级数据产品介绍如下表
官网产品数据介绍https://sentiwiki.copernicus.eu/web/sentinel-5p
Product typeParameterL2__O3____Ozone (O3) total column总柱含量L2__O3_TCLOzone (O3) tropospheric column对流层柱含量L2__O3__PROzone (O3) profile总剖面数据L2__NO2___Nitrogen Dioxide (NO2)total and tropospheric columns总柱和对流层柱含量L2__SO2___Sulfur Dioxide (SO2) total column总柱含量L2__CO____Carbon Monoxide (CO) total column总柱含量L2__CH4___Methane (CH4) total column总柱含量L2__HCHO__Formaldehyde (HCHO) total column总柱含量L2__CLOUD_Cloud fraction, albedo, top pressure云量、反照率、云顶大气压L2__AER_AIUV Aerosol Index紫外区域气溶胶指数L2__AER_LHAerosol Layer Height (mid-level pressure)气溶胶层高度中等气压UV productSurface Irradiance/erythemal doseL2__NP_BDx, x3, 6, 7Suomi-NPP VIIRS Clouds云产品辅助数据
L2级数据产品根据生产时间/质量分为三种数据流
近实时数据流near-real-timeNRTI卫星成像3小时后即可获取数据可能不完整或存在质量缺陷。离线数据流OfflineOFFL一般成像后几天即可获得。再次处理数据流ReprocessingRPRO有些数据可能经过了多次处理获得的最新的质量最佳的版本。
一般情况下长期的时序变化研究不可以混用不同级别的数据流推荐使用最新的RPRO数据以保证数据质量。
数据下载
下载教程Sentinel–哨兵系列数据下载地址及流程更新【202310】
⛄Sentinel-5P数据预处理
本博客以L2级HCHO产品近实时数据流HCHO总柱含量为例 S5P_NRTI_L2__HCHO___20241015T052816_20241015T053316_36303_03_020601_20241015T061049.nc 方法1ENVI处理 GLT几何校正法利用输入的几何文件生成一个地理位置查找表文件Geographic Lookup Table, GLT从该文件中可以了解到某个初始像元在最终输出结果中实际的地理位置。 地理位置查找表文件是一个二维图像文件文件中包含两个波段:地理校正影像的行和列。文件对应的灰度值表示原始影像每个像素对应的地理位置坐标信息用有符号整型储存它的符号说明输出像元是对应于真实的输入像元是由邻近像元生成的填实像元infill pixel。符号为正时说明使用了真实的像元位置值符号为负时说明使用了邻近像元的位置值值为0说明周围7个象元内没有邻近像元位置值。 GLT文件包含初始影像每个像元的地理定位信息它的校正精度很高。避免了通过地面控制点利用二次多项式几何校正法对低分辨率影像数据的处理。 1打开步骤Open as→Generic Formasts→HDF5/NetCDF-4或者直接将数据文件拖进ENVI。
2地理位置数据和浓度产品数据分开读取将经纬度数据读取成一个图层将HCHO柱浓度数据读取为另一个图层。 3制作GLT文件在ENVI工具箱中搜索Build GLT在弹出的对话框中Input X Geometry Band选择经度数据Input Y Geometry Band选择纬度数据。 在接下来的对话框中投影信息选择WGS-84Geographic Lat/Lon。 4在Build Geometry Lookup File Parameters对话框中像元大小选择默认旋转角度Rotation为0正上方为北。 5利用GLT文件对HCHO数据进行校正工具箱中搜索Georeference from GLT在Input Geometry Lookup File对话框中选择GLT文件在Input Data File对话框中选择待校正的HCHO数据文件。 存在问题针对ENVI方式的处理结果小编存在诸多疑问。我们通过ENVI加载数据或者Panoply软件加载数据获取数据的行列358×450但是经过ENVI处理之后分辨率Pixel: 0.044004 DegreesDims: 849 x 527 x 1。ENVI中列在前行在后出现行列数不对应的问题 方法2Panoply可视化展示
Panoply是一款跨平台的数据可视化软件工具专门用于绘制地理参考以及其他数组数据包括netCDF、HDF、GRIB等格式。Panoply运行需要JAVA环境配置JAVA环境后解压程序包双击Panoply.exe即可。 官网下载https://www.giss.nasa.gov/tools/panoply/ 博客下载遥感数据处理及可视化-第三方辅助工具集 方法3S5Processor处理
为了把S5P的NC数据正确转为TIFF数据试了多种方法在GitHub上面发现了一个R包S5Processor S5Processor包地址https://github.com/MBalthasar/S5Processor R及RStudio下载安装教程参考博客①、参考博客②、参考博客③ R包镜像地址https://cran.r-project.org/src/contrib/Archive/
1首先这个包需要借助Rtools下载Rtools根据自己运行环境选择合适的版本安装即可。
Rtools下载地址https://cran.r-project.org/bin/windows/Rtools/
2提示没有devtools则直接安装该包即可。 3安装S5Processor包
devtools::install_github(MBalthasar/S5Processor)在安装过程中可能会遇到依赖包缺少的问题比如maptools、rgdal、rgeos包需要在R包镜像地址中检索并下载对应版本的包进行手动安装。 4实例代码
基于R语言使用S5Processor包对S5P数据进行TIFF格式转换便于在GIS中进行分析。
devtools::install_github(MBalthasar/S5Processor)
library(S5Processor)
library(ncdf4)
library(ggplot2)
library(dismo)
library(maptools)
library(raster)
library(geosphere)
library(rgdal)
library(rgeos)
library(sp)
S5P_HCHO - S5P_process(input D:\\Desktop\\data\\S5P_NRTI_L2__HCHO\\S5P_NRTI_L2__HCHO___20241015T052816_20241015T053316_36303_03_020601_20241015T061049.nc, product 6)
writeRaster(S5P_HCHO, D:\\Desktop\\data\\S5P_NRTI_L2__HCHO\\S5P_HCHO_R.tif, format GTiff, overwrite TRUE)针对代码中product参数值相信很多朋友疑惑如何获取呢 大部分博客中NO2实例的取值都是一样的比如39但是对于其它产品数据如何确定呢 小编开始也很疑惑最后通过不断调试代码能够输出产品列表判断出对应编号。 5转换结果 存在问题针对R包方式的处理结果小编仍然存在诸多疑问。我们通过ENVI加载数据或者Panoply软件加载数据获取数据的行列358×450但是通过R包方式处理之后分辨率Pixel: 0.240748 x 0.17972 DegreesDims: 155 x 129 x 1。ENVI中列在前行在后依然出现行列数不对应的问题并且和ENVI处理结果的行列数也不相同但是两者结果的整体趋势保持一致。 多谢多谢 笔者不才请多交流
欢迎大家关注预览我的博客BlogHeartLoveLife 能力有限敬请谅解
