操作记录文档：WUDAPT 进行LCZ分类

WUDAPT

网站链接：http://www.wudapt.org/
已有LCZ图下载链接：https://wudapt.cs.purdue.edu/

LCZ分类具体步骤

WUDAPT流程是一种监督分类的方法，需要勾画训练样本，WUDAPT提供了40个城市的样本数据，如果我们的研究区域不包括在其中，就需要自己进行训练区域勾画。

准备工作

首先，需要下载所需要的Landsat影像，USGS下载即可
之后，需要下载Google Earth和SAGA，网站里有提供下载链接
还有一些数据文件的下载
整个分类文件夹的下载：https://drive.google.com/drive/folders/0B83nwq2eGktbUGxSV196a2tpYTA
打开之后下载city_template.zip，解压

数据下载

我们以武汉为例，下载了2020年的Landsat8的一景影像。

定义研究区域

大概就是选取影像的一部分进行研究
链接：http://www.wudapt.org/prepfeat_overview/path2step3b/
在Google Earth里进行操作，之后在SAGA里转换投影

数据准备

链接：http://www.wudapt.org/prepfeat_overview/path2step3c/
主要是裁剪、重采样，完成之后保存至feature文件夹
不需要第8、9波段

勾画训练样本

链接：http://www.wudapt.org/create-lcz-training-areas/step-4/
这步比较重要，也比较烦人

分类

网站似乎更新了，找不到旧版分类的过程，在这里记录一下。
打开SAGA GIS软件和处理过的landsat影像

Geoprocessing→Imagery→Classification→Local Climate Zone Classification

分类图另存为GeoTIFF

分类精度评价

参考：https://desktop.arcgis.com/zh-cn/arcmap/latest/manage-data/raster-and-images/accuracy-assessment-for-image-classification.htm
分类图需要转换为16位
参考：https://blog.csdn.net/weixin_40625478/article/details/106898712

在Google Earth标注GroundTruth区域，得到KML文件，打开ArcGIS，使用KMLTo Layer
然后，转换为raster

这里会出现一个问题，Raster文件的Value与Name不对应，需要使用Reclassify工具，value 重分类为1 2 3 4.。。。101 102 。。。

然后就OK了。
下一步 Create Accuracy Assessment Points


下一步Update Accuracy Assessment Points

附：
LCZ各类别颜色RGB值

	  # set color for each value
	  colors.SetColorEntry(1, (140,   0,       33))
	  colors.SetColorEntry(2, (209,   0,        0))
	  colors.SetColorEntry(3, (255,   0,        0))
	  colors.SetColorEntry(4, (191,   77,       0))
	  colors.SetColorEntry(5, (255,   102,      0))

	  colors.SetColorEntry(6, (255,   153,      85))
	  colors.SetColorEntry(7, (250,   238,      5))
	  colors.SetColorEntry(8, (188,   188,    188))
	  colors.SetColorEntry(9, (255,   204,    170))
	  colors.SetColorEntry(10, (85,    85,     85))

	  colors.SetColorEntry(11, (0,   106,      0))
	  colors.SetColorEntry(12, (0,   170,     0))
	  colors.SetColorEntry(13, (100,   133,      37))
	  colors.SetColorEntry(14, (185,   219,    121))
	  colors.SetColorEntry(15, (0,     0,    0))

	  colors.SetColorEntry(16, (251,   247,    174))
	  colors.SetColorEntry(17, (106,   106,    255))