IDL处理葵花8Himawari-8标准HSD数据——制作大气校正数据集（卫星角度数据）

IDL制作卫星天顶角和方位角数据集，对应NC文件中的SAA、SAZ

PRO tempCOMPILE_OPT idl2e=envi(/headless)locRaster=e.OpenRaster('**\loc.dat')locData=locRaster.GetData()pi=3.14159265358979323846d2r=0.01745329252;pi/180r2d=57.295779513;180/picols=(locdata.dim)[0]lines=(locdata.dim)[1]outdata=make_array(cols,lines,2,/FLOAT,VALUE=0)lon=locData[*,*,0]lat=locData[*,*,1]alt=locData[*,*,2]#单位：km;计算卫星天顶角、方位角satz=obslook_sat(lon, lat, alt);输出角度数据集outraster=ENVIRaster(outdata,URI='*\satz.dat')outraster.SavelocRaster.CLoseoutraster.Closee.Close
END

代码中的’**\loc.dat’是影像的像元对应的经纬度数据，可以查看https://blog.csdn.net/qq_33339770/article/details/102957857
其中有生成loc文件的过程。

角度数据集生存代码参看：https://download.csdn.net/download/urbancorbie/10343539?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522158769212919726869052219%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fall.57649%2522%257D&request_id=158769212919726869052219&biz_id=1&utm_source=distribute.pc_search_result.none-task-download-2_allfirst_rank_v2~rank_v25-1
生成卫星天顶角的函数如下：

;卫星天顶角、方位角
Function obslook_sat,lon,lat,altCOMPILE_OPT idl2pi=3.14159265358979323846d2r=0.01745329252;pi/180r2d=57.295779513;180/picols=(lon.dim)[0]rows=(lon.dim)[1]clat=atan(0.993305616*tan(lat*d2r));'Rpol^2 / Req^2'=0.993305616clon=make_array(cols,rows,2,/FLOAT,value=d2r*lon);clon=d2r*lonslon=make_array(cols,rows,2,/FLOAT,value=d2r*140.7);'sub_lon'=140.7,slon=d2r*140.7;计算卫星与地面距离;  'Earth's polar radius'=6356.7523,'(Req^2 / Rpol^2) / Req^2'=0.00669438444rp=6356.7523/(sqrt(1-0.00669438444*cos(clat)*cos(clat)))+alt;计算地面坐标x0=rp*cos(clat)*cos(clon)y0=rp*cos(clat)*sin(clon)z0=rp*sin(clat);计算卫星坐标;'Distance from Earth's center to virtual satellite'=42164kmxs=42164*cos(slon)ys=42164*sin(slon)zs=make_array(cols,rows,2,/FLOAT,value=0);zs=0;计算差向量rx=xs-x0ry=ys-y0rz=zs-z0;计算地心坐标系下的极坐标形式的差向量;南方,极坐标下：x=r*sin(el)*cos(az)ls=sin(lat*d2r)*cos(clon)*rx+sin(lat*d2r)*sin(clon)*ry-cos(lat*d2r)*rz;西方,极坐标下：y=r*sin(el)*sin(az),计算的是东方lse=-sin(clon)*rx+cos(clon)*ry;天顶,极坐标下：z=r*cos(el)lz=cos(lat*d2r)*cos(clon)*rx+cos(lat*d2r)*sin(clon)*ry+sin(lat*d2r)*rz;计算方位角：az=atan(y/x)az=atan(-lse/ls)az[where(ls gt 0)] += piaz[where(az lt 0)] += pi*2;计算高度角rg=sqrt(rx*rx+ry*ry+rz*rz)el=asin(lz/rg)el=el*r2daz=az*r2d;方位角修正为[-180.0,180.0]az[where(az GT 180.0)] -= 360.0;天顶角为高度角的余角return,[az,90.0-el]
END

IDL处理葵花8Himawari-8标准HSD数据——制作大气校正数据集（卫星角度数据）相关推荐

python大气校正_Python调用C2RCC实现批量数据的大气校正
目录 1.前言 2.配置环境 3.程序 1.前言 C2RCC大气校正方法在SNAP中作为一个工具出现,安装好SNAP后,就会有该工具.但是发现在SNAP中直接打开C2RCC对OLCI数据进行大气校正时 ...
MODISL1B数据FLAASH大气校正
MODIS数据是著名的大尺度低空间分辨率数据源,一般用于大范围监测.宏观信息的快速提取.在MODIS数据预处理时,最常做的是几何校正.对MODIS数据进行几何校正的工具很多,诸如MCTK ...
哨兵2数据预处理大气校正分为2017之前与2017年之后
哨兵2数据下载可以从USGS和欧空局下载,下载网站自行搜索. 由于2016年年底 , 欧空局改变了对哨兵2 数据的命名.哨兵2数据其他博客说2017年之前的数据,可以采用envi5.3打开,同时 en ...
StaMPS中使用GACOS数据进行大气校正，解缠、计算形变速率
StaMPS中使用GACOS数据进行大气校正 1 GACOS数据下载 2 查看GACOS用法 3 在StaMPS中使用GACOS进行大气校正 4 stamps6 相位解缠 stamps7时间形变速率 ...
Sentinel2（哨兵二号）1C数据转2A数据（大气校正）
1.将下载下来的1C数据进行解压,文件夹以.safe结尾. 2.Sen2Cor软件下载安装,路径不要有中文,参照论文(基于 Sentinel-2 的老挝北部刀耕火种农业遥感监测及其特征研究)数据预处理 ...
Worldclim网站CIMP6未来数据下载（校正处理好的数据）
CIMP6-2.1新版数据下载:(以下为校正过的数据,全球范围,精度2.5m,格式ASC,每个文件分别为处理好的BIO1-BIO19变量,19个处理好的ASC文件) CIMP6-BCC-CSM2-MR ...
matlab 导入/制作/导出 STK .tle卫星文件/数据（直接/间接，好几种方法）
导入/制作/导出卫星tle文件整理: 代码: 整理: 1.在STK help中的 inFrance Help里的core code 里能够导入的文件只有.a文件 2.导出.tle文件参考拜火老师的 ...
多通道时序数据读取matlab,科学网—卫星时序数据-平滑-插补-MATLAB - 李雷的博文...
%程序用于平滑处理数据,以及将间隔时间数据插值成每天 clc clear %光学卫星获取的数据存在受云影响的情况,通过滤波的方式将极低或极高等不合理值进行处理(以MODIS中的LAI为例) [a,b] ...
ENVI软件对Landsat-8数据进行辐射定标、大气校正、提取NDVI、估算植被覆盖度等操作【图说GIS】
目录一.前言二.数据下载第一步访问网站第二步登录账户第三步点击高级检索, 下载数据三.软件下载及安装四.辐射定标第一步打开数据第二步辐射定标五.计算区域平均高程第一步 ...
《大数据》第1期“研究”——大数据是数据、技术，还是应用
大数据是数据.技术,还是应用朱扬勇1,2,熊赟1,2 1. 复旦大学计算机科学技术学院上海 201203:2. 上海市数据科学重点实验室(复旦大学) 上海 201203 摘要:通常认为大数据是一 ...

IDL处理葵花8Himawari-8标准HSD数据——制作大气校正数据集（卫星角度数据）

IDL处理葵花8Himawari-8标准HSD数据——制作大气校正数据集（卫星角度数据）相关推荐

最新文章

热门文章