CN116644497A - 基于数字表面模型的屋顶光伏日照阴影分析方法及系统 - Google Patents

基于数字表面模型的屋顶光伏日照阴影分析方法及系统 Download PDF

Info

Publication number
CN116644497A
CN116644497A CN202310606885.9A CN202310606885A CN116644497A CN 116644497 A CN116644497 A CN 116644497A CN 202310606885 A CN202310606885 A CN 202310606885A CN 116644497 A CN116644497 A CN 116644497A
Authority
CN
China
Prior art keywords
shadow
solar
grid
angle
sun
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
Application number
CN202310606885.9A
Other languages
English (en)
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.)
Shanghai Energy Technology Development Co ltd
Original Assignee
Shanghai Energy Technology Development Co ltd
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 Shanghai Energy Technology Development Co ltd filed Critical Shanghai Energy Technology Development Co ltd
Priority to CN202310606885.9A priority Critical patent/CN116644497A/zh
Publication of CN116644497A publication Critical patent/CN116644497A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Structural Engineering (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Remote Sensing (AREA)
  • Computational Mathematics (AREA)
  • Computer Graphics (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

本发明提供一种基于数字表面模型的屋顶光伏日照阴影分析方法及系统，基于数字表面模型的屋顶光伏日照阴影分析方法包括：确定日照阴影分析目标区域的太阳方位角和太阳高度角；在各时间点进行日照阴影分析，并得到预设时间的日照阴影分布栅格；在所述日照阴影分布栅格中筛选所述目标区域的可开发区域。本发明通过计算太阳方位角建立变换矩阵旋转DSM数据，计算太阳高度角并使用可调整计算步长的改进阴影堆积算法计算区域各时刻的阴影遮挡情况，将各时刻的阴影遮挡栅格数据采用累加的栅格镶嵌方式获得目标区域日照时长栅格，最后根据目标区域内各地点的日照市场筛选适用于屋顶光伏铺设的区域，为屋顶光伏资源评估提供参考。

Description

基于数字表面模型的屋顶光伏日照阴影分析方法及系统
技术领域
本发明涉及光伏技术领域，具体涉及一种基于数字表面模型的屋顶光伏日照阴影分析方法及系统。
背景技术
日照阴影是影响光伏排布设计及实际发电量的重要因素，对感兴趣区域日照阴影的模拟分析有其现实实用意义。目前大部分日照阴影分析工具都依托于一些CAD软件或三维建模软件，主要用于待建住宅小区的屋顶日照时长模拟。对于CAD，由于CAD相关工具受限于CAD软件本身数据格式的局限性需要对区域内各个建筑及其他设施逐个设置建筑高度，并通过手动设置参数进行日照分析，无法对建筑屋顶区域的日照阴影情况进行分析，并且忽视了植被的遮挡，无法应用于屋顶光伏资源评估场景。
而对于SketchUp等三维建模软件能够实现对建筑屋顶区域的日照模拟，但是需要事先完成区域内建筑设施及植被的三维建模工作，前期建模工作量大，后期分析耗时久，且分析结果只能以模型文件的格式输出，无法对日照阴影分析结果进行后处理及空间分析。
基于此，本申请发明人提出一种基于数字表面模型的屋顶光伏日照阴影分析方法及系统，以期解决上述技术问题。
发明内容
本发明要解决的技术问题是为了克服现有技术中CAD软件无法应用于屋顶光伏资源评估、三维建模软件无法对日照阴影分析结果进行后处理及空间分析的缺陷，提供一种基于数字表面模型的屋顶光伏日照阴影分析方法及系统。
本发明是通过下述技术方案来解决上述技术问题的：
本发明提供一种基于数字表面模型的屋顶光伏日照阴影分析方法，其特点在于，包括如下步骤：
S1：确定日照阴影分析目标区域的太阳方位角和太阳高度角；
S2：在各时间点进行日照阴影分析，并得到预设时间的日照阴影分布栅格；
S3：在所述日照阴影分布栅格中筛选所述目标区域的可开发区域。
根据本发明的一个实施例，所述步骤S1包括：
基于DSM数据的经纬度与时间计算各时间点太阳方位角和高度角；其中，所述DSM数据由无人机正射获取。
根据本发明的一个实施例，所述步骤S1还包括：
根据DSM数据的投影坐标系，建立栅格平面坐标和经纬度对应关系，通过计算栅格中心点坐标得到所述目标区域中心的经纬度。
根据本发明的一个实施例，所述步骤S1还包括：
根据需要分析的日期计算太阳赤纬角；根据所述目标区域中心的经度和时间计算真太阳时与时角；根据所述目标区域中心的纬度、所述太阳赤纬角与所述时角计算太阳高度角；根据所述太阳高度角、所述目标区域中心的纬度以及所述太阳赤纬角计算太阳方位角。
根据本发明的一个实施例，所述步骤S2包括：
S21：先对所述DSM数据预处理，并在所述太阳高度角非90度时对所述DSM数据进行栅格旋转，使得旋转后的所述DSM栅格行方向与太阳入射方向平行，并且太阳入射光线由栅格右侧射入；其中，旋转角度及方向计算公式如下：
Degrotate＝azimuth-90；
其中，Degrotate为栅格旋转角度，azimuth为太阳方位角；
S22：对完成旋转的所述DSM数据进行归一化处理，计算阴影最大影响范围作为日照阴影分析的边界条件，计算公式如下：
其中，Shadowmax为阴影最大影响范围，nDSMmax为归一化DSM最大值，scale为DSM数据分辨率，tan(altitude)为太阳高度角正切值；
S23：根据栅格的序列号生成阶梯式栅格矩阵，与归一化后的所述DSM数据叠加生成堆积栅格，计算方式如下：
其中，Shadowcast为堆积栅格，nDSM为归一化栅格，n为nDSM栅格矩阵最大列序号，tan(altitude)为太阳高度角正切值；
S24：以步长值作为偏移量比较所述堆积栅格列方向的不同切片，得到所述日照阴影分布栅格。
根据本发明的一个实施例，所述步长数列表达式如下：
step＝{2x}，x∈[0,log2Shadowmax+1)；
其中，step为步长；Shadowmax为阴影最大影响范围。
根据本发明的一个实施例，所述步骤S2还包括：
在得到所述日照阴影分布栅格后使用累加的方式对各时间点阴影分布栅格进行栅格镶嵌，获得全天日照阴影分布栅格。
本发明还提供一种基于数字表面模型的屋顶光伏日照阴影分析系统，所述基于数字表面模型的屋顶光伏日照阴影分析系统采用如上所述的基于数字表面模型的屋顶光伏日照阴影分析方法，其特点在于，所述分析系统包括：
计算模块，用于确定日照阴影分析目标区域的太阳方位角和太阳高度角；
分析模块，用于在各时间点进行日照阴影分析，并得到预设时间的日照阴影分布栅格；
输出模块，用于在所述日照阴影分布栅格中筛选所述目标区域的可开发区域。
本发明还提供一种电子设备，包括：处理器和存储器，所述存储器存储可在所述处理器上运行的程序或指令，所述程序或指令被所述处理器执行实现如上所述的基于数字表面模型的屋顶光伏日照阴影分析方法。
本发明还提供一种可读存储介质，所述可读存储介质上存储程序或指令，所述程序或指令被处理器执行时实现如上所述的基于数字表面模型的屋顶光伏日照阴影分析方法。
本发明的积极进步效果在于：
本发明提供的基于数字表面模型的屋顶光伏日照阴影分析方法及系统，通过计算太阳方位角建立变换矩阵旋转DSM数据，计算太阳高度角并使用可调整计算步长的改进阴影堆积算法计算区域各时刻的阴影遮挡情况，将各时刻的阴影遮挡栅格数据采用累加的栅格镶嵌方式获得目标区域日照时长栅格，最后根据目标区域内各地点的日照市场筛选适用于屋顶光伏铺设的区域，为屋顶光伏资源评估提供参考。
附图说明
本发明的上述的以及其他的特征、性质和优势将通过下面结合附图和实施例的描述而变得更加明显，其中：
图1为本发明基于数字表面模型的屋顶光伏日照阴影分析方法的流程图；
图2为本发明基于数字表面模型的屋顶光伏日照阴影分析方法的一实施方式中的DSM数据图；
图3为本发明基于数字表面模型的屋顶光伏日照阴影分析方法的一实施方式个的日照阴影分布栅格图；
图4为本发明电子设备的结构示意图。
具体实施方式
为让本发明的上述目的、特征和优点能更明显易懂，以下结合附图对本发明的具体实施方式作详细说明。
现在将详细参考附图描述本发明的实施例。现在将详细参考本发明的优选实施例，其示例在附图中示出。在任何可能的情况下，在所有附图中将使用相同的标记来表示相同或相似的部分。此外，尽管本发明中所使用的术语是从公知公用的术语中选择的，但是本发明说明书中所提及的一些术语可能是申请人按他或她的判断来选择的，其详细含义在本文的描述的相关部分中说明。此外，要求不仅仅通过所使用的实际术语，而是还要通过每个术语所蕴含的意义来理解本发明。
参照图1，本发明提供一种基于数字表面模型的屋顶光伏日照阴影分析方法，包括如下步骤：
S1：确定日照阴影分析目标区域的太阳方位角和太阳高度角。
在一个实施方式中，本发明基于DSM数据的经纬度与时间按各时间点太阳方位角和高度角。其中，DSM数据由无人机正射获取。
具体地，无人机可以实现半小时完成1.5km2数据采集的数字表面模型(digitalsurface model，DSM)，如此可以解决现有技术中CAD日照工具参数设置复杂、三维建模软件前期工作量大的问题。
根据DSM数据的投影坐标系，建立栅格平面坐标和经纬度对应关系，通过计算栅格中心点坐标得到目标区域中心的经纬度。经纬度包括目标区域的经度和目标区域的纬度。
在得到目标区域的经纬度后可以计算太阳方位角和太阳高度角。具体地，根据需要分析的日期计算太阳赤纬角；根据所述目标区域中心的经度和时间计算真太阳时与时角；根据所述目标区域中心的纬度、所述太阳赤纬角与所述时角计算太阳高度角；根据所述太阳高度角、所述目标区域中心的纬度以及所述太阳赤纬角计算太阳方位角。
进一步地，对于太阳赤纬角，太阳赤纬又称赤纬角，是地球赤道平面与太阳和地球中心的连线之间的夹角。赤纬角以年为周期，在+23"27'与-23"27"的范围内移动，成为季节的标志。每年6月21日或22日赤纬达到最大值+23°27称为夏至，该日中午太阳位于地球北回归线正上空，是北半球日照时间最长、南半球日照时间最短的一天。随后赤纬角逐渐减少至9月21日或22日等于零时全球的昼夜时间均相等为秋分。至12月21日或22日赤纬减至最小值-23"27为冬至，此时阳光斜射北半球，昼短夜长而南半球则相反。当赤纬角又回到零度时为春分即3月21日或22日，如此周而复始形成四季。因赤纬值日变化很小，一年内任何一天的赤纬角6可用下式计算：sino＝0.39795cos[0.98563(N-173)]。式中N为日数，自1月1日开始计算。如此，当确定需要分析的日期后，便可以根据上述公式计算得到太阳赤纬角。
对于真太阳时和时角，从天文学上来说，时间可分为平太阳时和真太阳时。平太阳时就是平时参考的国家地区报时时间，如北京时间，其每天的时间间隔是相等的。真太阳时是以当地太阳位于正南向的瞬时为正午。由于太阳与地球之间的距离和相对位置随时间在变化，以及地球赤道与其绕太阳运行的轨道的处平面的不一致，因而真太阳时与钟表指示的时间(平太阳时)之间总会有所差异，它们的差值即为时差。最大时差可达16min。一年中只有4次时差为零。计算太阳位置时应采用真太阳时。
在中国地区，真太阳时的换算公式为：
真太阳时(t)＝北京时间+时差；
时差＝(当地经度-120°)/15°；
对于时角，太阳时角是指日面中心的时角，即从观测点天球子午圈沿天赤道量至太阳所在时圈的角距离。
计算公式为：ω＝15×(ST-12)；
其中，ST为真太阳时，以24小时计。
如此，通过目标区域中心的精度和时间便可以计算得出真太阳时与时角。
对于太阳高度角，太阳高度角随着地方时和太阳的赤纬的变化而变化。太阳赤纬角(与太阳直射点纬度相等)以δ表示，观测地地理纬度用表示(太阳赤纬与地理纬度都是北纬为正，南纬为负)，地方时以t表示，有太阳高度角的计算公式：
其中，δ为太阳赤纬角；为目标区域的纬度；t为地方时。
根据上述的太阳高度角、目标区域中心的维度以及太阳赤纬角可以计算得到太阳方位角。
S2：在各时间点进行日照阴影分析，并得到预设时间的日照阴影分布栅格。
其中，DSM获取的数据图可以参照图2，日照阴影分布栅格图可以参照图3。
具体地，步骤S2包括：
S21：先对所述DSM数据预处理，并在所述太阳高度角非90度时对所述DSM数据进行栅格旋转，使得旋转后的所述DSM栅格行方向与太阳入射方向平行，并且太阳入射光线由栅格右侧射入；其中，旋转角度及方向计算公式如下：
Degrotate＝azimuth-90；
其中，Degrotate为栅格旋转角度，azimuth为太阳方位角；
S22：对完成旋转的所述DSM数据进行归一化处理，计算阴影最大影响范围作为日照阴影分析的边界条件，计算公式如下：
其中，Shadowmax为阴影最大影响范围，nDSMmax为归一化DSM最大值，scale为DSM数据分辨率，tan(altitude)为太阳高度角正切值；
S23：根据栅格的序列号生成阶梯式栅格矩阵，与归一化后的所述DSM数据叠加生成堆积栅格，计算方式如下：
其中，Shadowcast为堆积栅格，nDSM为归一化栅格，n为nDSM栅格矩阵最大列序号，tan(altitude)为太阳高度角正切值；
S24：以步长值作为偏移量比较所述堆积栅格列方向的不同切片，得到所述日照阴影分布栅格。
进一步地，通过设置日照阴影分析计算步长，由于被距离越远的栅格遮挡的可能性越小，因此使用等比例步长加快计算速度。其中，最大步长不超过Shadowmax。步长数列表达式如下：
step＝{2x}，x∈[0,log2Shadowmax+1)；
其中，step为步长；Shadowmax为阴影最大影响范围。
需要说明的是，如不需要考虑时长或者需求最高精度的日照阴影分析，可将步长数列表达式修改为：
step＝{x}，x∈[0,Shadowmax+1)；
具体地，以步长值作为偏移量比较堆积栅格列方向的不同切片。例如，当步长值为1时，比较堆积栅格第1列至倒数第2列的切片与第2列至最后1列的切片，将较大值作为堆积栅格该行列号的新值。当步长值为32时，比较堆积栅格第1列至倒数第32列的切片与第32列至最后1列的切片，较大值作为堆积栅格该行列号的新值。完成所有步长的堆积栅格切片比较后，将最终的堆积栅格结果与原始堆积栅格做比较，若某行列号的最终堆积栅格值大于原始堆积栅格值，则说明该位置被阴影遮挡，赋值为0，反之赋值为1。最终得到该时间点的日照阴影分布栅格，1表明被太阳直射，0表明被阴影遮挡。
进一步地，在得到日照阴影分布栅格后，使用累加的方式对各时间点阴影分布栅格进行栅格镶嵌，获得全天日照阴影分布栅格。
其中，栅格值为该时段各位置受太阳直射小时数。
S3：在所述日照阴影分布栅格中筛选所述目标区域的可开发区域。
例如，本发明使用冬至日上午7时至下午17时的时间段内的日照小时数大于等于8小时作为筛选屋顶光伏可开发区域的阈值。在获得的全天日照阴影分布栅格中栅格值大于等于8的赋值1，其他栅格赋值Nodata，最终生成屋顶光伏可开发区域栅格。
本发明提供的基于数字表面模型的屋顶光伏日照阴影分析方法，使用易于大范围获取的DSM数据作为本发明的唯一输入数据，无参数化全自动实现数据的预处理，避免了目前各日照阴影分析中前期需要人工参与进行建筑高度设置或建筑三维建模等需要耗费大量人工的弊端。
而且，可以根据实际需求设置计算步长，适用于多尺度多应用场景的日照阴影分析。当使用等比例步长进行日照阴影分析计算，并筛选屋顶光伏可开发面积时，比标准步长step＝{x}(其中x∈[0,shadowmax+1))计算快19倍(Intel i7 6700hq，16G内存)。计算面积1.5km2，分辨率0.5m的DSM数据耗时29.99s，效率远高于三维软件中的日照阴影分析模块。
并且，可以灵活调整日照时长预制，对最终获取的全天日照阴影分布栅格进行屋顶光伏可开发区域提取筛选，并生成可开发区域的栅格或矢量数据，便于使用其他CAD软件或者GIS软件进行后处理。
本发明还提供一种基于数字表面模型的屋顶光伏日照阴影分析系统，采用如上述的基于数字表面模型的屋顶光伏日照阴影分析方法，分析系统包括计算模块、分析模块以及输出模块，计算模块用于确定日照阴影分析目标区域的太阳方位角和太阳高度角；分析模块用于在各时间点进行日照阴影分析，并得到预设时间的日照阴影分布栅格；输出模块用于在所述日照阴影分布栅格中筛选所述目标区域的可开发区域。
获取模块通过计算太阳方位角建立变换矩阵旋转DSM数据，分析模块计算太阳高度角并使用可调整计算步长的改进阴影堆积算法计算区域各时刻的阴影遮挡情况，将各时刻的阴影遮挡栅格数据采用累加的栅格镶嵌方式获得目标区域日照时长栅格，最后输出模块根据目标区域内各地点的日照市场筛选适用于屋顶光伏铺设的区域，为屋顶光伏资源评估提供参考。
参照图4，本发明还提供一种电子设备900，包括：处理器901和存储器902，存储器902存储可在处理器901上运行的程序或指令，程序或指令被处理器901执行上述的基于数字表面模型的屋顶光伏日照阴影分析方法。该程序或者指令被处理器901执行时，实现上述基于数字表面模型的屋顶光伏日照阴影分析方法实施方式的各个过程，且能达到相同的技术效果，为避免重复，在此不做赘述。
本发明还提供一种可读存储介质，可读存储介质上存储程序或指令，程序或指令被处理器执行时实现如上述的基于数字表面模型的屋顶光伏日照阴影分析方法。该程序或者指令被处理器执行时，实现上述基于数字表面模型的屋顶光伏日照阴影分析方法实施方式的各个过程，且能达到相同的技术效果，为避免重复，在此不做赘述。
对于本领域技术人员来说，上述发明披露仅仅作为示例，而并不构成对本申请的限定。虽然此处并没有明确说明，本领域技术人员可能会对本申请进行各种修改、改进和修正。该类修改、改进和修正在本申请中被建议，所以该类修改、改进、修正仍属于本申请示范实施方式的精神和范围。
同时，本申请使用了特定词语来描述本申请的实施方式。如“一个实施方式”、“一实施方式”、和/或“一些实施方式”意指与本申请至少一个实施方式相关的某一特征、结构或特点。因此，应强调并注意的是，本说明书中在不同位置两次或多次提及的“一实施方式”或“一个实施方式”或“一替代性实施方式”并不一定是指同一实施方式。此外，本申请的一个或多个实施方式中的某些特征、结构或特点可以进行适当的组合。
本申请的一些方面可以完全由硬件执行、可以完全由软件(包括固件、常驻软件、微码等)执行、也可以由硬件和软件组合执行。以上硬件或软件均可被称为“数据块”、“模块”、“引擎”、“单元”、“组件”或“系统”。处理器可以是一个或多个专用集成电路(ASIC)、数字信号处理器(DSP)、数字信号处理器件(DAPD)、可编程逻辑器件(PLD)、现场可编程门阵列(FPGA)、处理器、控制器、微控制器、微处理器或者其组合。此外，本申请的各方面可能表现为位于一个或多个计算机可读介质中的计算机产品，该产品包括计算机可读程序编码。例如，计算机可读介质可包括，但不限于，磁性存储设备(例如，硬盘、软盘、磁带……)、光盘(例如，压缩盘CD、数字多功能盘DVD……)、智能卡以及闪存设备(例如，卡、棒、键驱动器……)。
计算机可读介质可能包含一个内含有计算机程序编码的传播数据信号，例如在基带上或作为载波的一部分。该传播信号可能有多种表现形式，包括电磁形式、光形式等等、或合适的组合形式。计算机可读介质可以是除计算机可读存储介质之外的任何计算机可读介质，该介质可以通过连接至一个指令执行系统、装置或设备以实现通讯、传播或传输供使用的程序。位于计算机可读介质上的程序编码可以通过任何合适的介质进行传播，包括无线电、电缆、光纤电缆、射频信号、或类似介质、或任何上述介质的组合。
同理，应当注意的是，为了简化本申请披露的表述，从而帮助对一个或多个发明实施方式的理解，前文对本申请实施方式的描述中，有时会将多种特征归并至一个实施方式、附图或对其的描述中。但是，这种披露方法并不意味着本申请对象所需要的特征比权利要求中提及的特征多。实际上，实施方式的特征要少于上述披露的单个实施方式的全部特征。一些实施方式中使用了描述成分、属性数量的数字，应当理解的是，此类用于实施方式描述的数字，在一些示例中使用了修饰词“大约”、“近似”或“大体上”来修饰。除非另外说明，“大约”、“近似”或“大体上”表明数字允许有±20％的变化。相应地，在一些实施方式中，说明书和权利要求中使用的数值参数均为近似值，该近似值根据个别实施方式所需特点可以发生改变。在一些实施方式中，数值参数应考虑规定的有效数位并采用一般位数保留的方法。尽管本申请一些实施方式中用于确认其范围广度的数值域和参数为近似值，在具体实施方式中，此类数值的设定在可行范围内尽可能精确。
本发明虽然以较佳实施方式公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以做出可能的变动和修改。因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施方式所作的任何修改、等同变化及修饰，均落入本发明权利要求所界定的保护范围之内。

Claims (10)

1.一种基于数字表面模型的屋顶光伏日照阴影分析方法，其特征在于，包括如下步骤：
S1：确定日照阴影分析目标区域的太阳方位角和太阳高度角；
S2：在各时间点进行日照阴影分析，并得到预设时间的日照阴影分布栅格；
S3：在所述日照阴影分布栅格中筛选所述目标区域的可开发区域。
2.根据权利要求1所述的基于数字表面模型的屋顶光伏日照阴影分析方法，其特征在于，所述步骤S1包括：
基于DSM数据的经纬度与时间计算各时间点太阳方位角和高度角；其中，所述DSM数据由无人机正射获取。
3.根据权利要求2所述的基于数字表面模型的屋顶光伏日照阴影分析方法，其特征在于，所述步骤S1还包括：
根据DSM数据的投影坐标系，建立栅格平面坐标和经纬度对应关系，通过计算栅格中心点坐标得到所述目标区域中心的经纬度。
4.根据权利要求3所述的基于数字表面模型的屋顶光伏日照阴影分析方法，其特征在于，所述步骤S1还包括：
根据需要分析的日期计算太阳赤纬角；根据所述目标区域中心的经度和时间计算真太阳时与时角；根据所述目标区域中心的纬度、所述太阳赤纬角与所述时角计算太阳高度角；根据所述太阳高度角、所述目标区域中心的纬度以及所述太阳赤纬角计算太阳方位角。
5.根据权利要求2所述的基于数字表面模型的屋顶光伏日照阴影分析方法，其特征在于，所述步骤S2包括：
S21：先对所述DSM数据预处理，并在所述太阳高度角非90度时对所述DSM数据进行栅格旋转，使得旋转后的所述DSM栅格行方向与太阳入射方向平行，并且太阳入射光线由栅格右侧射入；其中，旋转角度及方向计算公式如下：
Degrotate＝azimuth-90；
其中，Degrotate为栅格旋转角度，azimuth为太阳方位角；
S22：对完成旋转的所述DSM数据进行归一化处理，计算阴影最大影响范围作为日照阴影分析的边界条件，计算公式如下：
其中，Shadowmax为阴影最大影响范围，nDSMmax为归一化DSM最大值，scale为DSM数据分辨率，tan(altitude)为太阳高度角正切值；
S23：根据栅格的序列号生成阶梯式栅格矩阵，与归一化后的所述DSM数据叠加生成堆积栅格，计算方式如下：
其中，Shadowcast为堆积栅格，nDSM为归一化栅格，n为nDSM栅格矩阵最大列序号，tan(altitude)为太阳高度角正切值；
S24：以步长值作为偏移量比较所述堆积栅格列方向的不同切片，得到所述日照阴影分布栅格。
6.根据权利要求5所述的基于数字表面模型的屋顶光伏日照阴影分析方法，其特征在于，所述步长数列表达式如下：
step＝{2x}，x∈[0,log2Shadowmax+1)；
其中，step为步长；Shadowmax为阴影最大影响范围。
7.根据权利要求5所述的基于数字表面模型的屋顶光伏日照阴影分析方法，其特征在于，所述步骤S2还包括：
在得到所述日照阴影分布栅格后使用累加的方式对各时间点阴影分布栅格进行栅格镶嵌，获得全天日照阴影分布栅格。
8.一种基于数字表面模型的屋顶光伏日照阴影分析系统，所述基于数字表面模型的屋顶光伏日照阴影分析系统采用如权利要求1-7任一项所述的基于数字表面模型的屋顶光伏日照阴影分析方法，其特征在于，所述分析系统包括：
计算模块，用于确定日照阴影分析目标区域的太阳方位角和太阳高度角；
分析模块，用于在各时间点进行日照阴影分析，并得到预设时间的日照阴影分布栅格；
输出模块，用于在所述日照阴影分布栅格中筛选所述目标区域的可开发区域。
9.一种电子设备，其特征在于，包括：处理器和存储器，所述存储器存储可在所述处理器上运行的程序或指令，所述程序或指令被所述处理器执行实现如权利要求1-7任一项所述的基于数字表面模型的屋顶光伏日照阴影分析方法。
10.一种可读存储介质，其特征在于，所述可读存储介质上存储程序或指令，所述程序或指令被处理器执行时实现如权利要求1-7任一项所述的基于数字表面模型的屋顶光伏日照阴影分析方法。
CN202310606885.9A 2023-05-26 2023-05-26 基于数字表面模型的屋顶光伏日照阴影分析方法及系统 Pending CN116644497A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310606885.9A CN116644497A (zh) 2023-05-26 2023-05-26 基于数字表面模型的屋顶光伏日照阴影分析方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310606885.9A CN116644497A (zh) 2023-05-26 2023-05-26 基于数字表面模型的屋顶光伏日照阴影分析方法及系统

Publications (1)

Publication Number Publication Date
CN116644497A true CN116644497A (zh) 2023-08-25

Family

ID=87622431

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310606885.9A Pending CN116644497A (zh) 2023-05-26 2023-05-26 基于数字表面模型的屋顶光伏日照阴影分析方法及系统

Country Status (1)

Country Link
CN (1) CN116644497A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117195016A (zh) * 2023-11-07 2023-12-08 长江三峡集团实业发展(北京)有限公司 污水处理模式确定方法、装置、计算机设备及存储介质

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117195016A (zh) * 2023-11-07 2023-12-08 长江三峡集团实业发展(北京)有限公司 污水处理模式确定方法、装置、计算机设备及存储介质
CN117195016B (zh) * 2023-11-07 2024-02-06 长江三峡集团实业发展(北京)有限公司 污水处理模式确定方法、装置、计算机设备及存储介质

Similar Documents

Publication Publication Date Title
Walch et al. Big data mining for the estimation of hourly rooftop photovoltaic potential and its uncertainty
Nguyen et al. Incorporating shading losses in solar photovoltaic potential assessment at the municipal scale
Freitas et al. Modelling solar potential in the urban environment: State-of-the-art review
Alam et al. Shadow effect on photovoltaic potentiality analysis using 3D city models
CN103823927A (zh) 跟踪式光伏设备的阵列式排布方法
CN114881399B (zh) 一种基于gf7遥感影像的光伏发电潜力与经济性评估方法
CN109344491A (zh) 一种考虑天气状态变化和云层遮挡的太阳辐照度建模方法
CN116644497A (zh) 基于数字表面模型的屋顶光伏日照阴影分析方法及系统
CN106777634A (zh) 基于Arcgis的光伏阵列阴影计算方法
CN103778331B (zh) 一种建筑光伏系统中太阳能资源的计算方法
Díaz et al. An adaptive solar radiation numerical model
CN103985060A (zh) 一种基于测光数据的光伏电站发电量计算方法
CN116187807A (zh) 光伏发电量评估方法、装置、计算机设备、介质和产品
Cellura et al. A photographic method to estimate the shading effect of obstructions
Devetaković et al. Design of solar systems for buildings and use of BIM tools: Overview of relevant geometric aspects
CN110595477A (zh) 一种基于遗传算法的根据视频中太阳影子进行定位的方法
CN112085260B (zh) 一种基于云图预测的光伏功率超短期预测方法
CN111914396B (zh) 基于高分辨率dem数据的次网格地形三维地表太阳辐射强迫效应快速参数化方法
Abramczyk Parametric building forms rationalizing the incident direct solar irradiation
Alam et al. Resolution in photovoltaic potential computation
KR20220086189A (ko) 3차원 기반의 태양광 발전량 산출 장치 및 그 방법
CN116049946B (zh) 基于三维模型的建筑群朝向模拟方法、系统和存储介质
Braun A Scalable Approach for Spatio-Temporal Assessment of Photovoltaic Electricity Potentials for Building Facades of Entire Cities
CN112395764A (zh) 新能源光伏出力的新型模拟计算方法
Kaynak et al. GPU-Based Dynamic Solar Potential Estimation Tool Using 3D Plans

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