後でフュージョン テーブルまたは他の GIS アプリケーションにインポートする R でシェープファイルを作成しようとしています。
まず、カナダのすべての国勢調査区を含む空のシェープファイルをインポートしました。CT の一意の ID に基づいてシェープファイルに他のデータ (表形式) を添付し、結果をマッピングしました。現時点では、Vancouver のものだけが必要で、Vancouver の CT と新しくアタッチされた属性データのみを含むシェープファイルをエクスポートしたいと考えています。
これが私のコードです(プライバシー上の理由により一部省略されています):
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
これを使用した他の例を見てきました: writePolyShapeでシェープファイルを作成します。試してみましたが、ある程度効きました。.shp、.dbf、および .shx ファイルが作成されました。.prj ファイルがありません。作成方法がわかりません。シェープファイルを作成するためのより良い方法はありますか?
この問題に関する助けをいただければ幸いです。