找回密码
 立即注册
查看: 3725|回复: 6

Openstreetmap的问题,报错Error in minY:maxY : NA/NaN argument

[复制链接]
发表于 2017-10-17 02:37:26 | 显示全部楼层 |阅读模式
求各位大神帮忙看看是什么问题~~~

load(file.choose())
crs.val<-CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")
## Load an OpenStreetMap backdrop
ul<-as.vector(cbind(bbox(en.county)[2,2], bbox(en.county)[1,1]))
lr<-as.vector(cbind(bbox(en.county)[2,1], bbox(en.county)[1,2]))
# download the map tile
MyMap<-openmap(ul,lr,9,'osm')


运行openmap这一步报错  Error in minY:maxY : NA/NaN argument
我输入的数据是一个Large SpatialPolygonDataFrame, ul/lr所获取的值分别为num[1:2] 416919 1506734及num[1:2] 525421 1441727

求问这是什么问题??
跪谢!
回复

使用道具 举报

发表于 2017-10-17 10:03:03 | 显示全部楼层
UL和LR是用经纬度表示的,经度的最大值是180,纬度的最大值是90。但是你输入的数据明显不在这个范围,所以报错。
回复

使用道具 举报

 楼主| 发表于 2017-10-17 10:31:45 | 显示全部楼层
chunfei_xu 发表于 2017-10-17 10:03
UL和LR是用经纬度表示的,经度的最大值是180,纬度的最大值是90。但是你输入的数据明显不在这个范围,所以 ...

感谢大神!然后问一下,现在我改成下边这样了,
> crs.val<-CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")
> ## Load an OpenStreetMap backdrop
> # download the map tile
> MyMap<-openmap(c(-78,42),
+                c(-80,44),
+                9,'osm')
> plot(MyMap, removeMargin=FALSE)


然而输出的结果只有一个蓝色的矩形,而没有streetmap,这是为什么?
12.png
谢谢

回复

使用道具 举报

发表于 2017-10-17 11:46:06 | 显示全部楼层
经纬度的值应该都是大于零的,然后78和80之间的差别太小,可以放大一点,42和44也一样。可以再试试看!至于经纬度如何唯一定位一个位置,我也还没有搞明白,理论上一个经纬度可以代表四个地理位置。
回复

使用道具 举报

发表于 2017-10-17 13:40:43 | 显示全部楼层
我似乎搞明白了,如果用(X1,Y1)代表UL,(X2,Y2)代表LR。那么当X1>X2是表示北纬,否则表示南纬;Y1<Y2是表示表示东经,否则表示西经。
回复

使用道具 举报

 楼主| 发表于 2017-10-17 23:00:54 | 显示全部楼层
感谢大神赐教!搞定了,谢了
回复

使用道具 举报

发表于 2017-10-18 09:05:29 | 显示全部楼层
有帮助就好,
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Archiver|手机版|小黑屋|R语言中文网

GMT+8, 2024-11-22 14:28 , Processed in 0.036351 second(s), 21 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表