Homer, C., C. Huang, L. Yang, B. Wylie and M. Coan, 2004. Development of a 2001 national land cover database for the United States. Photogrammetric Engineering and Remote Sensing.
The USGS acknowledges the support of MDA Federal Inc in development of data in this zone.
The NLCD 2001 was created by partitioning the U.S. into mapping zones. A total of 66 mapping zones were delineated within the conterminous U.S. based on ecoregion and geographical characteristics, edge matching features and the size requirement of Landsat mosaics. Mapping zone 31 encompasses whole or portions of several states including South Dakota and Nebraska. Questions about the NLCD mapping zone 31 can be directed to the NLCD 2001 land cover mapping team at the USGS/EROS, Sioux Falls, SD (605) 594-6151 or mrlc@usgs.gov.
By combining the training data with Landsat spectral data and other ancillary data, percent imperviousness prediction models were developed using a regression tree algorithm (named "Cubist"). All data layers used as input for running the Cubist in order to establish the prediction model for imperviousness are listed in the "Supplemental Information" section. Three ancillary data sets were used for urban/suburban delineation masking, which is essential to confine the imperviousness mapping within the developed areas.
After the best combination of possible input data layers had been determined and the final percent impervious layer produced, the urban mask created was used to eliminate those pixels from the final percent impervious file that fall outside identified impervious surface areas. Thus, the impervious pixels identified within the mask were retained, while all non-impervious surface pixels were removed from the final product. This corrected for any small areas of pixels that may have been included as part of the impervious surface layer because their spectral signature was not covered in the Cubist modeling process. Finally, visual inspection of imperviousness layer was made with limited manual editing to eliminate non-urban areas based on area of interest delineated by the mapping team.
Acquisition dates of Landsat ETM+ (TM) scenes used for land cover classification in zone 31 are as follows:
SPRING-
Index 1 for Path 29/Row 31 on 05/03/02 = Scene_ID 7029031000212350
Index 2 for Path 29/Row 32 on 05/13/00 = Scene_ID 7029032000013450
Index 3 for Path 30/Row 30 on 04/11/03 = Scene_ID 7030030000310150
Index 4 for Path 30/Row 31 on 05/04/00 = Scene_ID 7030031000012550
Index 4 for Path 30/Row 32 on 05/04/00 = Scene_ID 7030032000012550
Index 5 for Path 31/Row 29 on 05/01/99 = Scene_ID 5031029009912110
Index 6 for Path 31/Row 30 on 05/14/01 = Scene_ID 7031030000113450
Index 7 for Path 31/Row 31 on 04/28/01 = Scene_ID 7031031000111850
Index 8 for Path 31/Row 32 on 04/09/00 = Scene_ID 7031032000010050
Index 9 for Path 32/Row 29 on 05/02/00 = Scene_ID 7032029000012350
Index 9 for Path 32/Row 30 on 05/02/00 = Scene_ID 7032030000012350
Index 9 for Path 32/Row 31 on 05/02/00 = Scene_ID 7032031000012350
Index10 for Path 33/Row 29 on 05/12/01 = Scene_ID 7033029000113250
Index10 for Path 33/Row 30 on 05/12/01 = Scene_ID 7033030000113250
Index11 for Path 33/Row 31 on 04/13/02 = Scene_ID 7033031000210350
LEAF ON (Summer)-
Index 1 for Path 29/Row 31 on 06/30/00 = Scene_ID 7029031000018250
Index 2 for Path 29/Row 32 on 07/03/01 = Scene_ID 7029032000118450
Index 3 for Path 30/Row 30 on 08/08/00 = Scene_ID 7030030000022150
Index12 for Path 30/Row 30 on 08/24/00 = Scene_ID 7030030000023750
Index 4 for Path 30/Row 31 on 07/07/00 = Scene_ID 7030031000018950
Index 4 for Path 30/Row 32 on 07/07/00 = Scene_ID 7030032000018950
Index 5 for Path 31/Row 29 on 07/30/00 = Scene_ID 7031029000021250
Index 5 for Path 31/Row 30 on 07/30/00 = Scene_ID 7031030000021250
Index 6 for Path 31/Row 31 on 07/14/00 = Scene_ID 7031031000019650
Index 7 for Path 32/Row 29 on 07/03/02 = Scene_ID 5032029000218410
Index 8 for Path 32/Row 30 on 07/03/99 = Scene_ID 7032030009918450
Index 9 for Path 32/Row 31 on 08/06/00 = Scene_ID 7032031000021950
Index10 for Path 33/Row 29 on 07/12/00 = Scene_ID 7033029000019450
Index10 for Path 33/Row 30 on 07/12/00 = Scene_ID 7033030000019450
Index13 for Path 33/Row 30 on 09/20/02 = Scene_ID 7033030000226350
Index11 for Path 33/Row 31 on 07/28/00 = Scene_ID 7033031000021050
LEAF-OFF (Fall)-
Index 1 for Path 29/Row 31 on 09/21/01 = Scene_ID 7029031000126450
Index 1 for Path 29/Row 32 on 09/21/01 = Scene_ID 7029032000126450
Index 2 for Path 30/Row 30 on 10/25/99 = Scene_ID 7030030009929850
Index 3 for Path 30/Row 31 on 11/15/01 = Scene_ID 7030031000131950
Index 4 for Path 30/Row 32 on 09/12/01 = Scene_ID 7030032000125550
Index 5 for Path 31/Row 29 on 11/06/01 = Scene_ID 7031029000131050
Index 6 for Path 31/Row 30 on 10/18/00 = Scene_ID 7031030000029250
Index 6 for Path 31/Row 31 on 10/18/00 = Scene_ID 7031031000029250
Index 6 for Path 31/Row 32 on 10/18/00 = Scene_ID 7031032000029250
Index 7 for Path 32/Row 29 on 09/26/01 = Scene_ID 7032029000126950
Index 8 for Path 32/Row 30 on 10/09/00 = Scene_ID 7032030000028350
Index 9 for Path 32/Row 31 on 09/10/01 = Scene_ID 7032031000125350
Index10 for Path 33/Row 29 on 09/20/02 = Scene_ID 7033029000226350
Index11 for Path 33/Row 30 on 10/16/00 = Scene_ID 7033030000029050
Index11 for Path 33/Row 31 on 10/16/00 = Scene_ID 7033031000029050
Landsat data and ancillary data used for the imperviousness prediction -
Data Type of DEM composed of 1 band and a Continuous Variable Type.
Data Type of Slope composed of 1 band and a Continuous Variable Type.
Data Type of Aspect composed of 1 band of Categorical Variable Type.
Data type of Position Index composed of 1 band of Continuous Variable Type.
Ancillary data used for urban area masking include:
1. NOAA Defense Meteorological Satellite Program 1 km city light data (1 band)
2. MRLC 1992 classification (1 band)
3. TIGER 2000 roads (1 band)