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 10). samapriya/porder: porder: Simple CLI for Planet ordersV2 API (Version 0.4.9). Zenodo.
http://doi.org/10.5281/zenodo.3478850
- Prerequisites
- Installing porder
- Windows Setup
- Getting started
- porder Simple CLI for Planet ordersv2 API
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
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
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
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.
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
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.
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]
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.
This prints the tool version and escapes. Simple use would be
porder version
Just a simple tool to print your planet subscription quota quickly.
This does exactly as it sounds, it encodes your credential files to base64 for use with gcs.
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.
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.
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.
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.
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.
or without the cloud filter
This allows you to split your idlist into small csv files incase you wanted to created batches of orders.
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.
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"
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 |
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 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
The tool now allows you to estimate the total download size for a specific order.
An example setup look like the following
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
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.
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.
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
- Fixed issue with gdal import for pipwin windows.
- Fixed import issue with stats endpoint.
- Replaced concurrency check with stats endpoint to get queued and running orders.
- Change pipwin cache refresh time to two weeks.
- Fixed issue with queuing state for orders and downloads.
- Handles refreshing pipwin cache and better error handling
- Fixed issue with downloading unique manifest ID for zip files.
- Updated ReadMe with improved documentation.
- Handles installation of windows specific libraries using pipwin.
- General improvements
- 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
- Fixed issues with setup.py and pyproj version.
- Improved ReadMe instructions.
- 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
- Fixed issue with shapely self intersection using buffer(0).
- General improvements
- 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
- Removed deprecated bundles from bundles list.
- Improved parameter description
- Added capability to pass subscription id when submitting order.
- 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.
- Better integration for quota tool
- Updates information while waiting for idlist
- Updated requirements
- Added async downloader for python 3.4
- Checks for existing files before spawning processes
- Better handling of multiprocessing output
- Added a quick version tool
- Fixed issue with order name when no ops are used.
- Used file basename for splitting the idlist.
- 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
- 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.
- Enhances idlist to execute faster search and return using Planet CLI
- Included better error handling while placing order.
- Added tool to convert folder with shapefiles to GeoJSONs
- Improved overlap calculations for larger geometries
- Added a geometry simplification tool to reduce number of vertices
- Skysat area are calculated using EPSG:3857 to resolve metadata EPSG issue
- General improvements
- Fixed issue with area calculation estimates
- General improvements
- Now functions without limit on the number of assets in the idlist
- Parses possible asset types if only item type is supplied for idlist
- Now estimates area before and after clip when you run idlist
- General improvements
- Now exports only csv idlist
- Fixed count with concurrency check
- Fixed pysmartdl install issues
- Added concurrent orders check
- version and os resolve for shapely
- Added msavi
- Fixed issues with GeoJSON read
- Fixed issues with empty JSON append
- General improvements to the tool
- Added band math indices for PlanetScope item
- Fixed issues with retry for downloader
- General improvements to the tool
- Made fixes to have python 3.X compatability
- General improvements and bug fixes
- Fixed issue with Python 3 CSV write compatability
- Fixed issues with Shapely instance issue
- Fixed issue with clipboard access in headless setup
- Fixed issue and extension for multiprocessing downloader
- Overall general improvements to the tool
- Fixed issue and improved idlist and sort
- Fixed issue with clip tool
- Overall general improvements to the tool
- 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
- Now allows for all downloads or download using extension
- Polling for order to complete and automatically download
- General improvements
- Merged contribution by David Shean
- Fixed issues with op equals None
- Fixed issues with relative import
- Improved Py3 compatability
- General improvements
- Added exponential backoff for pydl
- Fixed issues with dependency
- General overall improvements
- Created strict geoinstersection to avoid orders to fail
- Improvements to overlap function
- General overall improvements
- Added overlap function to idlist
- Added multiprocessing downloader with rate limit and extension filter
- General overall improvements
- Fixed issues with import modules