A Regional Mapping Method for Oilseed Rape Based on HSV Transformation and Spectral Features

This study proposed a colorimetric transformation and spectral features-based oilseed rape extraction algorithm (CSRA) to map oilseed rape at the provincial scale as a first step towards country-scale coverage. Using a stepwise analysis strategy, our method gradually separates vegetation from non-vegetation, crop from non-crop, and oilseed rape from winter wheat. The wide-field view (WFV) images from Chinese Gaofen satellite no. 1 (GF-1) at six continuous flowering stages in Wuxue City, Hubei Province, China are used to extract the unique characteristics of oilseed rape during the flowering period and predict the parameter of the CSRA method. The oilseed rape maps of Hubei Province from 2014 to 2017 are obtained automatically based on the CSRA method using GF-1 WFV images. As a result, the CSRA-derived provincial oilseed rape maps achieved at least 85% overall accuracy of spatial consistency when comparing with local reference oilseed rape maps and lower than 20% absolute error of provincial planting areas when comparing with agricultural census data. The robustness of the CSRA method is also tested on other satellite images including one panchromatic and multispectral image from GF-2 and two RapidEye images. Moreover, the comparison between the CSRA and other previous methods is discussed using the six GF-1 WFV images of Wuxue City, showing the proposed method has better mapping accuracy than other tested methods. These results highlight the potential of our method for accurate extraction and regional mapping capacity for oilseed rape.


Introduction
Oilseed rape (Brassica napus L.) is a major cash crop, as well as one of the principal sources of edible oil, high energy, and protein meal plant species in countries with moderate climate (e.g., China, Europe, Canada, India, Australia, etc.) [1][2][3].Accurate maps for oilseed rape (OR) distribution are essential to guide macro decision-making which maintains balance in agricultural supply and demand.However, the OR sowing areas in China are mainly generated by the ground samples method at the county scale, which lacks detailed spatial distribution with a time lag [4,5].The lack of an effective and efficient approach to the mapping OR at the large regional scale (e.g., provincial, country, or global scale) impedes its efficient management.In the last decades, the emergence of remote sensing with a large spatial coverage and in near real-time has facilitated the monitoring of crop identification and of reflectance, the visual difference caused by color changes during the flowering period is another important feature for OR, which means that color space technology provides a new opportunity to identify and extract OR in the flowering period.A color space composed of specific bands has been commonly used to distinguish land cover types based on the colors [26,27].The HSV color space is widely used to perform the target identification and extraction.For instance, Pekel et al. [28] proposed an HSV-based algorithm to manage vegetation in arid and semi-arid regions.He also developed a real-time method to detect the water surface based on the HSV transformation of Moderate Resolution Imaging Spectroradiometer time series data [29].Lessel et al. [30] presented the basic framework for a semi-automated crop identification method that identified soil and vegetation patterns for various crop-cycle scenarios using H values.
To explore the phenology feature of OR during the flowering period, imagery with high temporal resolution is needed.High spatial resolution imagery is necessary for the small size of crop fields [31], especially in Southern China [23].The satellite imagery including Landsat, WorldView-2, and Satellite Pour l'Observation de la Terre (SPOT) are the main selections for crop discrimination [32,33].However, these datasets usually have low temporal resolution leading to difficulty in acquiring accurate crop information over large areas during the growing season [34].Currently, GF satellite imagery has become a popular selection for agricultural monitoring in China [35,36].GF is the program name of the Chinese high resolution land-observation system containing a series of civilian satellites.As the first satellite of the Chinese high-resolution land-observation system, GF-1 carries two imaging system: a panchromatic and multispectral (PMS) system which is used for common high resolution observation, and a WFV system which is used for large-scale observation [37].The characteristics of middle spatial resolution (16 m) coupled with high temporal resolution (two days), large swath width (800 km combining four cameras), and free acquisition promotes the application of GF-1 WFV images in crop monitoring [38][39][40].GF-2 is a sun-synchronous satellite with spatial resolution better than 1 m.The PMS imaging system on GF-2 equips with two cameras of 1 m panchromatic and 4 m multispectral resolutions with 23 km swath width [41], which is well satisfied for the demand of detailed scale crop mapping.Thus, the spatial and temporal resolution of the GF series imagery have a great advantage in agricultural monitoring both at regional and local scales and more frequent updates.
To partly solve the deficiency of existing methods, as well as to take advantage of the data from the newly-launched GF satellites, the objectives of this study are to: (1) develop an effective and generic extraction method for OR based on HSV transformation and spectral features; (2) produce annual OR planting maps of Hubei Province from 2014 to 2017 and evaluate the products as a first step towards country scale coverage; and (3) assess the robustness of the developed method.

Study Area
Hubei Province (Figure 1a) is selected as the study area based on its status as a major Chinese agricultural province with the largest OR production in China.Hubei Province is located in Central China and extends from 108 • 21 E to 116 • 07 E in longitude and 29 • 05 N to 33 • 20 N in latitude, with the Yangtze and Hanjiang Rivers running through.The topography is characterized by an inclination overall going down from west to east with various landforms.The climate is humid subtropical and is largely controlled by the East Asian monsoon.The advantageous natural conditions promoted the excellent agricultural environment since ancient China.The crop species growing during winter growth season (from late September every year to mid-May in next year) in Hubei Province are OR and WW, based on the records of agro-meteorological stations.Hubei Province contributes one-sixth of the OR planting areas and products of China according to the statistical data from the China Statistical Yearbook.There are three OR producing areas with 80% planting areas and production in Hubei Province.The corresponding names for the different OR producing areas are: the Jianghan Plain (blue ellipse in Figure 1a, mainly containing Jinzhou, Xiantao, Tianmen, Qianjiang and part of Wuhan, Xiaogan, Yichang), the Eastern E Plain (E is a traditional abbreviation for Hubei Province, purple ellipse in Figure 1a, mainly containing Ezhou, Huanggang, Huangshi), and the Central-North E (red ellipse in Figure 1a, mainly containing Jinmen, Xiangyang and part of Xiaogan, Suizhou, Yichang).In fact, except for the high mountain area of Western Hubei Province, the whole mainland of Hubei Province is suitable for the growth of OR.
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 4 of 24 E (red ellipse in Figure 1a, mainly containing Jinmen, Xiangyang and part of Xiaogan, Suizhou, Yichang).In fact, except for the high mountain area of Western Hubei Province, the whole mainland of Hubei Province is suitable for the growth of OR.
The OR experiment bases in Wuxue City and Shayang County of our research team were built in 2014 and 2016, respectively.Wuxue City (Figure 1b) is located in Eastern Hubei Province and famous for OR with "two low" (low erucic acid and low glucosinolate).Shayang County is located in Central Hubei Province and famous for OR with high production.OR is a traditional dominated crop species during winter growth season both in Wuxue City and Shayang County.During the OR planting and observing experiments, investigation of crop species in Wuxue City was conducted in 2015 and 2016, as well as in Shayang County in 2016.

