Ecological daily flow variation regulation and control method for river reach based on suitable habitat overlap of fishes

Abstract

The present invention provides an ecological daily flow variation regulation and control method for a river reach based on a suitable habitat overlap of fishes. The method comprises: drawing a distribution diagram of reference flows, regulated and controlled flow increments, and suitable overlaps of fishes; and obtaining a daily flow mean value of the current day as a reference flow, searching for the distribution diagram to obtain a regulated and controlled flow increment range of which a suitable overlap of fishes is greater than a preset value, and using the sum of the reference flow and the regulated and controlled flow increment range as a flow regulation and control range of the next day. According to the present invention, quantitative assessment of suitable overlaps of fishes in the spawning areas under different reference flows and different flow variations are carried out, and then the constraint value of the ecological flow variation of a river reach is obtained. The method is applicable to different research river reaches and fishes having different spawning properties, has good adaptability and high practicability, can effectively reduce the adverse effects of scheduling and operation of hydropower stations on aquatic ecology of river reaches and fish reproduction, and promotes continuous high-quality development of hydropower.

Description

A method for regulating ecological daily flow changes in river sections based on the overlap rate of suitable fish habitats Technical Field

[0001] The present invention belongs to the technical field of flow control, and in particular relates to a method for controlling ecological daily flow changes in river sections based on the overlap rate of suitable fish habitats. Background Art

[0002] Fish are widely used as indicator species for evaluating and protecting river ecosystems. Their reproductive behavior is significantly influenced by river flow processes. Hydropower station construction alters natural river flow processes, blocking fish migration pathways and disrupting fish habitats. Ecologically regulating hydropower stations to discharge appropriate ecological flows to meet the hydrodynamic conditions necessary for fish spawning is one effective way to maintain healthy river ecosystems.

[0003] Existing studies on suitable ecological flow during the fish breeding season mostly focus on qualitative research on suitable flow for fish breeding. There is a lack of quantitative research on the change in flow per unit time from the perspective of fish breeding needs, making it impossible to accurately describe the changes in ecological daily flow in river sections. Summary of the Invention

[0004] In view of the defects of the existing technology, the present invention provides a method for regulating the ecological daily flow changes of river sections based on the overlap rate of suitable fish habitats, which can effectively solve the above problems.

[0005] The technical solution adopted in the present invention is as follows:

[0006] The present invention provides a method for regulating the change of river section ecological daily flow based on the overlap rate of suitable fish habitats, comprising the following steps:

[0007] Step 1: Statistically analyze the historical flow of the study section during the fish breeding period to determine the upper limit of the baseline flow RF of the study section. max and the lower limit value of the reference flow RF min ;Preset reference flow interval △RF;

[0008] According to the flow control requirements, the upper limit value QF of the flow control increment is preset max And control the flow rate increment lower limit QF min ; Pre-set flow rate increment spacing △QF;

[0009] Step 2: At the upper limit of the reference flow rate RF max To the lower limit value of reference flow RF min Between, according to the reference flow interval △RF, traverse to obtain n reference flows, and represent each reference flow as reference flow RF i , where i = 1, 2, ..., n;

[0010] In regulating the flow increment upper limit QF max To the lower limit value QF of the control flow increment min Between, according to the control flow increment interval △QF, traverse to obtain m control flow increments, and express each control flow increment as control flow increment QF j , where j = 1, 2, …, m;

[0011] Step 3: For n base flows and m control flow increments, each base flow RF i And each control flow increment QF j The corresponding fish suitable overlap rate is calculated, and the distribution map of baseline flow-regulated flow increment-fish suitable overlap rate is drawn;

[0012] Step 4, obtain the average daily flow of the current day as the baseline flow, find the distribution map of the baseline flow-control flow increment-fish suitable overlap rate, obtain the control flow increment range where the fish suitable overlap rate is higher than the preset value, and use the baseline flow plus the control flow increment range as the flow control range for the next day's flow.

[0013] Preferably, in step 3, the following method is used to obtain the reference flow RF i and control flow increment QF j Corresponding fish suitable overlap rate:

[0014] Step 3.1: Determine multiple study flows based on the baseline flow and the range of the regulated flow increment for the study section. For each study flow, obtain the corresponding fish reproduction suitability cloud map for the study flow using the following method:

