Próbuję utworzyć shapefile w R, który później zaimportuję do Fusion Table lub innej aplikacji GIS.

Na początek zaimportowałem pusty shapefile zawierający wszystkie traktaty spisowe w Kanadzie. Załączam inne dane (w formacie tabelarycznym) do shapefile na podstawie unikalnego identyfikatora CT i mapowałem moje wyniki. W tej chwili potrzebuję tylko tych w Vancouver i chciałbym wyeksportować plik shapefile zawierający tylko CT Vancouver oraz moje nowo dołączone dane atrybutów.

Oto mój kod (niektóre części pominięto ze względu na prywatność):

shape <- readShapePoly('C:/TEST/blank_ct.shp') #Load blank shapefile
shape@data = data.frame(shape@data, data2[match(shape@data$CTUID, data2$CTUID),]) #data2 is my created attributes that I'm attaching to blank file
shape1 <-shape[shape$CMAUID == 933,] #selecting the Vancouver CTs

Widziałem inne przykłady używające tego:writePolyShape aby utworzyć plik kształtu. Próbowałem tego i działało to do pewnego stopnia. Utworzono pliki .shp, .dbf i .shx. Brakuje mi pliku .prj i nie wiem, jak go utworzyć. Czy istnieją lepsze metody tworzenia plików shapefiles?

Jakakolwiek pomoc w tej sprawie byłaby bardzo mile widziana.