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