Remote Sensing Data
The remote sensing dataset used in this study is composed of sixty-six GF-1 WFV images, one GF-2 PMS image, and two RapidEye images.Among them, six GF-1 WFV images covering Wuxue City and representing continuous flowering stages of OR are used to learn and predict the parameter of the CSRA method.One GF-1 WFV, one GF-2 PMS, and two RapidEye images covering part of Shayang County in Jinmen are used to explore the robustness of the CSRA method.A summary of information regarding the application, sensor, acquisition date, number, coverage area, and phenology stage for these images are listed in Table 1.The other GF-1 WFV images are used to produce OR maps of Hubei Province from 2014 to 2017, whose detailed information, like Table 1, are listed in Table A1 The OR experiment bases in Wuxue City and Shayang County of our research team were built in 2014 and 2016, respectively.Wuxue City (Figure 1b) is located in Eastern Hubei Province and famous for OR with "two low" (low erucic acid and low glucosinolate).Shayang County is located in Central Hubei Province and famous for OR with high production.OR is a traditional dominated crop species during winter growth season both in Wuxue City and Shayang County.During the OR planting and observing experiments, investigation of crop species in Wuxue City was conducted in 2015 and 2016, as well as in Shayang County in 2016.

Remote Sensing Data
The remote sensing dataset used in this study is composed of sixty-six GF-1 WFV images, one GF-2 PMS image, and two RapidEye images.Among them, six GF-1 WFV images covering Wuxue City and representing continuous flowering stages of OR are used to learn and predict the parameter of the CSRA method.One GF-1 WFV, one GF-2 PMS, and two RapidEye images covering part of Shayang County in Jinmen are used to explore the robustness of the CSRA method.A summary of information regarding the application, sensor, acquisition date, number, coverage area, and phenology stage for these images are listed in Table 1.The other GF-1 WFV images are used to produce OR maps of Hubei Province from 2014 to 2017, whose detailed information, like Table 1, are listed in Table A1.The GF-1 WFV images are downloaded from the China Center for Resources Satellite Date and Application (CRESDA).The GF-2 PMS image is obtained from the Hubei Data and Application Center of High-Resolution Land-Observation Systems.The RapidEye images were bought by the research team.The series of GF imagery have the same spectral bands (blue, green, red, and near-infrared), and the detailed information can be browsed on the website of CRESDA.Compared with GF imagery, RapidEye imagery with 5 m spatial resolution and one-day temporal resolution has similar visible and near-infrared bands, but one more red edge band.The preprocessing of GF-1 WFV images includes three protocols: geometric correction, radiometric calibration, and atmospheric correction.Geometric correction was conducted according to [37] with the assistance of ASTER GDEM V2 data.Radiometric calibration performs the conversion of image digital numbers to absolute at-sensor radiance.For each image, it was processed using the Environment for Visualizing Images (ENVI) 5.3 software (Harris Geospatial Solutions, Inc.: Broomfield, CO, USA) with the updating calibration parameters [42] published in CRESDA, obtained by a large number of calibration experiments in Chinese calibration fields.Atmospheric correction converts the radiance to reflectance.In this study, atmospheric correction was performed using the Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) model in ENVI and executed using the Interactive Data Language (IDL) (Harris Geospatial Solutions, Inc.: Redlands, CA, USA).Before correction, the cloud cover on the images were masked by visual interpretation.The related FLAASH parameters were obtained according to the acquisition time and imaging conditions.
The GF-2 PMS image is divided into the multispectral (MSS) image and panchromatic image.The preprocessing procedure of MSS image is same as GF-1 WFV images.The preprocessing of panchromatic image contains geometric correction and radiometric calibration.Geometric correction was conducted as same as MSS image.Radiometric calibration was conducted in ENVI by converting the digital numbers to the top of atmosphere reflectance.Afterwards, the Gram-Schmidt pan-sharpening method in ENVI was used to fuse the MSS and panchromatic images, which obtained a multispectral image with a 1 m resolution.The RapidEye images are L1C level products, which need preprocessing including geometric correction and atmospheric correction.The geometric correction and atmospheric correction procedures are same as GF-1 WFV images.
The fused GF-2 PMS image covering the overlap region among GF-2 PMS, RapidEye, and GF-1 WFV images was clipped (Figure 1c) and processed according to [25].The processed result was further refined in ArcGIS Desktop 10.5 software (Environmental Systems Research Institute, Inc.: Broomfield, CO, USA) by visual interpretation based on the field investigation of Shayang County in 2016.Then, the reference OR map (Figure 1d) in this region was obtained by converting the modified boundaries of OR fields to raster data with 1 m resolution.At last, the reference OR map were resampled to 4 m, 5 m, and 16 m resolution using pixel aggregate resampling method, which were used to assess the accuracy of CSRA-derived OR results on GF-2 PMS MSS, RapidEye, and GF-1 WFV images, respectively.

Ancillary Data
The ASTER GDEM V2 data, used for the geometric correction of satellite images, is a global digital elevation model (DEM) with 30 m resolution.The V2 data within the range of Hubei Province are downloaded from the Geospatial Data Cloud.
The phenology of OR in this study is obtained from the agro-meteorological stations in and near Hubei Province (purple hexagon in Figure 1a).The flowering period of OR divides into four stages: EFS (from the beginning of flowering to the flowers proportion less than 25% in the field), FBS (flowers proportion is more than 75%), BFS (flowers proportion is between 25% and 75%), and EnFS (flowers proportion is less than 25% and close to pod stage).Thus, a complete flowering period of OR are composed of EFS, BFS, FBS, BFS, and EnFS, successively.The flowering period in Hubei Province appears as an obvious gap in Yichang City.The left region of Yichang belonging to the upper reaches of the Yangtze River with a flowering period from early March to late March, whereas the right region of Yichang belongs to the middle reaches of the Yangtze River with a flowering period from mid-March to early April.The date gap between the two regions is nearly ten days.
The statistical planting areas of OR in Hubei Province from 2014 to 2016 (the data of 2017 will be published at the end of 2018) at provincial and municipal levels are used to validate the CSRA-derived OR planting maps.The statistical data are downloaded from the National Statistical Bureau of China and the Hubei Statistical Bureau.
Google Earth (GE) releases free images in high spatial resolution that may provide some potential for regional land use/cover mapping.Studies regarding land cover classification and validation using GE images achieved good effects [43,44].Thus, we employed the GE images to assist with the validation of OR classification results in this study.Historical GE images with 1 m resolution at FBS in part of Qianjiang City in 2015 (Figure 2a) and Wuxue City in 2016 (Figure 2c) were downloaded to validate the provincial OR planting maps.Registration of the GE images was performed based on GF-1 WFV images firstly, then object-oriented classification was conducted in eCognition 8.6 software to extract OR as the local reference OR map.The multi-resolution segmentation parameters for scale, shape, and compactness used in GE image processing were 30, 0.2, and 0.5, respectively.The nearest neighbor method classification method was used.
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 6 of 24 The phenology of OR in this study is obtained from the agro-meteorological stations in and near Hubei Province (purple hexagon in Figure 1a).The flowering period of OR divides into four stages: EFS (from the beginning of flowering to the flowers proportion less than 25% in the field), FBS (flowers proportion is more than 75%), BFS (flowers proportion is between 25% and 75%), and EnFS (flowers proportion is less than 25% and close to pod stage).Thus, a complete flowering period of OR are composed of EFS, BFS, FBS, BFS, and EnFS, successively.The flowering period in Hubei Province appears as an obvious gap in Yichang City.The left region of Yichang belonging to the upper reaches of the Yangtze River with a flowering period from early March to late March, whereas the right region of Yichang belongs to the middle reaches of the Yangtze River with a flowering period from mid-March to early April.The date gap between the two regions is nearly ten days.
The statistical planting areas of OR in Hubei Province from 2014 to 2016 (the data of 2017 will be published at the end of 2018) at provincial and municipal levels are used to validate the CSRAderived OR planting maps.The statistical data are downloaded from the National Statistical Bureau of China and the Hubei Statistical Bureau.
Google Earth (GE) releases free images in high spatial resolution that may provide some potential for regional land use/cover mapping.Studies regarding land cover classification and validation using GE images achieved good effects [43,44].Thus, we employed the GE images to assist with the validation of OR classification results in this study.Historical GE images with 1 m resolution at FBS in part of Qianjiang City in 2015 (Figure 2a) and Wuxue City in 2016 (Figure 2c) were downloaded to validate the provincial OR planting maps.Registration of the GE images was performed based on GF-1 WFV images firstly, then object-oriented classification was conducted in eCognition 8.6 software to extract OR as the local reference OR map.The multi-resolution segmentation parameters for scale, shape, and compactness used in GE image processing were 30, 0.2, and 0.5, respectively.The nearest neighbor method classification method was used.

