https://github.com/amacd31/geofabric_tools

Tools for working with the Australian Hydrological Geospatial Fabric (geofabric)
https://github.com/amacd31/geofabric_tools

australia geospatial geospatial-processing hydrology

Last synced: 10 days ago
JSON representation

Tools for working with the Australian Hydrological Geospatial Fabric (geofabric)

• Host: GitHub
• URL: https://github.com/amacd31/geofabric_tools
• Owner: amacd31
• License: other
• Created: 2017-02-04T12:28:44.000Z (almost 9 years ago)
• Default Branch: master
• Last Pushed: 2022-02-05T11:40:55.000Z (almost 4 years ago)
• Last Synced: 2024-04-24T11:11:12.947Z (almost 2 years ago)
• Topics: australia, geospatial, geospatial-processing, hydrology
• Language: Python
• Size: 33.2 KB
• Stars: 6
• Watchers: 3
• Forks: 1
• Open Issues: 2
• Metadata Files:
  • Readme: README.rst
  • License: LICENSE

Awesome Lists containing this project

README

          
Geofabric Tools

===============

Tools for working with the Australian Hydrological Geospatial Fabric (geofabric).

Data from the Bureau of Meterology: http://www.bom.gov.au/water/geofabric/index.shtml

Download and unzip: ftp://ftp.bom.gov.au/anon/home/geofabric/SH_Network_GDB_V3_2.zip

or the previous release: ftp://ftp.bom.gov.au/anon/home/geofabric/version2/SH_Network_GDB_V2_1_1.zip

Convert SH_Catchments.gdb into a spatialite database::

    $ time ogr2ogr -f SQLite -dsco SPATIALITE=YES geofabric.sqlite SH_Network_GDB/SH_Network.gdb

    

    real    9m51.489s

    user    7m38.091s

    sys     0m7.016s

Add some indexes to the spatialite database::

    $ sqlite3 geofabric.sqlite "CREATE INDEX idx_ahgfnetworkconnectivitydown_to_id ON ahgfnetworkconnectivitydown (to_id);"

    $ sqlite3 geofabric.sqlite "CREATE INDEX idx_ahgfcatchment_hydroid ON ahgfcatchment (hydroid);"

    $ sqlite3 geofabric.sqlite "CREATE INDEX idx_ahgfcatchment_netnodeid ON ahgfcatchment (netnodeid);"

Extract catchment boundary upstream from sub-catchment containing a latitude and longitude::

    $ time python geofabric_tools/catchment_tools.py 148.821,-35.592

    INFO:__main__:Extracting: 148.821,-35.592

    INFO:__main__:Saving as: 1170184.json

    real    0m4.067s

    user    0m2.904s

    sys     0m1.126s

Optionally include an ID to use when saving the output file::

    $ time python geofabric_tools/catchment_tools.py 148.821,-35.592:410730

    INFO:__main__:Extracting: 148.821,-35.592

    INFO:__main__:Saving as: 410730.json

    real    0m5.810s

    user    0m4.645s

    sys     0m0.885s

By defaults all sinks inside the extracted boundary are assumed to be part of the catchment.

A catchment boundary can be extracted with the sinks excluded from the resulting polygon by

passing the --exclude-sinks option. For example::

    $ time python geofabric_tools/catchment_tools.py --exclude-sinks 148.821,-35.592:excluded_sinks

    INFO:__main__:Extracting: 148.821,-35.592:excluded_sinks

    INFO:__main__:Saving as: excluded_sinks.json

    real    0m6.915s

    user    0m5.936s

    sys     0m0.963s