Urbgeo/gisdb ver 0.0

Fra hbrgeo
Revisjon per 22. nov. 2016 kl. 18:42 av Hhiso (diskusjon | bidrag)

Etablering av UrbGeo GISdatabase

Versjon 0.0, 22. november 2016.

Dette dokumentet viser steg-for-steg dokumentasjon av GIS-databasen urbgeo med DBMS PostgreSQL og overføring til ArcGIS server og etablering av enkel web-applikasjon basert på data og geodata. Dersom du berre er interessert i å laste ned geodata gå til siste hovedpunkt i dokumentasjonen,"Import av data".

Forusetninger

Urbgeo er laget ved hjelp av Postgresql. Forutsetningen er at PostgreSQL 9.x med tilhørende GIS-ekstensjon postgis er installert på ein linux-maskin. Denne versjonen av dokumentasjonen er basert på Ubuntu Linux 14.04 LTS, postgreSQL 9.5 og postgis 2.2.2. I tillegg er det nødvendig å installere pg_similarity for å kunne bruke meir avanserte streng-samanlikning. Det er også en forutsetning at database-brukeren har privilegier til å opprette databaser. Videre blir alle filer lagret i katalogen /home/[brukernamn]/urbgeo med underkataloger.


Denne versjonen av urbgeo bruker også geonorge.no og BAAT-konto er nødvendig både for å laste ned adressedata og for å gjere oppslag mot adressedata hos geonorge. Det er mulig å lage opne alternativ til både nedlasting og adressesøk, f.eks. å bruke de opne tenestene hos geonorge eller andre leverandører av adressedata for Norge.

Lage database

Opprette database frå linux prompt
~$cd
~$mkdir urbgeo
~$cd urbgeo
~$createdb urbgeo

Logge på database urbgeo frå linux prompt
~$psql -d urbgeo

Prompt i psql er databasenamn og =#

Laste postgis ekstensjon

urbgeo=#create extension postgis;


For å sjekke at alt er i orden kan du logge på databasen med psql -d urbgeo vil prompten vere urbgeo=#

Logg av med \q

Lasting av hjelpetabellar

Første steg blir å laste databasen med hjelpetabellar som i all hovudsak inneheld data om den administrative inndelinga, namn på byar og kjeldeidentifikasjon. Utgangspunkt er data frå NSD (Norsk Samfunnsvitenskaplig Datatjenste), Arkivverket og Hallstein Myklebost, Norges tettbygde steder.

Nedanfor er script for å opprette databaser og laste inn data frå dei ulike formata

Opprette database

-- Om du må lage helt ny database, bruk de to neste kommandoene
--
--DROP DATABASE IF EXISTS "urbgeo";
--CREATE DATABASE "urbgeo";

\c "urbgeo";

Tettsted og by

Første hjelpetabell er definisjon og data for tettsteder. Data er henta frå Hallstein Myklebust, Tettsteder i Norge 1875-1950 og lagret i et csv-dokument med ";" som kolonne-skille. Myklebosts tettstadnummer er nytta som primærnøkkel for byane i urbgeo. HMID=Hallstein Myklebust ID. Merk at encoding er sett til Latin1 sidan csv-fila er laga med MS-Excel under Windows 10.


--
-- Table structure for table 'hallmykl'
--


DROP TABLE IF EXISTS "hallmykl";

CREATE TABLE "hallmykl" (
  "hmid" INTEGER NOT NULL, 
  "fylke" VARCHAR(25), 
  "tettstad" VARCHAR(25), 
  "kommune" VARCHAR(25), 
  "f1875" INTEGER, 
  "f1890" INTEGER, 
  "f1900" INTEGER, 
  "f1910" INTEGER, 
  "f1920" INTEGER, 
  "f1930" INTEGER, 
  "f1946" INTEGER, 
  "f1959" INTEGER, 
  "hmtype" VARCHAR(50), 
  "hmkode" VARCHAR(50), 
  "kommuner" VARCHAR(255), 
  "asciiname" VARCHAR(25), 
  currknr integer,
  PRIMARY KEY (hmid)
);

