基于AI Earth的威宁草海水质反演

威宁草海水质反演

贵州威宁草海国家级自然保护区, 位于贵州省西部威宁县县城西南面,保护区面积120平方公里, 其中水域面积46.5平方公里。是一个完整、典型的高原湿地生态系统、是黑颈鹤等228种鸟类的重要越冬地和迁徙中转站;是中国著名三大高原湖泊(草海、滇池,青海湖)之一、贵州最大的高原天然淡水湖泊、中国Ⅰ级重要湿地、国家4A级旅游景区;是世界十大观鸟基地,被美国国家地理杂志评选为世界上最受欢迎的旅游胜地。

初始化环境

import aie
aie.Authenticate()
aie.Initialize()

检索数据

# 导入roi
shp = aie.FeatureCollection('user/17bf545043ad434aa03f516c8bb311eb')
geometry = shp.geometry()
# geometry = geometry.buffer(1500)
map = aie.Map(
    center=geometry.getCenter(),
    height=800,
    zoom=11
)
vis_params = {
    'color': '#00FF00'
}
map.addLayer(
    geometry,
    vis_params,
    'shapefile',
    bounds=geometry.getBounds()
)
map

[图片上传失败...(image-30c9b-1669469460447)]


# 指定检索数据集,可设置检索的空间和时间范围,以及属性过滤条件(如云量过滤等)
dataset = aie.ImageCollection('SENTINEL_MSIL2A') \
             .filterBounds(geometry) \
             .filterDate('2020-05-13', '2021-05-15') \
             .filter(aie.Filter.lte('eo:cloud_cover',10.0))
# 显示影像
# vis_params = {
#     'bands': ['B4', 'B3', 'B2'],
#     'min': 0,
#     'max': 2000,
# }
# map.addLayer(
#     dataset,
#     vis_params,
#     'True Color (432)',
#     bounds=dataset.getBounds()
# )
# map

水体反演指数计算

image = dataset.select(["B2", "B3", "B4"]).mosaic().clip(geometry)

# date = image.get("date").getInfo();  
b2 = image.select("B2") 
b3 = image.select("B3") 
b4 = image.select("B4") 

# 叶绿素α浓度
chl1 = ((b4.divide(b3)).pow(aie.Image(2))).multiply(aie.Image(4.089))
chl2 = (b4.divide(b3)).multiply(aie.Image(0.476))                                                   
chlorophyll = chl1.subtract(chl2).add(aie.Image(29.733)) 

# 悬浮物浓度
suspension = ((b3.divide(b2)).pow(aie.Image(6.0823))).multiply(aie.Image(119.62))
# 透明度
transparency = suspension.pow(aie.Image(-0.67)).multiply(aie.Image(284.15))

# 基于叶绿素α浓度的营养状态指数
TLIchia = (chlorophyll.multiply(aie.Image(10.86))).add(aie.Image(2.5))
# 基于透明度的营养状态指数
TLIsd = aie.Image(51.18).subtract(transparency.log().multiply(aie.Image(19.4)))
chlorophyll_vis  = {
    'min': 30,
    'max': 38,
    'palette': ['#ff2b18','#fe9a0c','#f9e71c','#d2fe79',
               '#8afccb','#39dfff','#6b85fe','#5813fc']
}


map.addLayer(
    chlorophyll,
    chlorophyll_vis,
    'Chlorophyll',
    bounds=chlorophyll.getBounds()
)

map
叶绿素α浓度
transparency_vis  = {
    'min': 0,
    'max': 10,
    'palette': ['#c2523c','#e18321','#f5c60b','#d6fa00',
               '#20e000','#16b569','#198090','#0b2c7a']
}


map.addLayer(
    transparency,
    transparency_vis,
    'Transparency',
    bounds=transparency.getBounds()
)

map
透明度
# unitScale操作,归一化操作
TLIchia_unit = TLIchia.unitScale(333, 408)
TLIsd_unit = TLIsd.unitScale(-9,63)
# max操作
Water_quality = TLIchia_unit.max(TLIsd_unit)
quality_vis  = {
    'min': 0,
    'max': 0.4,
    'palette': ['#5813fc','#6a70fd','#5eb7fe','#08fcfd','#95fcc3',
                '#ccfd83','#f2fe2a','#fdc521','#ff860e','#ff2b18']
}


map.addLayer(
    Water_quality,
    quality_vis,
    'WaterQuality',
    bounds=Water_quality.getBounds()
)

map
水体质量分布

水体质量分类

# Water_quality是水体污染指数,意思就是都是反过来,Water_quality三分类,
# 将quality小于等于0.39的区域设置为1,
# 将quality大于0.39小于等于0.52的区域设置为2,将quality大于0.52设置为3

Water_class = Water_quality.where(Water_quality.lte(aie.Image(0.39)),aie.Image(1))\
                           .where(Water_quality.gt(aie.Image(0.39)).And(Water_quality\
                                                      .lte(aie.Image(0.52))),aie.Image(2))\
                           .where(Water_quality.gt(aie.Image(0.52)),aie.Image(3))\


# map = aie.Map(
#     center=geometry.getCenter(),
#     height=800,
#     zoom=11
# )


# 正常水体区域为蓝色
vis_params = {
    'min': 1.0,
    'max': 3.0,
    'palette' : [
        '#97dbf2', '#ffffbe', '#fc3019'
    ]
}


map.addLayer(
    Water_class,
    vis_params,
    'WaterClass',
    bounds=Water_class.getBounds()
)
map
水体质量分级空间分布
# 导出数据
# task = aie.Export.image.toAsset(chlorophyll,'chlorophyll',10)
# task.start()

总结

利用AI Earth可用很快的进行遥感计算,利用哨兵2号影像进行对威宁草海进行水质反演,当然我这里没有具体的数据,也不知道做的效果怎么样,大家有数据的,可用帮我去看看,速度特别快的。后面附上ArcGIS Pro绘图的成果。
主要引用
[1]孙微,曹宇佳,马新文,魏丹.基于Sentinel-2和ZH-1遥感影像的松北区水质监测[J].地理信息世界,2021,28(01):118-123.


水体透明度空间分布图.jpg
水体叶绿素α浓度空间分布图.jpg
水体综合污染等级空间分布图.jpg
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,039评论 6 498
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,223评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 161,916评论 0 351
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,009评论 1 291
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,030评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,011评论 1 295
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,934评论 3 416
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,754评论 0 271
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,202评论 1 309
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,433评论 2 331
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,590评论 1 346
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,321评论 5 342
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,917评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,568评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,738评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,583评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,482评论 2 352

推荐阅读更多精彩内容