Как получить долготы и широты точек сетки в десятичных градусах вместо метра - PullRequest
Новая
       19

Как получить долготы и широты точек сетки в десятичных градусах вместо метра

0 голосов
/ 11 октября 2019

Я создаю сетку над провинцией Цзинсу с функцией make-сетки r-package (sp). Я хочу рассчитать расстояние между каждым центром сетки и местом отбора проб водного бентоса, но единицей долготы и широты сетки являются метры, а единицей долготы и широты водного бентоса являются десятичные градусы. Вопрос: как я могу перевести метры в десятичную степень? поэтому я могу рассчитать их расстояние.

Я пытался искать в интернете несколько дней, но до сих пор не могу решить эту проблему. Ваш ответ очень важен! 

library(sp)
library(rgdal)
Yixing <- readOGR(dsn = "D:/Data_R_afu/yixing", layer = "boundary")<br/>
grid <- makegrid(Yixing, cellsize = 1000) # cellsize in map units!<br/>
grid <- SpatialPoints(grid, proj4string = CRS(proj4string(Yixing)))<br/>
grid@coords<br/>
#       x1      x2<br/>
#[1,] 13311000 3627100<br/>
#[2,] 13311000 3628100<br/>
#[3,] 13312000 3628100<br/>
#[4,] 13313000 3628100<br/>
#...<br/>
#...<br/>
#...<br/>

waterbenthos<- readOGR(dsn = "D:/Data_R_afu/yixing", layer = 
                             "afu_waterbenthos_r_ln(Rt+1Dg+1)")<br/>
waterbenthos@coords<br/>
#      coords.x1 coords.x2<br/>
#[1,]  119.8853  31.37628<br/>
#[2,]  119.8853  31.37628<br/>
#[3,]  119.8853  31.37628<br/>
#[4,]  119.9273  31.40567<br/>
#...<br/>
#...<br/>
#...<br/>
...