本次內容為landsat-8遙感影像中水體的提取，用到了new water index(nwi)指數，效果較好。
文獻（丁鳳，2009）中還花一定篇幅對大氣校正進行闡述，意在體現其重要性。但個人認為如果只是對分類或提取，大氣校正其實是雞肋的（真香警告），所以本次辨證地對toa（大氣層頂反射率）計算結果nwi(toa)和sr（地表反射率）計算結果nwi(sr)進行了對比：在案例區域（長江黃岡段）中，nwi(toa)只要大於0就可以很好地提取出水體，而nwi(sr)則要自己調整閾值。但在後面的普適性驗證中，發現nwi(toa)>0這樣的條件雖然很簡單，實則並不具備普適性，容易將山體陰影也劃分為水體，所以最後還是在案例區域找到了適合nwi(sr)的閾值，並在我國大部分區域（華東，華南，華北，東北，西北，西南，台灣北，台灣南）做了驗證，發現-0.6這個閾值很適合landsat-8地表反射率產品計算的nwi指數做水體提取，對於快速提取水體有很大幫助。有興趣的可以自行驗證。（注：earth engine中的sr產品需要除10000才是小於1的反射率值）
以東北地區驗證為例，結果如下，藍色即為水體：
不足：提取結果中包含水田或蓄水池，且對於單像元的水體也可以進行提取，所以提取結果存在碎斑化現象。
/*
* extrction of water
* using landsat-8 imagery
* * author: deserts tsung, nov 20，2018
* chengdu university of information technology
**/// function of water extraction
// following new water index(nwi) was proposed by ding feng in 2009
var waterextraction = function(inputimage) ;
var sr = ee.image('landsat/lc08/c01/t1_sr/lc08_122039_20181026').divide(10000);
var wtrimg = waterextraction(sr);
map.addlayer(sr,,'sr');
map.addlayer(wtrimg.select('b5'),,'water');
// east p.r.china verified
// south p.r.china verified
// north p.r.china verified
// north-east area verified
// south-west area verified
// north-west area verified
// north taiwan area verified
// south taiwan area verified
下水道體現城市智慧型
近些年來，隨著急速的城市化，一些城市建設中忽視了與大自然和諧相處，怠慢了基礎設施尤其是城市排水和儲水設施的建設，遇上暴雨連綿，一些城市即成澤國。隨著不斷的城市化，大量水泥 玻璃 瀝青之類的非生態建築材料，以及人口加劇向城市集中，城市面積越來越大，使得城市中心密集區的環境問題越來越嚴重。大城市洪澇災害...
小評一下 水果豆漿機
啊哦，本來要睡覺了，突然想起今天sxj和qj同學的blog。幾個鳥人在發揮他們那齷齪的想象力，褻瀆著。哈哈，來吧出力的用biscuit 砸我吧。水果豆漿機的故事，過幾天再貼。個人思考 1 豆漿機精神已經過時了，宣告不是不需要了。豆漿機精神是上世紀末的精神 2 悶頭做出來的豆漿機 泛指，鬱悶總有人只認...
下水道井蓋為什麼是圓的？
事情是這樣子的，理查德 範曼同學在微軟找工作，下面是面試官同他的對話 面試官 現在我們要問乙個問題，看看你的創造性思維能力。不要想得太多，運用日常生活中的常識，描述一下你的想法。這個問題是，下水道的井蓋為什麼是圓的？範曼 它們並不都是圓的，有些是方的，的確有些圓井蓋，但我也看過方的，長方的。面試官 ...