Samples Collection
The land cover types in winter growth season of Hubei Province contain OR, WW, forest land, built-up land, water body, and bare land.Among them, crop species (OR and WW) and forest land make up the vegetation types.Built-up land, water bodies, and bare land compose non-vegetation types.In this study, sample collection for these land cover types is focused on Wuxue City, which is divided in two ways: field survey and visual interpretation.The field investigation with contents of the crop (OR and WW) was conducted in three experimental regions around Wuxue City in 2015 and 2016.The sampling method in each region was a simple random sample.The sample points (Figure 1b) of OR and WW were 2292 and 1669, respectively.Other land cover types, including forest land, built-up land, water bodies, and bare land, were collected on GF-1 WFV images by visual interpretation.These land cover types show obvious visual differences among each other, as well as when compared with crop species during the whole flowering period.The region of interest tool in ENVI was applied to collect samples of these land cover types with simple random sampling method.The sample points were evenly distributed over the entire image, and the number of forest land, built-up land, water bodeis, and bare land were 1662, 1922, 1810, and 621, respectively.For accuracy assessment of CSRA-derived results of six GF-1 WFV images of Wuxue City, the sample points of each type were divided into training samples and validation samples.Noticeably, the CSRA-derived results are a map with two classes (OR and NOR).Thus, WW, forest land, built-up land, water bodies, and bare land of the validation sample were merged into the NOR class.As a result, the pixel number of training samples of OR, WW, forest land, built-up land, water bodies, and bare land were 1278, 1069, 1512, 1802, 1710, and 577, respectively.The pixel number of the validation sample of OR and NOR were 1044 and 1014, respectively.

NDVI Model, HSV Transformation, and Normalization
VIs are quantitative measurements indicating the vigor of vegetation.The interest of these indices lies in their usefulness in the interpretation of remote sensing images and the improvement of classifications [45].Among them, NDVI is sensitive to the presence of green vegetation, the prediction of agricultural crops, and precipitation in semi-arid areas [45].In this study, NDVI is used to analyze the difference between vegetation and non-vegetation types.Equation ( 1) is the formula for calculating NDVI: where N IR, R represent the reflectance at near-infrared and red bands, respectively.Case studies have shown that color features are effective for crop discrimination [26,27].Compared with the RGB color space, the HSV color space has better performance to some extent in terms of decoupling chromaticity and luminance [28][29][30].Figure 3a describes the meaning of HSV components: H is defined as an angle between 0 • and 360 • representing a pure color (e.g., pure yellow, orange, or red, etc.); S refers to how much a pure color is diluted with white light; and V is the brightness of the color.The value of both S and V are between 0 and 1. Figure 3b shows the geometric relationship between the RGB and HSV systems.The transformation from the RGB color space to the HSV color space keeps the "black" vertex fixed at the origin, and the diagonal of the RGB color cube (named achromatic axis) is vertical so that the red, green, and blue axes are symmetrically arranged around the vertical axis that corresponds to the V axis of the HSV space [29].Equation ( 2) is used to perform the HSV transformation, which is a standardized transformation and a faster version as it does not contain any trigonometric functions or other functions requiring high-performance computation [29].Thus, it supports the efficiency of pixel-level calculation to a large extent of satellite imagery.
where R, G, B represent the reflectance at red, green, and blue bands, respectively.As the perceived color is determined by the shape of the (corrected) spectral reflectance, changes of HSV components are equivalent to the changes in land cover types.Especially, the change from one color to another color is interpreted as a land cover change [29].Thus, it is possible to associate a land cover type to a range of H, as well as the auxiliary decision by S and V.
where , ,  represent the reflectance at red, green, and blue bands, respectively.As the perceived color is determined by the shape of the (corrected) spectral reflectance, changes of HSV components are equivalent to the changes in land cover types.Especially, the change from one color to another color is interpreted as a land cover change [29].Thus, it is possible to associate a land cover type to a range of H, as well as the auxiliary decision by S and V.
The value of NDVI is between −1 and 1.The reflectance at blue, green, red, and near-infrared bands, S, and V are all between 0 and 1, whereas H is between 0 and 360.Moreover, H and V will be used to establish new indicators for extracting OR from other crop species.Thus, for further analysis, normalization for H is needed.And equation ( 3) is used to normalize H between 0 and 1: The value of NDVI is between −1 and 1.The reflectance at blue, green, red, and near-infrared bands, S, and V are all between 0 and 1, whereas H is between 0 and 360.Moreover, H and V will be used to establish new indicators for extracting OR from other crop species.Thus, for further analysis, normalization for H is needed.And equation ( 3) is used to normalize H between 0 and 1: where H means variable H with a value between 0 and 360, H norm means the normalized variable H with a value between 0 and 1.

