[ncl-talk] Extracting bit information using 'dim_gbits'..
Dennis Shea
shea at ucar.edu
Wed May 3 10:08:49 MDT 2017
I am not sure I understand but ...
Let us(226,4350) be of type unsigned short
[1] step-by-step
us = f->UnsignedShort
printVarSummary(us)
dim_us = *dimsizes*(us)
nleft = dim_us(0) ; 226
nright = dim_us(1) ; 4350
ibit = 9 ; skip 9 bits
nbit = 3 ; 10-12
nskip = 16 ; all 'short' are 16 bits
N = nleft*nright
us_1d = *ndtooned*(us),
ubit_1d = *dim_gbits*(us_1d, ibit, nbit, nskip, N)
ubit = *reshape*(ubit_1d, (/nleft,nright/) )
*delete*(us_1d, ubit_1d) ; clean up; not really necessary
[2} slicker but less clear
ubit = reshape( *dim_gbits(**ndtooned*(us), ibit, nbit, nskip, N),
*product*(*dimsizes*(us)) )
[3] As always
* printVarSummary*(ubit)
*printMinMax*(ubit, 0)
Good Luck
On Wed, May 3, 2017 at 2:19 AM, qudsia zafar <missquaddus at gmail.com> wrote:
> Hi,
>
> I am working with a 2-D array [226 x 4350] of type 'ushort' consisting of
> 16-bit integers. I want to extract bits 1-3 in correspondence with bits
> 10-12 from the array. i have seen the function 'dim_gbits' but am confused
> in setting the function parameters such as *ibit*, *nbits *and *nskip. *
>
> Would very much appreciate the help in this regard.
>
> Sincerely,
> qudsia
>
> _______________________________________________
> 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/20170503/25fe051f/attachment.html
More information about the ncl-talk
mailing list