单站点的值对场的相关性分析

二子搞研究

我这边有一个站点的实测值，但时间序列是不连续的，想要对再分析数据中的某些变量场做相关性分析，也就是点对场的相关性分析，查看之前的贴子可以使用tcorr函数，但因为我的是txt文件是不是需要：转换成grd数据→转换成格点数据→相关性分析，这么一个思路？但时间是不连续的，是否可以挑出这些时间或者将其余时间补上赋值为undef的值？@river

river

思路差不多，只有一个地方有问题，你要做的就是点对场的相关，干嘛要把点的数据转换成格点的呢。那不是成了场和场的相关了，那要用的函数就不同了

二子搞研究

river 发表于 2016-2-27 21:25
思路差不多，只有一个地方有问题，你要做的就是点对场的相关，干嘛要把点的数据转换成格点的呢。那不是成了 ...

哦哦，可是因为时间不连续，是需要补齐的吗，默认缺测值？如果设置好了，系统是不是就会自动跳过缺测值部分去求其他实测值的相关性？

river

二子搞研究 发表于 2016-2-29 09:24
哦哦，可是因为时间不连续，是需要补齐的吗，默认缺测值？如果设置好了，系统是不是就会自动跳过缺测值部 ...

设置缺测也行，但是一个时间序列里面有缺测，然后这个点和一个场相关，感觉结果会不会很挺奇怪······

二子搞研究

river 发表于 2016-2-29 21:03
设置缺测也行，但是一个时间序列里面有缺测，然后这个点和一个场相关，感觉结果会不会很挺奇怪····· ...

其实要是能把那些仅有数据的部分挑选出来做相关也行，可是这样就没有办法生成grd文件了吧，不知道大神您知不知道有啥好办法不呢

river

二子搞研究 发表于 2016-2-29 21:24
其实要是能把那些仅有数据的部分挑选出来做相关也行，可是这样就没有办法生成grd文件了吧，不知道大神您 ...

你仔细看下这个求相关函数的用法，我不太理解你说的时间不连续是怎样不连续
tcorr()

tcorr (expr1, expr2, tdim1, tdim2)

This function produces a spatial map of temporal correlation coefficients.

expr1   - a valid GrADS expression that varies in time
expr2   - a valid GrADS expression that varies in time and may also vary in X and Y
tdim1   - starting time dimension expression
tdim2   - ending time dimension expression
The expr1 time series is correlated to the time series at each grid point in expr2. The result is a grid of correlation coefficients that matches the X and Y dimensions of expr2.

Usage Notes

If both expr1 and expr2 vary only in time, the output is a single value.
Use the scorr function to do correlation over the spatial domain.
Example

This example calculates the temporal correlation between sea level pressure and the defined variable elnino, an areal average of surface temperature in the equatorial Pacific.
set x 1
set y 1
set z 1
set t 1 100
define elnino = aave(ts,lon=-160,lon=-80,lat=-10,lat=10)
set lon -180 180
set lat -90 90
set z 1
set t 1
d tcorr(elnino, slp, t=1, t=100)

二子搞研究

river 发表于 2016-3-1 07:09
你仔细看下这个求相关函数的用法，我不太理解你说的时间不连续是怎样不连续
tcorr()

嗯嗯好的，时间不连续，举个列子：2013.10.28 2014.1.2,2014.2.3,2014.3.4，就只有这些时间有数据，其他时间全无数据，当然实际比这个时间段长，但意思就是这样的

river

二子搞研究 发表于 2016-3-1 09:32
嗯嗯好的，时间不连续，举个列子：2013.10.28 2014.1.2,2014.2.3,2014.3.4，就只有这些时间有数据，其他 ...

那我觉得可以利用趋势分析之类的方法把缺测的数据补上，不然的话我觉得结果不一定好

二子搞研究

river 发表于 2016-3-1 20:14
那我觉得可以利用趋势分析之类的方法把缺测的数据补上，不然的话我觉得结果不一定好

可是缺测值好多啊，我再啃啃吧，多谢大神

喃音

水经验来着{:5_231:}