Я застрял в том, что я думаю, что проблема не в том, чтобы найти правильный проекционный вход для моего государственного самолета (штат Аляска, плоскость плоскости 4 NAD83 Feet) до использование с функцией sp::CRS.Преобразование координат состояния плоскости (Аляска 4) в широту и долготу
Я следую инструкциям для преобразования из состояния самолета в: Convert latitude/longitude to state plane coordinates
и я посмотрел в АСБ, который взял меня в http://trac.osgeo.org/proj/ из Ссылки, но я даже не могу сказать, если этот сайт? у меня будет то, что мне нужно.
Я использую учебник Hadley Уикхемы для построения шейпфайлов на: https://github.com/hadley/ggplot2/wiki/plotting-polygon-shapefiles
И я могу построить любого из шейпфайлы мне нравится из: http://munimaps.muni.org/moagis/download.htm
Моего код:
require(rgdal)
require(maptools)
require(ggplot2)
require(plyr)
my_dsn <- "directory where the shapefile is"
Assembly = readOGR(dsn=my_dsn, layer="assembly")
[email protected]$id = rownames([email protected])
Assembly.points = fortify(Assembly, region="id")
Assembly.df = join(Assembly.points, [email protected], by="id")
[email protected]$id = rownames([email protected])
Assembly.points = fortify(Assembly, region="id")
Assembly.df = join(Assembly.points, [email protected], by="id")
#Assembly.df$DISTRICT <- factor(Assembly.df$DISTRICT)
ggplot(Assembly.df) +
aes(long,lat,group=group) +
geom_path(color="black") +
coord_equal()
Пожалуйста, помогите мне найти входные данные, необходимые для преобразования прогнозов. Я новичок в работе с картографическими проекциями.
Что proj4string (Assembly) непосредственно после строки readOGR? Если нет NA, вы можете пойти с? SpTransform, все, что другие вещи бесполезны – mdsumner
Выход «+ proj = tmerc + lat_0 = 54 + lon_0 = -150 + k = 0.9999 + x_0 = 500000.0000000001 + y_0 = 0 + datum = NAD83 + units = us-ft + no_defs + ellps = GRS80 + towgs84 = 0,0,0 " – cylondude
Я предполагаю, что googlemaps использует проекцию Меркатора. Как я могу найти вход для CRS() для этой проекции? – cylondude