地理空间数据云数据覆盖不了一个市怎么办

这次的主要任务是学会下载哨兵二号的遥感影像数据，学长推荐了几个常用的网站如GEE、EOS、USGS、地理空间数据云，我计划把几个网站都尝试一下，下面就按顺序进行尝试：

一、GEE

网址：https://code.earthengine.google.com

GEE(Google Earth Engine)，即谷歌地球引擎,据说是世界上使用量最大的云平台（如果真是这样的话，那说明遥感和地信的分量确实很重了）。

登录这个网站需要科学上网，首先需要用一个谷歌账号登录(关于谷歌账号的注册门道也有点多，这里就不细说了），然后用该账号提出申请，申请界面如下：

虽然这里的邮箱使用的是qq邮箱，但还是很不建议使用，最好使用学校的官方邮箱，申请成果的可能性会大一些。

机构就选择所在的大学，类型选择acadamic即可，国家选择China，卜拆尺至于最后的“申请账号型高来干什么”，我当时是先用百度翻译把我的申请理由给翻译成英文后再复制上去的，具体理由就是说一下要做研究之类的。我当时这样申请之后可以说是秒通过了（邮箱里会收到通知），这也是在之前失败了n次之后的结果。

来到主界面，是这个样子的：

首先导入武汉的矢量文件，注意是压缩包（zip）格式：

通过如下代码可实现导入：

导入结果：

在资源搜索栏中寻找哨兵二号的数据集，并选择导入（IMPORT）：

可以看到代码编辑区为我们准备好了代码：

接下来使用如下代码下载影像：

var district = table

var dsize = district.size()

var district_geometry = district.geometry()

print(dsize)

Map.centerObject(district_geometry,7)

//去云

function maskS2clouds(image) {

var qa = image.select('QA60')

// Bits 10和11分别是云和卷云

var cloudBitMask = 1 <<10

var cirrusBitMask = 1 <<御丛 11

// Both flags should be set to zero, indicating clear conditions.

var mask = qa.bitwiseAnd(cloudBitMask).eq(0)

.and(qa.bitwiseAnd(cirrusBitMask).eq(0))

return image.updateMask(mask).divide(10000)

}

//处理数据，生成数据集

var dataset = ee.ImageCollection('COPERNICUS/S2')

.filterBounds(district_geometry)//区域筛选

.filterDate('2020-01-01', '2020-12-30')//时间范围

// 选择云量

.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 1))

.map(maskS2clouds)

.select(['B8', 'B4', 'B3'])

//可视化参数

var rgbVis = {

min: 0.0,

max: 0.3,

bands: ['B8', 'B4', 'B3'],//按照8波段、4波段、3波段放入通道合成

}

var IMAGE = dataset.mosaic()

//下载数据

Export.image.toDrive({

image:IMAGE,

description:'image_really_none',//设置下载任务tasks的名称

scale:10,//分辨率是10m

maxPixels: 1e13,//单幅影像输出的最大像元数

region:district_geometry,//下载范围

fileFormat: 'GeoTIFF',//设置影像导出格式，注意GeoTIFF格式需要所有波段存储类型一致（如不能同时存Int16和Int32）

formatOptions: {

cloudOptimized: true

}

})

可以在右上角的Tasks栏中看到加载的任务：

点击“Open in Drive”,可以从云端下载影像：

将下载后的影像加载到ENVI中进行查看：

二、USGS

网址：EarthExplorer (usgs.gov)

USGS是美国地质调查局的官方网站，同样需要科学上网，我这里使用的是学校邮箱注册的官网账号，登陆后就可以进行数据下载了。

首先，我还是通过想通过矢量数据直接定位武汉，所以就选择"KML/Shapefile Upload",再选择Shapefile格式，加载文件（注意这里加载的文件必须是zip压缩文件形式）：

但是很遗憾的是出现了如下状况：

这里显示我的武汉市范围矢量文件不只一个多边形，我有些奇怪，就把文件加载到arcgis里看了下，果然在右上角发现了一个小范围的多边形：

所以我选择从行政区划的角度来选择，在左侧栏中输入关键字“wuhan”即可得到结果：

为了找出下载的具体范围，我们需要进行选定，这里方法比较多，比如说下图中的画圆法，先点击一下可得圆心，再点击一下可得半径：

也可以选择多边形法，这样画的范围可能更为准确:

下面还可以选择影像时间，默认为全部时间段的影像：

接下来选择数据源，这里是哨兵二号：

接下来可以选择一些具体的参数，如云量：

接下来就可以检索到符合条件的一些影像数据了，可以选择第一个像脚印一样的符号来把影像加载到地图上来查看范围（自动会用不同颜色标识）：

选择第五个图标即可进行下载：

三、地理空间数据云

网址：地理空间数据云 (gscloud.cn)

在这个上面下数据是比较简单的了，之前做遥感实习时的数据也主要是从这上面下的。打开官网后选择“高级检索”，可以看到界面是非常简洁明了的：

选择数据集为哨兵二号，选择用矢量文件来划定区域：

比较奇怪的是，这里用矢量数据划定区域后找不到合适的影像（可能和之前在USGS下载时的问题一样，由于矢量文件不只一个多边形所导致）：

所以这里就采取另一种方法：先根据行政区划找到覆盖武汉市影像，然后再用武汉市范围的矢量数据进行剪裁，这就需要找到能够覆盖武汉市全部范围的图像，总共是3张（之前做实习下的是Landsat8的数据，也是三张即可覆盖武汉全境），注意要检索能够下载的数据，因为很多数据能被检索到但我们却不能下载：

总结

哨兵二号的数据可以从很多网站上下载得到，除了上面介绍的这些还有很多，担上面这三种应当可以满足大部分时候的需求了。

1下载《将shp文件导入到mysql的命令》程序，解压后打开，

2、将shp以弯亏塌及shp的相关文件和DOShere的doc文件放在同一目空春录

3、doc命令为：

shp2mysql shp文件名.shp 表名 数据库名 >生成的sql文件名埋圆.sql

4、修改生成的sql文件：

ALTER TABLE 表名 ADD the_geom GEOMETRY

add必须为：the_geom GEOMETRY

替换所有的,-1) )为) )

5、到mysql数据库导入sql文件即可

---