小編研究了一下之前學習的java，打開似曾相識的idea，想到在學校老師叫我們實現(xiàn)的一個編程作業(yè)，讀取遙感波段圖像的特征值（最值、中值、眾數(shù)、平均值、標準差之類的）。眾所周知，遙感波段數(shù)據(jù)就是一個二維數(shù)組，里面存儲的是每一個像元的DN值（灰度值也叫像元值），讀取圖像的特征值就是統(tǒng)計這些DN值的特征值。
為了讓運算量少一點，小編提取了Landsat5 TM數(shù)據(jù)的3波段和4波段的一個矩形柵格，像元數(shù)量是64*30，那么最棘手的問題就是怎么將這些像元值讀取到一個數(shù)組里面，我們想到使用arcmap里面的【轉(zhuǎn)換工具】|【從柵格轉(zhuǎn)出】|【柵格轉(zhuǎn)Ascll】，我們打開轉(zhuǎn)換出來的其中一個Ascll.TXT文件如下圖所示：
我們主要的操作對象是第7行開始的這64*30個整數(shù)。
package 實習;
import java.io.*;
import java.lang.annotation.ElementType;
import java.nio.charset.Charset;
import java.time.temporal.Temporal;
import java.util.ArrayList;
import java.util.Arrays;
import java.util.List;
public class 讀取遙感圖像 {
public static void main(String[] args) throws IOException {//創(chuàng)建一個文件對象 int[] a = 讀入像元值(new File("D:\\桌面\\3波段.txt"));//存放波段3的ASCLL碼值路徑 int[] b = 讀入像元值(new File("D:\\桌面\\4波段.txt"));//存放波段4的ASCLL碼值路徑 int[] C = 復制數(shù)組(a); int[] D = 復制數(shù)組(b); NDVI(C, D, new File("D:\\桌面\\4波段.txt")); double[] A = 統(tǒng)計特征值(a, "D:\\桌面\\3波段.txt");//對數(shù)組a進行的計算,并且傳回了平均數(shù)和標準差 double[] B = 統(tǒng)計特征值(b, "D:\\桌面\\4波段.txt");//對數(shù)組b進行了計算并且傳回了平均數(shù)和標準差 多波段圖像統(tǒng)計(C, D, A, B);}public static int[] 復制數(shù)組(int[] a) { int[] C = new int[a.length]; for (int i = 0; i < a.length; i++) { C[i] = a[i]; } return C;}static int w = 0;//用來控制是否追加字符用的public static int[] 讀入像元值(File q) throws IOException { BufferedReader b = new BufferedReader(new FileReader(q)); StringBuffer buffer = new StringBuffer();//Stringbuffer是動態(tài)字符串數(shù)組,可以放很多字符串 String value; int a = -7;//統(tǒng)計行數(shù) do { value = b.readLine();//循環(huán)一次就保存一行字符串（代表著一行像元值） a++;//保存柵格的行數(shù) if (a >= 0 && value != null)//第7行開始讀取 buffer.append(value);//append()是往動態(tài)字符串添加新的字符串，每循環(huán)一次存一行，兩行之間有空格 } while (value != null); String[] temp = buffer.toString().split(" "); int[] number = new int[temp.length]; //temp里面的N個字符串代表所有的像元值 if (w == 0) { FileWriter writer = new FileWriter("D:\\桌面\\遙感圖像統(tǒng)計值.txt");//第一次寫出直接清空之前的文件內(nèi)容 writer.write(q + "有" + a + "行," + temp.length / a + "列像元,一共" + temp.length + "個像元。" + "\n"); writer.close(); } else { FileWriter writer = new FileWriter("D:\\桌面\\遙感圖像統(tǒng)計值.txt", true);//第二次寫出則追加到之前的文本 writer.write(q + "有" + a + "行," + temp.length / a + "列像元,一共" + temp.length + "個像元。" + "\n"); writer.close(); } w++;//控制是否是第一次輸入 for (int i = 0; i < temp.length; i++) { if (Integer.parseInt(temp[i]) >= 0) number[i] = Integer.parseInt(temp[i]);//字符串數(shù)組轉(zhuǎn)換成整數(shù)數(shù)組
// if (i % ((temp.length / a)) == 0)
// System.out.print("\n" + “第” + (i / (temp.length / a) + 1) + “行：”);
// System.out.print(number[i] + " ");
} return number;}public static double[] 統(tǒng)計特征值(int[] b, String a) throws IOException { int e;//計算中位數(shù)用的 int k = b.length;//所有大于0的像元個數(shù) int t = 0;//用于輔助比較各個像元值出現(xiàn)的次數(shù) int s = 0;//存放出現(xiàn)次數(shù)最多的像元值（眾數(shù)） int sum = 0;//統(tǒng)計所有享元年值總和 for (int i = 0; i < b.length; i++) {//這里將值為-9999都篩選掉了 sum += b[i]; } for (int i = 0; i < k; i++) {//比較n-1次就夠了 for (int j = 0; j < k - 1; j++) {//每一個數(shù)的比較都比上一次少一次 if (b[j] > b[j + 1]) {//如果前面一個比后面一個大，那么大的都往后退 int temp = b[j]; b[j] = b[j + 1]; b[j + 1] = temp;//大的值都給了后面 } } } int[] y = new int[(b[k - 1]) + 1];//存放出現(xiàn)的灰度級 int[] m = new int[(b[k - 1]) + 1];//存放各個灰度級出現(xiàn)的次數(shù) FileWriter writer = new FileWriter("D:\\桌面\\遙感圖像統(tǒng)計值.txt", true); double[] f = new double[2]; f[0] = sum / k;//將平均值保存到浮點數(shù)組中 writer.write("\n" + a + "該遙感圖像DN值最大為:" + b[k - 1] + "，最小為：" + b[0] + ",平均數(shù)是：" + f[0] + "\n"); //計算均方根誤差(標準差) for (int i = 0; i < k; i++) { f[1] += ((sum / k) - b[i]) * ((sum / k) - b[i]); } f[1] = Math.sqrt(f[1] / k); writer.write("均方根誤差(標準差)是：" + f[1] + "\n"); //計算灰度級中位數(shù) if (k % 2 == 0) e = (b[k / 2] + b[k / 2 + 1]) / 2; else e = b[k / 2 + 1]; writer.write("中位數(shù)是：" + e + "\n"); //計算眾數(shù) for (int i = b[0]; i <= b[k - 1]; i++) {//范圍是最小值和最大值之間的唯一值 y[i] = i;//由于排序了，所以第一個肯定是最小的灰度級 for (int j = 0; j < k; j++) {//統(tǒng)計對象是所有的像元值 if (y[i] == b[j])//如果該灰度級等于其中一個像元值 m[i]++;//這個灰度級中的像元數(shù)量加一 } if (t < m[i]) {//不斷比較m t = m[i];//t是頻數(shù)最大的那個 s = y[i];//把頻數(shù)最大的灰度級賦給s } } writer.write("眾數(shù)是：" + s + "\n"); //直方圖 for (int i = b[0]; i <= b[k - 1]; i++) {//范圍是最小值和最大值之間的唯一值 y[i] = i;//由于排序了，所以第一個肯定是最小的灰度級 if (y[i] < 10) {//如果是個位數(shù)的灰度級就加兩個空格 writer.write(" "); } if (y[i] > 9 && y[i] < 100) {//如果是十位數(shù)的灰度級就加一個空格 writer.write(" "); } writer.write(y[i] + ":"); for (int j = 0; j < m[i]; j++) { writer.write("■"); } writer.write(m[i] + "\n");//m是每個灰度級出現(xiàn)的像元數(shù)量（頻數(shù)） } writer.close(); return f;}public static void NDVI(int[] c, int[] d, File m) throws IOException { float[] e = new float[c.length]; BufferedReader br = new BufferedReader(new FileReader(m)); BufferedWriter bw = new BufferedWriter(new FileWriter("D:\\桌面\\NDVI計算結(jié)果.txt")); //建立兩個緩沖區(qū)提高讀寫數(shù)據(jù)的效率 //可以換格式復制，很厲害 String str;//用來臨時存放數(shù)據(jù) for (int i = 0; i < 6; i++) { str = br.readLine(); bw.write(str);//把讀取的內(nèi)容寫進去 bw.newLine();//讀完一行就開辟新的一行 } for (int i = 0; i < c.length; i++) { e[i] = (float) (d[i] - c[i]) / (float) (d[i] + c[i]); bw.write(e[i] + " "); } bw.close();}public static void 多波段圖像統(tǒng)計(int[] a, int[] b, double[] A, double[] B) throws IOException { FileWriter writer = new FileWriter("D:\\桌面\\遙感圖像統(tǒng)計值.txt", true); double o = 0; for (int i = 0; i < a.length; i++) { o += (a[i] - A[0]) * (b[i] - B[0]); } o = o / a.length; writer.write("兩波段的協(xié)方差為：" + o + "\n"); o = o / (A[1] * B[1]); writer.write("兩波段的相關(guān)系數(shù)為：" + o); writer.close();}
}成都途遠GIS是一家專業(yè)致力于無人機航空攝影測繪、航空數(shù)據(jù)處理、GIS地理信息系統(tǒng)研發(fā)、數(shù)字孿生城市制作、數(shù)字沙盤模型等業(yè)務的創(chuàng)新型科技公司，為您提供一站式地理信息服務。
本文鏈接：http://www.dzac.com.cn/blog/465.html
本文標簽：遙感