[0015] Step 3.1.1: Grid the study river section into multiple grids; each grid is represented as Grid k , using a two-dimensional hydrodynamic numerical model, each grid is obtained k water depth and flow velocity at the study flow rate;

[0016] Step 3.1.2, based on each grid k The water depth and flow velocity under the study flow are used to find the pre-established water depth suitability curve and flow velocity suitability curve for protecting fish in the study river section; the water depth suitability curve is the relationship curve between water depth and the water depth suitability index for fish reproduction; the flow velocity suitability curve is the relationship curve between flow velocity and the flow velocity suitability index for fish reproduction, and the data of each grid are obtained. k The water depth suitability index D for fish reproduction under the study flow k and fish reproduction flow suitability index V k ;

[0017] Step 3.1.3, use the following formula to get each grid k Fish reproduction suitability index CSF under the study flow k : CSF k =V k ×D k ×C k

[0018] Where: C k Represents the grid k The suitability index of the bottom sediment factor under the study flow ranges from 0 to 1. The suitability index of the bottom sediment factor is determined based on the bottom sediment conditions of the fish spawning grounds collected on site. When the bottom sediment conditions of the fish spawning grounds are met, C k When the value is 1, it does not meet the spawning substrate conditions. k The value is 0;

[0019] Step 3.1.4, use the following formula to get each grid Grid k Suitable fish habitat area (WUA) under the study flow rate k : WUA k =CSF k ×A k

[0020] Among them: A k Grid k Projected area on the horizontal plane;

[0021] Step 3.1.5, based on each grid k Suitable fish habitat area (WUA) under the study flow rate k , obtain the fish reproduction suitability cloud map under the research flow;

[0022] In step 3.2, use the following formula to obtain the controlled flow rate F: F=RF i +QF j

[0023] Step 3.3, respectively, the reference flow RF i With the regulated flow F as the research flow, find the fish reproduction suitability cloud map obtained in step 3.1 and obtain the reference flow RF i The corresponding first fish reproduction suitability cloud map and the second fish reproduction suitability cloud map corresponding to the regulated flow F;

[0024] Step 3.4, pre-set the fish reproduction suitability threshold, and calculate the fish suitability overlap rate of the first fish reproduction suitability cloud map and the second fish reproduction suitability cloud map when the fish reproduction suitability is higher than the fish reproduction suitability threshold.

[0025] Preferably, step 3.4 is specifically as follows:

[0026] The first fish breeding suitability cloud map is binarized, and areas where the fish breeding suitability is greater than 0.5 are marked in black, and other areas are marked in white, thereby obtaining the first fish breeding suitability cloud map after binarization processing;

[0027] The second fish breeding suitability cloud map is binarized, and areas with fish breeding suitability greater than 0.5 are marked in black, and other areas are marked in white, thereby obtaining the second fish breeding suitability cloud map after binarization.

[0028] The fish suitability overlap rate of the first fish reproduction suitability cloud map after binarization and the second fish reproduction suitability cloud map after binarization was obtained using the following formula:

[0029] in:

[0030] R overlap It is the fish suitability overlap rate of the first fish reproduction suitability cloud map after binarization processing and the second fish reproduction suitability cloud map after binarization processing;

[0031] A first is the suitable distribution area of ​​the first fish reproduction suitability cloud map after binarization processing;

[0032] A second It is the suitable distribution area of ​​the second fish reproduction suitability cloud map after binarization processing.

[0033] The method for regulating river section ecological daily flow changes based on the overlap rate of suitable fish habitats provided by the present invention has the following advantages:

[0034] The method provided by this paper regulates daily ecological flow variation in river sections based on the overlap ratio of suitable fish habitats. By quantitatively evaluating the overlap ratio of suitable fish spawning areas at different baseline flow rates and varying flow variations, it derives a constraint value for the ecological flow variation in the river section. This method is applicable to different study river sections and fish species with different spawning characteristics. It has good adaptability and practicality, effectively reducing the adverse effects of hydropower station operation on river aquatic ecology and fish reproduction, and promoting the sustainable and high-quality development of hydropower. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] FIG1 is a flow chart of a method for regulating daily flow changes in river sections based on the overlap rate of suitable fish habitats provided by the present invention;

