classにハイホンがあるとうっとうしいので，Iris-以降を改めてclassとしてiris2を作っておく。 data iris2; set iris; class=SCAN(class,2,"-"); put class=; proc print data=iris2; run; setosa(class=set)だけのクラスを作りたければ，こうすればいいが，…

Irisデータ*1を読み込む data iris; infile 'I:\SAS\iris.csv' dlm=','; input sepalL sepalW petalL petalW class$; proc print data=iris; run; *1:http://archive.ics.uci.edu/ml/datasets/Iris

create table transkanto (sec int, aac character(5) primary key); copy transkanto from '/var/lib/postgresql/pgis/transkanto_sort_city_wi.csv' with csv; create table kantoswi as select * from kantowi4 inner join transkanto using (aac); creat…

create table p08 as select aac, ST_Union(the_geom) from "08" group by aac order by aac; create table p09 as select aac, ST_Union(the_geom) from "09" group by aac order by aac; create table p10 as select aac, ST_Union(the_geom) from "10" gr…

shp2pgsql -s 4612 -D -i -I -W cp932 N03-070401_08_EC01.shp "08" > hoge.sql shp2pgsql -s 4612 -D -i -I -W cp932 N03-070401_09_EC01.shp "09" >> hoge.sql shp2pgsql -s 4612 -D -i -I -W cp932 N03-070401_10_EC01.shp "10" >> hoge.sql shp2pgsql -s…

ncの計算がおかしかったので修正 ０始まりの市区町村コードがあるので，sec[i][j]の型をto_iからto_sに変更 lines = nc = nctmp = 0 sec = open("seckantowi2.csv") {|file| while l = file.gets sec[lines] = l.split(',') lines += 1 nctmp = l.split(',')…

seckantowi2.csv

読み込むテーブル名（７，１０，２４行目）と，アウトプットのファイル名（２カ所）を変更しただけ． #! /usr/bin/ruby1.9.1 # -*- coding: utf-8 -*- require "pg" #hogehoge = username #barbar = pass word of hogehoge #conn = PGconn.connect( "localho…

create table p13secwi as select sec as "aac", ST_Union(st_union) from p13swi group by sec order by sec; create table code13secwi as select aac from p13secwi;

create table trans13 (sec int, aac character(5) primary key); copy trans13 from '/var/lib/postgresql/pgis/trans13_sort_city_wi.csv' with csv; create table p13swi as select * from p13 inner join trans13 using (aac);

あまりに工夫の少ないRubyのコード(nb_c13.cgi)で恐縮です． PostGISの設定は例えばここ*1を参照． #! /usr/bin/ruby1.9.1 # -*- coding: utf-8 -*- require "pg" #hogehoge = username #barbar = pass word of hogehoge #conn = PGconn.connect( "localhost…

ST_Intersects*1 — ジオメトリ/ジオグラフィが「空間的にインタセクトする」(空間に共有部分がある)場合には，TRUEを返します．そうでない(非接続)場合はFALSEを返します．ジオグラフィ -- 許容量は0.00001メートルです(どの閉じた点でもインタセクトを考慮…

ST_Union*1 — ジオメトリ結合の点集合を表現するジオメトリを返します． テーブル"13"からaacでグルーピングして，the_geomを結合する．それのうち，aac（市区町村コード）とST_Union（ジオメトリ）を選んで，aacで並べ替えて新しいテーブルp13を作る． crea…

島しょ地域は， ・距離がずいぶん離れている（集積性を出すので） ・傾向が違っているだろうと思われる ので，別途解析するとして，以下のように，最後の番号は外しておこう lines = nc = 0 sec = open("sec13.csv") {|file| while l = file.gets sec[lines]…

本記事では，「国土数値情報（行政区域（面）） 国土交通省」のデータを利用しています．国土数値情報ダウンロードサービス*1より， （JPGIS準拠）データのうち， 国土骨格の ● ベクトルデータ 行政区域（面）をダウンロードする．二次医療圏と合わせたいの…

lines = nc = 0 sec = open("sec13.csv") {|file| while l = file.gets sec[lines] = l.split(',') lines += 1 nc = l.split(',').size end } foo = File.open("trans13.csv",'w') ii = 0 out = for i in 0...(lines) for j in 1...nc if (sec[i][j]=="")||(…

平成22年3月末時点での二次医療圏 (PDF)*1 より，二次医療圏コード，2列目以降は市区町村コード のCVSファイル（sec13.csv）を作成しておく 1301,13101,13102,13103,13105,13106,,,, 1302,13109,13111,,,,,,, 1303,13110,13112,13113,,,,,, 1304,13104,13114…

初めての投稿です。 #コードのテスト for(i in 1:1000000000){cat("hoge\n")} HTMLのテスト <pre>#コードのテスト for(i in 1:1000000000){cat("hoge\n")} </pre>を記事で書くと，#コードのテスト for(i in 1:1000000000){cat("hoge\n")}