genericmappingtools / remotes.jl Goto Github PK

Remote sensing data processing

License: Other

Julia 100.00%

landsat-8 sentinel-2 modis gmt sst geospatial gdal julia

remotes.jl's Introduction

RemoteS.jl

Package to perform operations with satellite data. Easy to use in computing true color images with automatic contrast stretch, many spectral indices and processing of MODIS L2 files. It depends heavily, and so far only, on GMT.jl. Since GMT relies on GDAL to read many raster formats most of satellite and other remote sensing data can be read and processed by this package. Plotting satellite orbits and finding AQUA and TERRA scene names is also possible via the SatelliteToolbox package(s).

Example, a true color image from Landsat8

Load the Red, Green & Blue bands from a Landsat8 scene from AWS (it takes a little time)

using RemoteS, GMT

Ir = gmtread("/vsicurl/http://landsat-pds.s3.amazonaws.com/c1/L8/037/034/LC08_L1TP_037034_20160712_20170221_01_T1/LC08_L1TP_037034_20160712_20170221_01_T1_B4.TIF");
Ig = gmtread("/vsicurl/http://landsat-pds.s3.amazonaws.com/c1/L8/037/034/LC08_L1TP_037034_20160712_20170221_01_T1/LC08_L1TP_037034_20160712_20170221_01_T1_B3.TIF");
Ib = gmtread("/vsicurl/http://landsat-pds.s3.amazonaws.com/c1/L8/037/034/LC08_L1TP_037034_20160712_20170221_01_T1/LC08_L1TP_037034_20160712_20170221_01_T1_B2.TIF");
Irgb = truecolor(Ir, Ig, Ib);
imshow(Irgb)

See more examples ath the Gallery ection of the documentation

Install

] add RemoteS

remotes.jl's People

Contributors

Stargazers

Watchers

Forkers

hdfeos yangleir genericaltimetrytools whigg

remotes.jl's Issues

Histogram stretch and overlay of other rasters

Continuing the discussion from this post.

The bright dots in the following image are targets, the low Rayleigh corrected reflectances values are expected since land was masked out and water absorbs most of the radiation.

According to this tutorial, I used histogram(Ir_img, auto=true, bin=10, show=true) to find good stretch points

I used imshow(Ir_img, stretch=[0 100]), didn't notice any difference on the Ir_img plot and receive this warning the following options were not consumed in grdimage => [:stretch].

I am trying to reduce the darkness of the water in the RGB image.

Also, is it possible to overlay a fourth raster on the RGB plot? If not, is this a future feature?

If you'd like for me to do this for you, comment TagBot fix on this issue.
I'll open a PR within a few hours, please be patient!