Справка, отображающая Географические Данные на графике в R, использующем PBSMapping и Файлы форм
Если Вы находитесь на среде Windows тогда существует, вероятно, мало, можно сделать, поскольку направляющие, кажется, запускаются действительно медленно в соответствии с Windows. Я имел тот же опыт в Windows и должен был переместить мою установку в VM Linux для создания его действительно быстрым (я также использовал автотест).
2 ответа
Вы также можете посмотреть на эти связанные вопросы, особенно на ответы Эдуардо:
Похоже, что PBSmapping использует грубую эвристику для расчета проекции из файла .prj. (см. справку (importShapefile)). Я лично не понимаю всего, что находится внутри файла prj, но, используя этот веб-сайт www.spatialreference.org, я считаю, что ваша карта соответствует
http://www.spatialreference.org/ref/epsg/26912/
Каждый раз, когда Я получаю новый файл формы, нахожу его систему проецирования на этом веб-сайте, а затем ищу строку proj4, которая в данном случае "+ proj = utm + zone = 12 + ellps = GRS80 + datum = NAD83 + units = m + no_defs"
(Как я уже сказал, я не знаю PBSmapping, но вы можете прочитать это, используя maptools следующим образом)
library(maptools)
sf=readShapeSpatial("SGID93_DEMOGRAPHIC_CensusTracts2000.shp",proj4string=CRS("+proj=utm +zone=12 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"))
, а затем преобразовать в латлонги с использованием
library(rgdal)
sftransformed=spTransform(sf,CRS("+proj=longlat"))
и
plot (sftransformed, axes = T)
дает график с правильными единицами измерения по осям.
Не уверены, понимает ли PBSmapping строку proj4? Похоже, это неправда.