[0036] FIG2 is a topographic map of the river section studied by the present invention;

[0037] FIG3 is a schematic diagram of a water depth suitability curve provided by the present invention;

[0038] FIG4 is a schematic diagram of a flow rate suitability curve provided by the present invention;

[0039] FIG5 is a cloud diagram of fish reproduction suitability under various research flow rates provided by the present invention;

[0040] FIG6 is a diagram showing the calculation process of the suitable overlap rate of fish provided by the present invention;

[0041] FIG7 is a schematic diagram of a distribution diagram of a baseline flow rate, a regulated flow rate increment, and a fish suitable overlap rate provided by the present invention. DETAILED DESCRIPTION

[0042] In order to make the technical problems, technical solutions and beneficial effects solved by the present invention more clearly understood, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0043] The technical problem to be solved by this invention is to overcome the shortcomings of existing technologies in quantifying flow rate changes per unit time during the fish breeding season. This method proposes a quantitative control method for ecological flow rate changes that takes into account the reproductive suitability of fish in the study section. By quantitatively evaluating the overlap rate of suitable fish spawning areas at different baseline flow rates and different flow rate changes, the constraint value of the ecological flow rate change in the river section is derived. This method is applicable to different study sections and fish with different spawning characteristics. It has good adaptability and strong practicality. It can effectively reduce the adverse effects of hydropower station scheduling and operation on the aquatic ecology and fish reproduction of the river section, and promote the sustainable and high-quality development of hydropower.

[0044] The present invention provides a method for regulating the ecological daily flow change of a river section based on the overlap rate of suitable fish habitats. Referring to FIG1 , the method comprises the following steps:

[0045] Step 1: Statistically analyze the historical flow of the study section during the fish breeding period to determine the upper limit of the baseline flow RF of the study section. max and the lower limit value of the reference flow RF min ;Preset reference flow interval △RF;

[0046] According to the flow control requirements, the upper limit value QF of the flow control increment is preset max And control the flow rate increment lower limit QF min ; Pre-set flow rate increment spacing △QF;

[0047] As an example, the fish breeding period is from April to July. Based on the historical flow monitoring data of the hydrological station in the studied river section, the distribution probability of the flow during the fish breeding period is calculated, and 10%-90% is taken as the upper and lower limits of the baseline flow in the studied river section, which are 200m 3 / s and 1400m 3 / s, per 2.5m 3 / s one level as the reference flow interval △RF, a total of 481 groups of reference flow. The flow increment range is -500m 3 / s~500m 3 / s, the flow rate increment interval △QF is set to 10m 3 / s.

[0048] Step 2: At the upper limit of the reference flow rate RF max To the lower limit value of reference flow RF min Between, according to the reference flow interval △RF, traverse to obtain n reference flows, and represent each reference flow as reference flow RF i , where i = 1, 2, ..., n;

[0049] In regulating the flow increment upper limit QF max To the lower limit value QF of the control flow increment min Between, according to the control flow increment interval △QF, traverse to obtain m control flow increments, and express each control flow increment as control flow increment QF j , where j = 1, 2, …, m;

[0050] Step 3: For n base flows and m control flow increments, each base flow RF i And each control flow increment QF j The corresponding fish suitable overlap rates were calculated, and a distribution map of baseline flow-regulated flow increment-fish suitable overlap rate was drawn.

[0051] As shown in Figure 7, this is a schematic diagram of the distribution map of baseline flow-regulated flow increment-fish suitability overlap rate. In practical applications, the fish suitability overlap rates of all possible situations within the range of 10%-90% of the water frequency in the study river section over the years can be obtained, and a schematic diagram of the distribution map of baseline flow-regulated flow increment-fish suitability overlap rate can be drawn.

[0052] In step 3, the following method is used to obtain the reference flow RF i and control flow increment QF j Corresponding fish suitable overlap rate:

