Skip to content

Instantly share code, notes, and snippets.

@mhiley
mhiley / cutline.geojson
Created June 13, 2019 16:03
Gdal issue file attachment 2/2
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@mhiley
mhiley / all_aoi_files.vrt
Created June 13, 2019 16:02
Gdal issue file attachment 1/2
<VRTDataset rasterXSize="34116" rasterYSize="33467">
<SRS>PROJCS["NAD_1983_HARN_WISCRS_Langlade_County_Feet",GEOGCS["GCS_North_American_1983_HARN",DATUM["NAD83_High_Accuracy_Reference_Network",SPHEROID["GRS 1980",6378137,298.2572221010042,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6152"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic_1SP"],PARAMETER["latitude_of_origin",45.15423710527778],PARAMETER["central_meridian",-89.03333333333333],PARAMETER["scale_factor",1.0000627024],PARAMETER["false_easting",651000],PARAMETER["false_northing",345405.421],UNIT["US survey foot",0.3048006096012192,AUTHORITY["EPSG","9003"]]]</SRS>
<GeoTransform> 6.0804913399999996e+05, 5.0000000000000000e-01, 0.0000000000000000e+00, 3.5421392200000002e+05, 0.0000000000000000e+00, -5.0000000000000000e-01</GeoTransform>
<VRTRasterBand dataType="Byte" band="1">
<ColorInterp>Red</ColorInterp>
<SimpleSource>
<SourceFilename relativeToVRT="0">/vsicurl/https://s3.us-east-2
@mhiley
mhiley / test.py
Last active February 10, 2016 18:08
access to NetCDF attributes with McIDAS-V loadGrid
boomstick()
path = '/path/to/data/'
f = path+'DR_ABI-L2-CMIPC-M3C01_G16_s20152322000151_e20152322002542_c20152322002596.nc'
cmi = loadGrid(filename=f, field='CMI')
dqf = loadGrid(filename=f, field='DQF')
l = activeDisplay().createLayer('Color-Shaded Plan View', cmi)
""" Example McIDAS-V script to load Lab 5 SEVIRI images
See below for instructions on how to use this script.
"""
# get a reference to the Local ADDE Dataset
dataset = getLocalADDEEntry(
dataset = 'MSG',
imageType = 'Channels 1-11'
)