Spectral and Color Characteristics of Land Cover Types on GF-1 WFV Images in Wuxue City
All the satellite images are conducted to execute the NDVI calculation, HSV transformation, and normalization.Then, the spectral feature (reflectance at each band, NDVI) and color feature (H norm , S, and V) values of the training sample on GF-1 WFV images of Wuxue City are exported.Afterwards, the average and standard deviation values of each land cover type are calculated to analyze the multi-temporal features.Finally, we found that NDVI, reflectance at near-infrared band (NIR), H norm , and V are sensitive to separate OR from other land cover types based on stepwise analysis.Figure 4 shows the multi-temporal features regarding spectral and color features applied to CSRA method.
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 9 of 24 V are sensitive to separate OR from other land cover types based on stepwise analysis.Figure 4 shows the multi-temporal features regarding spectral and color features applied to CSRA method.The multi-temporal NDVI in Figure 4a reveal that the NDVI values of vegetation types (OR, WW, and forest land) are always much larger than non-vegetation types (built-up land, water body, and bare land) even through the NDVI values of each type change over time.The rationale is that all these images are in plant growth season that vegetation species have large biomass, which can be reflected by NDVI.Thus, NDVI is used to separate non-vegetation types from vegetation types firstly.
Afterwards, we found there is a large difference in NIR among the vegetation types.As shown in Figure 4b, the NIR of crop types (OR and WW) are much larger than non-crop types at six flowering stages.The reason is that the crop species growing in Southern China in this period are in the reproductive growth stage, whose canopy, with dense leaf blade strongly reflecting the solar energy, while the forest type in this period are in the nutritional growth stage, with leaf blades starting growth.Thus, the reflectance in the near-infrared band of the forest type are small.Therefore, NIR is used to distinguish crop and non-crop types.
In the rest of the crop types, the reflectance at each band between OR and WW of images at EFS and EnFS have large similarity, caused by the flower proportion influencing the reflectance, as described in Section 1.However, we notice the difference between OR and WW at EFS and EnFS are still large in the HSV color space.In Figure 4c, H norm of WW is larger than OR at each stage, whereas V of WW is slightly smaller than OR.Thus, the ratio of V and H norm can be used to enlarge this difference.Additionally, the threshold of the ratio of V and H norm can be used to classify OR and WW.However, the differences at each flowering stage are not uniform, which would cause varying The multi-temporal NDVI in Figure 4a reveal that the NDVI values of vegetation types (OR, WW, and forest land) are always much larger than non-vegetation types (built-up land, water body, and bare land) even through the NDVI values of each type change over time.The rationale is that all these images are in plant growth season that vegetation species have large biomass, which can be reflected by NDVI.Thus, NDVI is used to separate non-vegetation types from vegetation types firstly.
Afterwards, we found there is a large difference in NIR among the vegetation types.As shown in Figure 4b, the NIR of crop types (OR and WW) are much larger than non-crop types at six flowering stages.The reason is that the crop species growing in Southern China in this period are in the reproductive growth stage, whose canopy, with dense leaf blade strongly reflecting the solar energy, while the forest type in this period are in the nutritional growth stage, with leaf blades starting growth.Thus, the reflectance in the near-infrared band of the forest type are small.Therefore, NIR is used to distinguish crop and non-crop types.
In the rest of the crop types, the reflectance at each band between OR and WW of images at EFS and EnFS have large similarity, caused by the flower proportion influencing the reflectance, as described in Section 1.However, we notice the difference between OR and WW at EFS and EnFS are still large in the HSV color space.In Figure 4c, H norm of WW is larger than OR at each stage, whereas V of WW is slightly smaller than OR.Thus, the ratio of V and H norm can be used to enlarge this difference.Additionally, the threshold of the ratio of V and H norm can be used to classify OR and WW.However, the differences at each flowering stage are not uniform, which would cause varying thresholds at each flowering stage and additional mixing (Figure 4d) when analyzing all the samples together.This mix can be solved by dividing the H norm -V space into three parts (part 1: V ≥ 0.07 and H norm ≤ 0.25; part 2: V ≥ 0.12 and 0.25 < H norm ≤ 0.42; and part 3: 0.07 ≤ V < 0.12 and 0.25 < H norm ≤ 0.42).The separability of OR and WW in each part is large.Thus, the ratio of V and H norm , named RRCI (ratio oilseed rape colorimetric index), is defined as Equation ( 4) to classify OR and WW in this study:

Workflow of Extracting Oilseed Rape
The workflow of applying CSRA method on remote sensing data to extract OR in this study is shown in Figure 5.It includes: (1) Data collecting and preprocessing, which has been introduced in Section 2.2.
(2) Developing the CSRA method based on GF-1 WFV-derived samples of Wuxue City.Section 2.3.2 has introduced the general idea of CSRA, which is a stepwise exclusive approach to achieve the goal of extracting OR.NDVI, NIR, and RRCI are successively used to separate vegetation from non-vegetation types, crop from non-crop types, and OR from WW.The threshold value in each step is obtained by the histogram thresholding method, which is the refined overlap point beyond two groups [3].Finally, the CSRA approach is performed to define the classification rules of the decision tree as its effectiveness and efficiency having been proved in the field of remote sensing classification [5,46].The CSRA method is executed using IDL.The acquisition of threshold values and detailed DT of the CSRA method will be reported in Section 3.1.
(3) Applying the CSRA method on remote sensing images to produce OR planting maps automatically.The OR maps includes: six GF-1 WFV-derived OR planting maps of Wuxue City, abbreviated as Result 1; GF-1 WFV-derived OR planting maps of Hubei Province from 2014 to 2017, abbreviated as Result 2; one GF-1 WFV-derived, one GF-2 PMS MSS-derived, and two RapidEye-derived OR planting maps of the overlap region of these four images, abbreviated as Result 3. (4) Validating the OR maps.The detailed validation contents regarding validating methods and dataset will be introduced in Section 2.3.4.

Results Validation
Result 1 and Result 3 are validated using the confusion matrix.Overall accuracy, kappa coefficient, user accuracy, and producer accuracy are used to assess the accuracy of the CSRA method.The validating samples of Result 1 are the validation samples introduced in Section 2.2.3.The validating sample of Result 3 are the OR and NOR pixels of the reference OR map in that region (Figure 1d).
The validation for Result 2 includes: local spatial consistency validation and area comparison against agricultural census data at provincial and municipal levels.The local spatial consistency is validated using confusion matrix basing on the GE image-derived OR reference maps (Figure 2b,d).The provincial level area comparison is validated using relative error between the CSRA-derived OR planting areas and the statistical areas of agricultural census data.Equation ( 5) is the formula for the relative error: RE = (Extraction − Reference)/Reference × 100% (5) where RE means the relative error, Extraction means the estimated areas, and Reference means the statistical areas.The municipal level area comparison is validated using decision coefficients (R 2 ) between the CSRA-derived OR planting areas and the statistical areas of agricultural census data.
beyond two groups [3].Finally, the CSRA approach is performed to define the classification rules of the decision tree as its effectiveness and efficiency having been proved in the field of remote sensing classification [5,46].The CSRA method is executed using IDL.The acquisition of threshold values and detailed DT of the CSRA method will be reported in Section 3.

Classification Rules of the CSRA Method
(1) Vegetation and non-vegetation types.Figure 6a-f shows the histogram results of NDVI between vegetation and non-vegetation types corresponding to 03/12, 03/17, 03/25, 03/28, 04/02, and 04/10, respectively.The threshold of NDVI for classifying vegetation from non-vegetation types is 0.3 at each flowering stage.Thus, the first step of the CSRA method is classifying vegetation pixels from non-vegetation pixels using the rule: NDVI ≥ 0.3.(2) Crop and non-crop types.Figure 7, regarding the histogram groups of NIR for crop and forest land, reveals that NIR is useful to distinguish them with a threshold of 0.23. Figure 7f shows some overlap with part of the FL samples having NIR values larger than 0.23, which is caused by the beginning of tree's nutritional growth stage since early April every year in Yangtze River Basin, China.These samples were gathered and analyzed in the next step.Thus, the second step of the CSRA method is separating the crop pixels from FL using the rule: NIR ≥ 0.23.(3) OR and WW. Figure 8 shows the performance of RRCI for distinguishing OR and WW, which has good separability in each part of the H norm -V space.The thresholds for parts 1-3 are 0.36, 0.43, and 0.25, respectively.In addition, the remaining forest land samples in step (2) can be divided into two groups: the first group with V values smaller than 0.07 and the second group located in part 3 of the H norm -V space.The RRCI values of the second group are all smaller than 0.2.Thus, the FL samples with NIR values larger than 0.23 have no confusion with OR.Therefore, the third step of the CSRA method is separating OR from WW using the rules: RRCI ≥ 0.36 when V ≥ 0.07 and H norm ≤ 0.25; RRCI ≥ 0.43 when V ≥ 0.12 and 0.25 < H norm ≤ 0.42; and RRCI ≥ 0.25 when 0.07 ≤ V < 0.12 and 0.25 < H norm ≤ 0.42.
between the CSRA-derived OR planting areas and the statistical areas of agricultural census data.

Classification Rules of the CSRA Method
(1) Vegetation and non-vegetation types.Figure 6a-f shows the histogram results of NDVI between vegetation and non-vegetation types corresponding to 03/12, 03/17, 03/25, 03/28, 04/02, and 04/10, respectively.The threshold of NDVI for classifying vegetation from non-vegetation types is 0.3 at each flowering stage.Thus, the first step of the CSRA method is classifying vegetation pixels from non-vegetation pixels using the rule: NDVI ≥ 0.3.(2) Crop and non-crop types.Figure 7, regarding the histogram groups of NIR for crop and forest land, reveals that NIR is useful to distinguish them with a threshold of 0.23. Figure 7f shows some overlap with part of the FL samples having NIR values larger than 0.23, which is caused by the beginning of tree's nutritional growth stage since early April every year in Yangtze River Basin, China.These samples were gathered and analyzed in the next step.Thus, the second step of the CSRA method is separating the crop pixels from FL using the rule: NIR ≥ 0.23.(3) OR and WW. Figure 8 shows the performance of RRCI for distinguishing OR and WW, which has good separability in each part of the H norm -V space.The thresholds for parts 1-3 are 0.36, 0.43, and 0.25, respectively.In addition, the remaining forest land samples in step (2) can be divided into two groups: the first group with V values smaller than 0.07 and the second group located in part 3 of the H norm -V space.The RRCI values of the second group are all smaller than 0.2.Thus, the FL samples with NIR values larger than 0.23 have no confusion with OR.Therefore, the third step of the CSRA method is separating OR from WW using the rules: RRCI ≥ 0.36 when V ≥ 0.07 and   ≤ 0.25; RRCI ≥ 0.43 when V ≥ 0.12 and 0.25 <   ≤ 0.42; and RRCI ≥ 0.25 when 0.07 ≤ V < 0.12 and 0.25 <   ≤ 0.42.According to [30], the H value of vegetation types is between 60 and 180 because of the "green peak" of reflectance for vegetation types.Thus, the rules of H norm ≥ 0.167 is needed.Finally, all the rules are used to construct the decision tree of the CSRA method for extracting OR, which is shown in Figure 9.
divided into two groups: the first group with V values smaller than 0.07 and the second group located in part 3 of the H norm -V space.The RRCI values of the second group are all smaller than 0.2.Thus, the FL samples with NIR values larger than 0.23 have no confusion with OR.Therefore, the third step of the CSRA method is separating OR from WW using the rules: RRCI ≥ 0.36 when V ≥ 0.07 and   ≤ 0.25; RRCI ≥ 0.43 when V ≥ 0.12 and 0.25 <   ≤ 0.42; and RRCI ≥ 0.25 when 0.07 ≤ V < 0.12 and 0.25 <   ≤ 0.42.According to [30], the H value of vegetation types is between 60 and 180 because of the "green peak" of reflectance for vegetation types.Thus, the rules of   ≥ 0.167 is needed.Finally, all the rules are used to construct the decision tree of the CSRA method for extracting OR, which is shown in Figure 9.

Mapping Oilseed Rape in Wuxue City Using the CSRA Method and Accuracy Assessment
The presented CSRA approach is firstly applied to generate the spatial distributions of OR in Wuxue City using the six GF-1 WFV images.Figure 10a-f, ranged chronologically for a more intuitive contrast, are the results corresponding to 03/17, 04/02, 04/10, 03/12, 03/25, and 03/28, respectively.As shown in Figure 10, OR is mainly planted in Southern Wuxue City because of the excellent growth conditions composing of the running Yangtze River and the flat topography, while the Northern part is mainly scattered and small fields due to hilly terrain.Figure 10 also shows that the planting areas of OR in Wuxue City from 2014 to 2015 are basically stable, whereas from 2015 to 2016 there appears to be an obvious decrease.In fact, the planting areas of OR in Wuxue City from 2014 to 2016 are 30, 28.96, and 20.7 Kha, respectively.Thus, the tendency reflecting OR planting areas are coincidental between the CSRA-derived results and the agricultural census data.