[0053] Step 3.1: Determine multiple study flows based on the range of the baseline flow and the control flow increment of the study section; for each study flow, obtain the corresponding study flow fish reproduction suitability cloud map; for example, the baseline flow range is 200m 3 / s and 1400m 3 / s, the reference flow spacing △RF is 2.5m 3 / s; the flow rate increment range is -500m 3 / s~500m 3 / s, the flow rate increment interval △QF is set to 10m 3 / s. Therefore, when the reference flow rate is 200m 3 / s, the flow rate increment is -500m 3 / s, the research flow is 700m 3 / s; when the reference flow rate is 200m 3 / s, the flow rate increment is -490m 3 / s, the research flow is 690m 3 / s; in this way, each regulated flow increment and baseline flow are traversed to obtain multiple research flows.

[0054] Step 3.1 is as follows:

[0055] Step 3.1.1, as shown in Figure 2, is to create a topographic map of the study river section; grid the study river section into multiple grids; represent each grid as Grid k , using a two-dimensional hydrodynamic numerical model, each grid is obtained k water depth and flow velocity at the study flow rate;

[0056] Specifically, the topographic characteristics of the river section are collected and studied, and a two-dimensional hydrodynamic numerical model is established. Based on the two-dimensional hydrodynamic numerical model, the grid k The water depth and flow velocity under the study flow rate are as follows: The two-dimensional hydrodynamic numerical model can be the MIKE21 numerical model constructed using the shallow water equation.

[0057] Step 3.1.2, based on each grid k The water depth and flow velocity under the study flow are used to find the pre-established water depth suitability curve and flow velocity suitability curve for protecting fish in the study river section. The water depth suitability curve is a relationship curve between water depth and the water depth suitability index for fish reproduction, as shown in Figure 3, which is a schematic diagram of a water depth suitability curve. The flow velocity suitability curve is a relationship curve between flow velocity and the flow velocity suitability index for fish reproduction, as shown in Figure 4, which is a schematic diagram of a flow velocity suitability curve. The grid of each grid is obtained. k The water depth suitability index D for fish reproduction under the study flow k and fish reproduction flow suitability index V k ;

[0058] In this step, the water depth suitability curve and flow velocity suitability curve for protecting fish in the study river section are obtained by analyzing the fish habitat model established for the study river section.

[0059] Step 3.1.3, use the following formula to get each grid kFish reproduction suitability index CSF under the study flow k : CSF k =V k ×D k ×C k

[0060] Where: C k Represents the grid k The suitability index of the bottom sediment factor under the study flow ranges from 0 to 1. The suitability index of the bottom sediment factor is determined based on the bottom sediment conditions of the fish spawning grounds collected on site. When the bottom sediment conditions of the fish spawning grounds are met, C k When the value is 1, it does not meet the spawning substrate conditions. k The value is 0;

[0061] The substrate types required for spawning by broodstock of indigenous fish in the study river section are mainly pebbles (particle size > 64mm) and gravel (4mm < particle size < 64mm).

[0062] Therefore, in this invention, the fish breeding suitability is considered by considering three factors: flow velocity, water depth and bottom quality. The product method is used to calculate the flow rate of each grid. k Comprehensive suitability index CSF k .

[0063] Step 3.1.4, use the following formula to get each grid Grid k Suitable fish habitat area (WUA) under the study flow rate k : WUA k =CSF k ×A k

[0064] Among them: A k Grid k Projected area on the horizontal plane;

[0065] Step 3.1.5, based on each grid k Suitable fish habitat area (WUA) under the study flow rate k , a cloud map of fish reproduction suitability under the research flow rate is obtained; as shown in FIG5 , a cloud map of fish reproduction suitability under various research flow rates is obtained;

[0066] In step 3.2, use the following formula to obtain the controlled flow rate F: F=RF i +QF j

[0067] Step 3.3, respectively, the reference flow RF iWith the regulated flow F as the research flow, find the fish reproduction suitability cloud map obtained in step 3.1 and obtain the reference flow RF i The corresponding first fish reproduction suitability cloud map and the second fish reproduction suitability cloud map corresponding to the regulated flow F;

[0068] Step 3.4, pre-set the fish reproduction suitability threshold, and calculate the fish suitability overlap rate of the first fish reproduction suitability cloud map and the second fish reproduction suitability cloud map when the fish reproduction suitability is higher than the fish reproduction suitability threshold.

[0069] Refer to Figure 6, which shows the calculation process of the fish suitable overlap rate. Step 3.4 is as follows:

