odc docs pages add
parent
341b2d3c48
commit
010259e8e4
|
@ -0,0 +1,71 @@
|
|||
:Author: Brian Hamlin
|
||||
:Reviewer: Alex Leith
|
||||
:Version: osgeolive 16
|
||||
:License: Creative Commons Attribution 3.0 Unported (CC BY 3.0)
|
||||
|
||||
@LOGO_odc@
|
||||
@OSGEO_KIND_odc@
|
||||
|
||||
@NAME_odc@
|
||||
================================================================================
|
||||
|
||||
Open Data Cube
|
||||
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
||||
|
||||
The Open Data Cube (ODC) is a community of people and organisations building capability
|
||||
for working with earth observation data. At the core, the ODC is a software library and
|
||||
set of command line tools. And keeping the core together we have a group of organisations
|
||||
collaborating on strategy and who may maintain specific implementations of the ODC.
|
||||
|
||||
The objective of the ODC is to increase the impact of satellite data by providing an open
|
||||
and freely accessible exploitation tool, and to foster a community to develop, sustain,
|
||||
and grow the breadth and depth of applications.
|
||||
|
||||
|
||||
At a technical level, an implementation of the ODC is made up of three things:
|
||||
data, an index and software:
|
||||
|
||||
* Data is usually file based, either in local directories of GeoTIFFs or NetCDF files,
|
||||
but data can be anything that GDAL can read, including Cloud Optimised GeoTIFFs stored on AWS’ S3.
|
||||
|
||||
* For the Index, the ODC uses PostgreSQL as a database to store a list of Products
|
||||
(a specific data type, like Landsat 8 Analysis Ready Data) and Datasets (a single instance
|
||||
of a Product, for example, a single Landsat 8 scene). The index enables a user to ask for
|
||||
data at a time and location, without needing to know specifically where the required files
|
||||
are stored and how to access them.
|
||||
|
||||
* The Software at the core of the ODC is a Python library that enables a user to index data
|
||||
(add records to the Index), to ingest data (optimise indexed data for performance), to query data
|
||||
(returning data in a standard data format) and a wide range of other functions related
|
||||
to managing data.
|
||||
|
||||
|
||||
@SCREENSHOT_odc@
|
||||
|
||||
Core Features
|
||||
--------------------------------------------------------------------------------
|
||||
|
||||
* Journal quality cartography.
|
||||
* Gridding and contouring functions.
|
||||
|
||||
Details
|
||||
--------------------------------------------------------------------------------
|
||||
|
||||
**Website:** @WEB_odc@
|
||||
|
||||
**Licence:** Apache License Version 2
|
||||
|
||||
**Software Version:** |version-odc|
|
||||
|
||||
**Supported Platforms:** GNU/Linux, Mac OSX, MS Windows
|
||||
|
||||
**API Interfaces:** Command line
|
||||
|
||||
**Support:** https:///
|
||||
|
||||
|
||||
@VMDK_odc@
|
||||
@QUICKSTART_odc@
|
||||
|
||||
.. presentation-note
|
||||
Generic Mapping Tools, or GMT is a collection of tools that allow users to manipulate (x,y) and (x,y,z) data sets
|
|
@ -0,0 +1,127 @@
|
|||
:Author: Brian M Hamlin
|
||||
:Reviewer: Alex Leith
|
||||
:Version: osgeolive16.0
|
||||
:License: Creative Commons Attribution-ShareAlike 3.0 Unported (CC BY-SA 3.0)
|
||||
:Copyright: 2023 by The OSGeo Foundation
|
||||
|
||||
@LOGO_odc@
|
||||
@OSGEO_KIND_odc@
|
||||
@VMDK_odc@
|
||||
|
||||
|
||||
********************************************************************************
|
||||
@NAME_odc@ Quickstart
|
||||
********************************************************************************
|
||||
|
||||
The Open Data Cube Core provides an integrated gridded data
|
||||
analysis environment for decades of analysis ready earth observation
|
||||
satellite and related data from multiple satellite and other acquisition
|
||||
systems.
|
||||
|
||||
In this quick Open Data Cube tutorial we will use tools.
|
||||
|
||||
.. contents:: Contents
|
||||
:local:
|
||||
|
||||
Create a map
|
||||
============
|
||||
|
||||
As our input dataset we will use TODO from TODO dataset, which is accessible by GDAL using the GDAL-GRASS driver.
|
||||
|
||||
From a quick look at the DTM metadata, open a terminal and use the GDAL command `gdalinfo`:
|
||||
|
||||
::
|
||||
|
||||
gdalinfo /home/user/grassdata/nc_basic_spm_grass7/PERMANENT/cellhd/elevation
|
||||
|
||||
.. Cameron Review Comment:
|
||||
For each command described, we should show the output of the command, either as text
|
||||
or screenshot.
|
||||
A user should be able to look at this quickstart, and work out what to expect, without
|
||||
running the quickstart.
|
||||
|
||||
We can see the DTM is in a projected coordinate system and the data range is between 50 and 160 meters.
|
||||
|
||||
We can convert the DTM in a GMT Compatible netCDF file format with `gdal_translate` as follow:
|
||||
|
||||
::
|
||||
|
||||
gdal_translate -of GMT /home/user/grassdata/nc_basic_spm_grass7/PERMANENT/cellhd/elevation elevation.gmt
|
||||
|
||||
For our GMT example we will first transform the dataset into geographic coordinates Longitude-Latitude (EPSG:4326) using the `gdalwarp` command:
|
||||
|
||||
::
|
||||
|
||||
gdalwarp elevation.gmt geo_elevation.gmt -t_srs EPSG:4326
|
||||
|
||||
Make a color palette
|
||||
====================
|
||||
|
||||
To assign a color map to our grid we can make a GMT color palette tables
|
||||
with the `makecpt` command. We will set the range and number of colors with the flag `-T`
|
||||
and assign a colormap with the flag `-C`.
|
||||
|
||||
::
|
||||
|
||||
gmt makecpt -Chaxby -T50/160/10 -Z > elevation.cpt
|
||||
|
||||
You can experiment with other colortables, see `man gmt makecpt` for a complete list.
|
||||
|
||||
Define the gradient
|
||||
===================
|
||||
|
||||
For our shaded relief example we need to specify show how the color intensity is distributed along the grid.
|
||||
To do so, we will compute the directional derivative, or gradient, from our grid using the `grdgradient` command.
|
||||
In `grdgradient` we will use the `-Ne` flag to normalize the output using a cumulative Laplace distribution and the `-A` flag for Azimuthal direction for a directional derivative.
|
||||
|
||||
::
|
||||
|
||||
gmt grdgradient geo_elevation.gmt -Ne0.8 -A100 -fg -Ggradient.nc
|
||||
|
||||
The output is a netCDF file which can then be used as color intensity by the `grdimage` command as follow:
|
||||
|
||||
::
|
||||
|
||||
gmt grdimage geo_elevation.gmt -Igradient.nc -JM6i -P -Ba -Celevation.cpt -K > GMT_tut_OSGeoLive.ps
|
||||
|
||||
Display the color bar
|
||||
=====================
|
||||
|
||||
Finally we can add a colorbar with `psscale`:
|
||||
|
||||
::
|
||||
|
||||
gmt psscale -DjTC+w5i/0.25i+h+o0/-1i -Rgeo_elevation.gmt -J -Celevation.cpt -I0.4 -By+lm -O >> GMT_tut_OSGeoLive.ps
|
||||
|
||||
|
||||
Here the flag `-D` is used to select the colorbar position `TC: Top Center`, the size `5x0.25` expressed in inches `i` the orientation by `h: horizontal`
|
||||
and the flag `-By` to set a label `m` on the `y axis`.
|
||||
|
||||
The output ps document should look like:
|
||||
|
||||
.. image:: /images/projects/gmt/GMT_tut_OSGeoLive.png
|
||||
:scale: 70 %
|
||||
:alt: North Caroline - sample elevation data
|
||||
:align: center
|
||||
|
||||
|
||||
Things to try
|
||||
================================================================================
|
||||
|
||||
.. TODO do these needs to be installed to go thru the tutorial?
|
||||
https://readthedocs.org
|
||||
|
||||
* Install the following packages:
|
||||
|
||||
* datacube-core
|
||||
* python3-pystac
|
||||
|
||||
* `TODO tutorial <https://docs.generic-mapping-tools.org/latest/tutorial>`
|
||||
|
||||
What next?
|
||||
==========
|
||||
|
||||
* Documentation: http://datacube-core.readthedocs.io/en/latest/
|
||||
* Website: @WEB_odc@
|
||||
* Support: Slack <http://slack.opendatacube.org>
|
||||
|
Loading…
Reference in New Issue