本文介绍了在 r 中将 Lambert 共形圆锥投影转换为 wgs84的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!

问题描述

我有 Lambert 共形圆锥投影 x,y 信息.
我需要 WGS84 坐标.但我不知道什么是 lcc.
我在下面提供了 lcc 信息.
有没有办法在 r 中将 lcc 转换为 WGS84?

I have Lambert conformal conic projection x,y information.
I need the WGS84 coordinate. But I don't know what is lcc exactly.
I have provided the lcc information below.
Is there a way to convert lcc to WGS84 in r?

example lcc x,y : xy <- cbind(c(509535.7, 514535.7),c(201098.6, 201098.6))

lcc 信息:
第一标准平行纬度:30.0
第二标准纬线纬度:60.0
原点纬度:38.0原点经度:126.0
计算点东东:43
计算点北距:136
4 边经纬度:左上(43.3935, 123.3102), 左下(31.7944, 123.7613),
右上(43.2175, 132.7750), 右下(31.6518, 131.6423)

lcc information :
Latitude of first standard parallel : 30.0
Latitude of second standard parallel : 60.0
Origin latitude : 38.0 Origin longitude : 126.0
Easting of computation point : 43
Northing of computation point : 136
4 edge lon,lat point : left-upper(43.3935, 123.3102), left-lower(31.7944, 123.7613),
right-upper(43.2175, 132.7750), right-lower(31.6518, 131.6423)

推荐答案

您所要求的实际上是不可能的,因为 Lambert 等角圆锥曲线是地图投影,而 WGS84 是基准.很可能您的 LCC 数据已经与 WGS84 相关.我假设您想将 LCC 转换为经度/纬度.(而且我还假设您的输入数据是 WGS84.)

What you are asking for is really not possible because Lambert conformal conic is a map projection while WGS84 is a datum. It may well be that your LCC data already are relative to the WGS84. I assume you want to transform LCC to longitude/latitude. (And I am also assuming that your input data is WGS84.)

xy <- cbind(c(509535.7, 514535.7),c(201098.6, 201098.6))

library(sp)
library(rgdal)
crs <- CRS("+proj=lcc +lat_1=30 +lat_2=60 +lat_0=38 +lon_0=126 +datum=WGS84")
p <- SpatialPoints(xy, proj4string=crs)
g <- spTransform(p, CRS("+proj=longlat +datum=WGS84"))
coordinates(g)

对于栅格数据(因为您将 'raster' 作为关键字列出),请参阅 raster::projectRaster

For raster data (since you list 'raster' as a keyword), see raster::projectRaster

这篇关于在 r 中将 Lambert 共形圆锥投影转换为 wgs84的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持!

07-13 19:27