CN114080540A - 信息处理装置、信息处理方法、程序和感测系统 - Google Patents
信息处理装置、信息处理方法、程序和感测系统 Download PDFInfo
- Publication number
- CN114080540A CN114080540A CN202080046901.8A CN202080046901A CN114080540A CN 114080540 A CN114080540 A CN 114080540A CN 202080046901 A CN202080046901 A CN 202080046901A CN 114080540 A CN114080540 A CN 114080540A
- Authority
- CN
- China
- Prior art keywords
- measurement
- macro
- section
- information processing
- processing apparatus
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000010365 information processing Effects 0.000 title claims abstract description 128
- 238000003672 processing method Methods 0.000 title claims description 6
- 238000005259 measurement Methods 0.000 claims abstract description 566
- 238000004364 calculation method Methods 0.000 claims abstract description 151
- 238000004458 analytical method Methods 0.000 claims abstract description 139
- 238000012821 model calculation Methods 0.000 claims abstract description 119
- 238000001514 detection method Methods 0.000 claims abstract description 103
- 238000012545 processing Methods 0.000 claims description 39
- 238000000034 method Methods 0.000 claims description 32
- 230000007613 environmental effect Effects 0.000 claims description 25
- 230000008569 process Effects 0.000 claims description 24
- 238000013507 mapping Methods 0.000 claims description 23
- 230000004044 response Effects 0.000 claims description 21
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 claims description 12
- 229930002875 chlorophyll Natural products 0.000 claims description 12
- 235000019804 chlorophyll Nutrition 0.000 claims description 12
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims description 12
- 239000002689 soil Substances 0.000 claims description 9
- 230000003595 spectral effect Effects 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 230000010287 polarization Effects 0.000 claims description 7
- 238000003786 synthesis reaction Methods 0.000 claims description 7
- 238000003384 imaging method Methods 0.000 description 37
- 241000196324 Embryophyta Species 0.000 description 35
- 238000005516 engineering process Methods 0.000 description 24
- 230000006870 function Effects 0.000 description 20
- 239000000872 buffer Substances 0.000 description 16
- 238000004891 communication Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 14
- 239000003086 colorant Substances 0.000 description 10
- 230000000704 physical effect Effects 0.000 description 8
- 230000029553 photosynthesis Effects 0.000 description 7
- 238000010672 photosynthesis Methods 0.000 description 7
- 230000011218 segmentation Effects 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 7
- 238000009434 installation Methods 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 230000002194 synthesizing effect Effects 0.000 description 5
- 230000012010 growth Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 208000036855 Left sided atrial isomerism Diseases 0.000 description 2
- 238000004497 NIR spectroscopy Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000035790 physiological processes and functions Effects 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 238000001429 visible spectrum Methods 0.000 description 2
- 208000003643 Callosities Diseases 0.000 description 1
- 206010020649 Hyperkeratosis Diseases 0.000 description 1
- 235000007688 Lycopersicon esculentum Nutrition 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 240000003768 Solanum lycopersicum Species 0.000 description 1
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005401 electroluminescence Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000007773 growth pattern Effects 0.000 description 1
- 230000012447 hatching Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C39/00—Aircraft not otherwise provided for
- B64C39/02—Aircraft not otherwise provided for characterised by special use
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/66—Arrangements or adaptations of apparatus or instruments, not otherwise provided for
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N2021/1793—Remote sensing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N2021/1793—Remote sensing
- G01N2021/1797—Remote sensing in landscape, e.g. crops
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N2021/3595—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N2021/8466—Investigation of vegetal material, e.g. leaves, plants, fruits
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Quality & Reliability (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Radiology & Medical Imaging (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Image Processing (AREA)
Abstract
提供了一种信息处理装置,包括宏观测量分析计算部件,被配置为计算来自宏观测量部件的检测数据,宏观测量部件适于以第一空间分辨率对测量对象的第一测量范围进行感测;微观测量分析计算部件,被配置为计算来自微观测量部件的检测数据,微观测量部件适于以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及逆模型计算部件,被配置为基于由微观测量分析计算部件确定的来自微观测量部件的检测数据,获取用于使用来自宏观测量分析计算部件的计算结果的逆模型计算的模型参数。
Description
相关申请的交叉引用
本申请要求2019年7月3日提交的日本优先权专利申请JP2019-124764的权益,其全部内容通过引用并入本文。
技术领域
本技术涉及信息处理装置、信息处理方法、程序和感测系统,特别地,涉及适合于生成植被状态等的测量结果的技术。
背景技术
例如,已经做出了很多努力来通过使用安装在小型飞行体(如无人机)中的成像装置,在飞行体在农田上空飞行时对植物的植被状态进行成像,来遥感植被状态。
PTL 1公开了一种通过对农田成像来执行遥感的技术。
引文列表
专利文献
PTL 1
日本专利号5162890
发明内容
技术问题
除了基于可见光(R(红色)、G(绿色)和B(蓝色))测量的形状测量之外,这种遥感还允许使用各种光波长和方法测量对象的物理属性、生理状态等。然而,可安装在小飞行体中的感测设备通常在尺寸、重量等方面受到限制。
例如,可以获取大量波长并执行分量分析等的超谱相机通常可能需要配置为获取二维图像的扫描机构,并且尺寸较大。因此,超谱相机很难安装在小尺寸的无人机或类似设备上。
另一方面,为了使用人造卫星进行感测,安装了能够执行高级感测的先进设备的人造卫星正在运行。然而,就空间分辨率而言,这种感测是不够的。
由于在一个空间分辨率单位中混合了各种对象,低空间分辨率不仅排除了形状的确定,而且还妨碍了仅对待检查对象的测量。
为了覆盖低空间分辨率,特定对象的测量值是使用包含关于测量对象形式的信息的模型(例如,将“形状”与“特征/环境响应”关联的模型等,具体而言,例如辐射传输特性模型)基于逆模型计算(逆计算)确定的。然而,在现实中,测量对象和模型之间的形状等的差异使得精确测量变得困难。
因此,期望提供一种系统以及用于该系统的信息处理装置,该系统能够在基于使用例如人造卫星的高性能传感器的遥感中获得更精确的测量结果。
对问题的解决方案
根据本技术的一个实施例的信息处理装置包括宏观测量分析计算部件,被配置为计算来自宏观测量部件的检测数据,宏观测量部件适于以第一空间分辨率对测量对象的第一测量范围进行感测;微观测量分析计算部件,被配置为计算来自微观测量部件的检测数据,微观测量部件适于以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及逆模型计算部件,被配置为基于由微观测量分析计算部件确定的来自微观测量部件的检测数据,获取用于使用来自宏观测量分析计算部件的计算结果的逆模型计算的模型参数。
用于逆计算的模型参数基于第二空间分解的测量结果生成。
另外,在根据上述本技术的实施例的信息处理装置中,逆模型计算部件可以使用基于由微观测量分析计算部件确定的来自微观测量部件的检测数据的模型参数,作为使用来自宏观测量分析计算部件的计算结果的逆模型计算中的逆模型的参数。
通过使用能够实现高性能感测的宏观测量,可以在逆模型计算的基础上实现高级测量。在这种情况下,基于实际测量并通过高空间分辨率的感测获得的参数用作逆模型的参数。
在根据上述本技术的实施例的信息处理装置中,逆模型计算部件可以使用基于第二测量范围的检测数据的模型参数来按第一空间分辨率的单位确定计算结果。
逆模型计算部件确定例如特征或环境响应的形式的计算结果。在这种情况下,使用与微观测量部件相关的第二测量范围的模型参数,以与宏观测量部件相关的第一空间分辨率的单位确定计算结果。
在根据上述本技术的实施例的信息处理装置中,逆模型计算部件可以按第一空间分辨率的单位确定作为计算结果的测量对象的特征或环境响应。
特征是指测量对象的静态形状或特性。环境响应是指测量对象的动态形状或特性。
在根据上述本技术的实施例的信息处理装置中,宏观测量部件可以在距离测量对象比微观测量部件更长的距离处进行感测。
与微观测量部件相比,宏观测量部件在距离测量对象更长距离处在更宽测量范围上执行测量。另一方面,与宏观测量部件相比,微观测量部件在距离测量对象较短的距离处在相对窄的测量范围内执行测量。
在根据上述本技术的实施例的信息处理装置中,逆模型计算部件可以获取进行聚类的每个测量区域的代表性个体的模型参数,为每个聚类测量区域切换模型参数,并将所得到的模型参数应用于逆模型。
通过计算来自微观测量部件的检测数据来确定由聚类产生的每个测量区域的模型参数,逆模型计算部件可以为由聚类产生的各个测量区域使用不同的模型参数。
在根据上述本技术的实施例的信息处理装置中,可以基于指定区域的用户输入进行所述聚类。
例如,农田被分割为例如种植不同作物的区域。例如,诸如农民之类的用户可以输入这样的信息。
在根据上述本技术的实施例的信息处理装置中,可以基于来自宏观测量部件的检测数据或来自微观测量部件的检测数据进行所述聚类。
例如,使用来自宏观测量部件或微观测量部件的检测数据允许彼此区分具有明显不同形状或状态的区域,从而自动执行聚类计算。
在根据上述本技术的实施例的信息处理装置中,可以基于指定区域的用户输入以及来自宏观测量部件的检测数据或来自微观测量部件的检测数据进行所述聚类。
即,执行反映手动输入和自动辨别两者的聚类计算。
在根据上述本技术的实施例的信息处理装置中,模型参数可以包括植物的三维结构、植物高度、平均叶子角度(平均叶子倾角)、植物覆盖率、LAI、叶绿素浓度、土壤的谱特性或太阳叶子比率中的任一个。
请注意,LAI(叶子面积指数)是表示植物叶子数量的叶子面积指数,因此值越大表示叶子越多。
在根据上述本技术实施例的信息处理装置中,微观测量部件可以包括作为微观测量传感器的可见光图像传感器、立体相机、用于激光图像检测和测距的传感器、偏振传感器或ToF(飞行时间)传感器中的任一个。
请注意,用于激光图像检测和测距的传感器称为光雷达(光检测和测距)。
在根据上述本技术的实施例的信息处理装置中,宏观测量部件可以包括作为宏观测量传感器的多谱相机、超谱相机、FTIR(傅里叶变换红外光谱仪)或红外传感器中的任一个。
根据上述本技术的实施例的信息处理装置还可以包括输出部件,其被配置为基于来自逆模型计算部件的计算结果生成和输出输出图像数据。
也就是说,信息处理装置允许将关于来自逆模型计算部件的计算结果的信息成像并呈现给用户。
在根据上述本技术的实施例的信息处理装置中,输出部件可以生成从来自逆模型计算部件的计算结果的颜色映射得到的输出图像数据。
在针对多个区域中的每一个获得来自逆模型计算部件的计算结果的情况下,生成呈现给用户的图像,使得向每个区域分配不同的颜色。
在根据上述本技术的实施例的信息处理装置中,输出部件可以生成根据从来自逆模型计算部件的计算结果的颜色映射得到的图像与另一图像的合成而得到的输出图像数据。
通过例如叠加或盖写来合成将不同颜色分配给每个区域的图像和另一图像。
在根据上述本技术的实施例的信息处理装置中,宏观测量部件可以安装在人造卫星中。
宏观测量部件安装在人造卫星中,以从天空中的远处测量测量对象,如农田。
在根据上述本技术的实施例的信息处理装置中,微观测量部件可以安装在能够被无线电控制或自动控制的飞行体中。
能够无线电控制或自动控制的飞行体的示例包括所谓的无人机、小型无线电控制固定翼飞机、小型无线电控制直升机等。
根据本技术的另一实施例的信息处理方法包括:由信息处理装置执行计算来自宏观测量部件的检测数据的宏观测量分析处理,宏观测量部件被配置成以第一空间分辨率对测量对象的第一测量范围进行感测;由信息处理装置执行计算来自微观测量部件的检测数据的微观测量分析处理,微观测量部件被配置成以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及由信息处理装置执行基于在微观测量分析处理中确定的来自微观测量部件的检测数据,获取用于使用在宏观测量分析处理中的计算结果的逆模型计算的模型参数的逆模型计算处理。
因此,信息处理装置能够针对测量对象生成由宏观测量和微观测量的组合产生的先进且准确的测量结果信息。
根据本技术的另一实施例的程序是使信息处理装置执行上述方法的程序。这有助于实现产生高级测量结果的计算机装置。
根据本技术的另一实施例的感测系统包括:宏观测量部件,被配置为以第一空间分辨率对测量对象的第一测量范围进行感测;微观测量部件,被配置为以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及上述信息处理装置。
因此,可以构建执行宏观测量和微观测量并进一步使用这些测量的结果来生成测量结果的系统。
附图说明
图1是根据本技术实施例的感测系统中的宏观测量部件和微观测量部件的说明图。
图2是根据实施例的农田遥感示例的说明图。
图3是根据实施例的宏观测量部件和微观测量部件的测量的说明图。
图4示出了根据本实施例的宏观测量部件和微观测量部件的测量范围和分辨率的图。
图5示出了逆计算导致错误结果的情况的图。
图6是根据实施例的聚类的描述图。
图7是根据实施例的信息处理装置的硬件配置的框图。
图8是根据实施例的信息处理装置的功能配置的框图。
图9是根据实施例的处理示例的流程图。
图10是根据实施例的微观测量分析计算处理的流程图。
图11是根据实施例的用于微观测量分析计算的图像的说明图。
图12是根据实施例的聚类的流程图。
图13示出了根据实施例的聚类的图。
图14是根据实施例的逆模型计算的流程图。
图15示出了用于根据实施例的逆模型计算的宏观测量和微观测量部件模型参数的图。
图16示出了根据实施例的使用颜色映射获得的输出图像的图。
图17示出了根据实施例的颜色映射图像和另一图像的合成的图。
图18示出了根据实施例的颜色映射图像和另一图像的合成的图。
具体实施方式
将按以下顺序描述实施例。
<1.感测系统的配置>
<2.信息处理装置的配置>
<3.处理示例>
<4.各种示例>
<5.结论和修改示例>
<1.感测系统的配置>
首先,将描述根据实施例的感测系统。
图1示出了包括在感测系统中的宏观测量部件2和微观测量部件3。
微观测量部件3在相对靠近测量对象4的位置处执行感测。执行一个感测单位的测量范围是被表示为微观测量范围RZ3的相对窄的范围。注意,这里使用的一个单位例如在相机的情况下指的是捕获一帧图像的范围,但是该单位可以取决于传感器类型而变化。
相对照,宏观测量部件2在比微观测量部件3更远离测量对象4的位置处执行感测。执行一个感测单位的测量范围被表示为比微观测量范围RZ3宽的宏观测量范围RZ2。然而,由宏观测量部件2执行一个感测单位的测量范围可以与微观测量范围RZ3相同。
在本实施例中,微观测量范围RZ3与宏观测量范围RZ2相同或更窄。也就是说,测量对象4中的微观测量范围RZ3的区域也被宏观测量范围RZ2覆盖。换句话说,微观测量范围RZ3是在其内执行微观测量部件3的微观测量和宏观测量部件2的宏观测量两者的范围。
如上所述使用宏观测量部件2和微观测量部件3的感测系统的示例是感测例如农田300的植被状态的系统,如图2所示。
图2示出了农田300的状态。最近,已经做出了很多努力来使用安装在小型飞行体200(例如,如图2所示的无人机)中的成像装置250遥感植被状态。
飞行体200可通过例如由操作员无线电控制或自动控制而在农田300上方飞行。
成像装置250被设置在飞行体200中以成像例如飞行体200下方的区域。当飞行体200沿着预定路线在农田300上方飞行时,成像装置250例如定期地捕获静止图像。
如上所述安装在飞行体200中的成像装置250用作图1中的微观测量部件3。由成像装置250捕获的图像被用作微观测量的检测数据。成像装置250的成像范围对应于微观测量范围RZ3。
此外,图2示出了位于天空中的人造卫星210。成像装置220安装在人造卫星210中,并且能够感测地球的表面侧。
成像装置220可以感测(成像)农田300。即,成像装置220用作宏观测量部件2。由成像装置220捕获的图像被用作宏观测量的检测数据。成像装置220的成像范围对应于宏观测量范围RZ2。
用作安装在飞行体200中的微观测量部件3的成像装置250(换句话说,特定微观测量传感器)被假定为可见光图像传感器(以R(红色)、G(绿色)和B(蓝色)成像可见光的图像传感器)、立体相机、光雷达(用于激光图像检测和测距的传感器)、偏振传感器、ToF传感器、NIR(近红外)成像相机等。
此外,微观测量传感器可以是在多个波段中执行成像的多谱相机,例如捕获NIR图像和R(红色)图像,并且能够基于获得的图像计算NDVI(标准化差异植被指数),只要传感器具有能够可操作地将传感器安装在飞行体200中的设备尺寸即可。NDVI是表示植被分布状况和活动性的指标。
上述传感器期望适于分析例如测量对象的特征、环境响应、环境状态(范围、分布等)等。请注意,特征是测量对象的静态形状和特性。环境响应是测量对象的动态形状和特性。环境状态是测量对象所在环境的状态,包括测量对象所在范围、测量对象分布或环境特性。
此外,所述传感器期望尺寸相对小且重量轻,并且容易安装在飞行体200中。
另一方面,用作安装在人造卫星210中的宏观测量部件2的成像装置220(换句话说,特定宏观测量传感器)可以是在多个波段中捕获图像(例如,NIR图像和R图像)的多谱相机、超谱相机、FTIR(傅里叶变换红外光谱)、红外传感器等。在这种情况下,接受相对大型的感测设备,并假设能够进行精确的感测。
这些宏观测量传感器是适于分析各种物理属性值(例如,关于光合作用等的信息)的传感器。
此外,由于设备尺寸、重量等原因,这些传感器难以安装在小型飞行体200中。然而,在本示例中的感测系统中,这样的传感器安装在人造卫星210中。
另外,标签信息被附加到由成像装置220和250捕获和获得的图像。标签信息包括成像日期和时间信息、以GPS(全球定位系统)数据的形式的位置信息(纬度/经度信息)、成像装置信息(关于相机的个体标识信息和型号信息等)、关于每个图像数据的信息(例如图像大小、波长和成像参数等信息)等。
注意,位置信息以及成像日期和时间信息也被用作将来自成像装置220的图像(检测数据)与来自成像装置250的图像(检测数据)相关联的信息。
如上所述,由安装在飞行体200中的成像装置250和安装在人造卫星210中的成像装置220获得的图像数据和标签信息被发送到信息处理装置1。信息处理装置1使用图像数据和标签信息来生成使用农田300作为测量对象的分析信息。另外,信息处理装置1执行将分析结果作为图像呈现给用户的处理。
信息处理装置1被实现为例如PC(个人计算机)、FPGA(现场可编程门阵列)、诸如智能电话或平板电脑之类的终端装置等。
注意,图1示出了与成像装置250分离的信息处理装置1,但是,例如,可以在包括成像装置250的单元中设置用作信息处理装置1的计算装置(微型计算机等)。
参考图3,将描述宏观测量部件2和微观测量部件3的作用。
宏观测量部件2在测量范围RZ2内执行测量,并使用模型执行逆计算,以获得植被等的输出或测量结果。例如,通过使用准备好的模型来解释宏观测量值,可以从测量值的混合中获得正确的输出。
然而,如果用于逆模型计算的模型与实际测量对象之间存在差异,则无法获得准确的输出。特别地,植物和动物因生长而改变形状,或因品种而在形状方面显著不同,因此该模型常常不合适。
因此,代表性个体的测量值被用于更新模型。为此,使用来自微观测量部件的检测数据。
微观测量部件3可以测量每个待测量的个体。例如,图示了个体OBa1、OBa2、OBa3…OBan,并且微观测量部件3可以测量或确定每个个体的特征、环境响应和环境状态、基于特征、环境响应和环境状态的区域标识等。
来自微观测量部件3的检测数据被用于测量代表性个体的特征或环境响应,并且特征或环境响应被用于改变模型参数。然后,可以使用基于实际测量的模型进行逆模型计算。
使用感测系统的测量的主要目的是宏观趋势(平均、合计、分布等),其可通过使用基于来自宏观测量部件2的检测数据的信息的模型逆计算来获得。在这种情况下,通过基于来自微观测量部件3的检测数据更新模型参数来确保测量精度。
图4A二维地示出了宏观测量范围RZ2和微观测量范围RZ3,并且还示出了宏观测量部件2和微观测量部件3的分辨率。
大单元格均表示宏观测量分辨率,而小单元格均表示微观测量分辨率。使用分辨率获得的信息对应于例如捕获图像的一个像素的信息。
即,安装在宏观测量部件2中的宏观测量传感器是具有与大单元格对应的分辨率的传感器,而安装在微观测量部件3中的微观测量传感器是具有与小单元格对应的分辨率的传感器。
例如,给定如虚线所示的测量对象,可以使用与由细线所示的小单元格对应的分辨率来确定测量对象的特征、环境响应、区域等,可以使用与由粗线所示的大单元格对应的分辨率来测量物理属性值等。
将描述需要这种感测系统的背景。
如上所述,最近频繁地执行使用诸如无人机之类的飞行体200的感测。这除了允许基于可见光(RGB)的测量来测量特征之外,还允许使用各种光波长和技术来测量对象的物理属性、生理状态等。然而,可安装在小飞行体200中的感测设备常常在尺寸、重量等方面受到限制。
获取大量波长并且能够执行分量分析等的超谱相机可能通常需要扫描机构以获取二维图像,并且尺寸较大。因此,除非飞行体很大,否则超谱相机很难安装。
此外,扫描可能需要一定的时间,还可能需要悬停,从而延长测量时间。扫描还常常妨碍利用飞行体200中的电池容量来测量大陆地。
此外,在扫描期间,飞行体200的可能振动可能降低测量精度。
此外,使用更高谱分辨率的FTIR方案原则上涉及长度较大且难以安装在飞行体200中的装备。
为了精确感测,可以通过安装大型成像器或提供多次曝光来提高S/N(信噪比)。然而,大型成像器包括大型光学系统,并且不适合安装在飞行体200中。由于飞行体200的悬停,多次曝光导致测量时间延长,并且飞行体200的可能振动降低了精度。
另外,通常,飞行体200的壳体暴露于阳光下并且具有高于正常温度的温度。
在精确感测中,传感器保持低温,以使得能够降低热噪声。室内使用的一些传感器(例如分光光度计)通过使用珀耳帖元件等保持低温来保持精确。然而,珀耳帖元件消耗大量电力,因此不适合安装在电力有限的飞行体200中。
使用如空调中的压缩机的热泵型温度调节设备具有高功率效率,但在尺寸或重量方面不适合于安装在飞行体200中。
另一方面,在卫星感测中,安装了使得能够进行先进感测的装备的人造卫星正在运行。然而,卫星感测在空间分辨率方面是不够的。
对于人造卫星210来说,安装超谱相机、FTIR或大型成像器以及如上所述的低温控制等并不那么困难。
然而,由于在一个空间分辨率单位中混合了各种对象,低空间分辨率不仅排除了形状的确定,而且还妨碍了仅对要检查的对象的测量。特别地,在植被测量的示例中,不利地混合了土壤、阴影等。
更具体地说,在使用人造卫星210的遥感中,从远位置执行测量,因此,除了要测量的物体之外,在与测量仪器的空间分辨率相对应的一个像素中混合不同的物体。测量值根据像素中每个物体的比率进行加权,并被称为混素(mixel)(混合像素)。例如,当要测量植被时,测量植物和土壤的混合物,妨碍了确定植物本身的物理属性值(叶绿素浓度等)。
为了解决这一问题,提出了各种方法。
“SAVI(土壤调整植被指数)”是一种植被指数,用于校正土壤反射率引起的波动。当LAI表示为“L”时,L=0(等于NDVI)用于高LAI,L=1用于低LAI。
SAVI=〔(NIR-RED)/(NIR+RED+L)〕×(1+L)
然而,可能需要使用假设值作为LAI的值,在LAI可能随位置变化的情况下(如在农田中)无法进行精确校正。
此外,影响遥感的不利影响不限于物体的二维比率。来自植物叶子的反射不是朗伯(Lambertian)反射,并且受光源入射角、叶子角度依赖性和测量仪器角度的影响(这种影响通常称为BRDF)。
植物的三维形状导致阴影、群落中光束的多次反射等,从而使测量值波动。
为了处理这样的复杂机制,使用描述反射率的角度依赖性的叶子辐射传输模型(示例:PROSAIL模型)或包含阴影的群落的光学模型来模拟光学行为。
即,为了覆盖卫星遥感的低空间分辨率并且应对影响测量的复杂因素,使用包括关于测量对象形式的信息的“模型(辐射传输特性模型等)”的逆计算被用于确定特定对象的测量值。
然而,在测量对象的延伸和模型的形状之间没有差异的情况下(如在植被的情况下在热带雨林中),这是合适的,但例如,在农田300的测量(勘察)中,形状本身要被测量,且不能被标识或正确测量。例如,特征等是未知的,因为,例如,作物在生长过程中形状发生变化,或者由于某种原因,作物无法适当地生长或即将死亡。
例如,图5A示出了具有基本相同形状的植物个体被广泛地排列成线,并且例如,LAI在宏观测量范围RZ2内基本上是一致的。
另一方面,图5B示出了其中种植具有大叶子(长得很大)的植物的区域H1和其中种植具有小叶子的植物的区域H2。图示的宏观测量范围RZ2跨越区域H1和H2。在这种情况下,LAI或植物高度随位置而变化。因此,宏观测量范围RZ2和用于逆模型计算的模型之间可能存在较大差异。
因此,在根据本实施例的感测系统中,由微观测量部件3进行的高分辨率测量包括对实际田地中的植物状态(LAI、植物覆盖度和高度、叶子平均角度、叶绿素浓度等)、垄(ridge)的安装状态等的测量,并利用实际物体的参数进行卫星感测的逆模型计算。这使得即使在测量对象具有变化的形状或状态或具有与由模型指示的标准形状不同的形状的情况下也可以进行测量。
在特定示例中,基于来自微观测量部件3的检测数据,使用真实物体的参数执行逆模型计算。因此,即使在测量对象的形状和状态未知的情况下,也计算正确的SIF(太阳诱导叶绿素荧光),以获得有关光合作用速率的信息。
模型参数假定为植物的三维结构、植物高度、平均叶子角度、植物覆盖度、LAI、叶绿素浓度、土壤的谱特性、太阳叶子比率等。
此外,根据本实施例的感测系统还包括聚类的思想。换句话说,感测系统包括将测量对象区域划分为聚类、测量每个聚类的代表性个体以及逐聚类地切换模型的聚类思想。
图6示意性地示出了在特定区域H10和H20中的每一个中执行感测的情况。
假设个体OBa1,OBa2,OBa3…OBan存在于H10区域中,且具有不同于个体OBa1、OBa2、OBa3…OBan的特征的特征的个体OBb1、OBb2、OBb3…OBbn存在于区域H20中。
聚类处理用于将区域H10和H20识别为不同的聚类。
此外,对于每个聚类,测量代表性个体以确定模型参数。
不同的模型参数用于各自的聚类。
即,在使用宏观测量部件2和微观测量部件3测量区域H10的情况下,使用基于相关聚类的代表性个体的模型参数执行逆模型计算。
在使用宏观测量部件2和微观测量部件3测量区域H20的情况下,使用基于相关聚类的代表性个体的模型参数执行逆模型计算。
以这种方式,对从聚类产生的各个区域使用不同的模型,以获得准确的测量结果。
另外,在图4A和图4B中,用虚线示出了测量对象个体,但是在该图中,微观测量部件3覆盖了区域H1和H2两者。假设区域H1和H2包括例如如图5B所示的具有不同叶子大小(生长方式)的群体,作为不同植物等。
同样在这种情况下,以下配置是可能的:执行逆模型计算,使得对于宏观测量范围RZ2,基于区域H1的微观测量的模型被用于与区域H1相关的部分,并且基于区域H2的微观测量的模型被用于与区域H2相关的部分。图4C示意性地示出了用于逆模型计算的模型的切换。
<2.信息处理装置的配置>
将描述上述感测系统中的信息处理装置1,其从宏观测量部件2和微观测量部件3获取检测信息并执行诸如分析之类的处理。
图7示出了信息处理装置1的硬件配置。信息处理装置1包括CPU(中央处理器)51、ROM(只读存储器)52和RAM(随机存取存储器)53。
CPU 51根据存储在ROM 52中的程序或从存储部件59加载到RAM 53中的程序执行各种类型的处理。RAM 53还适当地存储用于CPU 51执行各种类型的处理的数据等。
CPU 51、ROM 52和RAM 53经由总线54连接在一起。输入/输出接口55也连接到总线54。
输入/输出接口55可以连接到包括液晶面板、有机EL(电致发光)面板等的显示部件56、包括键盘、鼠标等的输入部件57、扬声器58、存储部件59、通信部件60等。
显示部件56可以与信息处理装置1集成,或者可以是与信息处理装置1分离的装备。
显示部件56基于来自CPU 51的指令在显示屏上显示各种分析结果等。另外,基于来自CPU 51的指令,显示部件56显示各种操作菜单、图标、消息等,即,作为GUI(图形用户界面)操作以提供显示。
输入部件57表示使用信息处理装置1的用户所使用的输入设备。
例如,输入部件57被假定为诸如键盘、鼠标、键、拨盘、触摸面板、触摸板或遥控器之类的各种操作器和操作设备中的任何一个。
输入部件57检测用户的操作,并且由CPU 51解释与输入操作相对应的信号。
存储部件59包括存储介质,例如HDD(硬盘驱动器)、固态存储器等。存储部件59存储例如从宏观测量部件2和微观测量部件3接收的检测数据和分析结果以及各种其他类型的信息。另外,存储部件59用于存储用于分析处理等的程序数据。
通信部件60经由包括因特网的网络执行通信处理,并与外围部分中的装备通信。
例如,通信部件60可以是与微观测量部件3(成像装置250)和宏观测量部件2(成像装置220)通信的通信设备。
驱动器61也根据需要连接到输入/输出接口55,并且诸如存储器卡之类的存储设备6根据需要安装在输入/输出接口55中,以允许向存储设备6写入数据和从存储设备6读取数据。
例如,从存储设备6读取的计算机程序安装在存储部件59中,并且根据需要,由CPU51处理的数据存储在存储部件中。不用说,驱动器61可以是用于诸如磁盘、光盘或磁光盘之类的可移除存储介质的记录和再现驱动器。磁盘、光盘、磁光盘等也是存储设备6的模式。
注意,根据本实施例的信息处理装置1不限于具有如图7所示的硬件配置的单一信息处理装置(计算机装置)1,并且可以将多个计算机装置配置成系统。多个计算机装置可以被配置成使用LAN等的系统,或者通过利用因特网等的VPN(虚拟专用网络)等被布置在远程位置。该多个计算机装置可以包括通过云计算服务可用的计算机装置。
另外,图7中的信息处理装置1可以实现为诸如台式或笔记本类型之类的个人计算机,或者诸如平板终端或智能电话之类的便携式终端。此外,诸如测量装置、电视装置、监视器装置、成像装置或用作信息处理装置1的设施管理装置之类的电子装备可以配备根据本实施例的信息处理装置1。
例如,具有这种硬件配置的信息处理装置1包括CPU 51的计算功能、ROM 52、RAM53和存储部件59的存储功能、通信部件60和驱动器61的数据获取功能以及显示部件56的输出功能等。已安装的软件功能提供如图8所示的功能配置。
即,作为主要部件,信息处理装置1设置有图8所示的数据输入部件10、分析执行部件20以及数据保存和输出部件30。
这些处理功能由CPU 51启动的软件实现。
包括在软件中的程序从网络下载或从存储设备6(例如,可移除存储介质)读取,并安装在图7中的信息处理装置1中。或者,该程序可以预先存储在存储部件59等中。然后,CPU51启动程序以激活如上所述的每个部件的功能。
另外,使用例如RAM 53中的存储区域或存储部件59中的存储区域来实现各种缓冲器等的存储功能。
另外,图8示出了作为信息处理装置1的外部装置的宏观测量部件2、微观测量部件3、网络5、存储设备6和操作输入部件7。
如上所述,宏观测量部件2例如安装在人造卫星210中。宏观测量传感器2S是大型传感器,例如超谱相机或FTIR,并且可以安装在人造卫星210中,但是难以安装在飞行体200中。这些传感器通常在可见光谱之外工作,主要用于测量物理属性。
微观测量部件3安装在飞行体200中。微观测量传感器3S是诸如RGB相机或立体相机之类的小型传感器,并且易于安装在飞行体200中。一般地,微观测量传感器3S工作在可见光谱中,主要用于测量测量对象的特征和环境响应。
网络5被假定为例如因特网、家庭网络、LAN(局域网)、卫星通信网络或各种其他网络中的任何一个。
如上所述,存储设备6主要被假定为可移除存储介质,例如存储器卡或盘状记录介质。
操作输入部件7是用户可以通过其提供手动输入的设备,并且可以被认为是图7中的输入部件57。
操作输入部件7可以与信息处理装置1集成,或者可以是与信息处理装置1分离的装置。
信息处理装置1中的数据输入部件10对应于接受来自上述外部装置的数据输入的功能,并且包括传感器输入部件11和12以及程序和模型输入部件13。
传感器输入部件11输入来自宏观测量部件2的宏观测量传感器2S的检测信息。例如,通过宏观测量部件2和图7中的通信部件60之间的通信,可以直接接收来自宏观测量传感器2S的检测数据。
或者,可以通过通信部件60经由网络5接收来自宏观测量传感器2S的检测数据。
此外,可以经由存储设备6获取来自宏观测量传感器2S的检测数据。
传感器输入部件12输入来自微观测量部件3的微观测量传感器3S的检测信息。例如,可以通过微观测量部件3和通信部件60之间的通信直接接收来自微观测量传感器3S的检测数据,通过通信部件60经由网络5接收、经由存储设备6获取来自微观测量传感器3S的检测数据等。
注意,诸如光源的光谱校正之类的处理可以在传感器输入部件11和12中进行。
程序和模型输入部件13通过经由网络5从服务器下载程序或模型或从存储设备6读取程序或模型来获取合适的程序或模型。模型是指用于逆计算的模型参数。
分析执行部件20包括宏观测量分析计算部件21、宏观测量分析值缓冲器22、微观测量分析计算部件23、微观测量分析值缓冲器24、位置映射部件25、逆模型计算程序和模型保持部件26(以下称为“保持部件26”)、逆模型计算部件27和聚类计算部件28。
宏观测量分析计算部件21根据由传感器输入部件11获取的来自宏观测量传感器2S的检测数据,执行用于确定物质分量的量等的计算。
例如,宏观测量分析计算部件21根据来自超谱相机或FTIR的多波长数据,基于NIRS(近红外光谱)和FLD方法(夫琅和费线辨别)方法(太阳黑线:夫琅和费线)等,计算植被指数和SIF(叶绿素荧光)。
请注意,太阳暗线O2A的波长宽度约为1nm且较薄,因此使用诸如超谱相机或FTIR之类的传感器的该感测通常是适合的。在本示例中的感测系统中,由于此类装备易于安装在人造卫星210中,因此宏观测量分析计算部件21基于来自宏观测量部件2的检测数据执行计算。
宏观测量分析值缓冲器22临时保存由宏观测量分析计算部件21处理的数据。
例如,宏观测量分析值缓冲器22保存由宏观测量分析计算部件21计算的SIF、从宏观测量部件2通知的位置信息等。
微观测量分析计算部件23执行用于从由传感器输入部件12获取的来自微观测量传感器3S的检测数据中获得适当信息的计算。
例如,微观测量分析计算部件23计算LAI、平均叶子角度、太阳叶子比率等。太阳叶子是指暴露在阳光下的叶子或叶子的一部分。
微观测量分析值缓冲器24临时保存由微观测量分析计算部件23处理的数据。
例如,微观测量分析值缓冲器24保存由微观测量分析计算部件23确定的关于LAI、平均叶子角度、太阳叶子比率等的信息,以及从微观测量部件3通知的位置信息,并且进一步保存RGB图像、NDVI图像等。
位置映射部件25执行用于从具有不同分辨率或成像单位(测量范围RZ2和RZ3)的一组图像中提取相同点的计算。例如,GPS信息或正交镶嵌(orthomosaic)处理用于将由宏观测量分析计算部件21处理的信息与由微观测量分析计算部件23处理的信息对准。
逆模型计算部件27是使用来自宏观测量分析计算部件21和微观测量分析计算部件23的计算结果执行逆模型计算的功能。
即,逆模型计算部件27基于由微观测量分析计算部件23确定的来自微观测量部件3的检测数据生成模型参数。
然后,逆模型计算部件27使用模型参数作为逆模型的参数,以使用来自宏观测量分析计算部件21的计算结果执行逆模型计算。
注意,在基于来自微观测量部件3的检测数据生成模型参数之后,逆模型计算部件27可以将模型参数与关于日期、时间和位置的信息相关联,并且使得结果信息保存在预定存储部件中。日期、时间和位置指的是执行用于确定模型参数的原始微观测量的日期、时间和位置(例如,GPS信息)。
保持部件26保持由程序和模型输入部件13获取的逆模型计算程序和模型默认参数。逆模型计算部件27基于这些程序和模型执行逆模型计算。
聚类计算部件28执行聚类计算。例如,聚类计算部件28基于通过操作输入部件7的用户输入,执行与将要测量的农田300等的到区域的分割相对应的聚类。例如,用户指定田间边界,跨越边界种植不同作物或处于不同发育阶段的相同作物。这允许用户执行可选的聚类分割。
另外,聚类计算部件28可以经由网络5或经由存储设备6获取关于用于聚类的指示的信息。
此外,聚类计算部件28可以基于从微观测量分析计算部件23获得的信息或从宏观测量分析计算部件21获得的信息来执行自动聚类。
从微观测量分析计算部件23获得并用于聚类的信息可以是来自微观测量部件3的检测数据,或者由微观测量分析计算部件23根据检测数据计算的LAI、平均叶子角度、太阳叶子比率等。
从宏观测量分析计算部件21获得并用于聚类的信息可以是来自宏观测量部件2的检测数据(多谱光谱信息等),或者由宏观测量分析计算部件21根据检测数据计算的SIF等。
可以基于上述不同类型信息之间的差异来进行聚类分割。此外,假设执行图像纹理的分析、机器学习等。
此外,在生理特征或环境响应以及形式或特征方面存在任何差别的情况下,该差别可用于聚类。
另外,在手动输入和自动聚类确定的组合中,在用户指定区域并且用户指定的一个区域包括由宏观测量部件自动辨别的多个聚类的情况下,可以设定自动子分割的聚类。
来自聚类计算部件28的聚类信息被提供给逆模型计算部件27,并在逆模型计算期间被参考。
数据保存和输出部件30包括分析数据缓冲器31、颜色映射部件32、图像合成部件33、图形生成部件34、图像输出部件35和数据输出部件36。
关于来自逆模型计算部件27的计算结果的信息临时保存在分析数据缓冲器31中。
在逆模型计算部件27例如确定SIF的情况下,分析数据缓冲器31保持该信息。另外,分析数据缓冲器31可以保持RGB图像或NDVI图像。
为了可视化并显示作为来自逆模型计算部件27的计算结果而获得的物理值,颜色映射部件32执行例如使用每一级RGB原色来将特定范围的物理值转换为从蓝色到红色的颜色等级(gradation)的计算处理。
图像合成部件33执行将颜色映射的物理值数据粘贴到原始空间区域或将颜色映射的物理值数据叠加到RGB图像上的计算处理。
为了可视化和显示数据,图形生成部件34执行创建图形的计算处理,例如,以虚线表示物理值或将二维物理值转换为散布图。
图像输出部件35将通过颜色映射部件32、图像合成部件33和图形生成部件34的处理而生成的图像数据输出到外部显示部件56以进行显示。或者,图像输出部件35执行输出,以使用网络5将生成的图像数据发送到外部装置,或者将图像数据归档并将结果文件保存在存储设备6中。
数据输出部件36输出关于保存在分析数据缓冲器31中的来自逆模型计算部件27的计算结果的信息。例如,数据输出部件36执行输出,用于使用网络5将关于逆模型计算结果的信息发送到外部装置,或者将关于逆模型计算结果的信息归档并将结果文件保存在存储设备6中。
<3.处理示例>
将描述包括上述功能的信息处理装置1的处理示例。
图9示出了信息处理装置1的处理示例。
在步骤S101中,信息处理装置1使用传感器输入部件11的功能输入来自宏观测量部件2的测量值。
在步骤S102中,信息处理装置1使用宏观测量分析计算部件21的功能执行宏观测量分析计算。例如,对于有关光合作用的信息,执行SIF计算。对于SIF计算,已知基于太阳光谱中的暗线的FLD方法。
在步骤S103中,信息处理装置1使用传感器输入部件12的功能输入来自微观测量部件3的测量值。
在步骤S104中,信息处理装置1使用微观测量分析计算部件23的功能执行微观测量分析计算。
图10示出了步骤S104中微观测量分析计算中的处理示例。
注意,假设微观测量分析计算部件23已经获取了图11中所示的RGB图像、NIR图像、R图像和偏振传感器角度信息。
在图10的步骤S201中,微观测量分析计算部件23执行图像结构分析计算。具体地说,图像被分割。在简单情况下,如图11所示,可以将图像分割为单元格。
然而,可以通过使用图像识别等来根据对象等进行辨别来分割图像。例如,标识作为测量对象的个体,并根据标识出的个体对图像进行分割。
在步骤S202中,微观测量分析计算部件23基于分割单位计算LAI、平均叶子角度和太阳叶子比率。
LAI可根据植物覆盖率确定。植物覆盖率可通过将与特定值或更大的NDVI相对应的像素数除以相关分割单位中的测量点数(像素数)来确定。
注意,NDVI可以根据R图像和NIR图像确定。也就是说,NDVI的值由以下公式确定:
NDVI=(NIR-R)/(NIR+R)
其中,“R”是可见区域红色的反射率,“NIR”是近红外区域中的反射率。NDVI具有被标准化在“-1”到“1”的范围内的数值,NDVI的正值越大,表明植被越密集。
叶子角度可以通过从偏振传感器获取角度信息来确定。
太阳叶子比率是暴露在阳光下的叶子的比率。可以将太阳叶子提取为具有特定值或更大的NDVI和特定值或更大的NIR的像素。因此,可以通过将这样提取的像素的数量除以相关分割单位中的测量点的数量(像素的数量)来确定太阳叶子比率。
在步骤S203中,微观测量分析计算部件23确定是否针对所有分割单位完成LAI、平均叶子角度和太阳叶子比率的计算。在计算未完成的情况下,则在步骤S202中,为下一分割单位计算LAI、平均叶子角度和太阳叶子比率。
一旦对所有分割单位完成处理,微观测量分析计算就结束。
在图10中的上述处理中,图9中步骤S104的执行使得信息处理装置1在步骤S105中使用聚类计算部件28的功能执行聚类。
图12示出了聚类计算的示例。
在步骤S301中,聚类计算部件28为测量对象指定由用户手动输入的分割。
图13A示意性地示出了与测量对象相对应的农田300。请注意,区域AR1至AR6的定义是为了便于描述,并不一定对应于不同类型植被的区域。但是,假设阴影区域AR3是种植与其他区域中的作物不同的作物的区域。
在这种情况下,用户通过手动输入指定用粗线所示的边界。例如,假设指示了农田中种植植物的整个区域(田地边界),并指定了区分种植不同类型作物的区域AR3的边界。
聚类计算部件28执行粗线指示的聚类分割,反映基于用户初步知道的此类信息的输入。
在步骤S302中,聚类计算部件28使用从宏观测量分析计算部件21获得的信息和从微观测量分析计算部件23获得的信息执行自动聚类。使用例如SIF量、LAI、平均叶子角度、太阳叶子比率等来执行聚类。
图13B示出了作为与多个测量相关的微观测量范围RZ3的测量范围a、b、c和d。在这种情况下,测量范围a对应于用于测量区域AR3的微观测量范围RZ3,测量范围b对应于用于测量区域AR4的微观测量范围RZ3,测量范围c对应于用于测量区域AR5的微观测量范围RZ3,并且测量范围d对应于用于测量区域AR6的微观测量范围RZ3。
假设自动聚类涉及执行例如将区域分割为具有不同LAI的聚类的处理。假设LAI值在测量范围a、b和c之间变化,但在测量范围c和d之间基本相同。
区域AR3与区域AR4的作物和LAI不同。假设区域AR4、AR5和AR6的作物相同,但只有区域AR4涉及不同生长情况。那么,将区域AR4设定为单独的聚类是合适的。
在步骤S303中,聚类计算部件28基于用户输入对聚类设定进行子分割,使得该子分割反映步骤S304中的自动聚类确定。
在上述示例中,如图13C中的粗线所示,区域AR4被新设定为一个聚类。
例如,区域AR3和区域AR4分别被假定为聚类CL3和聚类CL4。注意,未提及区域AR1和AR2,但在区域AR1和AR2具有与区域AR5和AR6相同的LAI的情况下,区域AR1、AR2、AR5和AR6被共同指定为一个聚类CL1。
例如,如上所述,通过根据用户输入执行聚类并使用与宏观测量或微观测量相关的值执行用于子分割的自动聚类确定来执行聚类计算。
不用说,可以排他地根据用户输入或者通过使用与宏观测量或微观测量相关的值执行用于子分割的自动聚类确定来设定聚类。
在执行上述聚类计算之后,在图9的步骤S106中,信息处理装置1使用位置映射部件25的功能执行位置映射。
即,信息处理装置1将宏观测量分析计算结果与微观测量分析计算结果对准。
在步骤S107中,信息处理装置1使用逆模型计算部件27的功能执行逆模型计算。图14示出了逆模型计算的处理示例。
注意,图15A示意性地示出了基于宏观测量的SIF计算的区域。SIF以单元格为单位来确定,每个单元格被图示为宏观测量分辨率(宏观分辨率单位W1至Wn)。
在步骤S401中,逆模型计算部件27针对一个宏观分辨率单位读取由宏观测量分析计算部件21计算的SIF。例如,逆模型计算部件27首先读取宏观分辨率单位W1的SIF。
在步骤S402中,逆模型计算部件27针对与宏观分辨率单位相对应的聚类获取由微观测量分析计算部件23确定的参数,即LAI、平均叶子角度和太阳叶子比率。
图15B示出了上述测量范围a、b和c(=d)的LAI、平均叶子角度和太阳叶子比率。换句话说,如图13C所示,LAI、平均叶子角度和太阳叶子比率是区域AR3的聚类CL3、区域AR4的聚类CL4和区域AR1+AR2+AR5+A6的聚类CL1的模型参数。
例如,宏观分辨率单位W1对应于聚类CL1,因此获得图15B中的聚类CL1的模型参数。
在步骤S403中,逆模型计算部件27执行逆模型计算。也就是说,根据基于宏观测量获得的SIF确定期望的物理属性值(例如,测量对象的特征)。
在这种情况下,使用基于在步骤S402中获取的微观测量的参数作为模型参数。
因此,即使在仅仅通过由宏观测量部件2执行的感测不知道测量对象的形状或状态的情况下,也可以根据实际测量对象的形状和状态来确定正确的SIF。
通过从步骤S404返回到步骤S401来重复上述处理,直到对所有宏观分辨率单位W1到Wn执行处理。换句话说,从宏观分辨率单位W1到宏观分辨率单位Wn,使用基于对应聚类的微观测量的模型参数执行逆模型计算。
注意,已经描述了SIF的示例,但例如,在基于宏观测量获得NDVI的情况下,所需的物理属性值(例如,作为环境响应的叶绿素浓度)可通过根据基于宏观测量的NDVI的逆计算来确定。
一旦对所有宏观分辨率单位完成了处理,逆模型计算部件27进入步骤S405,以将计算结果写入分析数据缓冲器31中。在这种情况下,为每个宏观分辨率单位W1到Wn写入计算结果。
注意,在上述处理中,对应于每个宏观分辨率单位的聚类的模型参数被应用于宏观分辨率单位,但是可以不特别地执行聚类。例如,在相同植物在整个农田300中以类似生长条件种植的情况下,步骤S105中的聚类可能不导致聚类分割。另外,假设不执行步骤S105中的处理的处理示例。
在这些情况下,通过将通过测量单个微观测量范围RZ3而获得的模型参数应用于所有宏观分辨率单位来执行逆模型计算就足够了。
在上述处理中,在执行图9中的步骤S107之后,信息处理装置1使用数据保存和输出部件30的功能来执行步骤S108中的颜色映射、步骤S109中的图像合成和步骤S110中的图像输出。
因此,用户可以使用显示部件56等来检查计算结果。
将描述在这种情况下经受颜色映射等的输出图像的示例。
图16是通过对如上所述获得的每个宏观分辨率单位的逆模型计算结果执行颜色分配(颜色映射)来生成图像的示例。这里使用的“颜色分配”是指预设对应于每个数值范围的颜色,根据对象的值选择颜色,并将颜色分配给相关像素。
图16A示出了每个宏观分辨率单位的逆模型计算结果。如图16B所示,将颜色分配给这些值以生成颜色映射图像。在图像中,与计算结果的值相对应的颜色被分配给每个区域。
注意,附图基于阴影线、圆点等类型表示颜色差异。此外,对于不存在有效计算结果的宏观分辨率单位(例如,不存在太阳叶子的部分等),指示“无数据”。例如,为指示“无数据”的区域分配背景色(白色)。
在向用户呈现这样的颜色映射图像的情况下,每个宏观分辨率单位中的区域的逆模型计算结果用颜色表示,并且该图像允许用户容易地识别宏观测量结果。
接下来,图17是合成具有分配给植被处于特定状态的区域的颜色的图像的示例。图17A示出了每个宏观分辨率单位的逆模型计算结果的值。图17B是描绘提取的太阳叶子的图像NDVIp-pr(NDVI植物过滤Par过滤)。描绘提取的太阳叶子的图像NDVIp-pr是指描绘提取范围的图像,该提取范围具有特定值或更大的NDVI和特定值或更大的NIR。
然后,将颜色分配给每个宏观分辨率单位中的太阳叶子部分,以生成如图17C所示的颜色映射图像。只有太阳叶子部分具有与逆模型计算结果相应的颜色。因此,该图像允许用户容易地识别宏观测量结果以及每个区域中的太阳叶子分布。
接下来,图18是在可见光图像(RGB图像)上的叠加显示的示例。
图18A示出了每个宏观分辨率单位的逆模型计算结果的值。图18B是RGB图像。
如图18C所示,在RGB图像中,叠加根据逆模型计算结果的值而分配给每个宏观分辨率单位的颜色。该图说明了将颜色应用于相关像素部分。
换句话说,在图像中,在RGB图像上表示指示计算结果的颜色。因此,在用户通常观看的图像上容易识别宏观测量结果。
注意,可以执行使用分配给相关像素的颜色的盖写(overwrite),来代替叠加(overlay)。
例如,如上所述生成并在图16、图17和图18中示出了输出图像,并将其显示在显示部件56上,使用网络5发送到外部装置,或者与保存在存储设备6中的结果文件一起被归档。然后用户可以利用该分析结果。
<4.各种示例>
上述描述假设植被感测。然而,根据本公开的技术适用于各种领域。
例如,在诸如办公楼之类的建筑物中使用中央热源的情况下,可以通过宏观测量来识别整个建筑物的能量使用量。
此外,可以对建筑物的某一部分(例如,占据某一楼层的商务办公室)进行特定测量,作为微观测量。
在使用基于宏观测量的信息执行逆模型计算的情况下,可以基于微观测量设置模型参数。
此外,例如,在劳动统计等领域中,某一时段的失业率的转变可以作为宏观测量进行测量,而在每个季节中的失业人员的属性等可以作为微观测量进行测量。
在使用基于宏观测量的信息执行逆模型计算的情况下,可以基于微观测量设置模型参数。
<5.结论和修改示例>
上述实施例产生以下效果。
根据本实施例的信息处理装置1包括宏观测量分析计算部件21,该宏观测量分析计算部件21计算来自宏观测量部件2的检测数据,该宏观测量部件2以宏观测量分辨率(第一空间分辨率)执行针对测量对象的宏观测量范围RZ2(第一测量范围)的感测。信息处理装置1还包括微观测量分析计算部件23,该微观测量分析计算部件23计算来自微观测量部件3的检测数据,该微观测量部件3针对微观测量范围RZ3(第二测量范围)以微观测量分辨率(第二空间分辨率)执行感测,微观测量分辨率高于宏观测量分辨率,微观测量范围RZ3包括在宏观测量范围RZ2中。此外,信息处理装置1包括逆模型计算部件27,该逆模型计算部件27基于由微观测量分析计算部件23确定的来自微观测量部件3的检测数据,获取用于使用来自宏观测量分析计算部件21的计算结果的逆模型计算的模型参数。
通过如上所述基于微观测量计算模型参数,可以获得适用于使用来自宏观测量分析计算部件的计算结果的逆模型计算的逆模型的参数。
另外,在本实施例中,逆模型计算部件27使用基于由微观测量分析计算部件23确定的来自微观测量部件3的检测数据的模型参数,作为用于使用来自宏观测量分析计算部件21的计算结果的逆模型计算中的逆模型的参数。
在执行逆模型计算的情况下,以及在形状或状态变化等情况下,测量对象的形状或状态可能不同于由模型指示的标准形状或状态,导致测量精度降低。在本实施例中,基于基于微观测量的实际测量生成逆模型参数。因此,即使在形状或状态变化或测量对象没有由模型指示的标准形状的情况下,也可以实现精确测量。
特别地,在这种情况下,通过使用来自微观测量部件3(其可以以高空间分辨率执行感测)的检测数据,可以使逆模型参数变得更合适。这使得能够基于与测量对象的特征和环境响应相对应的逆模型计算进行测量,这不能仅通过宏观测量部件2来实现。
更具体地说,在使用飞行体200的高分辨率测量中,测量实际田地中的植物状态(植物覆盖率和高度、叶子的平均角度、叶绿素浓度等)或垄的设定状态,并且,将真实对象的参数用于基于人造卫星210的感测的逆模型计算。然后,即使在形状或状态变化或测量对象没有由模型指示的标准形状的情况下,也可以实现精确测量。
此外,例如,这种测量系统的使用允许飞行体200在不测量SIF的情况下获得光合作用信息。有利的是,人造卫星210也可以实现高分辨率信息的获取,这是迄今为止没有预料到的。
在该实施例中,逆模型计算部件27使用基于微观测量范围RZ3(第二测量范围)的检测数据的模型参数来确定宏观测量分辨率(第一空间分辨率)的单位的计算结果(参见图14和图15)。
因此,使用逆模型计算,可以在宏观测量范围RZ2内以宏观测量分辨率的单位确定反映微观测量的测量值。
特别地,来自人造卫星210的图像被包括在输出范围中,因此,作为覆盖大于飞行体200所覆盖的范围的范围的图像,获得逆模型计算结果。
在该实施例中,逆模型计算部件27以宏观测量分辨率的单位确定测量对象的特征或环境响应作为计算结果。
这实现了确定测量对象的静态形状或特性或测量对象的动态形状或特性的感测。
例如,通过获得植物的形状、植被指数、关于光合作用的信息等,实现适合农业的遥感。
例如,作为关于光合作用的信息,即使在测量对象的形状或状态未知的情况下,也可以获得正确的SIF(太阳诱导叶绿素荧光)和从正确的SIF计算的各种类型的信息。
在该实施例中,宏观测量部件2在距离测量对象4(例如农田300)比微观测量部件3更长的距离处执行感测。
在宏观测量部件2位于相对远离测量对象4的位置的情况下,相对大型的装置或设备容易实现为宏观测量部件2或配备有宏观测量部件2的装置。
注意,微观测量部件3安装在飞行体200中,而宏观测量部件2安装在人造卫星210中,但是宏观测量部件2可以安装在飞行体200中,例如无人机。例如,宏观测量部件2可以安装在飞行在天空中更高位置的飞行体200中,以感测宏观测量范围RZ2。
在实施例中提及的示例中,逆模型计算部件27获取从聚类产生的每个测量区域的代表性个体的模型参数,切换从聚类产生的每个测量区域的模型参数,并将所得到的模型参数应用于逆模型(见图14和图15)。
通过计算来自微观测量部件的检测数据以确定从聚类产生的每个测量区域的模型参数,逆模型计算部件27可以对从聚类产生的各个测量区域使用不同的模型参数。
因此,即使在宏观测量范围RZ2内,也可以针对各个聚类区域(例如,针对种植不同作物的各个区域)获得适当的不同测量结果。
请注意,通过为每个聚类指定作物名称,不仅可以自动切换模型参数(例如高度),还可以自动切换模型本身(反映品种(例如番茄和玉米)之间形状差异的模型)。
在实施例中提及的示例中,所述聚类基于指定区域的用户输入。
例如,在农田300等中,在各个区域中种植不同的作物。例如,诸如农民之类的用户可以输入这样的信息。
这允许清楚地获取关于种植不同作物的区域、不同时段开始种植作物的区域等的信息。因此,可以针对每个区域适当地确定作物的特征或环境响应等。
此外,反映用户输入的聚类允许针对用户期望的每个区域获得测量结果。
注意,图12示出了基于用户输入、宏观测量值和微观测量值的聚类,但是,例如,可以执行仅基于用户输入的聚类计算,其中仅执行图12中的步骤S301。
另外,在本实施例中提及的示例中,基于来自宏观测量部件2的检测数据或来自微观测量部件3的检测数据执行聚类。
例如,通过使用来自宏观测量部件2或微观测量部件3的检测数据(或从检测数据计算的信息),可以将具有确定地不同的形状或状态的区域彼此区分,从而允许自动聚类计算。
这允许对种植不同作物的区域、在不同时期开始种植作物的区域等进行自动聚类。因此,可以准确地确定每个区域的作物的特性或环境响应等,而无需用户付出很大的努力。
图12示出了基于用户输入、宏观测量值和微观测量值的聚类。然而,例如,可以执行仅基于宏观测量值和微观测量值的聚类计算,其中仅执行图12中的步骤S302。此外,仅基于宏观测量值的聚类计算和仅基于微观测量值的聚类计算是可能的。
在实施例中提到的示例中,所述聚类基于指定区域的用户输入以及来自宏观测量部件的检测数据或来自微观测量部件的检测数据。
也就是说,执行反映手动输入和自动辨别两者的聚类计算。
在图12所示的示例中,聚类计算部件28使用来自宏观测量部件2的检测数据、来自微观测量部件3的检测数据和来自操作输入部件7的输入数据执行聚类。在这种情况下,除了用户输入的准确性和请求的反映之外,还提供了具有不同植被状况的区域的自动辨别系数,从而允许获得每个更合适区域的测量结果。
在图12中的示例中,通过自动辨别对基于用户输入的聚类进行子分割。然而,相比之下,可以通过用户输入对基于自动辨别的聚类进行子分割。
本实施例中描述的模型参数可以是植物的三维结构、植物高度、平均叶子角度(平均叶子倾角)、植物覆盖率、LAI、叶绿素浓度、土壤谱特性或太阳叶子比率中的任意一个或多个。
图15B中提到的示例是LAI、平均叶子角度和太阳叶子比率。可根据测量目的应用其他模型参数。
因此,模型参数可适当地应用于适合于植被遥感的测量。
在本实施例中提到的示例中,微观测量部件3包括可见光图像传感器、立体相机、用于激光图像检测和测距的传感器、偏振传感器或ToF传感器中的任何一个,作为微观测量传感器3S。
这些传感器适合于分析测量对象的特征、环境响应、范围、分布等,例如,形状分析等。
此外,传感器相对容易安装在飞行体200中,并且适合于作为小型无人飞行体(例如无人机)的飞行体200的操作。
在本实施例中提及的示例中,宏观测量部件2包括多谱相机、超谱相机、傅里叶变换红外光谱仪或红外传感器中的任何一个,作为宏观测量传感器2S。
这些传感器适用于分析各种物理属性值,如有关光合作用的信息。
此外,传感器相对难以安装在飞行体200中。因此,例如,将传感器安装在人造卫星210中允许促进作为小型无人飞行体(例如无人机)的飞行体200的操作。
本实施例中的信息处理装置1包括数据保存和输出部件30,该数据保存和输出部件30基于来自逆模型计算部件27的计算结果生成和输出图像数据。
在不改变计算结果的情况下,来自逆模型计算部件27的计算结果可能不适合作为人类观看的图像(难以理解来自图像的评估结果)。因此,数据保存和输出部件30将计算结果转换为适合呈现给人类的图像,并将所得到的图像输出到显示部件56、网络5或存储设备6。因此,可以向用户提供计算结果容易理解的图像。
在本实施例中提及的示例中,数据保存和输出部件30生成由来自逆模型计算部件27的计算结果的颜色映射产生的输出图像(参见图16)。
即,在针对与宏观分辨率单位相对应的每个区域获得来自逆模型计算部件27的计算结果的情况下,呈现给用户的图像被生成为其中颜色被分配给各个区域的图像。
因此,可以向用户提供允许逐色识别分析结果的图像。
在本实施例中提及的示例中,数据保存和输出部件30生成由来自逆模型计算部件27的计算结果的颜色映射产生的图像和另一图像的合成产生的输出图像(参见图17和图18)。
通过叠加或盖写,通过合成另一图像和由颜色映射产生的图像,数据保存和输出部件30可以向用户提供图像,该图像允许逐个颜色地识别每个区域的评估结果,同时允许借助另一图像识别每个区域。
在本实施例中提及的示例中,宏观测量部件2安装在人造卫星210中。
人造卫星210包括相对高的功能,并且可以容易地将大型传感器安装在人造卫星210中。因此,人造卫星210适于安装执行高级感测的宏观测量部件2。
例如,通过允许大量农民、感测实施组织等共享人造卫星210的宏观测量部件2,可以降低操作成本,并且可以有效地利用宏观测量传感器2S。
注意,代替人造卫星210,飞行体200或相对较大的飞行体可以配备宏观测量部件2,并在高于微观测量部件3执行感测的位置的位置处执行感测。
在本实施例中提及的示例中,微观测量部件3安装在能够被无线电控制或自动控制的飞行体200中。
能够无线电控制或自动控制的飞行体200的示例包括所谓的无人机、小型无线电控制固定翼飞机、小型无线电控制直升机等。
小飞行体200在距离测量对象(例如农田300)相对较低的高度处执行感测。在这种情况下,小飞行体200适于以高空间分辨率进行感测。
另外,避免在飞行体200中安装宏观测量部件2有助于小飞行体200的操作,并且使得能够降低感测成本。
在上述示例中,根据本实施例的信息处理装置1包括保持部件26,其保持从外部装置输入的逆模型计算程序和逆模型。
即,信息处理装置1允许从外部装置获取定义插值分析计算部件的计算算法的程序。
例如,从诸如网络5或存储设备6之类的外部装置获取用于插值分析计算的程序,并将其保存在保持部件26中。逆模型计算部件是在程序的基础上进行的。这使得信息处理装置1能够执行各种逆模型计算。
根据本实施例的程序使得信息处理装置1执行宏观测量分析计算处理,计算来自宏观测量部件2的检测数据,该宏观测量部件2对测量对象的宏观测量范围RZ2以宏观测量分辨率进行感测。该程序还使信息处理装置1执行微观测量分析计算处理,计算来自微观测量部件3的检测数据,该微观测量部件3以微观测量分辨率对微观测量范围RZ3进行感测,微观测量分辨率高于宏观测量分辨率,微观测量范围RZ3包括在宏观测量范围RZ2中。该程序还使得信息处理装置1执行使用基于由微观测量分析计算部件23确定的来自微观测量部件3的检测数据的模型参数作为用于使用来自宏观测量分析计算部件21的计算结果的逆模型计算的逆模型的参数的逆模型计算处理。
即,该程序使信息处理装置执行图9、图10和图14中的处理。
这样的程序有助于实现根据本实施例的图像处理装置1。
这样的程序可以预先存储在例如内置于诸如计算机装置的设备中的记录介质、包括CPU的微型计算机中的ROM等中。或者,程序可以临时或持久地存储在可移除记录介质中,例如半导体存储器、存储卡、光盘、磁光盘或磁盘。另外,这种可移除记录介质可以作为所谓的软件包提供。
另外,除了从可移除记录介质安装到个人计算机等中之外,还可以经由诸如LAN或因特网的网络从下载站点下载这样的程序。
注意,本文描述的效果仅是说明性的,并不打算被限制,并且可以产生其他效果。
本领域技术人员应当理解,根据设计要求和其他因素,可以进行各种修改、组合、子组合和变更,只要它们在所附权利要求或其等同物的范围内。
注意,本技术也可以如下所述进行配置。
(1)一种信息处理装置,包括:
宏观测量分析计算部件,被配置为计算来自宏观测量部件的检测数据,宏观测量部件适于以第一空间分辨率对测量对象的第一测量范围进行感测;
微观测量分析计算部件,被配置为计算来自微观测量部件的检测数据,微观测量部件适于以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及
逆模型计算部件,被配置为基于由微观测量分析计算部件确定的来自微观测量部件的检测数据,获取用于使用来自宏观测量分析计算部件的计算结果的逆模型计算的模型参数。
(2)根据上述(1)的信息处理装置,其中
逆模型计算部件使用基于由微观测量分析计算部件确定的来自微观测量部件的检测数据的模型参数,作为使用来自宏观测量分析计算部件的计算结果的逆模型计算中的逆模型的参数。
(3)根据上述(2)的信息处理装置,其中
逆模型计算部件使用基于第二测量范围的检测数据的模型参数来按第一空间分辨率的单位确定计算结果。
(4)根据上述(2)或(3)的信息处理装置,其中
逆模型计算部件按第一空间分辨率的单位确定作为计算结果的测量对象的特征或环境响应。
(5)根据上述(1)至(4)中任一项的信息处理装置,其中
宏观测量部件在距离测量对象比微观测量部件更长的距离处进行感测。
(6)根据上述(1)至(5)中任一项的信息处理装置,其中
逆模型计算部件获取进行聚类的每个测量区域的代表性个体的模型参数,为每个聚类测量区域切换模型参数,并将所得到的模型参数应用于逆模型。
(7)根据上述(6)的信息处理装置,其中
基于指定区域的用户输入进行所述聚类。
(8)根据上述(6)或(7)的信息处理装置,其中
基于来自宏观测量部件的检测数据或来自微观测量部件的检测数据进行所述聚类。
(9)根据上述(6)至(8)中任一项的信息处理装置,其中
基于指定区域的用户输入以及来自宏观测量部件的检测数据或来自微观测量部件的检测数据进行所述聚类。
(10)根据上述(1)至(9)中任一项的信息处理装置,其中
模型参数包括植物的三维结构、植物高度、平均叶子角度、植物覆盖率、LAI、叶绿素浓度、土壤的谱特性或太阳叶子比率中的任一个。
(11)根据上述(1)至(10)中任一项的信息处理装置,其中
微观测量部件包括作为微观测量传感器的可见光图像传感器、立体相机、用于激光图像检测和测距的传感器、偏振传感器或ToF传感器中的任一个。
(12)根据上述(1)至(11)中任一项的信息处理装置,其中
宏观测量部件包括作为宏观测量传感器的多谱相机、超谱相机、傅里叶变换红外光谱仪或红外传感器中的任一个。
(13)根据上述(1)至(12)中任一项的信息处理装置,还包括:
输出部件,被配置为基于来自逆模型计算部件的计算结果生成和输出输出图像数据。
(14)根据上述(13)的信息处理装置,其中
输出部件生成从来自逆模型计算部件的计算结果的颜色映射得到的输出图像数据。
(15)根据上述(13)的信息处理装置,其中
输出部件生成根据从来自逆模型计算部件的计算结果的颜色映射得到的图像与另一图像的合成而得到的输出图像数据。
(16)根据上述(1)至(15)中任一项的信息处理装置,其中
宏观测量部件安装在人造卫星中。
(17)根据上述(1)至(16)中任一项的信息处理装置,其中
微观测量部件安装在能够被无线电控制或自动控制的飞行体中。
(18)一种信息处理方法,包括:
由信息处理装置执行计算来自宏观测量部件的检测数据的宏观测量分析处理,宏观测量部件被配置成以第一空间分辨率对测量对象的第一测量范围进行感测;
由信息处理装置执行计算来自微观测量部件的检测数据的微观测量分析处理,微观测量部件被配置成以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及
由所述信息处理装置,执行使用基于在所述微观测量分析处理中确定的来自所述微观测量部件的检测数据的模型参数,作为用于使用宏观测量分析处理中的计算结果的逆模型计算的逆模型的参数的逆模型计算处理。
(19)一种使信息处理装置执行以下处理的程序:
计算来自宏观测量部件的检测数据的宏观测量分析处理,宏观测量部件被配置成以第一空间分辨率对测量对象的第一测量范围进行感测;
计算来自微观测量部件的检测数据的微观测量分析处理,微观测量部件被配置成以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及
使用基于在所述微观测量分析处理中确定的来自所述微观测量部件的检测数据的模型参数,作为用于使用宏观测量分析处理中的计算结果的逆模型计算的逆模型的参数的逆模型计算处理。
(20)一种感测系统,包括:
宏观测量部件,被配置为以第一空间分辨率对测量对象的第一测量范围进行感测;
微观测量部件,被配置为以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;
宏观测量分析计算部件,被配置为计算来自宏观测量部件的检测数据;
微观测量分析计算部件,被配置为计算来自微观测量部件的检测数据;以及
逆模型计算部件,被配置为使用基于由所述微观测量分析计算部件确定的来自所述微观测量部件的检测数据的模型参数,作为用于使用来自宏观测量分析计算部件的计算结果的逆模型计算的逆模型的参数。
参考标志列表
1信息处理装置、2宏观测量部件、2S宏观测量传感器、3微观测量部件、3S微观测量传感器、4测量对象、5网络、6存储设备、7操作输入部件、10数据输入部件、11传感器输入部件、12传感器输入部件、13程序和模型输入部件,20分析执行部件、21宏观测量分析计算部件、22宏观测量分析值缓冲器、23微观测量分析计算部件、24微观测量分析值缓冲器、25位置映射部件、26反向计算程序和模型保持部件、27反向模型计算部件、28聚类计算部件、30数据保存和输出部件、31分析数据缓冲器、32颜色映射部件、33图像合成部件、34图形生成部件、35图像输出部件、36数据输出部件、51CPU、52ROM、53RAM、54总线、55输入/输出接口、56显示部件、57输入部件、58扬声器、59存储部件、60通信部件、61驱动器、200飞行体、210人造卫星、220成像装置、250成像装置、300农田
Claims (20)
1.一种信息处理装置,包括:
宏观测量分析计算部件,被配置为计算来自宏观测量部件的检测数据,宏观测量部件适于以第一空间分辨率对测量对象的第一测量范围进行感测;
微观测量分析计算部件,被配置为计算来自微观测量部件的检测数据,微观测量部件适于以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及
逆模型计算部件,被配置为基于由微观测量分析计算部件确定的来自微观测量部件的检测数据,获取用于使用来自宏观测量分析计算部件的计算结果的逆模型计算的模型参数。
2.根据权利要求1所述的信息处理装置,其中
逆模型计算部件使用基于由微观测量分析计算部件确定的来自微观测量部件的检测数据的模型参数,作为使用来自宏观测量分析计算部件的计算结果的逆模型计算中的逆模型的参数。
3.根据权利要求2所述的信息处理装置,其中
逆模型计算部件使用基于第二测量范围的检测数据的模型参数来按第一空间分辨率的单位确定计算结果。
4.根据权利要求2所述的信息处理装置,其中
逆模型计算部件按第一空间分辨率的单位确定作为计算结果的测量对象的特征或环境响应。
5.根据权利要求1所述的信息处理装置,其中
宏观测量部件在距离测量对象比微观测量部件更长的距离处进行感测。
6.根据权利要求1所述的信息处理装置,其中
逆模型计算部件获取进行聚类的每个测量区域的代表性个体的模型参数,为每个聚类测量区域切换模型参数,并将所得到的模型参数应用于逆模型。
7.根据权利要求6所述的信息处理装置,其中
基于指定区域的用户输入进行所述聚类。
8.根据权利要求6所述的信息处理装置,其中
基于来自宏观测量部件的检测数据或来自微观测量部件的检测数据进行所述聚类。
9.根据权利要求6所述的信息处理装置,其中
基于指定区域的用户输入以及来自宏观测量部件的检测数据或来自微观测量部件的检测数据进行所述聚类。
10.根据权利要求1所述的信息处理装置,其中
模型参数包括植物的三维结构、植物的高度、平均叶子角度、植物覆盖率、LAI、叶绿素浓度、土壤的谱特性或太阳叶子比率中的任一个。
11.根据权利要求1所述的信息处理装置,其中
微观测量部件包括作为微观测量传感器的可见光图像传感器、立体相机、用于激光图像检测和测距的传感器、偏振传感器或ToF传感器中的任一个。
12.根据权利要求1所述的信息处理装置,其中
宏观测量部件包括作为宏观测量传感器的多谱相机、超谱相机、傅里叶变换红外光谱仪或红外传感器中的任一个。
13.根据权利要求1所述的信息处理装置,还包括:
输出部件,被配置为基于来自逆模型计算部件的计算结果生成和输出输出图像数据。
14.根据权利要求13所述的信息处理装置,其中
输出部件生成从来自逆模型计算部件的计算结果的颜色映射得到的输出图像数据。
15.根据权利要求13所述的信息处理装置,其中
输出部件生成根据从来自逆模型计算部件的计算结果的颜色映射得到的图像与另一图像的合成而得到的输出图像数据。
16.根据权利要求1所述的信息处理装置,其中
宏观测量部件安装在人造卫星中。
17.根据权利要求1所述的信息处理装置,其中
微观测量部件安装在能够被无线电控制或自动控制的飞行体中。
18.一种信息处理方法,包括:
由信息处理装置执行计算来自宏观测量部件的检测数据的宏观测量分析处理,宏观测量部件被配置成以第一空间分辨率对测量对象的第一测量范围进行感测;
由信息处理装置执行计算来自微观测量部件的检测数据的微观测量分析处理,微观测量部件被配置成以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及
由信息处理装置执行基于在微观测量分析处理中确定的来自微观测量部件的检测数据,获取用于使用在宏观测量分析处理中的计算结果的逆模型计算的模型参数的逆模型计算处理。
19.一种使信息处理装置执行以下处理的程序:
计算来自宏观测量部件的检测数据的宏观测量分析处理,宏观测量部件被配置成以第一空间分辨率对测量对象的第一测量范围进行感测;
计算来自微观测量部件的检测数据的微观测量分析处理,微观测量部件被配置成以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;以及
基于在微观测量分析处理中确定的来自微观测量部件的检测数据,获取用于使用在宏观测量分析处理中的计算结果的逆模型计算的模型参数的逆模型计算处理。
20.一种感测系统,包括:
宏观测量部件,被配置为以第一空间分辨率对测量对象的第一测量范围进行感测;
微观测量部件,被配置为以第二空间分辨率对第二测量范围进行感测,所述第二空间分辨率高于所述第一空间分辨率,第二测量范围被包括在测量对象的第一测量范围中;
宏观测量分析计算部件,被配置为计算来自宏观测量部件的检测数据;
微观测量分析计算部件,被配置为计算来自微观测量部件的检测数据;以及
逆模型计算部件,被配置为基于由微观测量分析计算部件确定的来自微观测量部件的检测数据,获取用于使用来自宏观测量分析计算部件的计算结果的逆模型计算的模型参数。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019-124764 | 2019-07-03 | ||
JP2019124764A JP7415348B2 (ja) | 2019-07-03 | 2019-07-03 | 情報処理装置、情報処理方法、プログラム、センシングシステム |
PCT/JP2020/025083 WO2021002280A1 (en) | 2019-07-03 | 2020-06-25 | Information processing apparatus, information processing method, program, and sensing system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114080540A true CN114080540A (zh) | 2022-02-22 |
Family
ID=71575716
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202080046901.8A Pending CN114080540A (zh) | 2019-07-03 | 2020-06-25 | 信息处理装置、信息处理方法、程序和感测系统 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20220299433A1 (zh) |
EP (1) | EP3994659A1 (zh) |
JP (1) | JP7415348B2 (zh) |
CN (1) | CN114080540A (zh) |
WO (1) | WO2021002280A1 (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6970946B1 (ja) * | 2021-03-07 | 2021-11-24 | 西日本技術開発株式会社 | 分布図作成装置、分布図作成方法、及び、プログラム |
CN117042595A (zh) * | 2021-03-30 | 2023-11-10 | 索尼集团公司 | 信息处理装置、信息处理方法和程序 |
WO2022209418A1 (ja) * | 2021-03-31 | 2022-10-06 | ソニーグループ株式会社 | 情報処理装置、情報処理方法、プログラム |
CN117980959A (zh) | 2021-09-27 | 2024-05-03 | 索尼半导体解决方案公司 | 信息处理装置和信息处理方法 |
JP7316004B1 (ja) * | 2022-10-24 | 2023-07-27 | 国立研究開発法人農業・食品産業技術総合研究機構 | 情報処理装置、情報処理方法、およびプログラム |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4185075B2 (ja) | 2005-07-08 | 2008-11-19 | 株式会社エヌ・ティ・ティ・データ | 緑被率マップ作成装置、緑被率マップ作成方法、及び、プログラム。 |
JP5162890B2 (ja) | 2006-12-01 | 2013-03-13 | 株式会社サタケ | リモートセンシングにおける補正方法 |
JP5560157B2 (ja) | 2010-10-19 | 2014-07-23 | 株式会社日立製作所 | スペクトル情報抽出装置 |
JP6082162B2 (ja) | 2014-03-28 | 2017-02-15 | 株式会社日立製作所 | 画像生成システム及び画像生成方法 |
IL236606B (en) | 2015-01-11 | 2020-09-30 | Gornik Amihay | Standards and methods for agricultural monitoring |
JP6507927B2 (ja) | 2015-08-12 | 2019-05-08 | コニカミノルタ株式会社 | 植物生育指標測定装置、該方法および該プログラム |
EA039345B1 (ru) | 2015-11-08 | 2022-01-17 | Агровинг Лтд | Способ получения и анализа аэрофотоснимков |
US20180306905A1 (en) * | 2017-04-20 | 2018-10-25 | Analog Devices, Inc. | Method of Providing a Dynamic Region of interest in a LIDAR System |
US10638038B2 (en) * | 2017-07-27 | 2020-04-28 | Stmicroelectronics (Research & Development) Limited | System and method for enhancing the intrinsic spatial resolution of optical sensors |
-
2019
- 2019-07-03 JP JP2019124764A patent/JP7415348B2/ja active Active
-
2020
- 2020-06-25 CN CN202080046901.8A patent/CN114080540A/zh active Pending
- 2020-06-25 US US17/597,079 patent/US20220299433A1/en active Pending
- 2020-06-25 WO PCT/JP2020/025083 patent/WO2021002280A1/en unknown
- 2020-06-25 EP EP20739477.6A patent/EP3994659A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
EP3994659A1 (en) | 2022-05-11 |
US20220299433A1 (en) | 2022-09-22 |
WO2021002280A1 (en) | 2021-01-07 |
JP2021012433A (ja) | 2021-02-04 |
JP7415348B2 (ja) | 2024-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zheng et al. | Improved estimation of rice aboveground biomass combining textural and spectral analysis of UAV imagery | |
WO2021002280A1 (en) | Information processing apparatus, information processing method, program, and sensing system | |
Tu et al. | Optimising drone flight planning for measuring horticultural tree crop structure | |
Lu et al. | Improved estimation of aboveground biomass in wheat from RGB imagery and point cloud data acquired with a low-cost unmanned aerial vehicle system | |
Roth et al. | Predicting cover crop biomass by lightweight UAS-based RGB and NIR photography: an applied photogrammetric approach | |
Aasen et al. | Multi-temporal high-resolution imaging spectroscopy with hyperspectral 2D imagers–From theory to application | |
Chen et al. | Improving estimates of fractional vegetation cover based on UAV in alpine grassland on the Qinghai–Tibetan Plateau | |
Haghighattalab et al. | Application of unmanned aerial systems for high throughput phenotyping of large wheat breeding nurseries | |
Ahmed et al. | Hierarchical land cover and vegetation classification using multispectral data acquired from an unmanned aerial vehicle | |
Frazer et al. | Gap Light Analyzer (GLA), Version 2.0: Imaging software to extract canopy structure and gap light transmission indices from true-colour fisheye photographs, users manual and program documentation | |
JP7415347B2 (ja) | 情報処理装置、情報処理方法、プログラム、センシングシステム | |
Lu et al. | Mapping vegetation biophysical and biochemical properties using unmanned aerial vehicles-acquired imagery | |
Heiskanen et al. | Use of airborne lidar for estimating canopy gap fraction and leaf area index of tropical montane forests | |
Sharifi | Estimation of biophysical parameters in wheat crops in Golestan province using ultra-high resolution images | |
Mafanya et al. | Radiometric calibration framework for ultra-high-resolution UAV-derived orthomosaics for large-scale mapping of invasive alien plants in semi-arid woodlands: Harrisia pomanensis as a case study | |
Rasmussen et al. | The challenge of reproducing remote sensing data from satellites and unmanned aerial vehicles (UAVs) in the context of management zones and precision agriculture | |
EP4064182A1 (en) | Information processing device, information processing method, and program | |
CN108776106A (zh) | 一种基于无人机低空遥感的作物长势监测方法及系统 | |
Kurihara et al. | Unmanned Aerial Vehicle (UAV)-based hyperspectral imaging system for precision agriculture and forest management | |
CN115453555A (zh) | 一种草地生产力的无人机快速监测方法及系统 | |
Jafarbiglu et al. | Impact of sun-view geometry on canopy spectral reflectance variability | |
Komárek | The perspective of unmanned aerial systems in forest management: do we really need such details? | |
Panigada et al. | Chlorophyll concentration mapping with MIVIS data to assess crown discoloration in the Ticino Park oak forest | |
Wang et al. | The impact of variable illumination on vegetation indices and evaluation of illumination correction methods on chlorophyll content estimation using UAV imagery | |
Teixeira et al. | Evaluation of the radiometric accuracy of images obtained by a Sequoia multispectral camera |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |