[ncl-talk] CONUS-SOIL

[Rick Brownrigg]

Hi Michael,

It appears that that file is in big-endian order. Try adding:

setfileoption("bin","ReadByteOrder","BigEndian")

before the cbinread()


Hope that helps...
Rick


On Mon, Jul 15, 2019 at 10:46 AM Michael Notaro via ncl-talk <
ncl-talk at ucar.edu> wrote:

> I am now trying to read the CONUS-Soil dataset for bulk density.  I tried
> the following
>
> but ended up with values of around 320 rather then the expect values
> around 1.5.
>
> What might I be missing here?
>
>
> dims = (/ 11, 2984, 6936 /)
> codes = cbinread ("bd.bsq", dims, "ushort")
> soil=int2flt(ushorttoint(codes))
> soil=soil(:,::-1,:)/100.
>
> Data description
>
>
> http://www.soilinfo.psu.edu/index.cgi?soil_data&conus&data_cov&bd&datasets&geo
> The 30-arcsec gridded data are also available as a three-dimensional array
> of 16-bit unsigned binary integers, where each integer value represents
> the bulk density multiplied by 100. The data are organized as eleven
> two-dimensional layers, giving, respectively, mean BD for each of the 11
> standard soil layers. Within each layer, the values are arranged by
> horizontal rows. All values for the first layer are given first, starting
> with the value for the left (western) end of the top (northern-most) row,
> and continuing left to right along the row. The data for the last row of
> the first layer are then followed by the top row of the second layer, and
> so on.  The array contains 11 layers each containing 2984 rows of 6936
> values.
>
>
> Map of bulk density
>
> http://www.soilinfo.psu.edu/index.cgi?soil_data&conus&data_cov&bd&image
>
>
> Thanks, Michael
>
>
>
> Michael Notaro
> Associate Director
> Nelson Institute Center for Climatic Research
> University of Wisconsin-Madison
> Phone: (608) 261-1503
> Email: mnotaro at wisc.edu
>
>
> ------------------------------
> *From:* Dave Allured - NOAA Affiliate <dave.allured at noaa.gov>
> *Sent:* Tuesday, July 9, 2019 12:05 PM
> *To:* Michael Notaro
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] CONUS-SOIL
>
> The BSQ binary file is a flat binary 3-D array of unsigned bytes.  This is
> the simplest way to read in NCL.  I think I got the correct dimension
> order, but please verify it.
>
>     dims = (/ 11, 2984, 6936 /)
>     codes = cbinread ("layertext.bsq", dims, "ubyte")
>
>
> On Tue, Jul 9, 2019 at 10:15 AM Michael Notaro via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
> Does anyone have a recommendation of the best way in NCL to read in the
> CONUS-SOIL dataset
>
> for soil texture types at 11 soil levels?  It's available in
> band-sequential binary format or
>
> Arc/Info grid format or Arc/Info polygon format.
>
>
>
> http://www.soilinfo.psu.edu/index.cgi?soil_data&conus&data_cov&texture&datasets&geo
>
> Thanks, Michael
>
>
> Michael Notaro
> Associate Director
> Nelson Institute Center for Climatic Research
> University of Wisconsin-Madison
> Phone: (608) 261-1503
> Email: mnotaro at wisc.edu
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190715/7740c57c/attachment.html>