Question

我使用geosphere包创建了一个很好的lat / lon轨迹：

flightTraj  = greatCircle( c( originAptLon, originAptLat ), c( destinAptLon, destinAptLat ), n = nPts, sp = TRUE )

其属性为：

class       : SpatialLines 
features    : 1 
extent      : -180, 180, -52.74719, 52.74719  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=WGS84

我想投射到WGS84椭圆体的UTM区域14制图。我试过了：

projectedTrajectories  <- CRS("+proj=utm +zone=14 +datum=WGS84 +units=km +no_defs") %>% 
    spTransform( flightTraj, . )

但是，控制台显示：

non finite transformation detected:
  V1 lat  
  Error in .spTransform_Line(input[[i]], to_args = to_args, from_args = from_args,  : 
  failure in Lines 1 Line 1 points 
  In addition: Warning message:
  In .spTransform_Line(input[[i]], to_args = to_args, from_args =     from_args,  :
  6 projected point(s) not finite

请帮助。

Answer 1

“一切都很顺利......”

#!/bin/sh
sqlite3 test.db <<EOF
create table n (id INTEGER PRIMARY KEY,f TEXT,l TEXT);
insert into n (f,l) values ('john','smith');
select * from n;
EOF

将一个SpatialLines对象从lat / lon投影到utm

1 个答案: