Merge pull request #78 from minus34/202305

Updates for 202305 release

README.md

@@ -51,8 +51,8 @@ The behaviour of gnaf-loader can be controlled by specifying various command lin
 
 #### Optional Arguments
 * `--srid` Sets the coordinate system of the input data. Valid values are `4283` (the default: GDA94 lat/long) and `7844` (GDA2020 lat/long).
-* `--geoscape-version` Geoscape version number in YYYYMM format. Defaults to current year and last release month. e.g. `202302`.
-* `--previous-geoscape-version` Previous Geoscape release version number as YYYYMM; used for QA comparison. e.g. `202211`.
+* `--geoscape-version` Geoscape version number in YYYYMM format. Defaults to current year and last release month. e.g. `202305`.
+* `--previous-geoscape-version` Previous Geoscape release version number as YYYYMM; used for QA comparison. e.g. `202302`.
 * `--raw-gnaf-schema` schema name to store raw GNAF tables in. Defaults to `raw_gnaf_<geoscape_version>`.
 * `--raw-admin-schema` schema name to store raw admin boundary tables in. Defaults to `raw_admin_bdys_<geoscape_version>`.
 * `--gnaf-schema` destination schema name to store final GNAF tables in. Defaults to `gnaf_<geoscape_version>`.
@@ -69,7 +69,7 @@ The behaviour of gnaf-loader can be controlled by specifying various command lin
 * `--no-boundary-tag` DO NOT tag all addresses with some of the key admin boundary IDs for creating aggregates and choropleth maps.
 
 ### Example Command Line Arguments
-* Local Postgres server: `python load-gnaf.py --gnaf-tables-path="C:\temp\geoscape_202302\G-NAF" --admin-bdys-path="C:\temp\geoscape_202302\Administrative Boundaries"` Loads the GNAF tables to a Postgres server running locally. GNAF archives have been extracted to the folder `C:\temp\geoscape_202302\G-NAF`, and admin boundaries have been extracted to the `C:\temp\geoscape_202302\Administrative Boundaries` folder.
+* Local Postgres server: `python load-gnaf.py --gnaf-tables-path="C:\temp\geoscape_202305\G-NAF" --admin-bdys-path="C:\temp\geoscape_202305\Administrative Boundaries"` Loads the GNAF tables to a Postgres server running locally. GNAF archives have been extracted to the folder `C:\temp\geoscape_202305\G-NAF`, and admin boundaries have been extracted to the `C:\temp\geoscape_202305\Administrative Boundaries` folder.
 * Remote Postgres server: `python load-gnaf.py --gnaf-tables-path="\\svr\shared\gnaf" --local-server-dir="f:\shared\gnaf" --admin-bdys-path="c:\temp\unzipped\AdminBounds_ESRI"` Loads the GNAF tables which have been extracted to the shared folder `\\svr\shared\gnaf`. This shared folder corresponds to the local `f:\shared\gnaf` folder on the Postgres server. Admin boundaries have been extracted to the `c:\temp\unzipped\AdminBounds_ESRI` folder.
 * Loading only selected states: `python load-gnaf.py --states VIC TAS NT ...` Loads only the data for Victoria, Tasmania and Northern Territory
 
@@ -117,8 +117,8 @@ Should take 15-60 minutes.
 - A knowledge of [Postgres pg_restore parameters](https://www.postgresql.org/docs/14/app-pgrestore.html)
 
 ### Process
-1. Download the [GNAF dump file](https://minus34.com/opendata/geoscape-202302/gnaf-202302.dmp) or [GNAF GDA2020 dump file](https://minus34.com/opendata/geoscape-202302-gda2020/gnaf-202302.dmp) (~2.0Gb)
-2. Download the [Admin Bdys dump file](https://minus34.com/opendata/geoscape-202302/admin-bdys-202302.dmp) or [Admin Bdys GDA2020 dump file](https://minus34.com/opendata/geoscape-202302-gda2020/admin-bdys-202302.dmp) (~2.8Gb)
+1. Download the [GNAF dump file](https://minus34.com/opendata/geoscape-202305/gnaf-202305.dmp) or [GNAF GDA2020 dump file](https://minus34.com/opendata/geoscape-202305-gda2020/gnaf-202305.dmp) (~2.0Gb)
+2. Download the [Admin Bdys dump file](https://minus34.com/opendata/geoscape-202305/admin-bdys-202305.dmp) or [Admin Bdys GDA2020 dump file](https://minus34.com/opendata/geoscape-202305-gda2020/admin-bdys-202305.dmp) (~2.8Gb)
 3. Edit the _restore-gnaf-admin-bdys.bat_ or _.sh_ script in the supporting-files folder for your dump file names, database parameters and for the location of pg_restore
 5. Run the script, come back in 15-60 minutes and enjoy!
 
@@ -127,11 +127,11 @@ Parquet versions of all the tables are in a public S3 bucket for use directly in
 
 Geometries are stored as Well Known Text (WKT) strings with WGS84 lat/long coordinates (SRID/EPSG:4326). They can be queried using spatial extensions to analytical platforms, such as [Apache Sedona](https://sedona.apache.org/) running on [Apache Spark](https://spark.apache.org/).
 
-The files are here: `s3://minus34.com/opendata/geoscape-202302/parquet/` or `s3://minus34.com/opendata/geoscape-202302-gda2020/parquet/`
+The files are here: `s3://minus34.com/opendata/geoscape-202305/parquet/` or `s3://minus34.com/opendata/geoscape-202305-gda2020/parquet/`
 
 ### AWS CLI Examples:
-- List all datasets: `aws s3 ls s3://minus34.com/opendata/geoscape-202302/parquet/`
-- Copy all datasets: `aws s3 sync s3://minus34.com/opendata/geoscape-202302/parquet/ <my-local-folder>`
+- List all datasets: `aws s3 ls s3://minus34.com/opendata/geoscape-202305/parquet/`
+- Copy all datasets: `aws s3 sync s3://minus34.com/opendata/geoscape-202305/parquet/ <my-local-folder>`
 
 ## DATA LICENSES
 

docker/Dockerfile

@@ -1,7 +1,7 @@
 FROM debian:buster-slim
 
 # replaced the downloading of the Potgres dump files to use local files instead (for performance)
-# ARG BASE_URL="https://minus34.com/opendata/geoscape-202302"
+# ARG BASE_URL="https://minus34.com/opendata/geoscape-202305"
 # ENV BASE_URL ${BASE_URL}
 
 # Postgres user password - WARNING: change this to something a lot more secure
@@ -33,23 +33,23 @@ RUN echo "listen_addresses='*'" >> /etc/postgresql/15/main/postgresql.conf
 RUN mkdir -p /data
 WORKDIR /data
 
-ADD gnaf-202302.dmp .
-ADD admin-bdys-202302.dmp .
+ADD gnaf-202305.dmp .
+ADD admin-bdys-202305.dmp .
 
 # replace the add statements above if wanting to download Postgres dump files
 # RUN /data \
-#     && wget --quiet ${BASE_URL}/gnaf-202302.dmp \
-#     && wget --quiet ${BASE_URL}/admin-bdys-202302.dmp
+#     && wget --quiet ${BASE_URL}/gnaf-202305.dmp \
+#     && wget --quiet ${BASE_URL}/admin-bdys-202305.dmp
 
 RUN /etc/init.d/postgresql start \
-    && pg_restore -Fc -d postgres -h localhost -p 5432 -U postgres /data/gnaf-202302.dmp \
+    && pg_restore -Fc -d postgres -h localhost -p 5432 -U postgres /data/gnaf-202305.dmp \
     && /etc/init.d/postgresql stop \
-    && rm /data/gnaf-202302.dmp
+    && rm /data/gnaf-202305.dmp
 
 RUN /etc/init.d/postgresql start \
-    && pg_restore -Fc -d postgres -h localhost -p 5432 -U postgres /data/admin-bdys-202302.dmp \
+    && pg_restore -Fc -d postgres -h localhost -p 5432 -U postgres /data/admin-bdys-202305.dmp \
     && /etc/init.d/postgresql stop \
-    && rm /data/admin-bdys-202302.dmp
+    && rm /data/admin-bdys-202305.dmp
 
 EXPOSE 5432
 

docker/xx_code_snippets.sh

@@ -2,7 +2,7 @@
 cd /Users/$(whoami)/git/minus34/gnaf-loader/docker
 
 # build gnaf loader image
-docker build --squash --tag minus34/gnafloader:latest --tag minus34/gnafloader:202302 .
+docker build --squash --tag minus34/gnafloader:latest --tag minus34/gnafloader:202305 .
 
 # run gnaf loader container
 docker run --name=gnafloader --publish=5433:5432 minus34/gnafloader:latest

load-gnaf.py

@@ -214,7 +214,7 @@ def populate_raw_gnaf(pg_cur):
         # load all PSV files using multiprocessing
         geoscape.multiprocess_list("sql", sql_list, logger)
 
-        # fix missing geocodes (added due to missing data in 202302 release)
+        # fix missing geocodes (added due to missing data in 202305 release)
         sql = geoscape.open_sql_file("01-04-raw-gnaf-fix-missing-geocodes.sql")
         pg_cur.execute(sql)
 

postgres-scripts/01-04-raw-gnaf-fix-missing-geocodes.sql

@@ -1,4 +1,4 @@
--- workaround for missing default coordinates - 202302 release issue
+-- workaround for missing default coordinates - 202305 release issue
 with missing as (
     select address_detail_pid
     from raw_gnaf.address_default_geocode

postgres-scripts/02-02a-prep-admin-bdys-tables.sql

@@ -203,10 +203,10 @@ UPDATE admin_bdys.locality_bdys
 ;
 
 
--- -- add old locality_pids to unedited localities -- need to rollover old locality pids from GNAF 202302 release - not supplied in 202302 release
+-- -- add old locality_pids to unedited localities -- need to rollover old locality pids from GNAF 202305 release - not supplied in 202305 release
 -- UPDATE admin_bdys.locality_bdys as new
 --     SET old_locality_pid = old.old_locality_pid
--- FROM admin_bdys_202302.locality_bdys AS old
+-- FROM admin_bdys_202305.locality_bdys AS old
 -- WHERE new.locality_pid = old.locality_pid;
 
 

postgres-scripts/xx-04-02-manual-bdy-tags.sql

@@ -4,7 +4,7 @@
 
 
 -- fix 35 boatsheds
-update gnaf_202302.address_principal_admin_boundaries
+update gnaf_202305.address_principal_admin_boundaries
     set lga_pid = 'lgacbffb11990f2',
         lga_name = 'Hobart City'
 where locality_pid = 'loc0f7a581b85b7'

postgres-scripts/xx-add-elevation-to-gnaf.sql

@@ -43,7 +43,7 @@ DROP TABLE IF EXISTS temp_gnaf_100m_points;
 --
 -- SELECT ST_Value(dem.rast, gnaf.geom) as elevation,
 --        *
--- FROM gnaf_202302.address_principals as gnaf
--- INNER JOIN gnaf_202302.srtm_3s_dem as dem on ST_Intersects(gnaf.geom, dem.rast) limit 100;
+-- FROM gnaf_202305.address_principals as gnaf
+-- INNER JOIN gnaf_202305.srtm_3s_dem as dem on ST_Intersects(gnaf.geom, dem.rast) limit 100;
 
 

postgres-scripts/xx-alias-principals-with-different-coordinates.sql

@@ -6,9 +6,9 @@ SELECT als.gnaf_pid, als.street_locality_pid, als.locality_pid, als.alias_princi
 					ST_MakePoint(als.longitude, als.latitude)::geography,
 					ST_MakePoint(gnaf.longitude, gnaf.latitude)::geography
 				) as distance
-  FROM gnaf_202302.address_aliases as als
-  INNER JOIN gnaf_202302.address_alias_lookup as lkp on als.gnaf_pid = lkp.alias_pid
-  INNER JOIN gnaf_202302.address_principals as gnaf on lkp.principal_pid = gnaf.gnaf_pid
+  FROM gnaf_202305.address_aliases as als
+  INNER JOIN gnaf_202305.address_alias_lookup as lkp on als.gnaf_pid = lkp.alias_pid
+  INNER JOIN gnaf_202305.address_principals as gnaf on lkp.principal_pid = gnaf.gnaf_pid
   WHERE als.latitude <> gnaf.latitude
   OR als.longitude <> als.longitude
   order by ST_Distance(

postgres-scripts/xx-export-address-principals-to-csv.sql

@@ -6,5 +6,5 @@ COPY (
 				 address, locality_name, postcode, state, locality_postcode, confidence, 
 				 legal_parcel_id, mb_2016_code, mb_2021_code, latitude, longitude,
 				 geocode_type, reliability
-		FROM gnaf_202302.address_principals
+		FROM gnaf_202305.address_principals
 ) TO '/Users/hugh.saalmans/tmp/address_principals.psv' HEADER CSV;

postgres-scripts/xx-get-population-per-gnafpid.sql

@@ -22,7 +22,7 @@
 --WITH counts AS (
 --	SELECT mb_2016_code,
 --		   count(*) AS address_count
---	FROM gnaf_202302.address_principals
+--	FROM gnaf_202305.address_principals
 --	GROUP BY mb_2016_code
 --)
 --UPDATE testing.mb_2016_counts AS mb
@@ -35,7 +35,7 @@
 ---- add geoms
 --UPDATE testing.mb_2016_counts AS mb
 --  SET geom = bdys.geom
---  FROM admin_bdys_202302.abs_2016_mb as bdys
+--  FROM admin_bdys_202305.abs_2016_mb as bdys
 --  WHERE mb.mb_2016_code = bdys.mb_16code::bigint;
 --
 --ANALYSE testing.mb_2016_counts;
@@ -58,7 +58,7 @@ SELECT gnaf.gnaf_pid,
        mb.person,
        mb.address_count,
        gnaf.geom
-FROM gnaf_202302.address_principals as gnaf
+FROM gnaf_202305.address_principals as gnaf
 INNER JOIN testing.mb_2016_counts AS mb on gnaf.mb_2016_code = mb.mb_2016_code
 WHERE mb.address_count >= mb.dwelling
   AND mb.dwelling > 0
@@ -92,7 +92,7 @@ SELECT gnaf.gnaf_pid,
        mb.address_count,
        gnaf.geom,
        generate_series(1, ceiling(mb.dwelling::float / mb.address_count::float)::integer) as duplicate_number
-FROM gnaf_202302.address_principals as gnaf
+FROM gnaf_202305.address_principals as gnaf
 INNER JOIN testing.mb_2016_counts AS mb on gnaf.mb_2016_code = mb.mb_2016_code
 WHERE mb.address_count < mb.dwelling
     AND address_count > 0
@@ -219,7 +219,7 @@ WITH adr AS (
 	       mb.person,
            mb.address_count,
            gnaf.geom
-	FROM gnaf_202302.address_principals as gnaf
+	FROM gnaf_202305.address_principals as gnaf
 	INNER JOIN testing.mb_2016_counts AS mb on gnaf.mb_2016_code = mb.mb_2016_code
 	WHERE mb.address_count >= mb.person
 	  AND mb.dwelling = 0
@@ -253,7 +253,7 @@ WITH adr AS (
            mb.address_count,
            gnaf.geom,
 		   generate_series(1, ceiling(mb.person::float / mb.address_count::float)::integer) as duplicate_number
-	FROM gnaf_202302.address_principals as gnaf
+	FROM gnaf_202305.address_principals as gnaf
 	INNER JOIN testing.mb_2016_counts AS mb on gnaf.mb_2016_code = mb.mb_2016_code
 	WHERE mb.address_count < mb.person
 	  AND mb.address_count > 0

postgres-scripts/xx_calculate_partitions.sql

@@ -5,13 +5,13 @@ CREATE TABLE testing2.gnaf_partitions AS
 WITH parts AS(
     SELECT unnest((select array_agg(counter) from generate_series(1, 99, 1) AS counter)) as partition_id,
            unnest((select array_agg(fraction) from generate_series(0.01, 0.99, 0.01) AS fraction)) as percentile,
-           unnest((select percentile_cont((select array_agg(s) from generate_series(0.01, 0.99, 0.01) as s)) WITHIN GROUP (ORDER BY longitude) FROM gnaf_202302.address_principals)) as longitude
+           unnest((select percentile_cont((select array_agg(s) from generate_series(0.01, 0.99, 0.01) as s)) WITHIN GROUP (ORDER BY longitude) FROM gnaf_202305.address_principals)) as longitude
 ), parts2 AS (
-SELECT 0 AS partition_id, 0.0 AS percentile, min(longitude) - 0.0001 AS longitude FROM gnaf_202302.address_principals
+SELECT 0 AS partition_id, 0.0 AS percentile, min(longitude) - 0.0001 AS longitude FROM gnaf_202305.address_principals
 UNION ALL
 SELECT * FROM parts
 UNION ALL
-SELECT 100 AS partition_id, 1.0 AS percentile, max(longitude) - 0.0001 AS longitude FROM gnaf_202302.address_principals
+SELECT 100 AS partition_id, 1.0 AS percentile, max(longitude) - 0.0001 AS longitude FROM gnaf_202305.address_principals
 )
 SELECT partition_id,
        percentile,
@@ -43,7 +43,7 @@ WITH merge AS (
            name,
            state,
            st_intersection(bdy.geom, part.geom) AS geom
-    FROM admin_bdys_202302.commonwealth_electorates as bdy
+    FROM admin_bdys_202305.commonwealth_electorates as bdy
     INNER JOIN testing2.gnaf_partitions as part ON st_intersects(bdy.geom, part.geom)
 )
 INSERT INTO testing2.commonwealth_electorates_partitioned (partition_id, ce_pid, name, state, geom)
@@ -65,4 +65,4 @@ commit;
 
 select count(*) from testing2.commonwealth_electorates_partitioned;
 
-select count(*) from admin_bdys_202302.commonwealth_electorates_analysis;
+select count(*) from admin_bdys_202305.commonwealth_electorates_analysis;

postgres-scripts/xx_qa_table_counts.sql

@@ -3,14 +3,14 @@ SELECT new.table_name,
        new.aus - old.aus as difference,
 	   new.aus as new_aus,
        old.aus as old_aus
-	FROM gnaf_202302.qa as new
+	FROM gnaf_202305.qa as new
 	INNER JOIN gnaf_202102.qa as old ON new.table_name = old.table_name
 ;
 
 SELECT new.table_name,
        new.aus - old.aus as difference,
 	   new.aus as new_aus,
        old.aus as old_aus
-	FROM admin_bdys_202302.qa as new
+	FROM admin_bdys_202305.qa as new
 	INNER JOIN admin_bdys_202102.qa as old ON new.table_name = old.table_name
 ;

postgres-scripts/xx_test_state_electorates.sql

@@ -2,27 +2,27 @@
 
 
 
-DROP VIEW IF EXISTS raw_admin_bdys_202302.vw_tenp_state_electorates;
-CREATE VIEW raw_admin_bdys_202302.vw_tenp_state_electorates AS
+DROP VIEW IF EXISTS raw_admin_bdys_202305.vw_tenp_state_electorates;
+CREATE VIEW raw_admin_bdys_202305.vw_tenp_state_electorates AS
 SELECT dat.*,
 	   aut.name,
 	   bdy.se_ply_pid,
 	   bdy.geom
-	FROM raw_admin_bdys_202302.aus_state_electoral as dat
-	INNER JOIN raw_admin_bdys_202302.aus_state_electoral_class_aut as aut on dat.secl_code = aut.code
-	INNER JOIN raw_admin_bdys_202302.aus_state_electoral_polygon as bdy on dat.se_pid = bdy.se_pid
+	FROM raw_admin_bdys_202305.aus_state_electoral as dat
+	INNER JOIN raw_admin_bdys_202305.aus_state_electoral_class_aut as aut on dat.secl_code = aut.code
+	INNER JOIN raw_admin_bdys_202305.aus_state_electoral_polygon as bdy on dat.se_pid = bdy.se_pid
 -- 	where name = 'KEW'
 ;
 
-select * from raw_admin_bdys_202302.vw_tenp_state_electorates
+select * from raw_admin_bdys_202305.vw_tenp_state_electorates
 	where name = 'KEW'
 	order by se_pid,
 	         dt_create
 	;
 
 
 
-select * from raw_admin_bdys_202302.aus_state_electoral_polygon
+select * from raw_admin_bdys_202305.aus_state_electoral_polygon
 	where se_pid = 'VIC292'
 	order by se_pid,
 	         dt_create