如何从纬度和经度.tif栅格对象生成时区?

时间:2018-04-09 17:05:24

标签: r gis r-raster

我真的很感激任何帮助。

我需要按水密度和时区重新分类用水总量。

我是通过水密度来完成的,但我不知道如何按时区进行。

可以找到数据here

这是我的示例代码:

lowbound<- seq(0,cellStats(denty,"max", na.rm=TRUE),1000) #bands by 1k
upbound<-lowbound+1000
band<-upbound/1000
rclmat <- as.matrix(cbind(lowbound,upbound,band))
rclmat<-rclmat[1:length(rclmat[,1]),]
#here I need to reclassify by rclmat and time zone
#I would like some insights on how to get here time zone from lat/long to be used on reclassify function
water_density<- reclassify(denty, rclmat, include.lowest=TRUE, right=FALSE)

#stats per band
sumt<-zonal(total,water_density,"sum")
stats2<-cbind(sumt[,2])

非常感谢

0 个答案:

没有答案