porder: Simple CLI for Planet ordersV2 API  

Ordersv2 is the next iteration of Planet's API in getting Analysis Ready Data (ARD) delivered to you. Orders v2 allows you to improved functionality in this domain, including capability to submit an number of images in a batch order, and perform operations such as top of atmospheric reflectance, compression, coregistration and also enhanced notifications such as email and webhooks. Based on your access you can use this tool to chain together a sequence of operations. This tool is a command line interface that allows you to interact with the ordersv2 API along with place orders and download orders as needed. The tool also allows you to chain multiple processes together and additional functionalities will be added as needed. For exporting to cloud storages release 0.0.8 onwards has a configuration folder with config yml structures to be used with this tool. Simply replaces the fields as needed.

Please note: This tool is in no way an official tool or Planet offering, but is a personal project created and maintained by Samapriya Roy

If you use this tool to download data for your research, and find this tool useful, star and cite it as below

Samapriya Roy. (2019, October 24). samapriya/porder: porder: Simple CLI for Planet ordersV2 API (Version 0.5.2). Zenodo.
http://doi.org/10.5281/zenodo.3518295

Table of contents

• Prerequisites
• Installing porder
• Windows Setup
• Getting started
• porder Simple CLI for Planet ordersv2 API
  • porder version
  • porder quota
  • base64
  • shape to geojson
  • simplify
  • idlist
  • difflist
  • idsplit
  • idcheck
  • bundles
  • order
  • ordersize
  • stats
  • download
  • multipart download
  • multiprocessing download

Prerequisites

This assumes that you have native python & pip installed in your system, you can test this by going to the terminal (or windows command prompt) and trying. I recommend installation within virtual environment if you are worries about messing up your current environment.

python and then pip list

If you get no errors and you have python 2.7.14 or higher you should be good to go.

This command line tool is dependent on shapely and fiona and as such uses functionality from GDAL For installing GDAL in Ubuntu

sudo add-apt-repository ppa:ubuntugis/ppa && sudo apt-get update
sudo apt-get install gdal-bin
sudo apt-get install python-gdal

For Windows I found this guide from UCLA

Also for Ubuntu Linux I saw that this is necessary before the install

sudo apt install libcurl4-openssl-dev libssl-dev

Installing porder

Once you have shapely and the other libraries configured, to install porder: Simple CLI for Planet ordersv2 API you can install using two methods

pip install porder

For linux I found it helps to specify the pip type and use --user

pip install porder --user

or

pip3 install porder --user

or you can also try

git clone https://github.com/samapriya/porder.git
cd porder
python setup.py install

Windows Setup

Shapely and a few other libraries are notoriously difficult to install on windows machines so follow the steps mentioned here before installing porder. You can download and install shapely and other libraries from the Unofficial Wheel files from here download depending on the python version you have. Do this only once you have install GDAL. I would recommend the steps mentioned above to get the GDAL properly installed. However I am including instructions to using a precompiled version of GDAL similar to the other libraries on windows. You can test to see if you have gdal by simply running

gdalinfo

in your command prompt. If you get a read out and not an error message you are good to go. If you don't have gdal try Option 1,2 or 3 in that order and that will install gdal along with the other libraries

Option 1:

Starting from porder v0.4.5 onwards:

Simply run porder -h after installation. This should go fetch the extra libraries you need and install them. Once installation is complete, the porder help page will show up. This should save you from the few steps below.

Option 2:

If this does not work or you get an unexpected error try the following commands. You can also use these commands if you simply want to update these libraries.

pipwin refresh
pipwin install gdal
pipwin install pyproj
pipwin install shapely
pipwin install fiona
pipwing install geopandas

Option 3

For windows first thing you need to figure out is your Python version and whether it is 32 bit or 64 bit. You can do this by going to your command prompt and typing python.

For my windows machine, I have both 32-bit python 2.7.16 and 64-bit Python 3.6.6. You can get the python version at the beginning of the highlighted lines and the 32 or 64 bit within the Intel or AMD64 within the square brackets. Your default python is the one you get by just typing python in the command line. Then download the following packages based on the information we collect about our python type in the earlier step. We use unofficial binaries to install these. This step is only needed if you are on a windows machine if you are using a setup manager like anaconda you might be able to avoid this setup completely

At this stage if you were unable to install gdal then download the gdal binaries first, install that before everything else

gdal: https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal

Then follow along the following libraries * pyproj: https://www.lfd.uci.edu/~gohlke/pythonlibs/#pyproj * shapely: https://www.lfd.uci.edu/~gohlke/pythonlibs/#shapely * fiona: https://www.lfd.uci.edu/~gohlke/pythonlibs/#fiona * geopandas: https://www.lfd.uci.edu/~gohlke/pythonlibs/#geopandas

To choose the version that is correct for you use the python information you collected earlier For example for my python 3.6.6 and AMD 64 if I was installing shapely I would choose the following, here 36 means python 3.6 and amd64 refers to the 64bit we were talking about.

Shapely‑1.6.4.post2‑cp36‑cp36m‑win_amd64.whl

You will get a wheel file or a file ending with .whl. You can now simply browse to the folder or migrate to it in your command prompt. Once there if I am installing for my python 3.6 the command was. At this point we will make use of our trusted package installer that comes with python called pip. Note the choice of pip or pip3 depends on your python version usually you can get the pip to use with your python by typing

pip3 -V

you get a readout like this

pip 18.1 from c:\python3\lib\site-packages\pip (python 3.6)

if you have pip just replace that with pip -V

Then simply install the wheel files you downloaded using the following setup

pip3 install full path to Shapely‑1.6.4.post2‑cp36‑cp36m‑win_amd64.whl

in my case that would be

pip3 install "C:\Users\samapriya\Downloads\Shapely‑1.6.4.post2‑cp36‑cp36m‑win_amd64.whl"

Or you can use anaconda to install. Again, both of these options are mentioned on Shapely’s Official PyPI page.

Getting started

Make sure you initialized planet client by typing planet init. As usual, to print help:

To obtain help for a specific functionality, simply call it with help switch, e.g.: porder idlist -h. If you didn't install porder, then you can run it just by going to porder directory and running python porder.py [arguments go here]

porder Simple CLI for Planet ordersv2 API

The tool is designed to simplify using the ordersv2 API and allows the user to chain together tools and operations for multiple item and asset types and perform these operations and download the assets locally.

porder version

This prints the tool version and escapes. Simple use would be

porder version

porder quota

Just a simple tool to print your planet subscription quota quickly.

base64

This does exactly as it sounds, it encodes your credential files to base64 for use with gcs.

shape to geojson

This tool allows you to convert from a folder with multiple shapefiles to a folder with geojson that can then be used with the tool. It makes use of geopandas and reprojects your shapefile to make it compatible while passing onto the API for search and download.

simplify

This reduces the number of vertices for a multi vertex and complex GeoJSON. Extremely high vertex count (over 500) seem to fail and hence this tool allows you to export a new geojson with reduced vertices. It uses an implementation of the Visvalingam-Wyatt line simplification algorithm. This tool does work with and without Fiona, but Fiona installation is recommended.

idlist

Create an idlist for your geometry based on some basic filters,including geometry, start and end date and cloud cover. If no cloud cover is specified everything form 0 to 100% cloud cover is included. For now the tool can handle geojson,json and kml files. The output is a csv file with ids. The tool also allows you to make sure you get percentage overlap, when selecting image, for clip operations adjust it accordingly (usally --ovp 1 for orders not to fail during clip). The tool now also prints estimated area in Square kilometes for the download and estimated area if you clipped your area with the geometry you are searching (just estimates).

I have changed the setup to now do the following two things

• The number option is optional, so it can look for all images in the time range, but be careful if the area is too large, use at own risk. A better option is to supply the number.

• It is possible to often forget about the different asset types, so you can now not pass an item and the script will return every possible type of asset for each item type depending on the bundle.

A simple setup would be

To run an experiment to add additional filter, you can now pass an additional string or range filter or both flag for string and range filters, a setup would be. The additional filters are optional

porder idlist --input "Path to geojson file" --start "YYYY-MM-DD" --end "YYYY-MM-DD" --item "PSScene4Band" --asset "analytic" --outfile "Path to idlist.csv" --filters range:clear_percent:55:100 --number 20

porder idlist --input "Path to geojson file" --start "YYYY-MM-DD" --end "YYYY-MM-DD" --item "PSScene4Band" --asset "analytic" --outfile "Path to idlist.csv" --filters string:satellite_id:"1003,1006,1012,1020,1038" --number 20

porder idlist --input "Path to geojson file" --start "YYYY-MM-DD" --end "YYYY-MM-DD" --item "PSScene4Band" --asset "analytic" --outfile "Path to idlist.csv" --filters string:satellite_id:"1003,1006,1012,1020,1038" range:clear_percent:55:100 --number 20

The idlist tool can now use a multipolygon and iteratively look for scenes.

difflist

It is possible you already downloaded some images or metadata files, and your you want a difference idlist to create orders for only assets and item types you do not have. It takes in your local folder path, type image or metadata and some basic filters,including geometry, start and end date and cloud cover. If no cloud cover is specified everything form 0 to 100% cloud cover is included. For now the tool can handle geojson,json and kml files. The output is a csv file with ids.

A simple setup would be

or without the cloud filter

idsplit

This allows you to split your idlist into small csv files incase you wanted to created batches of orders.

