气象要素空间插值分析

基于气象站点

1、日地距离订正系数计算

代入公式计算日地距离订正系数（在EXCEL中完成）：

公式省略

其中（M2=2*pi()(Dn-1)/365）,Dn为计算日期，如此例计算2月15日，则为第36天，则M2=2*pi()(36-1)/365）=0.602497

2、日出日落时刻太阳时角
W0=arccos(-tan(a)*tan(b))其中a为赤纬值，如2月15日赤纬为-16.667°，b为每个站点的纬度，其中a,b都需要转换为十进制的小数度，方法为*pi()/180。结果应该与气象站具有共同数据列。

3、计算2月5日的日辐射量，依据下方公式

其中T=24*60，因为Io取经验值：0.0820单位为MJ.m^2/min，单位为每分钟，故2月5日当天24小时共有1440分钟
其中（1/p）^2为第一步中计算得到的日地距离订正系数，W0为日出日落时刻太阳时角在第二步中算处，并且该值在三角函数计算时不需要*pi()/180转换。该值即为日辐射量。

4、生成矢量数据
在ArcMap中通过Display XY生成基于经纬度的点坐标，且生成同时定义大地坐标系为WGS84坐标系。之后添加中国省级行政区划图，点面矢量图相叠合，为点矢量图进行投影转换，转换为面矢量图的投影系，效果如图

矢量点为每一个气象站点的代表，具有日照辐射量值属性（第3步计算得到），对该值进行IDW插值，进行地图整饰，出图之后如图

气象要素空间插值分析相关推荐

ArcGIS实验教程——实验十三：栅格空间插值分析
ArcGIS实验视频教程合集:<ArcGIS实验教程从入门到精通>(附配套实验数据) >>> 实验前必读:<完美解决ArcGIS10.x栅格空间插值报错无法进行和导 ...
你要的来了：ArcGIS空间插值分析方法权威解读
插值问题的提出??? 一.趋势面 Trend is a global polynomial interpolation that fits a smooth surface defined b ...
ArcGIS_空间插值分析
完美解决ArcGIS10.x栅格空间插值报错无法进行和导出插值栅格结果出错的问题
不少GISer,尤其是初学者,对栅格数据的操作,主要是栅格插值无法进行,或者是插值后的栅格结果无法导出,或者导出结果为空,或者是栅格格式无法转换,今天刘一哥GIS就帮您排忧解难吧. 其实,关于栅格(尤 ...
gstat | 空间插值（三）——克里金插值之泛克里金和简单克里金
本篇接着上篇继续介绍克里金插值.首先加载相关工具包和上篇使用的示例数据: library(gstat) library(sf) library(tidyverse) library(readxl) l ...
gstat | 空间插值（二）——克里金插值之普通克里金
说明:昨天的推文误把可吸入颗粒物当作PM2.5,实应该为PM10,这里修正后重发. 从本篇开始计划分三篇介绍克里金插值.与反距离权重插值不同,克里金插值是无偏估计,其中也涉及到模型估计.本篇先对普通克 ...
白话空间统计二十九：空间插值(二）
前文再续,书接上一回--上回书说到,空间插值可以对数据进行估算,这一张我们来具体说说空间插值的一些概念. 首先,来说说插值的输入与输出问题. 首先输入的肯定是用来进行插值的观测点数据,一般来说都是点要 ...
GEE：克里金 Kriging 空间插值（以陕西省2013年生物量为例）
作者:CSDN @ _养乐多_ 本文记录了在Google Earth Engine(GEE)平台上进行 Kriging 插值的介绍和代码案例.本文通过选取的2013年陕西省生物量样本点数据为例,利用 ...
arcgis中的插值分析了解
1.前言打开arcgis工具箱我们可以发现arcgis为我们提供了多种插值分析的工具,克里金插值.反距离权重插值.样条函数插值等,具体如图所示在工具箱中除地形转栅格和通过文件实现地形转栅格是用来绘 ...
第十五章栅格数据重分类、栅格计算器、插值分析
文章目录第十五章栅格数据分析第一章栅格数据重分类第一节栅格数据重分类第二节栅格重分类的使用第三节重分类的使用中的空值使用第四节重分类的案例:分类统计面积第五节坡度矢量分级图 ...

气象要素空间插值分析

基于气象站点

气象要素空间插值分析相关推荐

最新文章

热门文章