[0070] The first fish breeding suitability cloud map is binarized, and areas where the fish breeding suitability is greater than 0.5 are marked in black, and other areas are marked in white, thereby obtaining the first fish breeding suitability cloud map after binarization processing;

[0071] The second fish breeding suitability cloud map is binarized, and areas with fish breeding suitability greater than 0.5 are marked in black, and other areas are marked in white, thereby obtaining the second fish breeding suitability cloud map after binarization.

[0072] The fish suitability overlap rate of the first fish reproduction suitability cloud map after binarization and the second fish reproduction suitability cloud map after binarization was obtained using the following formula:

[0073] in:

[0074] R overlap It is the fish suitability overlap rate of the first fish reproduction suitability cloud map after binarization processing and the second fish reproduction suitability cloud map after binarization processing;

[0075] A first is the suitable distribution area of ​​the first fish reproduction suitability cloud map after binarization processing;

[0076] A second It is the suitable distribution area of ​​the second fish reproduction suitability cloud map after binarization processing.

[0077] Therefore, in the present invention, the fish suitable overlap rate adopts the reference flow RF i The overlapping area of ​​the suitability distribution area under the working condition and the suitability distribution area under the control flow F, divided by the reference flow RF i The suitability distribution area under the working condition, the obtained R overlap That is, the reference flow RF i As a reference to control the flow rate F-reference flow rate RF iThe difference between the two values ​​is taken as the fish suitability overlap rate of flow increment. overlap Also known as two flow states (reference flow RF i and the overlap ratio of spawning ground areas under regulated flow F).

[0078] Step 4, obtain the average daily flow of the current day as the baseline flow, find the distribution map of the baseline flow-control flow increment-fish suitable overlap rate, obtain the control flow increment range where the fish suitable overlap rate is higher than the preset value, and use the baseline flow plus the control flow increment range as the flow control range for the next day's flow.

[0079] For example, if the average daily flow rate is 1000m 3 / s, find the distribution diagram of baseline flow-regulated flow increment-fish suitable overlap rate shown in Figure 7, and obtain the range of regulated flow increment with a fish suitable overlap rate higher than 0.5, for example, 80-100m 3 / s, the flow control range of the next day's flow is 1080m 3 / s-1100m 3 / s.

[0080] The present invention provides a method for regulating the daily flow change of river sections based on the overlap rate of suitable fish habitats. For the first time, starting from the perspective of fish reproduction needs, by obtaining the overlap rate of suitable habitat areas under different flow conditions, the suitable incremental value of daily flow under different benchmark flows is quantitatively given as the flow regulation value for the next day. This method is applicable to different research river sections and fish with different spawning characteristics, and has good adaptability and strong practicality. This method fully considers the reproductive habits of fish, reduces the adverse effects of hydropower station scheduling and operation on the aquatic ecology of river sections and fish reproduction, and promotes the sustainable and high-quality development of hydropower.

[0081] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present invention. These improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A method for regulating the ecological daily flow variation of a river section based on the overlap rate of suitable habitats for fish, characterized in that, Including the following steps: Step 1: Conduct a statistical analysis of the historical flow during the fish spawning protection period in the study reach to determine the upper limit value RF of the reference flow and the lower limit value RF of the reference flow of the study reach max and the lower limit value RF of the reference flow min ; Preset the reference flow spacing △RF; Preset the upper limit value QF of the regulated flow rate increment and the lower limit value QF of the regulated flow rate increment according to the flow rate regulation requirements max and the lower limit value QF of the regulated flow rate increment min ; Preset the regulated flow increment spacing △QF; Step 2, between the reference flow rate upper limit value RF max and the reference flow rate lower limit value RF min traverse at the reference flow rate interval △RF to obtain n reference flow rates, and represent each reference flow rate as the reference flow rate RF i , where i = 1, 2, …, n; When adjusting the upper limit value QF of the flow rate increment max to the lower limit value QF of the flow rate increment min iterate through the flow rate increment spacing △QF, and obtain m flow rate increments for adjustment. Each flow rate increment for adjustment is represented as the flow rate increment QF j , where j = 1, 2, …, m; Step 3. For n reference flows and m regulation flow increments, for each reference flow RF i and each regulation flow increment QF j calculate their corresponding fish suitable overlap rates, and thus draw a reference flow - regulation flow increment - fish suitable overlap rate distribution map; Step 4: Obtain the daily flow mean of the current day as the reference flow, search the reference flow-regulated flow increment-fish suitable overlap rate distribution map, obtain the range of regulated flow increments with a fish suitable overlap rate higher than the preset value, and use the reference flow plus the range of regulated flow increments as the flow regulation range for the next day's flow.

