У меня есть сетка данных высот ascii в ячейках 100M.У меня есть широта / долгота для начальной точки, и мне нужно сложить эти 100-метровые приращения и получить новую координату широты / долготы, чтобы нарисовать цветную рамку.
Я использую R и буклет для рисования картыс данными высоты под серией маркеров.Я хочу, чтобы хороший слой с цветовой кодировкой на листовой карте отображал высоту.
coords <- function(cells_east, cells_north) {
cells_east =1000
cells_north=1000
start_lat <- 48.922582
start_lng <- -101.497821
kmeast<-cells_east/10 #each cell 100 m
kmnorth<-cells_north/10
#Latitude: 1 deg = 110.574 km
#Longitude: 1 deg = 111.320*cos(latitude) km
lat_east <- kmeast /110.574
lng_north <- kmnorth / (111.32 * cos(start_lat))
elevation <- dem$data[cells_east,cells_north] #elevation from ascii grid
ret <- c(lat_east+start_lat,lng_north+start_lng,elevation) #lat,lng,elev
}
Я перехожу к функциям преобразования из Интернета в нижней части, но кажется, что он отключен.Этот пример должен быть в 100 км к северу и востоку от начальной точки с юго-запада, но когда я пробиваю координаты в гугл-картах, это в некотором другом месте.Не другая сторона мира неправильна, но явно слишком неточна, чтобы быть полезной.