Mapping Oilseed Rape in Wuxue City Using the CSRA Method and Accuracy Assessment
The presented CSRA approach is firstly applied to generate the spatial distributions of OR in Wuxue City using the six GF-1 WFV images.Figure 10a-f, ranged chronologically for a more intuitive contrast, are the results corresponding to 03/17, 04/02, 04/10, 03/12, 03/25, and 03/28, respectively.As shown in Figure 10, OR is mainly planted in Southern Wuxue City because of the excellent growth conditions composing of the running Yangtze River and the flat topography, while the Northern part is mainly scattered and small fields due to hilly terrain.Figure 10 also shows that the planting areas of OR in Wuxue City from 2014 to 2015 are basically stable, whereas from 2015 to 2016 there appears to be an obvious decrease.In fact, the planting areas of OR in Wuxue City from 2014 to 2016 are 30, 28.96, and 20.7 Kha, respectively.Thus, the tendency reflecting OR planting areas are coincidental between the CSRA-derived results and the agricultural census data.
Table 2 lists the accuracy assessment results of these OR maps.The accuracy results have an increase from EFS to FBS whereas a decrease from FBS to EnFS.The rationale is that the RRCI separability between OR and other crop species also have an increase from EFS to FBS whereas there is a decrease from FBS to EnFS, since the color features are derived and affected by the reflectance at blue, green, and red bands.The overall accuracy of all results are larger than 89%, as well as kappa larger than 0.79.The producer accuracy of all results are larger than 91%.The user accuracy at EFS and EnFS are relatively worse because of the strong confusion between OR and WW at these phenology stages.However, the accuracy values are still larger than 85%.Thus, these accuracy results show good performance of the CSRA method on extracting OR.The results of Wuxue City demonstrated the effectiveness of the CSRA method on OR extracting at each flowering stages, which means the phenology problem is solved by the CSRA method.
The presented CSRA approach is firstly applied to generate the spatial distributions of OR in Wuxue City using the six GF-1 WFV images.Figure 10a-f, ranged chronologically for a more intuitive contrast, are the results corresponding to 03/17, 04/02, 04/10, 03/12, 03/25, and 03/28, respectively.As shown in Figure 10, OR is mainly planted in Southern Wuxue City because of the excellent growth conditions composing of the running Yangtze River and the flat topography, while the Northern part is mainly scattered and small fields due to hilly terrain.Figure 10 also shows that the planting areas of OR in Wuxue City from 2014 to 2015 are basically stable, whereas from 2015 to 2016 there appears to be an obvious decrease.In fact, the planting areas of OR in Wuxue City from 2014 to 2016 are 30, 28.96, and 20.7 Kha, respectively.Thus, the tendency reflecting OR planting areas are coincidental between the CSRA-derived results and the agricultural census data.Figure 11 shows the spatial distribution of OR in Hubei Province derived from GF-1 WFV images using the CSRA method.From 2014 to 2017, OR agriculture are concentrated in major regions near the Yangtze River and the Hanjiang River, whereas it occurs more sporadically in the western portion of the study area.The increased complexity of the topography in Western Hubei Province restricts the size of OR fields with most of the OR agriculture occurring in smaller river valleys or on terraced slopes.The largest planting area of OR are successively Jinzhou, Huanggang, and Jinmen, which corresponds to the major municipal City of Jianghan Plain, Eastern E Plain, and Central-North E mentioned in Section 2.1.The planting areas of OR estimated by the CSRA method from 2014 to 2017 are 1028.37,1003.39,965.23, and 811.83Kha, respectively.It has nearly no changes from 2014 to 2015, whereas there are successive decreases in 2016 and 2017.The decrease in 2016 mainly occurs in Tianmen, Xiantao, and Qianjiang.Cloud cover in Xiantao, Huangshi, and Jinmen in 2017 affect the accuracy of the results.However, Huanggang, Tianmen, and Suizhou have obviously decreased, except for those regions.The reason for the decrease are the increased planting costs and reduced market returns according to the agricultural announcement of the Hubei Provincial Department of Agriculture.
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 15 of 24 images is the second reason.Even in this situation, larger than 85% overall accuracy and nearly 90% producer accuracy still proves the CSRA-derived OR maps using GF-1 WFV images obtain high spatial consistency when compared with GE image-derived reference OR maps.The spatial matching of the OR planting maps is very important for practical applications.Local GE images are used to validate the spatial consistency of GF-1 WFV derived OR maps by the CSRA method.Figure 12a,d are the GE image-derived OR maps in part of Qianjiang City in 2015 and Wuxue City in 2016, respectively.Figure 12b,e shows the corresponding GF-1 WFV-derived OR maps of Figure 12a,d, respectively.The results show that the GF-1 WFV-derived OR maps have same spatial distribution compared with the GE image-derived OR maps.However, the GF-1 WFV-derived results have larger planting areas than the GE image-derived results due to the difference in spatial resolution.Furthermore, spatial analysis is conducted between the GF-1 WFV-derived and the resampled GE image-derived OR maps (Figure 12c,f).Additionally, the confusion matrix (Table 3) is performed to assess the accuracy using the spatial analysis results.The overall accuracy of Figure 12c is 88.68%, and Figure 12f is 85.31%.The user accuracy of GF-1 WFV-derived results are relatively low.The first reason for that is the difference of spatial resolution.For example, the GE image-derived OR maps can clearly extract OR fields and field borders, whereas the GF-1 WFV-derived results cannot distinguish fields borders from fields.The difficulty for co-registering GE images with GF-1 WFV images is the second reason.Even in this situation, larger than 85% overall accuracy and nearly 90% producer accuracy still proves the CSRA-derived OR maps using GF-1 WFV images obtain high spatial consistency when compared with GE image-derived reference OR maps.

Comparison of the Estimated Oilseed Rape Planting Areas with Agricultural Census Data
The comparison results (Table 4) between the GF-1 WFV-estimated OR planting areas and the agricultural census data at the provincial level confirm good accuracy.The absolute errors of the extracted annual total OR areas are less than 20%.The total OR planting areas of Jianghan Plain, Eastern E plain and Central-North E have provincial proportions of 78.47%, 76.24%, 77.67%, and 73.41% from 2014 to 2017, respectively, which are extremely well matched with the actual proportion of 80% published by the Hubei Provincial Department of Agriculture.The tendency of both estimated and statistical OR planting areas confirm a decrease occurring in 2016, which had the same conclusion with the analysis of Section 3.3.1.Furthermore, municipal-level validation of the OR planting areas extraction results is performed.The comparison results between the GF-1 WFV-estimated OR planting areas and the census data are shown in Figure 13.The solid line in the plot is the 1:1 line.The points in the plot are clustered near the 1:1 line, indicating that the GF-1 WFV-derived areas of OR is well correlated with the agricultural census data at the municipal level.The R 2 for 2014, 2015, and 2016 are 0.81, 0.73, and 0.76, respectively, which are all significant at the level of p < 0.01 (n = 17).The largest estimated mistakes of OR planting areas always appear in Shiyan, Enshi, and Xianning.One reason for that is the poor data availability caused by long-term cloud cover.Another reason is the omission error caused by the small fields in these mountain regions.The small fields are often processed as mixed pixels on GF-1 WFV images with uncertain reflectance.Generally, the GF-1 WFV-derived OR planting areas had good accuracy when comparing with agricultural census data.Furthermore, the tendency of the spatial distribution and local spatial consistency assessed in Section 3.3.1.also show a good effect.On the other hand, it is a time-consuming and labor-intensive task to implement a large-scale regional survey of OR planting areas and to obtain an annual spatial distribution map.Therefore, basing on the good spatial consistency, as well as the strong relationship between the GF-1 WFV-derived results and the agricultural census data from 2014 to 2017, the CSRA method can help to provide accurate OR area estimation and spatial distribution at the large-regional scale.Furthermore, municipal-level validation of the OR planting areas extraction results is performed.The comparison results between the GF-1 WFV-estimated OR planting areas and the census data are shown in Figure 13.The solid line in the plot is the 1:1 line.The points in the plot are clustered near the 1:1 line, indicating that the GF-1 WFV-derived areas of OR is well correlated with the agricultural census data at the municipal level.The R 2 for 2014, 2015, and 2016 are 0.81, 0.73, and 0.76, respectively, which are all significant at the level of p < 0.01 (n = 17).The largest estimated mistakes of OR planting areas always appear in Shiyan, Enshi, and Xianning.One reason for that is the poor data availability caused by long-term cloud cover.Another reason is the omission error caused by the small fields in these mountain regions.The small fields are often processed as mixed pixels on GF-1 WFV images with uncertain reflectance.Generally, the GF-1 WFV-derived OR planting areas had good accuracy when comparing with agricultural census data.Furthermore, the tendency of the spatial distribution and local spatial consistency assessed in Section 3.3.1.also show a good effect.On the other hand, it is a time-consuming and labor-intensive task to implement a large-scale regional survey of OR planting areas and to obtain an annual spatial distribution map.Therefore, basing on the good spatial consistency, as well as the strong relationship between the GF-1 WFV-derived results and the agricultural census data from 2014 to 2017, the CSRA method can help to provide accurate OR area estimation and spatial distribution at the large-regional scale.

