Add script for statsgo soil data processing

[jhamman]

I needed this asap so I just did a one off. Here's what I did:

!gdalwarp -t_srs PROJCS["Albers_Conical_Equal_Area",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,-0,-0,-0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["standard_parallel_1",29.5],PARAMETER["standard_parallel_2",45.5],PARAMETER["latitude_of_center",23],PARAMETER["longitude_of_center",-96],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["meters",1]] -te -2493045.0 177285.0 2342655.0 3310005.0 -tr 4000 4000 -srcnodata 255 -r average awc.tif awc_4000m.tif

from adlfs import AzureBlobFileSystem
fs = AzureBlobFileSystem('carbonplan', account_key=os.environ["BLOB_ACCOUNT_KEY"])

with fs.open('carbonplan-data/processed/statsgo/conus/4000m/raster.zarr', mode='wb') as target:
    with open('/home/jovyan/awc_4000m.tif', 'rb') as source:
        target.write(source.read())