A simple setup would be

idcheck

It is possible for you to modify the idlist, add or remove ids. Once done, this tool allows you to estimate the total area of images and area that intersect with your geometry or area if clipped.

A simple setup would be

bundles

Ordering using ordersv2 uses the concept of bundles. A bundle is a combination of multiple assets for an item that come together and are delivered as part of the overall fulfillment of the order. For example an analytic asset for PSScene4Band is a single tif file, however the analytic bundle for PSScene4Band includes analytic tiff file, the analytic_xml metadata and the udm data mask file as part of the bundle. You can find more information about bundles here. Thus the concept of bundles bring together single function to order and download multiple related assets. Since the list of bundles is long, this tool simply allows you to get every bundle type based on item type. The setup is simple

A simple setup would be

porder bundles --item "PSScene4Band"

order

This tool allows you to actually place the order using the idlist that you created earlier. the --op argument allows you to take operations, delivery and notifications in a sequence for example --op toar clip email performs Top of Atmospheric reflectance, followed by clipping to your geometry and send you an email notification once the order has completed, failed or had any any change of status. An important changes is the concept of passing bundles instead of using assets. The list of operations for the --op are below and ** the order of these operations is important**

clip|toar|comp osite|zip|zipall|compression|projection|kernel|aws|azu re|gcs|email : ndvi|gndvi|bndvi|ndwi|tvi|osavi|evi2|msavi2|sr

op	description
clip	Clip imagery can handle single and multi polygon verify or create geojson.io
toar	Top of Atmosphere Reflectance imagery generated for imagery
harmonize	Harmonize Dove R (instrument type PS2.SD) data to classic dove (instrument type PS)
composite	Composite number of images in a given order
zip	Zip bundles together and creates downloads (each asset has a single bundle so multiple zip files)
zipall	Create a single zip file containing all assets
compression	Use image compression
projection	Reproject before downloaing image
aws	Option called to specify delivery to AWS
azure	Option called to specify delivery to AZURE
gcs	Option called to specify delivery to GCS
email	Email notification to your planet registered email

You can now add some predefined indices for PlanetScope 4 band items with a maximum of 5 indices for a single setup . This is experimental. The list of indices include

Index	Source
Simple ratio (SR)	Jordan 1969
Normalized Difference Vegetation Index (NDVI)	Rouse et al 1973
Green Normalized Difference Index (GNDVI)	Gitelson et al 1996
Blue Normalized Difference Vegetation Index (BNDVI)	Wang et al 2007
Transformed Vegetation Index (TVI)	Broge and Leblanc 2000
Optimized Soil Adjusted Vegetation Index (OSAVI)	Rondeaux et al 1996
Enhanced Vegetation Index (EVI2)	Jian et al 2008
Normalized Difference Water Index (NDWI)	Gao 1996
Modified Soil-adjusted Vegetation Index v2 (MSAVI2)	Qi 1994

A simple setup with image clip with email notification would be

The same setup with delivery of each image, its metadata as a zip file would be. Note how we only added zip to the op list

A simple setup with Top of Atmospher reflectance and a few indices along with email notification would be

ordersize

The tool now allows you to estimate the total download size for a specific order.

An example setup look like the following

stats

The tool allows you to check on number of running and queued orders for both organization and user level. Using this is simple

porder stats

output should look like this:

Checking on all running orders...
Total queued order for organization: 0
Total running orders for organization: 1

Total queued orders for user: 0
Total running orders for user: 0

download

The allows you to download the files in your order, to a local folder. It uses the order url generated using the orders tool to access and download the files.

multipart download

The allows you to multipart download the files in your order, this uses a multiprocessing downloader to quickly download your files to a local folder. It uses the order url generated using the orders tool to access and download the files.

multiprocessing download

The uses the multiprocessing library to quickly download your files to a local folder. It uses the order url generated using the orders tool to access and download the files and includes an expotential rate limiting function to handle too many requests. To save on time it uses an extension filter so for example if you are using the zip operation you can use ".zip" and if you are downloading only images, udm and xml you can use ".tif" or ".xml" accordingly. For python 3.4 or higher, this switches to using an true async downloader instead of using multiprocessing.

A simple setup would be

Changelog

v0.5.2

• Added harmonization tool to harmonize PS2.SD to PS2.
• Improvements and error handling to quota tool
• Merged pull request 35 to keep download progress via enumerate.

v0.5.1

• Added utf-8 encoding for shapefile to geojson conversion
• Merged pull request 34 to refresh url once expired.

v0.4.9

• Fixed issue with gdal import for pipwin windows.
• Fixed import issue with stats endpoint.

v0.4.8

• Replaced concurrency check with stats endpoint to get queued and running orders.
• Change pipwin cache refresh time to two weeks.

