CN102105818A - 日照量的评价方法和评价装置 - Google Patents

日照量的评价方法和评价装置 Download PDF

Info

Publication number
CN102105818A
CN102105818A CN2009801293543A CN200980129354A CN102105818A CN 102105818 A CN102105818 A CN 102105818A CN 2009801293543 A CN2009801293543 A CN 2009801293543A CN 200980129354 A CN200980129354 A CN 200980129354A CN 102105818 A CN102105818 A CN 102105818A
Authority
CN
China
Prior art keywords
evaluation
mentioned
sunshine amount
unit
essential elements
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.)
Granted
Application number
CN2009801293543A
Other languages
English (en)
Other versions
CN102105818B (zh
Inventor
笹川正
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Pasco Corp
Original Assignee
Pasco Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Pasco Corp filed Critical Pasco Corp
Publication of CN102105818A publication Critical patent/CN102105818A/zh
Application granted granted Critical
Publication of CN102105818B publication Critical patent/CN102105818B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/12Sunshine duration recorders

Landscapes

  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Image Processing (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明的目的在于提供一种通过考虑地形条件而可提高精度的日照量的评价方法和评价装置。静止卫星图像数据(1)与三维地图模型(2)结合,根据静止卫星图像数据(1)的地面分辨率,设定基于平面坐标的网格(3),根据三维地图模型(2)和通过摄像日期时间确定的太阳位置,对背阴区域(4)进行运算,判断上述网格(3)的各单元(8)是否属于背阴区域(4),同时,参照表格(7)获得评价对象区域内各单元(8)的日照量评价要素(6),其中在该表格(7)中,静止卫星图像数据(1)中各像元的像元值(5)与各像元所对应的区域的日照量评价要素(6)相互对应,并且,当该单元(8)被判断为属于背阴区域(4)时,对该日照量评价要素(6)加入规定的校正计算,获得地形条件校正后的评价要素。

Description

日照量的评价方法和评价装置
技术领域
本发明涉及日照量的评价方法和评价装置。
背景技术
过去,对规定区域的日照量的评价采用静止卫星图像的技术,在专利文献1中所记载的技术被人们所知。该已有例,从静止卫星图像的像元值中获得与各像元对应区域的反射值与亮度温度,并与预先获得的过去同一时间段的晴天时的反射值等做比较,由此判断是晴天还是阴天,根据晴天或阴天各自所规定的计算式计算出日照量。计算式中,采用太阳天顶角、太阳常数、太阳-地球间的距离、臭氧吸收的透射率等各种透射率、云的日照吸收系数等,另外还采用上述反射值。又,上述反射值是像元的反射率除以太阳天顶角的余弦得到的。
另外,专利文献2中记载了如下方式,利用反射量与日照量的基于统计关系得到回归系数直接关联,根据反射量,也就是静止气象卫星图像的像元值来直接推算日照量。
现有技术文献
专利文献
专利文献1:日本特开平11-211560号公报
专利文献2:日本特开2004-118639号公报
发明内容
如以上已有例那样,静止卫星图像运用于日照量的评价时,仅仅考虑像元值,或者是再加上太阳和大气性质的天体、气象条件等,并没有考虑地形条件,因此,精度差,例如在农业收成或热岛现象分析的利用中,没有发挥出有效的评价机能。
本发明为解决上述缺陷,目的在于提供一种通过考虑地形条件而可提高精度的日照量的评价方法和评价装置。
根据本发明,上述目的通过如下方式实现,提供一种日照量的评价方法,该方法为:
在同时获得静止卫星图像数据1和摄像日期时间数据后,
将上述静止卫星图像数据1与三维地图模型2结合,
接着,根据结合数据中静止卫星图像数据1的地面分辨率,设定基于平面坐标的网格3,按照该网格3设定结合数据中的规定区域,并作为评价对象区域,
之后,根据上述三维地图模型2和通过上述摄像日期时间数据确定的太阳位置,对在上述评价对象区域内产生的背阴区域4进行计算,以规定的判断基准,判断上述网格3的各单元8是否属于背阴区域4,
同时,参照表格7获得上述评价对象区域内各单元8的日照量评价要素6,其中该表格7中,静止卫星图像数据1中各像元的像元值5与各像元所对应的区域的日照量评价要素6相互对应,并且,当该单元8被判断为属于背阴区域4时,取代上述日照量评价要素6而获得地形条件校正后的评价要素,对评价对象区域的日照量进行评价,上述地形条件校正后的评价要素为对上述日照量评价要素6加入规定的校正计算而算出。
根据本发明,能够考虑由地形产生的背阴,对日照量进行评价,因此,可以极大地提高因地形而导致的有背阴影响的区域的评价精度。为得到精确位置信息,采用三维地图数据21取代过去与静止卫星图像数据1结合的二维地图数据,根据三维地图模型2和太阳位置计算背阴区域4,从而可以掌握背阴的产生。
因此根据本发明,可对农业收成、农作物栽培的适当地点选择、精密农业的参数、热岛现象的进行分析等时以要求的精度,进行日照量评价。
另外,按照静止卫星图像数据1的地面分辨率,在静止卫星图像数据1与三维地图模型2结合形成的结合数据中设定基于平面坐标的网格3,并以该网格3的单元8为单位,可以在预计有较高分辨率的三维地图数据21和只能期待较低分辨率的静止卫星图像数据1之间,使分辨率不会徒然降低的情况下对背阴区域4进行判断。此时,以单元8的中心作为判断代表点28,在根据与太阳位置的关系对背阴区域4进行计算的三维地图模型2中,可将上述判断代表点28是否属于背阴区域4作为单元8单位的背阴判断基准。
从没有考虑产生背阴的已有技术中导出日照量评价要素6,并对其加入规定的校正计算,即可得到如上述那样判断为属于背阴区域4时的日照量的评价。该校正计算,可以考虑因背阴而损失的直达日照量,或根据经验法则的统计而设定的方式适当地进行。
另外,对于日照量评价的地形条件的考虑,除以上因地形产生的背阴之外,对于存在倾斜地形时的背阴也是有用的,此情况下,按下述方式构成日照量的评价方法。该方式为:
在同时获得静止卫星图像数据1和摄像日期时间数据后,
将上述静止卫星图像数据1与三维地图模型2结合,
接着,根据结合数据中静止卫星图像数据1的地面分辨率,设定基于平面坐标的网格3,按照该网格3设定结合数据中的规定区域,并作为评价对象区域,
之后,按照基于上述三维地图模型2所规定的判断基准,求得评价对象区域中网格3的各单元8的代表倾斜角以及代表倾斜方位,根据上述摄像日期时间确定的太阳位置、将假定各单元8为与日照方向垂直的平面(水平面)时接受的日照量评价,变换为获得上述代表倾斜角以及代表倾斜方位时接受的日照量评价时针对每个单元8计算倾斜校正系数,
并且,参照表格7获得上述评价对象区域中各单元8的日照量评价要素6,其中该表格7中,静止卫星图像数据1中各像元的像元值5与日照量评价要素6相互对应,该日照量评价要素6为各像元对应区域假定为水平面时的每个区域的要素。之后,对该日照量评价要素6按照上述校正系数做校正计算,得出单元8,并获得每个单元8的地形条件校正后的评价要素,对评价对象区域的日照量进行评价。
在本发明中,对日照量的评价,考虑接受日照的评价对象区域中地表面的倾斜,可以极大地提高倾斜地的日照量评价的精度。可以从三维地图数据21获得倾斜信息,比如,与上述同样,以单元8的中心作为判断代表点28,根据与周边单元8的判断代表点28的标高数据的差值,可以计算并判断各单元8的代表倾斜角以及代表倾斜方位。
因此根据本发明,与上述同样,也可对应农业收成、农作物栽培的适当地点选择、精密农业的参数、热岛现象的分析等时要求的精度,进行日照量评价。
可根据倾斜信息计算日照量评价的校正系数,也可对从没有考虑地表面倾斜的已有技术中导出的评价值,根据校正系数进行校正计算,即可反映出倾斜信息的日照量评价。可根据比如兰伯德余弦定律(Lambert’s cosine)等对校正系数作适当确定。
另外,通过由上述那样得到地表面的倾斜信息,根据太阳位置和地表面的倾斜信息可计算从地表面反射的太阳光的反射日照量,并添加到日照量的评价中。
以上2个发明中,作为与静止卫星图像数据1的像元值5一一对应而导出的日照量评价要素6,像上述已有例那样,可照原样采用日照量,或采用在掌握日照量时作为主要要素的云的日照透射系数等。另外,作为日照量评价的评价量,同样地,可照原样采用日照量,或采用合适的指数。在此情况下,采用作为日照量的评价指数可发挥机能的上述云的日照透射系数时,可以根据确定的太阳位置采用云的日照透射系数对日照量进行计算并评价。
以上描述的是通过采用所有因地形产生的背阴信息和地表面的倾斜信息,考虑了地形条件而进行日照量评价的情况,但是,也可将两者的信息合并后考虑,使得评价精度进一步提高。再加上,在上述中,对于作为评价要素的静止卫星图像数据1中,采用不同摄像时间的多个数据,从多个静止卫星图像数据1中获得多个日照量评价要素6和校正后的评价要素,并求出平均值,计算每个规定时间的日照量评价要素6等,可以把握日照量评价要素6等的时间的变化。
由以上的说明可知,本发明可提供一种可以通过考虑地形条件使精度提高的日照量的评价方法和评价装置,通过日照量的评价,可以对农业的收成、农作物栽培地点选定、精密农业的参数、以及热岛现象进行分析。
附图说明
图1是表示本发明的概略处理程序的流程图;
图2是本发明的计算机的方框图;
图3是表示日照量的计算步骤的详细处理程序的流程图;
图4是计算日照量的处理方案图,图4(a)是表示计算所对应的模型图,图4(b)是表示像元值·云的日照透射系数表格7的图。
具体实施方式
图1所表示的是本发明中日照量评价的处理程序的大致流程,图2表示的是实施该处理程序的计算机方框图。该计算机包括:具有获得机构9的输入模块20,存储三维地图数据21的存储模块22,运算模块23,将运算结果输出到图外的监视器等的输出模块24。日照量评价的处理首先由实施下述步骤(步骤S 1)开始,即,获得静止卫星图像数据1,该静止卫星图像数据1是最初计算机通过上述获得机构9,比如通过网上或合适的存储媒体从日本气象协会等获得。获得该静止卫星图像数据1时,同时获得由摄像日期时间信息和摄像区域的经度和纬度形成的坐标信息作为标注(annotation)数据。
上述静止卫星图像数据1是下述像元的集合,该像元的密度以地面分辨率为准。并且上述静止卫星图像数据1的像元值5为各像元的亮度数据。该像元值5是摄像对象的反射数据,比如反射率存在很大差异的陆地与海的像元值就有很大差异,因此,可以将像元值5中差异较大、呈线状配置的部位作为海岸线等来判断。
为了利用以上的静止卫星图像数据1对日照量进行评价,上述运算模块23,如图2所示,包括地图模型结合机构10、评价对象区域设定机构11、背阴判断机构12、倾斜校正运算机构14以及日照量评价机构13。
地图模型结合机构10将静止卫星图像数据1与上述三维地图数据21按如下方式相结合,先对静止卫星图像数据1做几何校正,之后将上述坐标信息作为媒介,与三维地图数据21配准并结合。又,该配准中,由于静止卫星图像数据1的坐标精度低,导致精度稍差,因此,在地图模型结合机构10中设置图中未表示的结合点设定部,提取出上述海岸线等作为结合点,并对在配准时根据结合点的匹配程度产生的误差进行修正。
评价对象区域设定机构11,设定日照量的评价区域,如图2所示,包括网格设定机构25、代表值设定机构26以及区域设定机构27。网格设定机构25在静止卫星图像数据1设定网格3,更准确的说,如上述那样,是将静止卫星图像数据1与三维地图数据21相结合而形成的合成数据上设定网格3。该网格3与静止卫星图像数据1的地面分辨率,也就是与静止卫星图像数据1上的一个像元对应的地面区域相对应。另外,代表值设定机构26,在形成的网格3的各单元8、8……中设定单一的代表坐标信息,在该实施方式中,将各单元8的中心点作为判断代表点28,并适用上述坐标信息。区域设定机构27,在合成数据内将指定的区域设定为评价对象区域,其中,该指定区域指的是通过上述输入模块20中包含的图中未标识的鼠标等输入而指定的区域,换言之仅仅是对应于指定区域的单个或多个单元8。另外,输入模块20中没有指定的评价对象区域时,也可进行将静止卫星图像数据1上的全部区域设定为评价对象区域的处理。
背阴判断机构12,具有背阴计算机构29和单元判断机构30,并利用上述合成数据中含有的三维地图数据21的三维地图模型2来判断背阴区域4。计算背阴区域4时所需的太阳位置可以根据上述摄像日期时间信息来确定,背阴计算机构29,参照存储模块22设置的日期时间·太阳坐标表格31,从摄像日期时间信息中获得太阳位置,并在三维地图模型2上对日照进行模拟,计算评价对象区域内的背阴区域4。单元判断机构30,在评价对象区域内以单元8为单位来判断是否属于背阴区域4,该实施方式中,将上述单元8的中心点作为判断代表点28,根据该判断代表点28是否属于背阴区域4来进行判断。
倾斜校正运算机构14包括倾斜提取机构32和倾斜校正系数设定机构33,其中,倾斜提取机构32根据上述三维地图模型2设定评价对象区域内的各单元8的倾斜信息,也就是设定各单元8的代表倾斜角和代表倾斜方位;该倾斜校正系数设定机构33在取得该代表倾斜角和代表倾斜方位时,计算各单元8接受的日照量的倾斜校正系数根据每个单元8设定,该各单元8接受的日照量是假设各单元8为水平面时接受的假定日照量。在本实施方式中,由位于该单元8的纵、横、斜周边8个方向的8个单元构成周边单元8,倾斜提取机构32根据周边单元8的标高信息来确定各单元8的倾斜信息。更具体地讲,例如,求得所有8个方向的周边单元8与位于中心位置的中心单元8的标高差,将所有8个方向的合计差值为最小时的中心单元8的倾斜角和倾斜方位确定为代表倾斜角和代表倾斜方位。另外,从假定各单元8为水平面时的日照量变换为取得代表倾斜角和代表倾斜方位时的日照量的倾斜校正系数以下述方式设定,即倾斜校正系数设定机构33将该倾斜校正系数设定在每个单元8上。该倾斜校正系数为相对于单元8的水平面倾斜信息的函数,可由例如逐次利用根据摄像日期时间确定的太阳位置而计算得出的方式形成。
日照量评价机构13,在本实施方式中,由计算评价对象区域内的自身日照量的机构构成,而且,本实施方式中的评价对象区域的日照量,是求出该评价对象区域内的各单元8的平均日照量而计算得出的。另外,通过求出摄像时间不同的多个静止卫星图像数据1、1……中的平均日照量,算出每小时、每天、每个月等任意时间单位的日照量。
假定计算对象区域为水平面,并且,对于通过现有的没有考虑由地形产生的背阴的日照量计算方法而推算的日照量,计算出的每个单元8的日照量可以反映出因上述背阴和倾斜而造成的日照量衰减。对于上述推算日照量的计算,具体地讲,参照存储模块22的像元值·云的日照透射系数表格7,从与各单元8对应的静止卫星图像数据1中的像元值5获得各单元8的云的日照透射系数(日照量评价要素6),另外,利用上述太阳位置导出太阳天顶角等,或从存储模块22读出等方式,酌情获得泛用的日照量计算式所需的各种计算参数,来计算上述日照量。
另外,上述像元值·云的日照透射系数表格7,将认定有一定关联性的像元值5与云的日照透射系数6相互对应,可以按如图4(b)所示的那样的方式构成。例如,可按如下方式形成,即,存储过去的像元值5和实测日照量,通过统计的方法推导出对应关系的方式。另外,对应于上述单元8的像元值5,一般地,是上述合成数据中该单元8的像元值5,但当对于地表面的太阳光入射角为极大的钝角时的太阳位置的情况下,也可认为是射入相关单元8的地表面的太阳光穿过其他单元8上面的云的情况,因此,也可以采用基于对应于入射角的其他单元8的像元值5的云的日照透射指数6。
另外,像这样得到推算的日照量后,当对单元判断机构30的相关单元8的判定为背阴时,将日照量调为0,或按照一定比例减少日照量,然后,根据上述倾斜校正系数重复进行校正计算,从而得出每个单元8的日照量(地形条件校正后的评价要素)。按上述方式求出属于评价对象区域的所有单元8的日照量,求出平均值,或者,从多个静止卫星图像数据1的计算结果求得平均值得到的日照量,由输出模块24输出,在图外的监视器中显示。例如,将与日照量对应而选定的颜色与评价对象区域重合,在地图数据上来显示日照量等进行显示。
根据图1,对以上的运算模块23的计算机操作流程进行说明。运算模块23实行如下步骤(步骤S2),通过获得机构9获得上述静止卫星图像数据1等,则从存储模块22读出三维地图数据21,通过地图模型结合机构10将静止卫星图像数据1与三维地图数据21相结合,接着,实施下述步骤(步骤S3),在像上述那样获得的结合数据上,通过网格设定机构25,按照平面坐标设定网格3,例如,该网格3是预先设定为与静止卫星图像数据1的地面分辨率相对应的网格。
之后,实施下述步骤(步骤S4),通过代表值设定机构26,利用上述三维地图数据21中的三维坐标信息,参照网格3中各单元8的判断代表点28的三维坐标信息,在各单元8中将该坐标信息作为与坐标相关的属性等而设定,同时,实施下述步骤(步骤S5),通过区域设定机构27,利用上述三维地图数据21中的地图和坐标信息,根据网格3的划分,对日照量的评价对象区域进行设定。
再之后,实施下述步骤(步骤S6),通过背阴运算机构29根据通过获得机构9获得的摄像日期时间信息,参照日期时间.太阳坐标表格31,确定太阳位置,通过该太阳位置与通过上述三维地图数据21能够构成的三维地图模型2,对评价对象区域中产生的背阴区域4进行计算后,实施下述步骤(步骤S7),通过单元判断机构30,判定评价对象区域内的各单元8是否属于背阴区域4,将判定结果作为各单元8的与背阴相关的属性而设定。于是,例如图4(a)所示,将静止卫星图像数据1面、三维地图数据21面、以及根据该三维地图数据21得到的三维地图模型2相互配准,像图中虚线所示,与像元的地面分辨率相对应设定网格3,对上述三维地图模型2,根据摄像日期时间信息确定位置的太阳34的阳光入射,产生图中阴影部分所示的背阴区域4的情况下,三维地图模型2的轮廓线上用白圈表示的各单元8的判断代表点28,如属于背阴区域4的单元8在三维地图数据21上用粗线表示的那样,被判断为背阴,剩余的单元8被判断为不是背阴。
另一方面,实施下述步骤(步骤S8),通过倾斜提取机构32,参照三维地图数据21从得到的与邻接的单元8的标高差值中提取各单元8的代表倾斜角和代表倾斜方位,然后,实施下述步骤(步骤S9),以射入根据该倾斜信息和太阳位置确定的各单元8的太阳光的入射角等为基准,对静止卫星图像数据1的像元值5和日照量的倾斜校正系数进行计算,并将倾斜校正系数作为各单元8与倾斜相关的属性而设定。
最后,实施下述步骤(步骤S10),通过日照量评价机构13,根据各单元8的像元值5、坐标信息、与背阴相关的属性、倾斜校正系数以及太阳位置,计算出评价对象区域的日照量。具体地说,如图3所示,首先,参照像元值·云的日照透射系数表格7,从各单元8的像元值5中获得云的透射系数6,不考虑因地形造成的背阴,从该云的透射系数6和太阳位置,以及各单元8的坐标信息等,计算假设各单元8为水平面时的推算日照量(步骤S10-1)。接着,参照各单元8的背阴的相关属性等(步骤S10-2),对于判断为不是背阴的单元8,进一步照原样参照倾斜校正系数,利用该倾斜校正系数,由推算日照量计算出日照量(步骤S10-4)。
另一方面,相反地,对于上述判断为背阴的单元8,设定从上述日照量变换为考虑了背阴的日照量时的背阴校正系数(步骤S10-3),之后,与上述倾斜校正系数结合,利用背阴校正系数从推算日照量中计算出日照量(步骤S10-4)。上述日照量的计算对应于评价对象区域内的所有单元8、8……并反复进行(步骤S10-5),得到评价对象区域内的所有单元8、8……的日照量后,求得平均值,计算评价对象区域的日照量(步骤S10-6)。另外,在利用摄像时间不同的多个静止卫星图像数据1,求得单位时间比较长的日照量的本发明中,对多个静止卫星图像数据1反复进行上述的评价对象区域的日照量计算(步骤S10-7),之后,对每个计算出的多个评价对象区域的日照量求平均值后,便可求得单位时间长的日照量(步骤S10-8)。
另外,上述中描述了在倾斜的计算以及取得对单元8的倾斜计算结果之前进行背阴的判定和取得对单元8的背阴判定结果的情况,但也可在背阴的判断等之前进行倾斜的计算等处理。
标号说明
标号1表示静止卫星图像数据;
标号2表示三维地图模型;
标号3表示网格;
标号4表示背阴区域;
标号5表示像元值;
标号6表示日照量评价要素;
标号7表示表格;
标号8表示单元;
标号9表示获得机构;
标号10表示地图模型结合机构;
标号11表示评价对象区域设定机构;
标号12表示背阴判断机构;
标号13表示日照量评价机构;
标号14表示倾斜校正运算机构。
权利要求书(按照条约第19条的修改)
1.一种日照量的评价方法,该方法中:
在同时获得静止卫星图像数据和摄像日期时间数据后,
将上述静止卫星图像数据与三维地图模型结合,
接着,根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域,
之后,根据上述三维地图模型和通过上述摄像日期时间数据确定的太阳位置,对在上述评价对象区域内产生的背阴区域进行计算,以规定的判断基准,判断上述网格的各单元是否属于背阴区域,
同时,参照表格获得上述评价对象区域内各单元的日照量评价要素,其中该表格中,静止卫星图像数据中各像元的像元值与各像元所对应的区域的日照量评价要素相互对应,并且,当该单元被判断为属于背阴区域时,取代上述日照量评价要素而获得地形条件校正后的评价要素,对评价对象区域的日照量进行评价,上述地形条件校正后的评价要素为通过对上述日照量评价要素加入规定的校正计算而算出。
2.(修改后)一种日照量的评价方法,该方法中:
在同时获得静止卫星图像数据和摄像日期时间数据后,
将上述静止卫星图像数据与三维地图模型结合,
接着,根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域,
之后,按照根据上述三维地图模型所规定的判断基准,求得评价对象区域中网格的各单元的代表倾斜角以及代表倾斜方位,根据上述摄像日期时间数据确定的太阳位置、将假定各单元为水平面时接受的日照量评价,变换为获得上述代表倾斜角以及代表倾斜方位时接受的日照量评价时针对每个单元计算倾斜校正系数,
并且,参照表格获得上述评价对象区域中各单元的日照量评价要素,其中该表格中,静止卫星图像数据中各像元的像元值与各像元对应区域假定为水平面时的每个区域的日照量评价要素相对应,之后,通过对该日照量评价要素按照上述倾斜校正系数做校正计算,并计算取得每个单元的地形条件校正后的评价要素,对评价对象区域的日照量进行评价。
3.权利要求1所述的日照量的评价方法,其中,在上述各单元中,分别设定以单元的中心点为基准的三维坐标数据,
并且,根据上述摄像日期时间数据确定的太阳位置以上述中心点的三维坐标在三维地图模型中是否属于背阴区域来判断各单元是否属于背阴区域。
4.权利要求2所述的日照量的评价方法,其中,在上述各单元中,分别设定以单元的中心点为基准的三维坐标数据,
并且,根据周边单元的中心点标高数据与中心单元的中心点标高数据的差值,计算各单元的代表倾斜角以及代表倾斜方位。
5.权利要求1所述的日照量的评价方法,其中,在上述表格中,设定云的日照透射系数作为日照量评价要素,
并且,根据通过上述摄像日期时间数据指定的太阳位置与上述日照透射系数计算日照量。
6.一种日照量的评价装置,该装置包括:
获得机构,该获得机构同时获得静止卫星图像数据和摄像日期时间数据;
地图模型结合机构,该地图模型结合机构将上述静止卫星图像数据与三维地图模型结合;
评价对象区域设定机构,该评价对象区域设定机构根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域;
背阴判断机构,该背阴判断机构根据上述三维地图模型和通过上述摄像日期时间确定的太阳位置,对在上述评价对象区域内产生的背阴区域进行运算,以规定的判断基准,判断上述网格的各单元是否属于背阴区域;
日照量评价机构,参照表格获得上述评价对象区域内各单元的日照量评价要素,其中,该表格中,上述静止卫星图像数据中各像元的像元值与各像元所对应的区域的日照量评价要素相互对应,并且,当该单元被判断为属于背阴区域时,取代上述日照量评价要素而获得地形条件校正后的评价要素,对评价对象区域的日照量进行评价,上述地形条件校正后的评价要素对上述日照量评价要素加入规定的校正计算而算出。
7.一种日照量的评价装置,该装置包括:
获得机构,该获得机构同时获得静止卫星图像数据和摄像日期时间数据;
地图模型结合机构,该地图模型结合机构将上述静止卫星图像数据与三维地图模型结合;
评价对象区域设定机构,该评价对象区域设定机构根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域;
倾斜校正运算机构,该倾斜校正运算机构按照根据上述三维地图模型所规定的判断基准,求得评价对象区域中网格的各单元的代表倾斜角以及代表倾斜方位,根据上述摄像日期时间数据确定的太阳位置、将假定各单元为水平面时接受的日照量评价,变换为获得上述代表倾斜角以及代表倾斜方位时接受的日照量评价时针对每个单元计算校正系数,
日照量评价机构,该日照量评价机构参照表格获得上述评价对象区域中各单元的日照量评价要素,其中该表格中,上述静止卫星图像数据中各像元的像元值与各像元对应区域假定为水平面时的每个区域的日照量评价要素相对应,之后,对该日照量评价要素按照上述校正系数做校正计算,并计算取得每个单元的地形条件校正后的评价要素,对评价对象区域的日照量进行评价。
8.一种日照量的评价方法,在该方法中:
在同时获得静止卫星图像数据和摄像日期时间数据后,
将上述静止卫星图像数据与三维地图模型结合,
接着,根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域,
之后,根据上述三维地图模型和通过上述摄像日期时间数据指定的太阳位置,对在上述评价对象区域内产生的背阴区域进行运算,以规定的判断基准判断上述网格的各单元是否属于背阴区域,
同时,按照根据上述三维地图模型所规定的判断基准,求得评价对象区域中网格的各单元的代表倾斜角以及代表倾斜方位,根据上述太阳位置、将假定各单元为水平面时接受的日照量评价,变换为取得上述代表倾斜角以及代表倾斜方位时接受的日照量评价时针对每个单元计算校正系数,
并且,参照表格获得上述评价对象区域中各单元的日照量评价要素,其中该表格中,静止卫星图像数据中各像元的像元值与各像元对应区域假定为水平面时的每个区域的日照量评价要素相对应,之后,当该单元被判断为属于背阴区域时,将对上述日照量评价要素进行校正计算,同时,
取代上述日照量评价要素而获得地形条件校正后的评价要素,对评价对象区域的日照量进行评价,上述地形条件校正后的评价要素为根据上述校正系数做进一步校正计算而算出。

Claims (8)

1.一种日照量的评价方法,该方法中:
在同时获得静止卫星图像数据和摄像日期时间数据后,
将上述静止卫星图像数据与三维地图模型结合,
接着,根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域,
之后,根据上述三维地图模型和通过上述摄像日期时间数据确定的太阳位置,对在上述评价对象区域内产生的背阴区域进行计算,以规定的判断基准,判断上述网格的各单元是否属于背阴区域,
同时,参照表格获得上述评价对象区域内各单元的日照量评价要素,其中该表格中,静止卫星图像数据中各像元的像元值与各像元所对应的区域的日照量评价要素相互对应,并且,当该单元被判断为属于背阴区域时,取代上述日照量评价要素而获得地形条件校正后的评价要素,对评价对象区域的日照量进行评价,上述地形条件校正后的评价要素为通过对上述日照量评价要素加入规定的校正计算而算出。
2.一种日照量的评价方法,该方法中:
在同时获得静止卫星图像数据和摄像日期时间数据后,
将上述静止卫星图像数据与三维地图模型结合,
接着,根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域,
之后,按照根据上述三维地图模型所规定的判断基准,求得评价对象区域中网格的各单元的代表倾斜角以及代表倾斜方位,根据上述摄像日期时间数据确定的太阳位置、将假定各单元为水平面时接受的日照量评价,变换为获得上述代表倾斜角以及代表倾斜方位时接受的日照量评价时针对每个单元计算倾斜校正系数,
并且,参照表格获得上述评价对象区域中各单元的日照量评价要素,其中该表格中,静止卫星图像数据中各像元的像元值与各像元对应区域假定为水平面时的每个区域的日照量评价要素相对应,之后,通过对该日照量评价要素按照上述校正系数做校正计算,并计算取得每个单元的地形条件校正后的评价要素,对评价对象区域的日照量进行评价。
3.权利要求1所述的日照量的评价方法,其中,在上述各单元中,分别设定以单元的中心点为基准的三维坐标数据,
并且,根据上述摄像日期时间数据确定的太阳位置以上述中心点的三维坐标在三维地图模型中是否属于背阴区域来判断各单元是否属于背阴区域。
4.权利要求2所述的日照量的评价方法,其中,在上述各单元中,分别设定以单元的中心点为基准的三维坐标数据,
并且,根据周边单元的中心点标高数据与中心单元的中心点标高数据的差值,计算各单元的代表倾斜角以及代表倾斜方位。
5.权利要求1所述的日照量的评价方法,其中,在上述表格中,设定云的日照透射系数作为日照量评价要素,
并且,根据通过上述摄像日期时间数据指定的太阳位置与上述日照透射系数计算日照量。
6.一种日照量的评价装置,该装置包括:
获得机构,该获得机构同时获得静止卫星图像数据和摄像日期时间数据;
地图模型结合机构,该地图模型结合机构将上述静止卫星图像数据与三维地图模型结合;
评价对象区域设定机构,该评价对象区域设定机构根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域;
背阴判断机构,该背阴判断机构根据上述三维地图模型和通过上述摄像日期时间确定的太阳位置,对在上述评价对象区域内产生的背阴区域进行运算,以规定的判断基准,判断上述网格的各单元是否属于背阴区域;
日照量评价机构,参照表格获得上述评价对象区域内各单元的日照量评价要素,其中,该表格中,上述静止卫星图像数据中各像元的像元值与各像元所对应的区域的日照量评价要素相互对应,并且,当该单元被判断为属于背阴区域时,取代上述日照量评价要素而获得地形条件校正后的评价要素,对评价对象区域的日照量进行评价,上述地形条件校正后的评价要素对上述日照量评价要素加入规定的校正计算而算出。
7.一种日照量的评价装置,该装置包括:
获得机构,该获得机构同时获得静止卫星图像数据和摄像日期时间数据;
地图模型结合机构,该地图模型结合机构将上述静止卫星图像数据与三维地图模型结合;
评价对象区域设定机构,该评价对象区域设定机构根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域;
倾斜校正运算机构,该倾斜校正运算机构按照根据上述三维地图模型所规定的判断基准,求得评价对象区域中网格的各单元的代表倾斜角以及代表倾斜方位,根据上述摄像日期时间数据确定的太阳位置、将假定各单元为水平面时接受的日照量评价,变换为获得上述代表倾斜角以及代表倾斜方位时接受的日照量评价时针对每个单元计算校正系数,
日照量评价机构,该日照量评价机构参照表格获得上述评价对象区域中各单元的日照量评价要素,其中该表格中,上述静止卫星图像数据中各像元的像元值与各像元对应区域假定为水平面时的每个区域的日照量评价要素相对应,之后,对该日照量评价要素按照上述校正系数做校正计算,并计算取得每个单元的地形条件校正后的评价要素,对评价对象区域的日照量进行评价。
8.一种日照量的评价方法,在该方法中:
在同时获得静止卫星图像数据和摄像日期时间数据后,
将上述静止卫星图像数据与三维地图模型结合,
接着,根据结合数据中静止卫星图像数据的地面分辨率,设定基于平面坐标的网格,按照该网格设定结合数据中的规定区域,并作为评价对象区域,
之后,根据上述三维地图模型和通过上述摄像日期时间数据指定的太阳位置,对在上述评价对象区域内产生的背阴区域进行运算,以规定的判断基准判断上述网格的各单元是否属于背阴区域,
同时,按照根据上述三维地图模型所规定的判断基准,求得评价对象区域中网格的各单元的代表倾斜角以及代表倾斜方位,根据上述太阳位置、将假定各单元为水平面时接受的日照量评价,变换为取得上述代表倾斜角以及代表倾斜方位时接受的日照量评价时针对每个单元计算校正系数,
并且,参照表格获得上述评价对象区域中各单元的日照量评价要素,其中该表格中,静止卫星图像数据中各像元的像元值与各像元对应区域假定为水平面时的每个区域的日照量评价要素相对应,之后,当该单元被判断为属于背阴区域时,将对上述日照量评价要素进行校正计算,同时,
取代上述日照量评价要素而获得地形条件校正后的评价要素,对评价对象区域的日照量进行评价,上述地形条件校正后的评价要素为根据上述校正系数做进一步校正计算而算出。
CN2009801293543A 2009-03-18 2009-06-09 日照量的评价方法和评价装置 Expired - Fee Related CN102105818B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2009-066770 2009-03-18
JP2009066770A JP4425983B1 (ja) 2009-03-18 2009-03-18 日射量の評価方法および評価装置
PCT/JP2009/002588 WO2010106582A1 (ja) 2009-03-18 2009-06-09 日射量の評価方法および評価装置

Publications (2)

Publication Number Publication Date
CN102105818A true CN102105818A (zh) 2011-06-22
CN102105818B CN102105818B (zh) 2013-06-05

Family

ID=42072707

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009801293543A Expired - Fee Related CN102105818B (zh) 2009-03-18 2009-06-09 日照量的评价方法和评价装置

Country Status (5)

Country Link
US (1) US8019117B2 (zh)
EP (1) EP2287638B1 (zh)
JP (1) JP4425983B1 (zh)
CN (1) CN102105818B (zh)
WO (1) WO2010106582A1 (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103136444A (zh) * 2013-01-29 2013-06-05 北京空间飞行器总体设计部 一种导航卫星不同控制模式下太阳帆板的受照分析方法
CN104580920A (zh) * 2013-10-21 2015-04-29 华为技术有限公司 一种成像处理的方法及用户终端
CN105911614A (zh) * 2016-07-04 2016-08-31 南通大学 一种基于数据缓存的日照时数计算方法
CN107439273A (zh) * 2013-05-13 2017-12-08 松下知识产权经营株式会社 温室用的日照量调整装置
CN110109196A (zh) * 2019-05-09 2019-08-09 洛阳众智软件科技股份有限公司 一种日照分析微区划分方法、装置及其设备
CN112106106A (zh) * 2018-05-09 2020-12-18 弗劳恩霍夫应用研究促进协会 关联卫星热图像数据以生成空间分辨率高的热量地图
TWI717796B (zh) * 2019-08-02 2021-02-01 慧景科技股份有限公司 以人工智慧估測地理位置日照量之系統、方法及儲存媒體
CN114501320A (zh) * 2021-12-29 2022-05-13 浙江百应科技有限公司 一种设备管控方法、装置以及系统

Families Citing this family (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011124287A (ja) * 2009-12-08 2011-06-23 Sony Corp 発電量予測装置、発電量予測システム、発電量予測方法及びコンピュータプログラム
JP5562895B2 (ja) * 2011-04-19 2014-07-30 中国電力株式会社 日射量推定装置、日射量推定方法及びプログラム
KR101890673B1 (ko) * 2011-09-21 2018-08-22 엘지전자 주식회사 일사량 추정을 위해 전천사진으로부터 운량을 계산하는 방법, 상기 계산한 운량을 이용한 태양광 발전량 예측 장치
US20130096832A1 (en) * 2011-10-18 2013-04-18 Dongaone Co., Ltd. Method and system for partitioning farming field
JP5840487B2 (ja) * 2011-12-28 2016-01-06 株式会社パスコ 日影領域検索装置、及び日影領域検索方法
JP2013242279A (ja) * 2012-05-23 2013-12-05 Eneres Corp 遮蔽係数算出装置、日射量予測装置、プログラム、および遮蔽係数算出方法
US20140039850A1 (en) * 2012-08-03 2014-02-06 Joon-Bum JEE Method and apparatus of forming solar radiation model applying topographical effects
US20140083413A1 (en) * 2012-09-24 2014-03-27 Brightsource Industries (Israel) Ltd. Method and apparatus for mapping cloud shading on the ground in a large area
US9091628B2 (en) 2012-12-21 2015-07-28 L-3 Communications Security And Detection Systems, Inc. 3D mapping with two orthogonal imaging views
US20160110663A1 (en) * 2014-10-17 2016-04-21 Ryan Craig Miller System and method for producing suitability score for energy management system on building rooftop
JP6325429B2 (ja) * 2014-12-16 2018-05-16 日本電信電話株式会社 環境変化検出システム、環境変化検出装置及び環境変化検出方法
US10082529B2 (en) 2015-02-16 2018-09-25 International Business Machines Corporation Estimating solar array power generation
JP6402932B2 (ja) * 2015-04-20 2018-10-10 清水建設株式会社 天候情報予測装置、天候情報予測方法、及びプログラム
ITUA20161393A1 (it) * 2016-03-07 2017-09-07 Flyby S R L Sistema di rilevazione della esposizione alla radiazione solare di un individuo.
WO2018148565A1 (en) * 2017-02-09 2018-08-16 Wove, Inc. Method for managing data, imaging, and information computing in smart devices
KR102404712B1 (ko) * 2017-11-30 2022-06-29 한국에너지기술연구원 인공신경망 기법을 이용한 일사량 산출을 위한 위성영상분석모델 생성 방법 및 시스템
ES2951153T3 (es) * 2019-01-11 2023-10-18 Cambras Gmbh Disposición de sensores y rastreo optimizado para sistemas CSP
KR102172967B1 (ko) * 2019-03-06 2020-11-02 한국에너지기술연구원 조도 정보를 이용한 위성 영상 보정 시스템 및 그 방법
JP7366556B2 (ja) * 2019-03-07 2023-10-23 Biprogy株式会社 太陽光線情報提供システム及び太陽光線情報提供プログラム
JP7407547B2 (ja) * 2019-09-25 2024-01-04 Biprogy株式会社 太陽光線情報提供システム、情報処理装置及び太陽光線情報提供プログラム
CN111597609B (zh) * 2020-04-28 2023-03-24 马辰龙 一种包含太阳辐射的基本运算单元及应用该单元的建筑能耗快速模拟方法
CN112258649A (zh) * 2020-10-27 2021-01-22 中国科学院空天信息创新研究院 一种基于光线投射的三维城市太阳辐射计算方法
CN112581595B (zh) * 2020-12-02 2023-12-19 中国人民解放军战略支援部队航天工程大学 一种多视角卫星影像一致性分析方法
KR102553456B1 (ko) * 2021-05-04 2023-07-11 한국에너지기술연구원 벡터 기반의 일사량 지도 생성 장치 및 그 방법
JP7129536B1 (ja) 2021-11-19 2022-09-01 株式会社パスコ 差分処理方法、差分処理装置及びプログラム
CN114257743A (zh) * 2021-12-16 2022-03-29 二十一世纪空间技术应用股份有限公司 调整相机成像参数方法、装置、终端及存储介质
JP7440686B1 (ja) 2023-03-29 2024-02-28 株式会社パスコ 差分処理方法、差分処理装置及びプログラム

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11211560A (ja) * 1998-01-21 1999-08-06 Fujitsu Ltd 衛星画像データを用いた日射量推定システム、方法及び記録媒体
JP2005051014A (ja) * 2003-07-28 2005-02-24 Atsushi Iga 太陽光発電量シミュレーション計算方法および計算プログラムを記録したコンピュータ読み取り可能なデータ記憶媒体
JP2006114838A (ja) * 2004-10-18 2006-04-27 Sharp Corp 太陽光発電設置診断システム、方法、プログラム
JP2006210750A (ja) * 2005-01-31 2006-08-10 Sanyo Electric Co Ltd 発電量予測システムおよび発電量予測方法
CN201087848Y (zh) * 2007-05-29 2008-07-16 南昌航空大学 建筑日照时数快速定量测量装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4202705B2 (ja) 2002-09-27 2008-12-24 株式会社パスコ レーザデータのノイズ除去方法
US7500391B2 (en) * 2006-04-27 2009-03-10 Ecometriks, Llc System and method for identifying the solar potential of rooftops

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11211560A (ja) * 1998-01-21 1999-08-06 Fujitsu Ltd 衛星画像データを用いた日射量推定システム、方法及び記録媒体
JP2005051014A (ja) * 2003-07-28 2005-02-24 Atsushi Iga 太陽光発電量シミュレーション計算方法および計算プログラムを記録したコンピュータ読み取り可能なデータ記憶媒体
JP2006114838A (ja) * 2004-10-18 2006-04-27 Sharp Corp 太陽光発電設置診断システム、方法、プログラム
JP2006210750A (ja) * 2005-01-31 2006-08-10 Sanyo Electric Co Ltd 発電量予測システムおよび発電量予測方法
CN201087848Y (zh) * 2007-05-29 2008-07-16 南昌航空大学 建筑日照时数快速定量测量装置

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103136444B (zh) * 2013-01-29 2015-10-21 北京空间飞行器总体设计部 一种导航卫星不同控制模式下太阳帆板的受照分析方法
CN103136444A (zh) * 2013-01-29 2013-06-05 北京空间飞行器总体设计部 一种导航卫星不同控制模式下太阳帆板的受照分析方法
CN107439273A (zh) * 2013-05-13 2017-12-08 松下知识产权经营株式会社 温室用的日照量调整装置
CN104580920B (zh) * 2013-10-21 2018-03-13 华为技术有限公司 一种成像处理的方法及用户终端
CN104580920A (zh) * 2013-10-21 2015-04-29 华为技术有限公司 一种成像处理的方法及用户终端
CN105911614B (zh) * 2016-07-04 2018-08-07 南通大学 一种基于数据缓存的日照时数计算方法
CN105911614A (zh) * 2016-07-04 2016-08-31 南通大学 一种基于数据缓存的日照时数计算方法
CN112106106A (zh) * 2018-05-09 2020-12-18 弗劳恩霍夫应用研究促进协会 关联卫星热图像数据以生成空间分辨率高的热量地图
CN112106106B (zh) * 2018-05-09 2023-11-07 弗劳恩霍夫应用研究促进协会 关联卫星热图像数据以生成空间分辨率高的热量地图
CN110109196A (zh) * 2019-05-09 2019-08-09 洛阳众智软件科技股份有限公司 一种日照分析微区划分方法、装置及其设备
CN110109196B (zh) * 2019-05-09 2021-07-30 洛阳众智软件科技股份有限公司 一种日照分析微区划分方法、装置及其设备
TWI717796B (zh) * 2019-08-02 2021-02-01 慧景科技股份有限公司 以人工智慧估測地理位置日照量之系統、方法及儲存媒體
CN114501320A (zh) * 2021-12-29 2022-05-13 浙江百应科技有限公司 一种设备管控方法、装置以及系统

Also Published As

Publication number Publication date
JP4425983B1 (ja) 2010-03-03
US20100310116A1 (en) 2010-12-09
EP2287638A1 (en) 2011-02-23
EP2287638A4 (en) 2011-10-12
JP2010217107A (ja) 2010-09-30
US8019117B2 (en) 2011-09-13
CN102105818B (zh) 2013-06-05
WO2010106582A1 (ja) 2010-09-23
EP2287638B1 (en) 2012-11-21

Similar Documents

Publication Publication Date Title
CN102105818B (zh) 日照量的评价方法和评价装置
Freitas et al. Modelling solar potential in the urban environment: State-of-the-art review
Matzarakis et al. Sky view factor as a parameter in applied climatology-rapid estimation by the SkyHelios model
CN107063197B (zh) 一种基于空间信息技术的水库特征曲线提取方法
Tovar-Pescador et al. On the use of the digital elevation model to estimate the solar radiation in areas of complex topography
CN112434617B (zh) 一种基于多源遥感数据的森林生物量变化监测方法及系统
US20050212794A1 (en) Method and apparatus for removing of shadows and shadings from texture images
US11983813B2 (en) Method and system for determining solar access of a structure
Quan et al. Probabilistic solar irradiance transposition models
CN109145494A (zh) 一种海面温度反演方法及系统
Lim et al. A land data assimilation system using the MODIS-derived land data and its application to numerical weather prediction in East Asia
Landier et al. Calibration of urban canopies albedo and 3D shortwave radiative budget using remote-sensing data and the DART model
CN114663786B (zh) 基于点云数据和计算机图形学的林分辐射通量计算方法
CN112595344A (zh) 一种遥感卫星高分相机在轨绝对辐射定标的方法及装置
Laue Using ray tracing to model agri-PV greenhouse energy production and PAR levels
Ilba Estimating the daily solar irradiation on building roofs and facades using Blender Cycles path tracing algorithm
Al Rammahi Specification of Actual Evapotranspiration in Iraq by GIS
Schneider et al. Normal water temperature and ice cover of the Laurentian Great Lakes a computer animation, data base, and analysis tool
Hurkmans Efficient solar potential estimation of 3D buildings: 3D BAG as use case
Smith et al. A system for estimating direct and diffuse photosynthetically active radiation from hemispherical photographs
Dobbs The Incorporation of atmospheric variability into DIRSIG
Koblet et al. Determination of length, area, and volume changes at Storglaciären, Sweden, from multi-temporal aerial images (1959–1999)
CN115982794A (zh) 一种利用山体等高线高效计算日照阴影区域的方法
Hussein et al. Quantifying the Reliability of Volumetric and Areal Calculation with UAV-Generated DEMs: A Comparative Study with Ground Truth Data
CN114692428A (zh) 一种基于云频率的光能资源分布估算方法和系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130605

Termination date: 20150609

EXPY Termination of patent right or utility model