Robustness Validation of the CSRA Method
Sections 3.2 and 3.3 have proved the effectiveness of the CSRA method on GF-1 WFV images.In this part, the CSRA method is used to explore its robustness using the GF-2 PMS MSS and RapidEye images.Figure 14 shows the spatial analysis results between the CSRA-derived OR maps and the reference OR map (Figure 1d). Figure 14a-d shows the results of the GF-2 PMS MSS image on 18 March 2016, the RapidEye images on 18 March 2016 (03/18) and 4 April 2016 (04/04), and the GF-1 WFV image on 28 March 2016, respectively.Figure 14e-h shows the corresponding results of the black rectanglar regions in Figure 14a-d.Table 5 shows the accuracy assessment results of these OR maps.These results show: (1) the GF-2 PMS MSS-derived result has the best spatial consistency, whereas GF-1 WFV has the worst; (2) all the misclassified pixels (omission and commission pixels) are distributed in the boundary area of OR fields; and (3) the RapidEye (04/04)-derived OR map has more omission and commission errors compared with the RapidEye (03/18)-derived result.The difference of the spatial consistency results between GF-2 PMS and GF-1 WFV are mainly caused by spatial resolution, which has been analyzed in Section 3.3.The GF-2 PMS MSS-derived OR map achieved 96.48% overall accuracy, which shows excellent spatial consistency with the reference OR map, and the GF-1 WFV-derived OR map also achieves 88.72% overall accuracy.The difference between the two RapidEye-derived results is caused by the acquisition date of the images.The images on 03/18 and 04/04 are at the FBS and EnFS, respectively, in which differences objectively exist when extracting OR.However, the accuracy assessment results (Table 5) show these differences are not large.The RapidEye-derived OR map on 04/04 achieves 91.8% overall accuracy, as well as 93.72% on 03/18.Both of these two images have good spatial consistency with the reference OR map.Thus, these results prove the good robustness of the CSRA method to some extent.spatial resolution, which has been analyzed in Section 3.3.The GF-2 PMS MSS-derived OR map achieved 96.48% overall accuracy, which shows excellent spatial consistency with the reference OR map, and the GF-1 WFV-derived OR map also achieves 88.72% overall accuracy.The difference between the two RapidEye-derived results is caused by the acquisition date of the images.The images on 03/18 and 04/04 are at the FBS and EnFS, respectively, in which differences objectively exist when extracting OR.However, the accuracy assessment results (Table 5) show these differences are not large.The RapidEye-derived OR map on 04/04 achieves 91.8% overall accuracy, as well as 93.72% on 03/18.Both of these two images have good spatial consistency with the reference OR map.Thus, these results prove the good robustness of the CSRA method to some extent.

Comparison with Previous Works
Table 6 lists the existing approaches for extracting OR and OR-liked crops using multispectral images.The previous studies have two characteristics: (1) based on spectral features (reflectance or VIs) during the flowering period, especially at the FBS; and (2) having a small study area.The limited study area and imagery acquisition data may generate maps that are difficult to compare for different regions and different flowering stages using these approaches [47].Moreover, large phenology differences exist caused by different cultivating dates and natural conditions, even on the same image [22,24].Thus, mapping OR using the existing method at the large-regional scale may suffer additional mistakes caused by phenology differences.In [21] and [25], provincial mapping approaches were developed.However, the images used in those studies were still at the FBS so that the phenology problem was still not solved.In short, regional mapping for OR using the existing methods may not achieve the goal.For the purpose of mapping OR at large-regional scale, we developed the CSRA method to solve the confusion of phenology difference and produced annual provincial OR planting maps.The results in the Sections 3.2 and 3.3 proved the effectiveness of the proposed method.The six GF-1 WFV images of Wuxue City are used to extract OR using the approaches in Table 6.The parameters of the SVM and MLC methods offer the best parameter association after multiple experiments.The parameters of VIs are obtained based on the methods applied in [24].The training and validation samples are same as Section 2.2.3.Moreover, NDVI ≥ 0.3 and NIR ≥ 0.23 are used to mask non-vegetation and non-crop pixels firstly so that the comparison conditions are fair.Finally, the methods including SVM, MLC, G, B/G, NDVI, and NGVI obtain better accuracy.Table 7 lists the accuracy results of these methods.The results show that: (1) all methods achieve the best accuracy on 03/25 and 03/28.The reason is that the spectral separability of OR and WW on images at the FBS are too large to classify; (2) The accuracy derived by VIs are not ideal at EFS and EnFS and change irregularly at different flowering stages.The rationale is that the VIs at EFS and EnFS show a large degree of confusion and the changing tendency of VIs are influenced by many factors; (3) Except for the CSRA method, the SVM method achieves better and stable accuracy.However, the SVM method suffers an accuracy decrease when using samples at FBS to train the images at EFS and EnFS, and using samples at EFS and EnFS to train the images at FBS and BFS.Therefore, if there were no images at the uniform flowering stage, the results might not be satisfactory when producing OR maps at the large-regional scale using the SVM, MLC, and VI methods.However, the fact is that the phenology stage of images for regional mapping are always not uniform.On the other hand, the CSRA method is proposed to solve this potential confusion caused by phenology differences.The overall accuracy and producer accuracy are excellent at each flowering stage.The user accuracy at EFS and EnFS are relatively smaller, but still larger than 85%.Therefore, the CSRA method has better performance on regional OR mapping than other previous methods.