COPY hallmykl FROM '/home/arne/urbgeo/data/hmexp.csv' DELIMITER ';' CSV HEADER QUOTE '"' ENCODING 'LATIN1' NULL '';

update "hallmykl" set kommune = 'Mo i Rana' where asciiname='Mo' or currknr=1833;
update "hallmykl" set asciiname = 'Mo i Rana' where asciiname='Mo' or currknr=1833;
alter table hallmykl add column url varchar(255);

Merk at noen av tettstedene i perioder er del av en kommune og var ikke selvstendig bykommune. Noen av tettstedene hadde også så liten befolkning at de vil falle bort fra de fleste definisjoner av "by". Og så er det noen av tettstedene som fikk bystatus etter 1910.

Utdrag/dump av tabellen vis hallmykl

Kommuner og fylker

Hovedtabell er nsd1910a

Neste datafil er kommunenefila for 1910. Denne er trekt ut frå kommunedatabasen i NSD. Kommunedatabasen har muligheit for å lage kommunefiler for eit tidspunkt (år). Den aktuelle fila er for år 1910. Lista frå NSD-kommunedatabase som inneheld nsd-kommunenr (nsdknr) og nsd-kommunenamn (nsdknam) er overført (copy-paste) til MS-Excel. I MS-Excel er kommunenr brukt for å leggje til fylkesnummer. Fylkesnummer er danna ved funksjonen avrund.ned(nsdknr/100;0). I tillegg er det konstruert ein fjerde kolonne som viser om kommuna er by eller landsbygd. Dette feltet er konstruert med funksjonen: Hvis(nsdknr-(fylkenr*100)<10;1;0). Dvs. dersom nsdknr - fylkesnr*100 er mindre enn 10, så er kommuna ein bykommune.

-- Oppretter kommunefil for Norge i 1910

DROP TABLE IF EXISTS "nsd1910a";

CREATE TABLE "nsd1910a" (
  "nsdknr" INTEGER NOT NULL, 
  "nsdknamn" VARCHAR(25), 
  "fylke" INTEGER, 
  "byland" INTEGER,  
   PRIMARY KEY (nsdknr) 
);

COPY nsd1910a FROM '/home/arne/urbgeo/data/nsdexp10.csv' DELIMITER ';' CSV HEADER QUOTE '"' ENCODING 'LATIN1' NULL '';

Alter table nsd1910a add column asciiname varchar(25);
Alter table nsd1910a add column hmid integer;
alter table nsd1910a add column url varchar(255);

-- Genererer asciiname uten nordiske teikn og standardiserer namn til gjeldande namn (Kristiania til Oslo)

update "nsd1910a" set asciiname = replace(nsdknamn,'å','aa');
update "nsd1910a" set asciiname = replace(asciiname,'Å','Aa');
update "nsd1910a" set asciiname = replace(asciiname,'æ','e');
update "nsd1910a" set asciiname = replace(asciiname,'ø','o');
update "nsd1910a" set asciiname = replace(asciiname,'Ø','O');
update "nsd1910a" set asciiname = replace(asciiname,'Æ','E');
update "nsd1910a" set asciiname = replace(asciiname,'Kristiania','Oslo');
update "nsd1910a" set asciiname = replace(asciiname,'Fredrikshald','Halden');
update "nsd1910a" set asciiname = 'Mo i Rana' where nsdknr=1833;

-- Sjekk av koblingene og standardisering

select count(*) from nsd1910a, hallmykl where nsd1910a.nsdknamn=hallmykl.tettstad;
select count(*) from nsd1910a, hallmykl where nsd1910a.asciiname=hallmykl.asciiname;

-- Setter hmid i nsd1910a fila 
update nsd1910a set hmid=hallmykl.hmid from hallmykl where nsd1910a.asciiname=hallmykl.asciiname;

-- Kontroll av tettsteder som ikke bystatus i 1910

select hallmykl.hmid, hallmykl.tettstad, nsd1910a.nsdknamn from hallmykl left join nsd1910a on hallmykl.asciiname=nsd1910a.asciiname where nsd1910a.hmid is null;

Utdrag/dump av tabellen vis nsd1910a

Tellekretser og bytabell i 1910


-- Oppretter tabell for teljekrets (tk)

DROP TABLE IF EXISTS "tk1910";

CREATE TABLE "tk1910" (
  "id" VARCHAR(25) NOT NULL, 
  "kildeid" INTEGER NOT NULL, 
  "kretsnr" VARCHAR(25) NOT NULL, 
  "tk_navn" VARCHAR(255), 
  "sogn" VARCHAR(50), 
  "prestegjeld" VARCHAR(50), 
  "herred_by" VARCHAR(50), 
  "merknader" VARCHAR(255),
  PRIMARY KEY (kildeid, kretsnr)
);

COPY tk1910 FROM '/home/arne/urbgeo/data/ft1910heleKrets.csv' DELIMITER ';' CSV HEADER QUOTE '"' ENCODING 'utf-8' NULL '';

-- Legg til ekstra felt id ed

Alter table tk1910 add column asciiname varchar(25);
alter table tk1910 add column byland INTEGER;
alter table tk1910 add column url varchar(255);

-- Standardiserer asciiname

update tk1910 set asciiname = replace(herred_by,'å','aa');
update tk1910 set asciiname = replace(asciiname,'Å','Aa');
update tk1910 set asciiname = replace(asciiname,'æ','e');
update tk1910 set asciiname = replace(asciiname,'ø','o');
update tk1910 set asciiname = replace(asciiname,'Ø','O');
update tk1910 set asciiname = replace(asciiname,'Æ','E');
update tk1910 set asciiname = replace(asciiname,'Kristiania','Oslo');
update tk1910 set asciiname = replace(asciiname,'Fredrikshald','Halden');
update tk1910 set asciiname = replace(asciiname,'Trondhjem','Trondheim');
update tk1910 set asciiname = replace(asciiname,'Porsgrund','Porsgrunn');
update tk1910 set asciiname = replace(asciiname,'Stenkjer','Steinkjer');
update tk1910 set asciiname = 'Sokndal' where prestegjeld is null and herred_by='Sogndal';

-- Merker teljekretser som er i by

update tk1910 set byland=1 from hallmykl  where tk1910.asciiname=hallmykl.asciiname;

-- Oppretter bytabell for 1900-teljinga

DROP TABLE IF EXISTS "by1910";

CREATE TABLE "by1910" (
  "hmid" INTEGER NOT NULL, 
  "bynamn" VARCHAR(25), 
  "nsdknr" INTEGER, 
  "asciiname" VARCHAR(25),
  "tkid" varchar(25)
);


-- laster bytabellen frå teljekretsar

insert into by1910 (hmid,bynamn,nsdknr,asciiname,tkid) select count(*), herred_by, count(*), min(asciiname), min(id) from tk1910 where byland=1 group by herred_by;

update by1910 set hmid=hallmykl.hmid from hallmykl where by1910.asciiname=hallmykl.asciiname;
update by1910 set nsdknr=nsd1910a.nsdknr from nsd1910a where by1910.asciiname=nsd1910a.asciiname;

alter table by1910 add  primary key (hmid);
alter table by1910 add column url varchar(255);


-- oppretter kildetabell frå teljekretser

DROP TABLE IF EXISTS "kilde1910";

CREATE TABLE "kilde1910" (
  "hmid" INTEGER, 
  "kildeid" INTEGER NOT NULL, 
  "asciiname" VARCHAR(25),
  PRIMARY KEY (kildeid) 
);

insert into kilde1910 (hmid,kildeid,asciiname) select count(*), kildeid, min(asciiname) from tk1910 where byland=1 group by kildeid;

-- Merker aller kilder med tettstadid (hmid)

update kilde1910 set hmid=hallmykl.hmid from hallmykl where kilde1910.asciiname=hallmykl.asciiname;


Folketellingen 1900

Utgår i denne versjon

Metadata

For framtidig bruk er det opprettet en metatabell for kilder.

DROP TABLE IF EXISTS "ub_kilde";
create table "ub_kilde" (
  id varchar(6) not null,
  startaar integer,
  sluttaar integer,
  dekning integer,
  filnamn varchar(64),
  url varchar(255),
  primary key (id) 
);

-- id: FTxxxx = folketelling, BGxxxx=Gravlagde, DP=Dopte, VI=Vigde
-- dekning 1=Norge, 2=fylke, 3=kommune, 4=tettstad

insert into ub_kilde values ('FT1910',1910,1910,1,'ft1910heleby.csv', 'http://digitalarkivet.arkivverket.no/ft/sok/1910'); 

Bosted og gatetabeller

Bosted og gatetabeller blir lastet fra FT1910. Hovedtabellen er bost1910.


DROP TABLE IF EXISTS "tmptab";
CREATE TABLE "tmptab" (
  autoid integer not null,
  id varchar(25),
  overid varchar(25),
  kildeid integer,
  kretsnr varchar(25),
  bostnr varchar(25),
  gateognr varchar(255),
  bydel varchar(50),
  gaardnr varchar(50),
  gardeier varchar(255),
  ant_tilst varchar(50),
  ant_hj varchar(50),
  ant_pp varchar(50),
  hustype varchar(255),
  sideb varchar(255),
  ubebodd varchar(25),
  ant_et varchar(25),
  ant_leil varchar(255),
  ant_naer varchar(255),
  ant_verk varchar(255),
  ant_ledig varchar(50),
  ant_kvist varchar(50),
  ant_kjeller varchar(50),
  merknad varchar(255), 
  primary key (autoid)
);

COPY tmptab FROM '/home/arne/urbgeo/data/ft1910heleby.csv' DELIMITER ';' CSV HEADER QUOTE '"' ENCODING 'utf-8' NULL '';

drop table if exists bost1910;
create table bost1910 (
    autoid integer not null,
    bostid varchar(25),
    overid varchar(25),
    kildeid integer,
    hmid integer, 
    gateognr varchar(255),
    bydel varchar(50),
    gaardnr varchar(50),
    gatenamn varchar(50),
    gatenr varchar(50),
    gatenum integer,
    gatelitr varchar(50),
    url varchar(255),
    primary key (autoid)
);

INSERT INTO bost1910 (autoid, bostid, overid, kildeid, gateognr, bydel, gaardnr) select autoid, id, overid, kildeid, gateognr, bydel, gaardnr from tmptab where kildeid is not null;

update bost1910 set hmid=kilde1910.hmid from kilde1910 where kilde1910.kildeid=bost1910.kildeid;
select * from bost1910 where hmid is null;

update bost1910 set gatenr = trim((regexp_split_to_array (gateognr, E'\\D+'))[2]);
update bost1910 set gatenamn = trim((regexp_split_to_array (gateognr, gatenr))[1]) where gatenr is not null;
update bost1910 set gatenamn = gateognr where gatenr is null;
update bost1910 set gatenamn = gateognr where gatenr = '';
update bost1910 set gatenr=null where gatenr='';
update bost1910 set gatelitr = trim((regexp_split_to_array (gateognr, gatenr))[2]) where gatenr is not null; 
update bost1910 set gatenum = to_number(gatenr,'999') where gatenr is not null;

select bost1910.hmid, hallmykl.tettstad, count(*) as antgrd from bost1910, hallmykl where ((gaardnr is not null) and (hallmykl.hmid=bost1910.hmid)) group by bost1910.hmid, tettstad ;


drop table if exists gatetab;
create table gatetab (
    hmid integer not null, 
    gatenamn varchar(50) not null,
    anthus integer,
    gateognr varchar(255),
    gatenr varchar(50),
    mingnum integer,
    maxgnum integer,
    aar integer,
    jwval double precision,
    std_gate varchar(50),
    primary key (hmid, gatenamn)
);

insert into gatetab (hmid, gatenamn, anthus, gateognr, gatenr, mingnum, maxgnum,aar) select hmid, gatenamn, count(*) as anthus, min(gateognr), min(gatenr), min(gatenum), max(gatenum), 1910 from bost1910 where gatenamn is not null group by hmid, gatenamn order by hmid, gatenamn;


