我希望使用SeaDAS 7.2对超过365个数据集进行像素提取。 这是我的数据集输入文件名:
Input Filename list:
2014-01-01_INDESO_BIO_1dAV-RT.nc
2014-01-02_INDESO_BIO_1dAV-RT.nc
2014-01-03_INDESO_BIO_1dAV-RT.nc
.....
如何定义文件名中的时间提取模式'因此,对于每个输入数据集,将生成具有相同文件名的分离输出数据:
Output Filename list:
2014-01-01_INDESO_BIO_1dAV-RT.txt
2014-01-02_INDESO_BIO_1dAV-RT.txt
2014-01-03_INDESO_BIO_1dAV-RT.txt
.....
答案 0 :(得分:0)
我已经使用R raster package
解决了我的问题。这是我的代码:
# define xy coordinate
# xy can be a single point or data frame
xy <- SpatialPoints(cbind(long, lat))
pixel <- extract(raster(raster_filename, varname='variabel_name'), xy)