BigQueryにShapeFileを取りこむ のバックアップソース(No.2)

• バックアップ一覧
• 差分 を表示
• 現在との差分 を表示
• バックアップ を表示
• BigQueryにShapeFileを取りこむ へ行く。
  • 1 (2021-05-03 (月) 10:16:19)
  • 2 (2021-05-03 (月) 10:51:39)
  • 3 (2021-05-03 (月) 14:08:58)
  • 4 (2021-05-04 (火) 10:26:00)
  • 5 (2021-05-04 (火) 14:46:09)
  • 6 (2021-05-04 (火) 17:20:06)
  • 7 (2021-05-09 (日) 09:57:12)

#author("2021-05-03T10:51:39+09:00","","")
*BigQueryにShpeFileを取り込む [#ma330452]

 ### ShpeFileからBigQuery用のWKTを作る
 # 目的:Google datapotalでの行政界ポリゴンがひどすぎる(超きたない）から自分で作る
 # ESRI市区町村界ShapeFileを使っているが、この方法では規定違反になるので取り扱い注意
 # 国土地理院の行政区データからのプログラムは後でやる
 library(sf)
 x <- st_read("japan_ver821/japan_ver821.shp")
 x2 <- x$geometry
 x3 <- st_as_text(x2)
 z2 <- x$JCODE
 z3 <- as.data.frame(z2)
 df1 <- cbind(z3,x3)
 df2 <- na.omit(df1)
 write.table(df2,"japan_ver821.csv",row.names = F,col.names = F, sep = ",")
 ## GCPでの後処理
 # csvをGCStregeで読み込んでからBigQueryでテーブル作成
 #
 # BigQueryへ読み込むスキーマ
 # JCODE	STRING	NULLABLE	
 # WKT_TEXT	STRING	NULLABLE
 # 一旦WKT_TEXTをSTRINGで読み込んでから、クエリでWKTに変換すること。直接やるとエラーになるため（後述）
 #
 # WKTの変換クエリ
 # SELECT *,ST_GEOGFROMTEXT(WKT_TEXT,make_valid => TRUE) AS WKT FROM `XXXX.XXXX` 
 # make_valid => TRUE　でポリゴンがループしている部分のエラーを修正

*コメント欄 [#t9b5bf47]
#comment(below)
- これを使うとデータポータルで美しい行政界地図が表示できます。もちろんJCODE付きです。 -- [[Rokinawa]] &new{2021-05-03 (月) 10:51:39};

&counter;