Significance, Uncertainty Analysis, and Implications for Extensive Applications
This work indicates the advances of the developed CSRA method on OR extraction and regional mapping.As shown in the results, the proposed CSRA method has good effectiveness and robustness, and the provincial OR maps have outstanding local spatial consistency with GE images and excellent fitting relation with census data.In addition to the improvement of accuracy and robustness, the unsupervised CSRA method greatly improved the mapping efficiency due to its automatic ability.OR is one of the most important winter crops of China, especially in the Yangtze images often have mixed pixels that may lead to overestimated OR planting areas.The limited availability of cloud-free images for the whole flowering period is another problem when mapping large-regional OR planting areas.The GF-1 WFV image with two days of temporal resolution still suffered from the absence of cloud-free images in 2017.Associating other remote sensing images, like SPOT 6, Landsat 8, Sentinel-2, HJ-CCD, etc., can increase the observation frequency, which will increase the efficiency of the methodology.The robustness of the CSRA method on GF-2 PMS and RapidEye images is tested in this study.Noticeably, the spectral bands of GF images and RapidEye images are similar.Thus, the robustness on images from other sensors are still needed to validate the reuslts.Finally, all the spectral and color features of the CSRA method are calculated using surface reflectance.The reflectance may be influenced by cloud cover and other factors [22].Therefore, to apply the proposed approach to other climatic and topography regions, a few points need to be considered carefully, such as accurate and uniform atmosphere correction is needs to be preprocessed; appropriate selection of the images' acquisition date on the flowering period needs to be considered; and, if possible, more high-resolution images with similar spectral bands of GF-1 WFV images should be used.
Future work could be conducted on: (1) combining machine learning methods to develop automatic classifying methods to map annual winter crop planting patterns in the Yangtze River Basin and tracking the spatial-temporal dynamics; (2) incorporating further imagery sources, such as synthetic aperture radar (SAR) satellite images, to avoid the limitation of clouds on optical satellite images; (3) applying a fusion strategy on GF-1 WFV images with higher spatial resolution images to achieve a more accurate estimation for planting areas; and (4) other possible improvements and applications in various regions to evaluate and accommodate the wide variability of crop management and climate regions.

Conclusions
For the purpose of regional OR mapping, this research extended our previous work by considering the phenology difference at the large-regional scale and presented an innovative and generic method based on HSV transformation and spectral features.The proposed CSRA method was applied on a provincial study area, Hubei Province, to produce annual OR maps from 2014 to 2017 using GF-1 WFV images.Our produced OR maps achieved good spatial consistency with reference OR maps and accurate estimation of planting areas consistent with agricultural census data.Furthermore, the annual OR maps accurately presented the change tendency of OR planting in Hubei Province.The findings from the annual OR maps can help provide policy-makers with useful reference for precision agricultural management.Moreover, the good robustness of the CSRA method on GF-1 WFV-liked images were tested on GF-2 PMS and RapidEye images.Additionally, the comparison between the CSRA and other previous methods revealed the effectiveness of the CSRA method on extracting OR at the early-flowering stage and the end-flowering stage.This is important because of the limited availability of cloud-free images at the full-blooming stages and the phenology difference causing by different cultivating date and natural conditions when producing regional OR maps.In addition, the CSRA method is an unsupervised method, indicating that the efficiency of regional OR mapping can be improved.The results from this study highlight the ability of the CSRA method on extracting OR under complicated conditions and mapping regional OR distribution accurately, which is important for further research on regional mapping for crop planting patterns and for precise and smart agriculture.

Figure 1 .
Figure 1.(a) Hubei Province with OR experimental base, running river, agro-meteorological station, and topography; (b) Wuxue City with field survey points of OR and WW; (c) the resized GF-2 PMS fusion image of overlap region between GF-2 PMS image and RapidEye images in Shayang County; and (d) the corresponding reference OR map of (c).

Figure 1 .
Figure 1.(a) Hubei Province with OR experimental base, running river, agro-meteorological station, and topography; (b) Wuxue City with field survey points of OR and WW; (c) the resized GF-2 PMS fusion image of overlap region between GF-2 PMS image and RapidEye images in Shayang County; and (d) the corresponding reference OR map of (c).

Figure 2 .
Figure 2. Local (a) Google Earth image of part of Qianjiang City on 25 March 2015; (b) oilseed rape map corresponding to (a); (c) Google Earth image of part of Wuxue City on 19 March 2016; and (d) oilseed rape map corresponding to (c).The NOR in the legend is the abbreviation for non-oilseed rape.

Figure 2 .
Figure 2. Local (a) Google Earth image of part of Qianjiang City on 25 March 2015; (b) oilseed rape map corresponding to (a); (c) Google Earth image of part of Wuxue City on 19 March 2016; and (d) oilseed rape map corresponding to (c).The NOR in the legend is the abbreviation for non-oilseed rape.

24 Figure 3 .
Figure 3.The HSV color space about (a) H, S, and V components; (b) the geometric relationship of transformation from RGB to HSV.The reflectance at red, green, and blue bands of remote sensing images are used to make the R, G, and B color composite in this study.The values of R, G, and B are between 0 and 1.

Figure 3 .
Figure 3.The HSV color space about (a) H, S, and V components; (b) the geometric relationship of transformation from RGB to HSV.The reflectance at red, green, and blue bands of remote sensing images are used to make the R, G, and B color composite in this study.The values of R, G, and B are between 0 and 1.

Figure 9 .
Figure 9.The decision tree for extracting oilseed rape based on the CSRA approach.

Figure 12 .
Figure 12.(a) Google Earth-derived oilseed rape map, (b) CSRA-derived oilseed rape map, and (c) spatial comparison between (a) and (b) at part of Qianjiang City in 2015.(d) Google Earth-derived oilseed rape map, (e) CSRA-derived oilseed rape map, and (f) spatial comparison between (d) and (e) in part of Wuxue City in 2016.

Figure 12 .
Figure 12.(a) Google Earth-derived oilseed rape map, (b) CSRA-derived oilseed rape map, and (c) spatial comparison between (a) and (b) at part of Qianjiang City in 2015.(d) Google Earth-derived oilseed rape map, (e) CSRA-derived oilseed rape map, and (f) spatial comparison between (d) and (e) in part of Wuxue City in 2016.

Figure 12 .
Figure 12.(a) Google Earth-derived oilseed rape map, (b) CSRA-derived oilseed rape map, and (c) spatial comparison between (a) and (b) at part of Qianjiang City in 2015.(d) Google Earth-derived oilseed rape map, (e) CSRA-derived oilseed rape map, and (f) spatial comparison between (d) and (e) in part of Wuxue City in 2016.

Figure 13 .
Figure 13.Comparison of CSRA-estimated oilseed rape acreages and agricultural census data at the municipal level in (a) 2014; (b) 2015; and (c) 2016.

Table 1 .
The information of GF-1 WFV, GF-2 PMS, and RapidEye images used to build and verify the CSRA method.EFS, FBS, BFS, and EnFS in the table are the abbreviations of early flowering stage, full blooming stage, blossoming stage, and end flowering stage, respectively.

Table 2 .
Accuracy assessment for CSRA-derived oilseed rape maps in Wuxue City.

Table 3 .
Confusion matrix of local accuracy assessment using Google Earth images.

Table 4 .
Comparison of oilseed rape planting areas derived using the CSRA algorithm and from the agricultural census data.

Table 5 .
Accuracy assessment for the CSRA-derived oilseed rape maps of GF-2 PMS MSS, RapidEye, and GF-1 WFV images in part of Shayang County.

Table 6 .
Comparison of the existing OR extracting approaches using multispectral images.
Note: MDC: minimum distance classifier; MLC: maximum likelihood classifier; SVM: support vector machines; ARTMAP: adaptive resonance theory mappings; DVI: difference vegetation indices; B/G: the ratio of reflectance at blue and green bands; RVI: ratio vegetation indices.The phenology stage in China was obtained by local agro-meteorological stations.

Table A1 .
The information of GF-1 WFV images used to produce OR maps of Hubei Province.