python处理smap level2c 数据---根据 iqc_flag进行数据质量控制（二进制&十进制）

最近在处理SMAP_level2c_sss数据时，再进行偏差估计时发现数据存在问题，根据卫星图像对比，发现在数据筛选之前未进行质量控制，因此有较大的问题。再次，进行一定的总结：

一、二进制与十进制

在此之前，先让我们简单了解一下基础知识吧~

什么是二进制?

二进制是计算机中广泛采用的一种算法。二进制数据是用0和1两个数码来表示的数，所以显示出来的数只有0、1。进位规则是“逢二进一”。当前的计算机系统使用的基本上是二进制系统，数据在计算机中主要是以补码的形式存储的。计算机中的二进制则是一个非常微小的开关，用1来表示“开”，0来表示“关”。
计算机内部表示数的字节单位是定长的，如8位，16位，或32位。
八位就是：0000 0000
其他类推
假如按照4位二进制的数进行转换，先看一张表格对比一下：
不难看出其规律：按照从右向左排序的方式，十进制转为二进制时依次为

举个栗子，

所以转换后就应该是：0101

这样说是不是一下子就明白了呐~~

还是不理解的自己上网看看视频就行啦~~

二、位运算符

在python中，位运算按照数据在内存中的二进制位（Bit）进行操作，同时， 位运算符只能用来操作整数类型，按照整数在内存中的二进制形式进行计算。具体的位运算符如下表所示：
简单来说，将两个十进制数进行比较，假如有4位二进制数
0000

a&b的意思是：a、b中对应的二进制位数进行比较，两个位都为 1 时，结果才为 1，否则为 0。

举个例子：

a=4 #二进制为：0100
b=5 #二进制为：0101
c=a&b      #既二进制对应比较 0100 
                       #  0101 
#同时为1,则对应位数显示为1，否则为0
#对应结果应该是：0100

最终显示的是十进位的结果！
同理，其他的就好明白啦~
a|b：
两个二进制位有一个为 1 时，结果就为 1，两个都为 0 时结果才为 0

a^b:
参与运算的两个二进制位不同时，结果为 1，相同时结果为 0。

“~” :
对参与运算的二进制位取反 .如~ 1为0，~0为1

<<:
将各个二进制位全部左移若干位，高位丢弃，低位补 0。由 << 右边的数字指定了移动的位数。

“>>”:
将左边的运算数的各二进位全部右移若干位，>> 右边的数字指定了移动的位数

三、如何进行位运算

下面针对我的smap盐度数据进行处理，通过查阅手册得知我的qc_flag为32 bit（位数）

每一位都表示一定的误差，当位数为0时，表示无表中体现的误差，位数为1时，有表中体现的误差。
简单来说，qc_flag中有1就说明，32 bit中有一位对应着一定的问题，可能是：陆地污染、风俗太大、银河辐射等等、、、

知道每一位对应的问题，我们再进行质量控制时，即可按需对比，比如我们想要风速小的数据，只有筛选（表中第12位）的数据是0的数据。

ok，下面打开nc文件中的iqc_flag,它告诉了我们如何比较你想要的位数：
打开发现，这些数字都是2的几次方（十进制数），这就告诉我们：
如果你想比较数据的第一位（时，将数据与1进行 & 比较
如果想比较第 8 bit ，将数据与256 进行比较 & 即可

import numpy as np
import netCDF4 as nc
filename ='G:\\SMAP\\01\\RSS_SMAP_SSS_L2C_r26542_20200120T090520_2020020_FNL_V04.0.nc'
f1 = nc.Dataset(filename)#
sss_smap_40km = f1.variables['sss_smap_40km'][:].data
sss_flg=f1.variables['iqc_flag'][:].data
#==========================================================
sss=np.array(sss_smap_40km[:,:,0])
qc=np.array(sss_flg[:,:,0])
# 比较第一位 bit
result=(qc&1)
# 取出质量控制之后的sss
sss_result=sss[result]
# 比较第八位 bit
result1=(qc&256)

当然，一般来说，我们肯定只想要最好的数据（没有上述表中所说明的问题），既 32 bit 中每一位都是0的数据，这样子就很简单啦

直接令所有的qc数据值等于0，再取出需要的数据即可：

import numpy as np
import netCDF4 as nc
filename ='G:\\SMAP\\01\\RSS_SMAP_SSS_L2C_r26542_20200120T090520_2020020_FNL_V04.0.nc'
f1 = nc.Dataset(filename)#
sss_smap_40km = f1.variables['sss_smap_40km'][:].data
sss_flg=f1.variables['iqc_flag'][:].data
#==========================================================
sss=np.array(sss_smap_40km[:,:,0])
qc=np.array(sss_flg[:,:,0])
# 比较第一位 bit
result=(qc==0)
# 取出质量控制之后的sss
sss_result=sss[result]

下面再进行你需要的处理（取出无效值）就可以啦~~

		                            
                           一个努力学习python的海洋小白
                            水平有限，欢迎指正！！！
                            欢迎评论、收藏。