R 检查点是否与多边形重叠
我有数据可以绘制每天风暴的纬度,如下所示(我删除了day ID,所以它只有年份和月份)R 检查点是否与多边形重叠,r,intersection,sp,rgdal,R,Intersection,Sp,Rgdal,我有数据可以绘制每天风暴的纬度,如下所示(我删除了day ID,所以它只有年份和月份) lat_lonDatlong/lat坐标应交换: coordinates(lat_lonDat) <- ~ Longitude + Latitude 坐标(横向水平)尾部(上方(横向水平,温度shp),15) GID_0名称_0 GID_1名称_1 NL_名称_1 GID_2名称_2 印度49印度2_1安得拉邦2.5_1克里希纳 50印度印度印度印度安得拉邦印度冈图尔 51印度工业2_1安得拉邦工
lat_lonDatlong/lat坐标应交换:
coordinates(lat_lonDat) <- ~ Longitude + Latitude
坐标(横向水平)尾部(上方(横向水平,温度shp),15)
GID_0名称_0 GID_1名称_1 NL_名称_1 GID_2名称_2
印度49印度2_1安得拉邦2.5_1克里希纳
50印度印度印度印度安得拉邦印度冈图尔
51印度工业2_1安得拉邦工业2.4_1冈图尔
52印度印度印度安得拉邦印度冈图尔
53印度印度印度安得拉邦印度冈图尔
54印度工业公司32印度泰兰加纳工业公司32印度纳尔贡达工业公司
印度第55工业区第32工业区第32.7工业区第32.7工业区纳尔贡达
印度第56工业区第32工业区第32.7工业区第32.7工业区纳尔贡达
57印度工业公司32蒂兰加纳工业公司32瓦兰加工业公司10蒂兰加工业公司
58印度印度印度印度
59印度工业公司32印度泰兰加纳工业公司32印度阿迪拉巴德
印度第60工业区印度第20工业区马哈拉施特拉邦第20.8工业区钱德拉普尔
印度61印度20印度1马哈拉施特拉邦20印度8印度1钱德拉普尔
62印度印度印度印度印度马哈拉施特拉邦
印度63印度印度20印度马哈拉施特拉邦20印度11印度冈迪亚
VARNAME_2 NL_NAME_2 TYPE_2 ENGTYPE_2 CC_2 HASC_2
公元49年基斯特纳区
公元50年
公元前51区
公元前52区
公元53年
第54区
北卡罗来纳州第55区
塔格纳56区
西澳第57区
公元58年
59.TG.AD区
MH.CH昌达区60号
61.MH.CH昌达区
伯克希尔哈撒韦州62区
63.MH.GO区
library(raster)
library(sp)
library(ggplot2)
temp.shp <- getData('GADM', country = 'IND', level = 2)
ggplot() +
geom_polygon(data = temp.shp, aes(x = long, y = lat, group = group), fill = NA, colour = "black") +
geom_point(data = lat_lonDat, aes(x = Longitude, y = Latitude, colour = 'red')) +
coord_quickmap() + theme_classic()
coordinates(lat_lonDat) <- ~ Latitude + Longitude
proj4string(lat_lonDat) <- proj4string(temp.shp)
over(lat_lonDat, temp.shp)
coordinates(lat_lonDat) <- ~ Longitude + Latitude
R> tail(over(lat_lonDat, temp.shp), 15)
GID_0 NAME_0 GID_1 NAME_1 NL_NAME_1 GID_2 NAME_2
49 IND India IND.2_1 Andhra Pradesh <NA> IND.2.5_1 Krishna
50 IND India IND.2_1 Andhra Pradesh <NA> IND.2.4_1 Guntur
51 IND India IND.2_1 Andhra Pradesh <NA> IND.2.4_1 Guntur
52 IND India IND.2_1 Andhra Pradesh <NA> IND.2.4_1 Guntur
53 IND India IND.2_1 Andhra Pradesh <NA> IND.2.4_1 Guntur
54 IND India IND.32_1 Telangana <NA> IND.32.7_1 Nalgonda
55 IND India IND.32_1 Telangana <NA> IND.32.7_1 Nalgonda
56 IND India IND.32_1 Telangana <NA> IND.32.7_1 Nalgonda
57 IND India IND.32_1 Telangana <NA> IND.32.10_1 Warangal
58 IND India IND.32_1 Telangana <NA> IND.32.1_1 Adilabad
59 IND India IND.32_1 Telangana <NA> IND.32.1_1 Adilabad
60 IND India IND.20_1 Maharashtra <NA> IND.20.8_1 Chandrapur
61 IND India IND.20_1 Maharashtra <NA> IND.20.8_1 Chandrapur
62 IND India IND.20_1 Maharashtra <NA> IND.20.5_1 Bhandara
63 IND India IND.20_1 Maharashtra <NA> IND.20.11_1 Gondiya
VARNAME_2 NL_NAME_2 TYPE_2 ENGTYPE_2 CC_2 HASC_2
49 Kistna <NA> District District <NA> IN.AD.KR
50 <NA> <NA> District District <NA> IN.AD.GU
51 <NA> <NA> District District <NA> IN.AD.GU
52 <NA> <NA> District District <NA> IN.AD.GU
53 <NA> <NA> District District <NA> IN.AD.GU
54 <NA> <NA> District District <NA> IN.TG.NA
55 <NA> <NA> District District <NA> IN.TG.NA
56 <NA> <NA> District District <NA> IN.TG.NA
57 <NA> <NA> District District <NA> IN.TG.WA
58 <NA> <NA> District District <NA> IN.TG.AD
59 <NA> <NA> District District <NA> IN.TG.AD
60 Chanda <NA> District District <NA> IN.MH.CH
61 Chanda <NA> District District <NA> IN.MH.CH
62 <NA> <NA> District District <NA> IN.MH.BH
63 <NA> <NA> District District <NA> IN.MH.GO