土地利用根据县域矢量边界按掩膜提取

Question

中科院全国范围的土地利用类型图，按掩膜提取出一个县的土地利用类型图，可以用同一个矢量边界数据提取多个不同年份的吗？还是要对应年份的矢量边界数据提取？
比较准确的县级矢量边界数据在哪些网站上怎么获取？

Answer 1

让阿豪来帮你解答，本回答参考chatgpt3.5编写提供，如果还有疑问可以评论或留言

问题回答：

1. 可以用同一个矢量边界数据提取多个不同年份的土地利用类型图。因为矢量边界数据通常是按照行政区划，如县级行政区划，划分出来的，这些行政区划的边界一般都是不变的，故在同一个县内提取不同年份的土地利用类型图时，可以使用同一个矢量边界数据来进行提取。
2. 要获取准确的县级矢量边界数据，可以到国家测绘地理信息局（www.nsdic.org.cn）或者各省市测绘局的网站上进行下载。以北京市为例，可以登录北京市测绘局的网站（www.bjmap.gov.cn）进行获取。 具体代码如下（以提取2001年和2010年北京市昌平区的土地利用类型图为例）：

import rasterio
import fiona
from rasterio.mask import mask
# 打开矢量边界数据
with fiona.open('changping.shp', 'r') as shapefile:
    shapes = [feature['geometry'] for feature in shapefile]
# 打开2001年土地利用类型图
with rasterio.open('land_use_2001.tif') as src:
    out_image, out_transform = mask(src, shapes, crop=True)
# 输出提取的土地利用类型图
with rasterio.open('land_use_2001_changping.tif', 'w', **src.meta) as dst:
    dst.write(out_image)
# 打开2010年土地利用类型图
with rasterio.open('land_use_2010.tif') as src:
    out_image, out_transform = mask(src, shapes, crop=True)
# 输出提取的土地利用类型图
with rasterio.open('land_use_2010_changping.tif', 'w', **src.meta) as dst:
    dst.write(out_image)

其中，changping.shp为昌平区的矢量边界数据，land_use_2001.tif和land_use_2010.tif为对应年份的土地利用类型图。运行以上代码可分别提取2001年和2010年北京市昌平区的土地利用类型图，并保存在当前目录下。