R语言ampStata已知某个
为了让大家更好的理解这个方法,欢迎各位培训班会员参加明晚8点的直播课「已知某个地点的经纬度,如何判断其是在北方还是南方?」
最近有个小伙伴问了这样一个问题:
已知某个地点的经纬度,如何判断其是在北方还是南方呢?
在R语言中这个问题非常简单,之前我们也推出过课程讲解:如何判断工企业在南方还是北方计算每个工企业距离秦岭淮河线的距离?(在平台上搜索即可找到)。
但是很多小伙伴还不会使用R语言,所以我们今天就来针对Stata用户讲解这个问题的解决方案。
大家都知道南北方实际上是个很模糊的概念,虽说是以秦岭-淮河为界,但是这个界太宽了,所以我们就选择一条看起来差不多的线:
使用R语言判断某个点是位于北方还是南方是个很简单的事情,例如待转换.csv文件是一些坐标点:
library(tidyverse)library(sf)read_csv("待转换.csv")%%mutate(lon=as.numeric(lon),lat=as.numeric(lat))%%st_as_sf(coords=c("lon","lat"),crs=)%%mutate(id=as.numeric(row.names(.)))-dfnread_sf("秦岭淮河以南.geojson")-southread_sf("秦岭淮河以北.geojson")-northbind_rows(dfn%%st_intersection(south)%%mutate(lon=st_coordinates(geometry)[,1],lat=st_coordinates(geometry)[,2])%%st_drop_geometry(),dfn%%st_intersection(north)%%mutate(lon=st_coordinates(geometry)[,1],lat=st_coordinates(geometry)[,2])%%st_drop_geometry())
为了让大家不使用R语言也能根据经纬度坐标判断南北方,我编写了一个shiny应用:northorsouth:
在线使用:
转载请注明:http://www.abuoumao.com/hykh/654.html