Utdrag/dump av tabellen vis gtab1910

Geokoding

Geokoding av stedsdata skjer på to måter, a) Stedsnavn for byer, fylker og kommuner blir geokodet ved oppslag mot ein eldre rutine hos kartverket og b) Gateadresser blir geokodet mot data lasta ned frå geonorge.

Geokoding av steder

Landtabellen, fylkestabell og bytabellen geokoda for å få punkt for landa i Norden, dei norske fylka, kommunene i 1910 og tettstadane i urbgeo. Geokoding av land og fylker utgår fra denne versjonen. Geokoding av tettsteder (byene) og kommuner skjer ved bruk av ein enkel php-rutine med oppslag mot geonorge.

Tenesten dokumentert hos geonorge.no

Capabilities for tjenesten https://ws.geonorge.no/SKWS2/services/SSR?wsdl

For å transformere koordinatsystema brukt av kartverket er biblioteket proj4 benyttet. Implementasjonen for php er lastet ned fra github og installert under lokal apache2 server (/var/www/html/vendor/proj4)

Kartverkets dokumentsjon av koder for koordinatsystem brukt av rutinene: SOSI koordinatsystem

<?php

$klient = new SOAPClient('https://ws.geonorge.no/SKWS2/services/SSR?wsdl');
$stedsnavn = "Hvitsten"; 
$koord="23";
$EPSG='EPSG:25833';

//Her treng du din BAAT-bruker og passord.
$data = array(
    'brukerid' => "[BAAT-BRUKER]",
    'passord' => "[PASSORD]",
    'aliasId' => "",
    'navn' => $stedsnavn,
    'fylkeKommuneListe' => "",
    'nordLL' => 0,
    'ostLL' => 0,
    'nordUR' => 0,
    'ostUR' => 0,
    'maxAnt' => 100,
    'sosiKoordSyst' => $koord
);

include("proj4php/vendor/autoload.php");

use proj4php\Proj4php;
use proj4php\Proj;
use proj4php\Point;

// Initialise Proj4
$proj4 = new Proj4php();

// Create two different projections.
$projutm    = new Proj($EPSG, $proj4);
$projWGS84  = new Proj('EPSG:4326', $proj4);

header("Content-Type: text/html;charset=utf-8");
$table="hallmykl";
$col=2;
$level=5;
if (isset($_GET['tabname'])) {
   $table = $_GET['tabname'];
}
if (isset($_GET['colno'])) {
   $col = $_GET['colno'];
}
if ($table == 'hallmykl') {
   $col=2;
   $level=4;
   $colnr=16; 
   $colk=3; 
} else {
   $col=1;
   $colk=1;
   $level=4;
   $colnr=0; 
}
include 'htmlhead.php';
print ("<body><pre>");
require('../config.php');
include 'db/opendb.php';

