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 39 encompasses whole or portions of several states including South Dakota, North Dakota, Minnesota, and Iowa. Questions about the NLCD mapping zone 39 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 39 are as follows:
SPRING-
Index 1 for Path 27/Row 29 on 04/22/03 = Scene_ID 7027029000311250
Index 2 for Path 28/Row 29 on 03/03/00 = Scene_ID 7028029000006350
Index 3 for Path 28/Row 30 on 04/21/03 = Scene_ID 5028030000311110
Index11 for Path 28/Row 30 on 09/30/01 = Scene_ID 7028030000127350
Index 4 for Path 29/Row 28 on 05/03/02 = Scene_ID 7029028000212350
Index 4 for Path 29/Row 29 on 05/03/02 = Scene_ID 7029029000212350
Index 4 for Path 29/Row 30 on 05/03/02 = Scene_ID 7029030000212350
Index12 for Path 29/Row 30 on 10/20/00 = Scene_ID 7029030000029450
Index 5 for Path 30/Row 28 on 05/18/02 = Scene_ID 5030028000213810
Index13 for Path 30/Row 28 on 08/27/01 = Scene_ID 7030028000123950
Index 6 for Path 30/Row 29 on 05/04/00 = Scene_ID 7030029000012550
Index14 for Path 30/Row 29 on 05/18/02 = Scene_ID 5030029000213810
Index 7 for Path 30/Row 30 on 04/11/03 = Scene_ID 7030030000310150
Index 8 for Path 31/Row 28 on 05/19/00 = Scene_ID 5031028000014010
Index 9 for Path 31/Row 29 on 05/01/99 = Scene_ID 5031029009912110
Index10 for Path 32/Row 28 on 05/02/00 = Scene_ID 7032028000012350
LEAF ON (Summer)-
Index 1 for Path 27/Row 29 on 07/05/01 = Scene_ID 7027029000118650
Index 2 for Path 28/Row 29 on 08/10/00 = Scene_ID 7028029000022350
Index 3 for Path 28/Row 30 on 08/13/01 = Scene_ID 7028030000122550
Index 4 for Path 29/Row 28 on 08/01/00 = Scene_ID 7029028000021450
Index 5 for Path 29/Row 29 on 08/04/01 = Scene_ID 7029029000121650
Index 5 for Path 29/Row 30 on 08/04/01 = Scene_ID 7029030000121650
Index 6 for Path 30/Row 28 on 08/24/00 = Scene_ID 7030028000023750
Index11 for Path 30/Row 28 on 07/13/02 = Scene_ID 7030028000219450
Index 7 for Path 30/Row 29 on 08/27/01 = Scene_ID 7030029000123950
Index12 for Path 30/Row 29 on 06/16/01 = Scene_ID 5030029000116710
Index 7 for Path 30/Row 30 on 08/27/01 = Scene_ID 7030030000123950
Index13 for Path 30/Row 30 on 06/16/01 = Scene_ID 5030030000116710
Index 8 for Path 31/Row 28 on 07/14/00 = Scene_ID 7031028000019650
Index 9 for Path 31/Row 29 on 07/30/00 = Scene_ID 7031029000021250
Index10 for Path 32/Row 28 on 07/03/02 = Scene_ID 5032028000218410
LEAF-OFF (Fall)-
Index 1 for Path 27/Row 29 on 11/10/01 = Scene_ID 7027029000131450
Index 2 for Path 28/Row 29 on 10/27/99 = Scene_ID 7028029009930050
Index 3 for Path 28/Row 30 on 09/30/01 = Scene_ID 7028030000127350
Index 4 for Path 29/Row 28 on 10/20/00 = Scene_ID 7029028000029450
Index 4 for Path 29/Row 29 on 10/20/00 = Scene_ID 7029029000029450
Index 4 for Path 29/Row 30 on 10/20/00 = Scene_ID 7029030000029450
Index 5 for Path 30/Row 28 on 09/25/00 = Scene_ID 7030028000026950
Index 9 for Path 30/Row 28 on 07/13/02 = Scene_ID 7030028000219450
Index 5 for Path 30/Row 29 on 09/25/00 = Scene_ID 7030029000026950
Index 5 for Path 30/Row 30 on 09/25/00 = Scene_ID 7030030000026950
Index 6 for Path 31/Row 28 on 09/27/01 = Scene_ID 5031028000127010
Index 7 for Path 31/Row 29 on 11/06/01 = Scene_ID 7031029000131050
Index10 for Path 31/Row 29 on 07/30/00 = Scene_ID 7031029000021250
Index 8 for Path 32/Row 28 on 09/26/01 = Scene_ID 7032028000126950
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)