首页 > 解决方案 > R代码将调查数据(包括GPS)与最近的经纬度的温度数据合并,以及调查月份前12个月(光栅包?)

问题描述

在 R 中,我想将包含受访者GPS 位置的dhs 调查数据温度数据air.mon.mean.v501.nc(陆地气温:1900-2017 网格化每月时间序列)(https://psl. noaa.gov/data/gridded/data.UDel_AirT_Precip.html)。


到目前为止:我的调查数据集(称为kids.dta”包含每位受访者的一项观察结果。除此之外,一个变量是“ lon ”,一个是“ lat ”,一个是“ interview_month ”,“ interview_year ”。纬度值范围从 -90 到 90;经度值范围从 -180 到 180;

例如,纬度 27.41974,经度 87.60922,interview_month 12,interview_year 2001

使用以下 R 代码,温度数据遵循与调查数据相同的坐标网格。

open_air_brick <- 砖(“air.mon.mean.v501.nc”,var="air")

open_air <- 旋转(open_air_brick)

请注意:温度数据的纬度和经度四舍五入:例如 lat 27.5, lon 87.5 (!)



我想要什么:我想要一个最终数据集,其中包含调查数据以及调查月份(调查月份和之前的 12 个月)之前的一年中的相应月度温度,基于它们的位置(lon / 纬度)。



我想我需要做的:

  1. 为调查数据中的每个观测值在温度数据中找到最接近的 lon/lat(例如 lat 27.41974 最接近 lat 27.5),并将其分配给调查数据

  2. 根据先前分配的最近经度和纬度将调查数据与温度数据合并,并且仅提取调查当月和调查前 12 个月的温度数据。


老实说,我不知道该怎么做。到目前为止,我一直在使用 R Raster 包

如果有人有任何想法或提示,我将不胜感激!

标签: rmergegeospatiallatitude-longituderaster

解决方案


推荐阅读