$strSQL="SELECT * FROM " . $table . ";";
//$strSQL="SELECT * FROM " . $table . ";";
$results = pg_query($con, $strSQL);
if ($results)  {
     while ($row = pg_fetch_row($results)) {
//        print ($row[$col] . "\t" . $row[$colnr] . "\t" . $row[$colk] . "\n");
        $data['navn']=$row[$col];
        $resultat = $klient->__soapCall('ssrSok', $data);
        if (is_null($resultat->ssrRecords[0])) {
//          print_r ($resultat); 
        } else {
         foreach ($resultat->ssrRecords as $adr) {
           if (intval($adr->for_komm) == intval($row[$colnr]) or strcmp(trim($adr->kom_navn),trim($row[$colk]))==0) {
   	      $svar=$adr->for_snavn; 
 	      $snord=$adr->for_nord1; 
 	      $sost=$adr->for_aust1; 
		// Create a point.
   	      $pointSrc = new Point($sost, $snord, $projutm);
              $pointDest = $proj4->transform($projWGS84, $pointSrc);
              $dlng=$pointDest->x;
              $dlat=$pointDest->y;
              $strISQL = "insert into ub_point_0 (loc_name,loc_type,nord,ost, lng,lat) values ('" . $row[$col] . "',$level," . $snord . "," . $sost . "," . $dlng . "," . $dlat .");\n";
              print   $strISQL;
              break;
           } else {
//              print ("\t" . $adr->for_komm .  "\t" . $adr->kom_navn . "\t" . $row[$col] . "\n");
           } 
         }
        } 
     }	
}
print ("

");

include 'htmlend.php'; ?>


Geokoding av gateadresser

Gatetabellen bost1910 blir "vasket" mot data lastet ned fra geonorge.

Koding av data med geodata

Generering av punkt-tabell

Oppretting av temporær tabell for punkt

drop table if exists ub_point_0;
create table ub_point_0 (gid serial,
"loc_name" varchar(50),
"loc_type" integer,
"x" numeric,
"y" numeric,
hbrgeoid varchar(18));

ALTER TABLE "ub_point_0" ADD PRIMARY KEY (gid);
SELECT AddGeometryColumn('','ub_point_0','geom','0','POINT',2);


Import av punkt fra land, fylke og gateadresser fra geokoding county - punkt for norske fylker country - punkt for land i Norden geo1910u - punkt for gateadresser geokhm - punkt for tettstader

drop table county;
drop table country;
drop table geo1910u;
drop table geokhm;

\i utm33/country.sql;
\i utm33/county.sql;
\i utm33/geo1910u.sql;
\i utm33/geokhm.sql;

Innsetting av punkt-data frå geo-koding

Type-koder; 0=ingen 1=land 2=fylke 3=kommune 4=tettstad (by) 5= 6= 7=gateadresse 8=hus 9=Anna


-- Setter inn fylke i temporær punkt-tabell med loc_type = 2 for fylke

insert into ub_point_0 (loc_name, loc_type, hmid, x, y, geom) select match_addr, 2, 0, x, y, geom from county;

-- Setter inn land i temporær punkt-tabell med loc_type = 1 for land

insert into ub_point_0 (loc_name, loc_type, hmid, x, y, geom) select match_addr, 1, 0, x, y, geom from country;

-- Setter inn tettstad i temporær punkt-tabell med loc_type = 4 for tettstad

insert into ub_point_0 (loc_name, loc_type, hmid, x, y, geom) select tettstad, 4, hmid, x, y, geom from geokhm;

-- Setter inn gater i temporær punkt-tabell med loc_type = 7 for gate

insert into ub_point_0 (loc_name, loc_type, hmid, x, y, geom) select gateognr, 7, hmid, x, y, geom from geo1910u;

Tildeling av hbrgeo-id


-- oppdaterer hbrgeoid i ub_point_0 på formatet 'nnExxxxxxxNyyyyyyy'.
-- Merk at string concat i postgresql er ||

-- select trim(to_char(loc_type,'9')) || 'E' || trim(to_char(x*100000,'0000000')) || 'N' || trim(to_char(y*100000,'0000000')) from ub_point_0 where x<32 and y<80 and x>5 and y>58;

update ub_point_0 set hbrgeoid='1' || trim(to_char(loc_type,'9')) || 'E' || trim(to_char(x*100000,'0000000')) || 'N' || trim(to_char(y*100000,'0000000'))  where x<32 and y<80 and x>5 and y>58;

Etablering av unik tabell for hbrgeoid med oppretting av tabell for unike punkt

drop table if exists ub_point;
create table ub_point (gid serial,
"loc_name" varchar(50),
"loc_type" integer,
"hmid" integer,
"x" numeric,
"y" numeric,
hbrgeoid varchar(18));

ALTER TABLE "ub_point" ADD PRIMARY KEY (hbrgeoid);
SELECT AddGeometryColumn('','ub_point','geom','0','POINT',2);

Etablering av unike punkt

Initiering av tabell med unike lokasjoner for hbgrgeoid



insert into ub_point (loc_name, loc_type, hmid, x, y, hbrgeoid, geom) select min(loc_name), min(loc_type), min (hmid), min(x), min(y), hbrgeoid, min(geom) from ub_point_0 where hbrgeoid is not null group by hbrgeoid order by hbrgeoid;

Koding av GIS data

GIS data består av to ulike datasett. Beskrivande data koda med geoid og sjølve gisdata: Punkt, linjer og polygon. Koplinga mellom data og gisdata er gjort med geoid.

Koding av data med geoid


-- Sette inn geokode i tabellen for bustader

Alter table gtab1910 add column hbrgeoid varchar(18);

update gtab1910 set hbrgeoid = ub_point.hbrgeoid from ub_point where ub_point.hmid=gtab1910.hmid and ub_point.loc_name=gtab1910.gateognr;

Export av Geo data

GIS data består av to ulike datasett. Beskrivande data (Analysedata) koda med geoid og sjølve gisdata: Punkt, linjer og polygon. Data koda med geoid (analysedata) og geodata blir eksportert kvar for seg.

Eksport av 1900 bosteder

For å kople bostedene i folketeljinga 1910 til geoid blir ein fil med hbrgeoid og primærnøkkelen for ft1910heleby exportert til en csv-fil. Feltet autoid er primærnøkkel for ft1910heleby.csv

COPY gtab1910 (autoid, hbrgeoid, hmid) to '/home/arne/urbgeo/export/bost1910.csv' DELIMITER ';' CSV HEADER QUOTE '"' ENCODING 'Latin1' NULL '';

Denne filen kan også lastast ned her: bost1910.csv Det kan vere lurt å laste denne fila i MS-Excel og konvertere den til MS-Excel 2003-format (*.xls) før du importerer den i ArcGIS.

Eksport av geodata

Punkt-data kan eksporterast enten som dump av tabellen ub_point eller som shapefile. Dump-fila kan importerast til ein annan postgresql (og andre SQL-databaser med postgis-ekstensjon.

Shapefila kan importerast både til desktop program som ArcGIS desktop eller QGIS. Den kan også importerast til ArcGIS Server

Eksport ved shapefile

Eksport av geodata med shapefile blir utført med konverteringsprogrammet pgsql2shp (Windows, Linux).

linux>pgsql2shp urbgeo ub_point
linux>zip -r ub_pointz.zip ub_point.* 

Programmet pgsql2shp vil lage filene ub_point.shp, ub_point.dbf, ub_point.shx, ub_point.prj og ub_point.cfg. Alle fem filene må flyttast til ArcGIS-området for shapefiler.

Last ned ub_pointz.zip ub_pointz.zip

Eksport ved database dump

Eksport av geodata (punkt, linje, polygon) for import i postgresql. ub_point har unike punkt og har primærnøkkel hbrgeoid

linux>pg_dump -t ub_point urbgeo -f urbgeo1.sql

Eksport av hjelpetabeller og datatabeller for import i postgresql

pg_dump -t hallmykl -t nsd1910a -t gtab1910 urbgeo -f urbgeo2.sql

Import av data

Import til ArcGIS

Vektordata

Last ned ub_pointz.zip og pakk ut zip-filen til ein katalog under ArcGIS Desktop eller QGIS desktop.

Primærnøkkel er hbrgeoid

Tilgjengelige analysedata er: Bosteder i folketellingen 1910 bost1910.csv Bruk autoid for å kople til bostedsfila for folketellinga 1910

Døde i Trondheim 1840-1911 [Ta_kontakt_med_hilde trhbegr.csv]

Analysedata inneheld nøkkelen hbrgeoid. Dette er ein framand-nøkkel hbrgeoid som må koplast til primærnøkkelen i geodata ub_point.

Rasterdata

Last ned Trondheim 1868 trh1868r.jpg og trh1868r.jpw

Last ned Trondheim 1883 (Del 1) trh1883r.jpg og trh1883r.jpw Last ned Trondheim 1883 (Del 2) trh1883vr.jpg og trh1883vr.jpw

Andre bakgrunnsdata

Kartverket WMS server: Åpne WMS tjenester

Import til annen Postgres-database

Slett eksisterande tabellar

drop table hallmykl;
drop table nsd1910a;
drop table gtab1910;

\i urbgeo1.sql
\i urbgeo2.sql