v0.4.7

• Fixed issue with queuing state for orders and downloads.

v0.4.6

• Handles refreshing pipwin cache and better error handling
• Fixed issue with downloading unique manifest ID for zip files.
• Updated ReadMe with improved documentation.

v0.4.5

• Handles installation of windows specific libraries using pipwin.
• General improvements

v0.4.4

• Manifest files for each asset is now written in format ItemID_manifest.json to avoid skipping manifest.json common file name.
• Simple and multipart downloader now show number of items remaining during the download.
• General improvements, bundles tool now prints Bundles:Name followed by assets included in the bundle

v0.4.3

• Fixed issues with setup.py and pyproj version.
• Improved ReadMe instructions.

v0.4.2

• Added geometry check functionality to multipolygon with shapely self intersection Issue 30.
• For multipolygons this also performs a vertex count check and simplifies polygon to fit under 500 vertices.
• General improvements

v0.4.1

• Fixed issue with shapely self intersection using buffer(0).
• General improvements

v0.4.0

• Fixed issue with placing reprojection request.
• Downloader can now download partial as well as completely successful orders.
• Added retry method for rate limit during downloading
• General improvements

v0.3.9

• Removed deprecated bundles from bundles list.
• Improved parameter description

v0.3.7

• Added capability to pass subscription id when submitting order.

v0.3.6

• Replaced asset in order tool with bundles.
• Created a new bundles tool to generate bundle list for an item type
• Improvements to the idlist tool now prints output as it makes progress.

v0.3.5

• Better integration for quota tool
• Updates information while waiting for idlist
• Updated requirements

v0.3.4

• Added async downloader for python 3.4
• Checks for existing files before spawning processes
• Better handling of multiprocessing output
• Added a quick version tool

v0.3.3

• Fixed issue with order name when no ops are used.
• Used file basename for splitting the idlist.

v0.3.2

• idlist tool can no use a multipolygon and iteratively look for scenes
• Orders clip tool can now handle multipolygon clip
• Added new tool zipall to handle single archive download in format ordername_date.zip

v0.3.1

• Can now support an additional string and range filter
• Check total area and clip area estimates from any idlist using idcheck.
• General improvements to the tool.

v0.3.0

• Enhances idlist to execute faster search and return using Planet CLI
• Included better error handling while placing order.

v0.2.8

• Added tool to convert folder with shapefiles to GeoJSONs

v0.2.7

• Improved overlap calculations for larger geometries
• Added a geometry simplification tool to reduce number of vertices

v0.2.6

• Skysat area are calculated using EPSG:3857 to resolve metadata EPSG issue
• General improvements

v0.2.5

• Fixed issue with area calculation estimates
• General improvements

v0.2.4

• Now functions without limit on the number of assets in the idlist
• Parses possible asset types if only item type is supplied for idlist

v0.2.3

• Now estimates area before and after clip when you run idlist
• General improvements

v0.2.1

• Now exports only csv idlist
• Fixed count with concurrency check

v0.2.0

• Fixed pysmartdl install issues
• Added concurrent orders check
• version and os resolve for shapely

v0.1.9

• Added msavi
• Fixed issues with GeoJSON read

v0.1.8

• Fixed issues with empty JSON append
• General improvements to the tool

v0.1.7

• Added band math indices for PlanetScope item
• Fixed issues with retry for downloader
• General improvements to the tool

v0.1.6

• Made fixes to have python 3.X compatability

v0.1.5

• General improvements and bug fixes

v0.1.4

• Fixed issue with Python 3 CSV write compatability
• Fixed issues with Shapely instance issue

v0.1.3

• Fixed issue with clipboard access in headless setup

v0.1.2

• Fixed issue and extension for multiprocessing downloader
• Overall general improvements to the tool

v0.1.0

• Fixed issue and improved idlist and sort
• Fixed issue with clip tool
• Overall general improvements to the tool

v0.0.8

• Improvements to operations in order tool
• Now supports export to gcs/azure/aws along with kernel, projection and compression
• base64 encoding tool for encoding gcs credentials
• Overall general improvements to the tool

v0.0.7

• Now allows for all downloads or download using extension
• Polling for order to complete and automatically download
• General improvements

v0.0.6

• Merged contribution by David Shean
• Fixed issues with op equals None
• Fixed issues with relative import
• Improved Py3 compatability
• General improvements

v0.0.5

• Added exponential backoff for pydl
• Fixed issues with dependency
• General overall improvements

v0.0.4

• Created strict geoinstersection to avoid orders to fail
• Improvements to overlap function
• General overall improvements

v0.0.3

• Added overlap function to idlist
• Added multiprocessing downloader with rate limit and extension filter
• General overall improvements

v0.0.2

• Fixed issues with import modules