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 Space Imaging LLC 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 66 encompasses whole or portions of several states including Maine, New Hampshire and Vermont. Questions about the NLCD mapping zone 66 can be directed to the NLCD 2001 land cover mapping team at the National Center, 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 66 are as follows:
SPRING-
Index 1 for Path 10/Row 29 on 05/17/03 = Scene_ID 7010029000313750
Index 2 for Path 11/Row 27 on 05/16/03 = Scene_ID 5011027000313610
Index 2 for Path 11/Row 28 on 05/16/03 = Scene_ID 5011028000313610
Index 3 for Path 11/Row 29 on 04/29/00 = Scene_ID 7011029000012050
Index10 for Path 11/Row 30 on 04/14/03 = Scene_ID 5011030000310410
Index 6 for Path 12/Row 27 on 05/12/02 = Scene_ID 7012027000213250
Index 8 for Path 12/Row 28 on 05/25/01 = Scene_ID 7012028000114550
Index 7 for Path 12/Row 29 on 05/04/02 = Scene_ID 5012029000212410
Index 9 for Path 12/Row 30 on 05/09/01 = Scene_ID 7012030000112950
Index 4 for Path 13/Row 29 on 05/08/01 = Scene_ID 5013029000112810
Index 4 for Path 13/Row 30 on 05/08/01 = Scene_ID 5013030000112810
Index 5 for Path 14/Row 29 on 05/07/01 = Scene_ID 7014029000112750
LEAF ON (Summer)-
Index 1 for Path 10/Row 29 on 06/28/01 = Scene_ID 7010029000117950
Index 2 for Path 11/Row 27 on 08/09/02 = Scene_ID 7011027000222150
Index 2 for Path 11/Row 28 on 08/09/02 = Scene_ID 7011028000222150
Index 2 for Path 11/Row 29 on 08/09/02 = Scene_ID 7011029000222150
Index 2 for Path 11/Row 30 on 08/09/02 = Scene_ID 7011030000222150
Index 3 for Path 12/Row 27 on 07/20/01 = Scene_ID 5012027000120110
Index 3 for Path 12/Row 28 on 07/20/01 = Scene_ID 5012028000120110
Index 4 for Path 12/Row 29 on 05/25/01 = Scene_ID 7012029000114550
Index 5 for Path 12/Row 30 on 06/21/02 = Scene_ID 5012030000217210
Index 6 for Path 13/Row 29 on 06/04/02 = Scene_ID 7013029000215550
Index 7 for Path 13/Row 30 on 08/31/99 = Scene_ID 7013030009924350
Index 8 for Path 14/Row 29 on 06/08/01 = Scene_ID 7014029000115950
LEAF-OFF (Fall)-
Index 1 for Path 10/Row 29 on 10/02/01 = Scene_ID 7010029000127550
Index 2 for Path 11/Row 27 on 10/04/02 = Scene_ID 5011027000227710
Index 3 for Path 11/Row 28 on 10/01/01 = Scene_ID 5011028000127410
Index 4 for Path 11/Row 29 on 10/22/00 = Scene_ID 7011029000029650
Index 5 for Path 12/Row 27 on 09/30/01 = Scene_ID 7012027000127350
Index 5 for Path 12/Row 28 on 09/30/01 = Scene_ID 7012028000127350
Index 6 for Path 12/Row 29 on 10/16/01 = Scene_ID 7012029000128950
Index 7 for Path 12/Row 30 on 10/16/01 = Scene_ID 7012030000128950
Index 8 for Path 13/Row 29 on 10/12/00 = Scene_ID 5013029000028610
Index 9 for Path 13/Row 30 on 10/20/00 = Scene_ID 7013030000029450
Index 9 for Path 14/Row 29 on 10/12/03 = Scene_ID 5014029000328510
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)