2. The method for regulating the ecological daily flow variation of a river section based on the overlap rate of suitable habitats for fish according to claim 1, wherein In step 3, the following method is used to obtain the reference flow rate RF i and the regulated flow rate increment QF j The corresponding suitable overlap rate of fish: Step 3.1: Determine multiple research flows according to the value ranges of the reference flow and the regulated flow increment in the research reach; For each research flow, a corresponding fish reproduction suitability cloud map is obtained. The method is as follows: Step 3.1.1, grid the study reach into multiple grids; represent each grid as Grid k , use a two-dimensional hydrodynamic numerical model to obtain the water depth and flow velocity of each grid Grid k at the study discharge; Step 3.1.2, based on each grid Grid k At the water depth and flow velocity under the study flow rate, look up the pre-established water depth suitability curve and flow velocity suitability curve for fish protection in the study river reach; wherein, the water depth suitability curve is the relationship curve between the water depth and the fish reproduction water depth suitability index; the flow velocity suitability curve is the relationship curve between the flow velocity and the fish reproduction flow velocity suitability index, and obtain each grid Grid k The fish reproduction water depth suitability index D at the study flow rate k and the fish reproduction flow velocity suitability index V k ; Step 3.1.3, use the following formula to obtain each grid Grid k The fish reproduction suitability index CSF at the research flow rate k : CSF k = V k × D k × C k Where: C k represents the grid Grid k The substrate factor suitability index at the said research flow rate, with a value range of 0 to 1; the substrate factor suitability index is determined according to the substrate conditions of the fish spawning ground collected on-site. When it meets the substrate conditions for spawning, C k takes a value of 1. When it does not meet the substrate conditions for spawning, C k takes a value of 0; Step 3.1.4, use the following formula to obtain each grid Grid k The suitable habitat area of fish WUA at the research flow rate k : WUA k = CSF k × A k Where: A k is the projected area of the grid Grid k on the horizontal plane; Step 3.1.5, based on each grid Grid k the suitable habitat area WUA of fish under the research flow rate k , to obtain the suitability cloud map of fish reproduction under the research flow rate; Step 3.2: Use the following formula to obtain the regulated flow F: F = RF i + QF j Step 3.3, using the reference flow rate RF i and the regulated flow rate F as the research flow rates respectively, search the fish breeding suitability cloud diagram obtained in Step 3.1 to obtain the first fish breeding suitability cloud diagram corresponding to the reference flow rate RF i and the second fish breeding suitability cloud diagram corresponding to the regulated flow rate F; Step 3.4: Preset the fish reproduction suitability threshold, and calculate the fish suitable overlap rate of the first fish reproduction suitability cloud map and the second fish reproduction suitability cloud map when the fish reproduction suitability is higher than the fish reproduction suitability threshold.

3. The method for regulating the ecological daily flow variation of a river section based on the overlap rate of suitable habitats for fish according to claim 2, characterized in that, Step 3.4 is specifically: Perform binarization processing on the first fish reproduction suitability cloud map, mark the area where the fish reproduction suitability is greater than 0.5 as black, and other areas as white, to obtain the binarized first fish reproduction suitability cloud map; Perform binarization processing on the second fish reproduction suitability cloud map, mark the area where the fish reproduction suitability is greater than 0.5 as black, and other areas as white, to obtain the binarized second fish reproduction suitability cloud map; The fish suitability overlap rate of the first fish reproduction suitability cloud map after binarization processing and the second fish reproduction suitability cloud map after binarization processing is obtained by using the following formula: Wherein: R overlap is the fish suitability overlap rate of the first fish breeding suitability cloud map after binarization processing and the second fish breeding suitability cloud map after binarization processing; A first is the suitable distribution area of the first fish breeding suitability cloud map after binarization processing; A second is the suitable distribution area of the second fish breeding suitability cloud map after binarization processing.

Images