[发明专利]历史遥感产品数据支持下的作物种植面积测量高效抽样方法

权利要求书

1.一种历史遥感产品数据支持下的作物种植面积测量高效抽样方法，其特征在于，所述方法包括以下几个步骤：

步骤一：对遥感数据标准化预处理；

步骤二：以辅助变量作为分层依据，采用两阶段分层抽样；

首先对历史中分辨率遥感影像作物种植面积进行提取，结合耕地地块数据，建立入样总体；

其次，采用不同的辅助变量作为分层标志，基于历史产品数据的分层标志的计算；

然后，将历史遥感产品用作入样总体和抽样框的合理性分析；

最后，对分层标志进行有效性研究。

2.根据权利要求1所述的高效抽样方法，所述的辅助变量为农作物种植面积多年变化率和面积规模，所述的两阶段分层抽样包括以下两个阶段：

第一阶段采用农作物种植面积多年变化率作为分层标志将入样总体分割成不同的子总体；

第二阶段在第一阶段划分的子总体前提下以面积规模为分层标志进一步分层，在各子总体中利用随机抽样方法抽选样本；

所述入样总体为区域内的耕地范围，所述两阶段分层抽样还包括抽样单元的建立，抽样单元的类型包括由自然边界构成的地块，乡、村行政边界或者规则的格网数据。所述分层标志为多年农作物种植面积变化率CACVM，定义为：

CACVM=SDMean=Σi=nm(S0i-S‾)2/t/S‾=Σi=nm(S0i-S‾)2/tS‾]]>

其中，Mean=Σi=nmS0i/N=Σi=nmS0i/t=S‾]]>

S_oi代表抽样单元内第i年的农作物种植面积，N代表历史年份，t代表历史总年数，Mean代表抽样单元内多年农作物种植面积平均值，SD代表单个抽样单元内多年农作物种植面积的标准差，CACVM代表抽样单元内种植面积年际间的变化率。R₁、R₂、…R_L代表一阶段分层后的各层数据集合，L为一阶段分层的层数，总体R为一阶段分层后各层的并集；

R＝R₁∪R₂∪...∪R_l

在一阶段分层内采用面积规模指标进行二阶段分层，R_i1，R_i2，…，R_im分别代表一阶段分层结果中第i层的二阶段分层结果，M为一阶段分层内二阶段分层的层数，一阶段每个分层集合为内部二阶段分层结果的并集，则：

R_i＝R_i1∪R_i2∪...∪R_ij(j＝1，2…M)

3.根据权利要求2所述的高效抽样方法，其特征在于：利用历史遥感产品数据建立两阶段分层标志(农作物种植面积多年变化率，面积规模)进行两阶段分层抽样，采用累计平方根法确定最优分层的划分点，即根据分层标志y的密度函数f(y)的开方累计计算分层。