浏览代码

i.atcorr: Worldview3 support added (filter functions with permission from apollomapping.com)

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@69012 15284696-431f-4ddb-bdfa-cd5b030d7da7
Markus Neteler 9 年之前
父节点
当前提交
24cfa0f754
共有 4 个文件被更改,包括 4057 次插入3 次删除
  1. 3 2
      imagery/i.atcorr/create_iwave.py
  2. 1934 1
      imagery/i.atcorr/iwave.cpp
  3. 18 0
      imagery/i.atcorr/iwave.h
  4. 2102 0
      imagery/i.atcorr/sensors_csv/worldview3.csv

+ 3 - 2
imagery/i.atcorr/create_iwave.py

@@ -51,8 +51,9 @@ def read_input(csvfile):
     should be a .csv file with the values
     should be a .csv file with the values
     of the filter function for each band in the sensor
     of the filter function for each band in the sensor
     one column for band
     one column for band
-    first line must have a header with sensor band name
+    first line must have a header with sensor band names
     first column is wavelength
     first column is wavelength
+    values are those of the discrete band filter functions
     """
     """
     infile = open(csvfile, 'r')
     infile = open(csvfile, 'r')
         
         
@@ -60,7 +61,7 @@ def read_input(csvfile):
     bands = infile.readline().split(',')
     bands = infile.readline().split(',')
     bands.remove(bands[0])
     bands.remove(bands[0])
     bands[-1] = bands[-1].strip()
     bands[-1] = bands[-1].strip()
-    
+    print "Number of bands found: %d" % len(bands)
     infile.close()
     infile.close()
     
     
     # create converter dictionary for import
     # create converter dictionary for import

文件差异内容过多而无法显示
+ 1934 - 1
imagery/i.atcorr/iwave.cpp


+ 18 - 0
imagery/i.atcorr/iwave.h

@@ -170,6 +170,23 @@ c        130 Green    "               ( 0.530-0.590 )                  c
 c        131 Red      "               ( 0.625-0.695 )                  c 
 c        131 Red      "               ( 0.625-0.695 )                  c 
 c        132 NIR      "               ( 0.760-0.890 )                  c 
 c        132 NIR      "               ( 0.760-0.890 )                  c 
 c        133 Pan      "               ( 0.480-0.830 )                  c 
 c        133 Pan      "               ( 0.480-0.830 )                  c 
+c        134 Pan  band of Worldview 3
+c        135 Coastal_MS7  "           
+c        136 Blue_MS4 "
+c        137 Green_MS3
+c        138 Yellow_MS6
+c        139 Red_MS2  "
+c        140 Red Edge_MS5
+c        141 NIR1_MS1 "
+c        142 NIR2_MS8 "
+c        143 SWIR1    "
+c        144 SWIR2    "
+c        145 SWIR3    "
+c        146 SWIR4    "
+c        147 SWIR5    "
+c        148 SWIR6    "
+c        149 SWIR7    "
+c        150 SWIR8    "
 c                                                                      c
 c                                                                      c
 c  note: wl has to be in micrometer                                    c
 c  note: wl has to be in micrometer                                    c
 c**********************************************************************/
 c**********************************************************************/
@@ -220,6 +237,7 @@ private:
 	void spot7(int iwa);
 	void spot7(int iwa);
 	void pleiades1a(int iwa);
 	void pleiades1a(int iwa);
 	void pleiades1b(int iwa);
 	void pleiades1b(int iwa);
+	void worldview3(int iwa);
 
 
 public:
 public:
 	/* To compute the equivalent wavelength needed for the calculation of the
 	/* To compute the equivalent wavelength needed for the calculation of the

文件差异内容过多而无法显示
+ 2102 - 0
imagery/i.atcorr/sensors_csv/worldview3.csv