JP5021795B2 - Method and apparatus for calculating carbon dioxide absorption effect - Google Patents

Description

The present invention relates to a calculation method and a calculation apparatus for a carbon dioxide absorption effect.

  In recent years, attention has been focused on the ability of trees to absorb carbon dioxide in relation to the reduction of carbon dioxide, which is a major cause of global warming. For example, in the Kyoto Protocol, which is a protocol of recent international conferences, the amount of CO2 reduction is captured using the managed forest as an absorption source.

Patent Document 1 describes a technique for calculating the carbon dioxide absorption effect of such forests. This Te is conventional odor, made to understand the uptake by forest of carbon dioxide directly. Specifically, the amount of carbon absorption is calculated by multiplying the annual growth amount of the forest by the carbon content, and this is converted into carbon dioxide to grasp the annual amount of carbon dioxide absorption. The growth amount is calculated by acquiring the biomass accumulation amount obtained by multiplying the trunk tree volume by the biomass coefficient in two periods before and after a predetermined measurement period, and converting this into the year. For the above-mentioned stand trunk volume, select a standard area of a certain area for each forest phase, measure the tree height in this standard area, calculate the breast height diameter from this tree height, and calculate the tree height from these tree height and breast height diameter per unit area. It is calculated by calculating the trunk volume of trees and multiplying by the forest area measured based on aerial photographs.

JP 2008-46837 A

  However, the above-mentioned conventional example limits the carbon dioxide absorption source to forests in the first place, so that it is possible to conduct analysis suitable for carbon dioxide reduction in areas where there are unavoidable forests such as urban areas. Is not always easy to say.

The present invention has been made to solve the above-described drawbacks, and an object of the present invention is to provide a method for calculating a carbon dioxide absorption effect that enables analysis suitable for carbon dioxide reduction in an area with few forests. Another object of the present invention is to provide an apparatus for calculating a carbon dioxide absorption effect that enables an analysis suitable for the reduction of carbon dioxide in an area with few forests.

It can be said that the absorption effect of carbon dioxide is not only expected in the forest as described above, but it is also meaningful to recognize plants in general if photosynthesis is considered. Meanwhile, the region of calculating the effect, that is, when calculating the target area is wide, although the use of such conventional Similarly aerial photograph is efficient, in this way results in the target plants generally, the calculation material The problem is how to grasp. In this regard, it is desirable that certain consideration can be given to the stability in which the calculated material is not easily lost, as required by the above-mentioned Kyoto Protocol.

The present invention has been made in consideration of the above,
Computer
Extracting the plant region 4 from the calculation target region 2 according to the distribution of pixels in the observed image data 3 obtained by observing the calculation target region 2 from the sky in a predetermined wavelength range in which the plant group 1 can be identified;
Based on the numerical terrain model 5 and the numerical surface layer model 6 coordinated to the observed image data 3, the difference Δh between the ground surface elevation ht and the surface surface elevation hs of the plant region 4 is set, and the plant region 4 is set to a predetermined area. Calculating the height of the plant group 1 for each division unit 7 divided by the unit;
The height Δh per unit area of the plant group 1 and the carbon dioxide absorption by the plant group 1 are set in each of the plurality of height range sections 8 of the plant group 1 set within a range of height less than the forest equivalent. Based on the carbon dioxide absorption amount e assigned in advance according to the causal relationship with the amount e, the carbon dioxide absorption amount e of the height range section 8 corresponding to the height Δh of the plant group 1 is set as the area in the whole area of the plant region 4. Accordingly, the steps of calculating the carbon dioxide absorption effect in the non-forest region in the calculation target region 2 are executed in accordance with the respective steps to configure the carbon dioxide absorption effect calculation method.

Plants as a calculation material generally have a tendency to form a group such that weeds suitable for the local environment are clustered. By paying attention to such trends, it is possible to grasp using satellite photographs etc. by grasping the presence or absence of a group, that is, vegetation, etc. instead of individual units when plants are generally accepted as calculation materials Thus, the calculation for a wide area can be advanced efficiently. Moreover, it can also be said that it has the stability which is hard to lose | disappear easily compared with an individual | organism | solid by having comprised the group.

  In addition, crops as a general plant account for a certain proportion, but this can be expected to form a group because it follows economic activities.

Vegetation 1 such as more vegetation, the calculation target region 2 by using the observed image data 3 of satellite photographs or aerial photographs of observing the calculation target region 2 from the sky by the predetermined wavelength range that can determine plant group 1 It can be extracted efficiently. The predetermined wavelength range in which the plant group 1 can be identified is, for example, a visible wavelength range in which the plant can be identified by green, which is a general plant color, or a near-infrared wavelength in which the plant can be identified by reflecting the plant very strongly In addition to a single wavelength region such as a region, it can also be configured as a plurality of wavelength regions. For example, a normalized vegetation index (a typical vegetation index composed of a near infrared region and a visible region (red) ( It is also possible to determine the presence or absence of a plant using NDVI (Normalized Difference Vegetation Index).

In addition, the plant group 1 extracted in this way may contain forests, but the plant group 1 other than forests useful for analysis in an area where there are few forests is specified using the height. can do. Specifically, a digital terrain model 5 (DTM: Digital Terrain Model) and a numerical surface model 6 (DSM: Digital Surface Model) coordinated with respect to the observation image data 3 described above are prepared, and the elevation between both models is prepared. By calculating the difference Δh, the height of the plant group 1 can be obtained efficiently. If this height is narrowed down to less than that of the forest, the plant group 1 other than the forest can be specified. In determining the height, the plant group 1 is improved by performing each division unit 7 obtained by dividing the plant region 4 obtained by extracting the plurality of plant groups 1 from the calculation target region 2 by a predetermined area unit. Can be identified.

  That is, for the plant region 4 from which a plurality of plant groups 1 are extracted and larger than the unit area of the plant group 1, for example, a division unit 7 is set which is divided using the area indicated by the single pixel of the observation image data 3 described above as a unit. In such a case, the height of each plant group 1 can be determined. The area of the division unit 7 is preferably suitable for determining whether or not the plant group 1 is a forest according to its height. For example, the area is smaller than the ground area corresponding to the single pixel described above, Although it is possible to increase the size, it is desirable to consider the resolution, accuracy, etc. in the numerical terrain model 5 described above.

The forests excluded from the calculation material as described above generally have a better carbon dioxide absorption effect than the plant groups 1 other than the forest, and thus the plant group 1 other than the forest is identified and used as the calculation material. As a result, it is possible to prevent a reduction in calculation accuracy due to erroneously including forests. In addition, when it is desired to calculate the carbon dioxide absorption effect including the forest, it is only necessary to make a comprehensive judgment by combining conventional calculations for the forest in addition to the above calculation .

If the plant group 1 other than the forest can be identified as described above, the calculation can be performed according to the causal relationship with the carbon dioxide absorption. For the calculation , it is sufficient to determine the carbon dioxide absorption amount e per unit area by the plant group 1 other than the forest by conducting a field survey in advance, and multiply this by the area. Thus, calculation according to the area can be performed. The carbon dioxide absorption amount e can be configured as an appropriate index in addition to the annual expected absorption amount of carbon dioxide as in the conventional example.

Moreover, the area of the plant group 1 for performing the above calculations can be specified by forming the division unit 7 described above in a predetermined area unit. Further, the unit area for determining the carbon dioxide absorption amount e described above by field survey or the like is preferably suitable for the determination of the carbon dioxide absorption amount e, but is the same as the unit area of the division unit 7 described above. Then, the processing load can be reduced. In this case, further, if the area occupied by the single pixel of the observed image data 3 is the same, the processing load can be further reduced. That is, for example, the observation image data 3, the numerical terrain model 5, and the numerical surface layer model 6 described above are coordinated with each other, and the information included in these data is managed and processed on the basis of the plane position to obtain a geographic information system (GIS). : Geographic Information System), the carbon dioxide absorption effect can be calculated based on this GIS data.

Moreover, when determining the height of the vegetation group as described above, a plurality of height range sections 8 are set within a range of heights less than that of the forest, and carbon dioxide is added to each height range section 8. If the carbon dioxide absorption amount e per unit area set individually according to the causal relationship with the absorption amount is assigned, the calculation accuracy can be further increased by utilizing the height of the plant closely related to the carbon dioxide absorption amount. . In this case, the height range category 8 can be determined in consideration of improvement in calculation accuracy. However, if the height range category 8 corresponding to a so-called land use category such as grassland or farmland is set, as described above, GIS. When the land use classification is set as the attribute data in the case of building, it can be confirmed and corrected by using it as it is. For example, when a field confirmation survey is performed in order to set a land use classification as attribute data, the calculation accuracy can be greatly increased.

Furthermore, although the case where the carbon dioxide absorption amount e is assigned to the plant group 1 has been described above, instead of this, a conversion rate for converting the carbon dioxide absorption amount to the forest is assigned, and the conversion rate is within the calculation target region 2. It is also possible to calculate the carbon dioxide absorption effect by the area equivalent to the forest obtained by multiplying the area of the division unit 7 having the plant group 1 in this case. In this case, the calculation result is compared with the carbon dioxide absorption capacity of the forest standard This makes it easier to achieve consistency with forest standards such as the aforementioned protocol.

In addition, considering the above-mentioned GIS, the calculation method of the carbon dioxide absorption effect is
Computer
After setting the calculation target area 2,
With reference to the observation image data storage unit 10 that stores the observation image data 3 that is coordinated to the map data 9 and observed from above the ground surface in a predetermined wavelength range that the plant group 1 can discriminate, on the map of the calculation target region 2 The plant region 4 is extracted from the calculation target region 2 by specifying the position on the map according to the pixel distribution in the observed image data 3 for the corresponding region of
Next, refer to the numerical terrain model storage unit 11 that stores the numerical terrain model 5 coordinated in the map data 9 and the numerical surface layer model storage unit 12 that stores the numerical surface model 6 coordinated in the map data 9. The difference Δh between the ground surface elevation ht and the surface elevation hs of the plant region 4 is calculated as the height of the plant group 1 for each division unit 7 obtained by dividing the plant region 4 by a predetermined area unit,
Thereafter, the height Δh per predetermined unit area of the plant group 1 and the plant group 1 are set in each of the plurality of height range sections 8 of the plant group 1 set within the range of the height less than the forest equivalent. Based on the carbon dioxide absorption amount e allocated in advance according to the causal relationship with the carbon dioxide absorption amount e, the carbon dioxide absorption amount e of the height range section 8 corresponding to the height Δh of the plant group 1 Can be configured by executing each process of calculating the carbon dioxide absorption effect in the non-forest region in the calculation target region 2 by accumulating according to the area, and in this case, the degree of freedom in selecting the calculation target region 2 Data processing with improved performance is established.

In addition, the calculation of the carbon dioxide absorption effect described above is
A map data storage unit 13 for storing the map data 9,
An observation image data storage unit 10 for storing observation image data 3 obtained by observing the ground surface from the sky in a predetermined wavelength range that is coordinated to the map data 9 and that the plant group 1 can distinguish;
A numerical terrain model storage unit 11 for storing the numerical terrain model 5 of the calculation target region 2 coordinated in the map data 9;
A numerical surface layer model storage unit 12 for storing the numerical surface layer model 6 of the calculation target region 2 coordinated in the map data 9;
The carbon dioxide absorption amount e per predetermined unit area of the plant group 1, on the basis of a causal relationship between the height of the plant group 1, the height range of the constant Mori Te Lin considerable range smell less than the height Tokoro Plant / carbon dioxide absorption table 14 stored separately for each category 8,
Target area setting means 15 for setting the calculation target area 2;
With reference to the map data storage unit 13 and the observation image data storage unit 10, the plant region 4 is calculated from the calculation target region 2 according to the distribution of pixels in the observation image data 3 for the corresponding region on the map of the calculation evaluation target region 2. Plant region extraction means 16 for extraction,
With reference to the numerical terrain model storage unit 11 and the numerical surface layer model storage unit 12, the difference Δh between the ground surface elevation ht and the surface elevation hs of the plant region 4 is divided into a predetermined area unit in the plant region 4 Plant height calculating means 17 for calculating the height of the plant group 1 for each divided unit 7,
With reference to the plant / carbon dioxide absorption table 14, it is determined to which height range section 8 the height Δh of the plant group 1 corresponds, and carbon dioxide per predetermined unit area of each divided unit 7 determining the absorption e, and effect calculating means 18 for calculating a carbon dioxide absorption effect in non-forest area in the calculation target region 2 the carbon dioxide absorption amount e by integrating according to the area in the entire plant area 4 This can also be realized by using a carbon dioxide absorption effect calculation device configured to include

  As is apparent from the above description, according to the present invention, analysis suitable for the reduction of carbon dioxide in an area where there are few forests is possible, which can contribute to global warming countermeasures.

BRIEF DESCRIPTION OF THE DRAWINGS It is a figure which shows this invention, (a) is a flowchart which shows the flow of the whole process, (b) is a content explanatory view of a plant and carbon dioxide absorption amount table. It is a block diagram which shows this invention. BRIEF DESCRIPTION OF THE DRAWINGS It is a figure explaining the structure of the geographic information system which concerns on this invention, (a) is a figure which shows the image which superimposed several information on the basis of the geographical position, (b) is a numerical terrain model and a numerical surface layer model It is a figure explaining the difference in the adoption standard of an altitude, etc. It is a figure which shows a calculation object area | region, (a) is a figure which shows the image specified in the wider area | region, (b) is a figure which shows the image on map data, (c) is the figure on observation image data It is a figure which shows an image. It is a figure which shows a calculation object area | region, (a) is a figure which shows the image which set the mesh, (b) is a figure which shows the image on a numerical terrain model, (c) is a figure which shows the image on a numerical surface layer model (D) is a figure which shows the image in the state which computed the height of the plant group. It is a flowchart of an effect calculation process.

This embodiment shows the case where the annual carbon dioxide absorption amount by the plant excluding the forest in the calculation target region 2 is directly obtained, and is shown in FIG. 1A by a computer configured by the block diagram shown in FIG. Processing is performed along the flow. As shown in FIG. 2, GIS data 20 is stored in the computer prior to processing, and this GIS data 20 is based on the plane position specified on the two-dimensional map data 9, and the observed image data 3, It is configured by managing the numerical terrain model 5 and the numerical surface layer model 6. FIG. 3A shows an overlay image of these data in GIS20.

The observation image data 3 is obtained by processing and analyzing an observation image obtained by a multispectral sensor mounted on an artificial satellite using a normalized vegetation index (NDVI). This NDVI makes it possible to grasp the distribution of the plant group 1 in the imaging region in units of pixels by using the reflection characteristics of the plant in a predetermined wavelength range. The normalized vegetation index is specifically calculated by NDVI = (NIR-R) / (NIR + R), where NIR is the reflectance in the near infrared wavelength region, and R is the reflectance in the visible red wavelength region. It is. The resolution of the observation image can be determined in consideration of the size of the calculation target region 2, the data processing efficiency, and the desired size of the plant group 1. The plant group 1 is grasped by using various vegetation indices other than the normalized vegetation index, using a hyperspectral sensor instead of the multispectral sensor, and using an aerial photographed image instead of a satellite photographed image. It is enough.

  The numerical terrain model 5 is generated by stereo matching of the above-described satellite-captured images. In particular, the ground surface of an area where plants are settled can be easily obtained with the photographing time around winter. Similarly, the numerical surface layer model 6 is also generated by stereo matching of satellite images, and in this case, it is set to a photographing time around summer when the plant is growing. Further, these models 5 and 6 are generated with high accuracy by complementing them using a plurality of shooting periods. Note that these models 5 and 6 are specifically composed of, for example, a triangulated irregular network (TIN), a voxel (volume cell) model, or the like, and are not a stereo matching of the above-described satellite-captured images. It may be generated by using a distance measuring device or using an aircraft or the like as a photographing platform.

  The observed image data 3 described above is aligned with the map data 9 described above by geometrically correcting the coastline and the like by using annotation data such as the latitude and longitude of the imaging region. Further, the numerical terrain model 5 and the numerical surface layer model 6 are photographed so that the ground reference point whose position coordinates are specified are included in the photographed image, and geometric correction is performed, for example, by using the ground reference point. To be combined. As shown in FIG. 2, the GIS data 20 described above includes map data 9 in the map data storage unit 13, observation image data 3 in the observation image data storage unit 10, and numerical terrain model 5 in the numerical terrain model storage unit 11. The numerical surface layer model 6 is configured to be stored in the numerical surface layer model storage unit 12.

The calculation of the carbon dioxide absorption amount using the above GIS data 20 is first advanced by setting the calculation target region 2 as the calculation target (step S1). FIG. 4A shows how the calculation target area 2 is set. In this embodiment, for example, the calculation target area 2 is specified by operating the mouse (not shown) and specifying the area on the map data 9. The area specification is specified from the mouse or the like via the input unit 21. The target area setting means 15 received in this way sets the area according to the area designation. In FIG. 4A, reference numeral 22 denotes an administrative division included in the map data 9 described above, and the setting of the calculation target region 2 is performed by specifying an administrative division, for example, in addition to the above-described region designation using a mouse or the like. You may do it.

In the example of this embodiment, the calculation target area 2 includes a forest 23, shrub land 24, grassland 25, farmland 26, settlement 27, road 28, and bare land 29 as shown in FIG. Configured.

When the calculation target region 2 is set as described above, next, the plant region 4 in which the plant group 1 exists is extracted (step S2). The extraction of the plant region 4 is performed on the map data 9 aligned with the observation image data 3 with respect to the position coordinates of the region corresponding to all the pixels indicating the presence of the plant group 1 in the observation image data 3 described above. It is done in search. The plant region extraction means 16 for extracting the plant region 4 reads the observation image data 3 of the region corresponding to the calculation target region 2 whose plane coordinates are specified by the map data 9 as described above, and distributes it on the observation image data 3. An area occupied by all of the pixels to be performed is identified as a plant area 4 using the plane coordinates of the pixels. FIG. 4C shows the observation image data 3 related to the calculation target region 2 in which the pixel region indicating the existence of the plant group 1 is hatched. As shown in this figure, regions corresponding to the forest 23, shrubland 24, grassland 25, and farmland 26 are all extracted as plant regions 4.

  Next, the height of the plant group 1 in the plant region 4 is calculated (step S3). This height is calculated using the numerical terrain model 5 and the numerical surface layer model 6 described above. Prior to this, the mesh 30 is first set to these models in order to ensure good accuracy in calculating the height. Is done. This mesh 30 corresponds to the fact that the plant group 1 is grasped with the resolution of the observation image data 3 in pixel units as described above, and enables the height of the plant group 1 to be calculated according to this resolution. Further, the cell 7 (division unit) position of the mesh 30 is set to correspond to the pixel position of the observation image data 3 with the same mesh size as the pixel size of the observation image data 3.

FIG. 5A shows the observed image data 3 expressed so that each pixel can be identified, and FIGS. 5B and 5C show images obtained by setting the mesh 30 to each of the numerical terrain model 5 and the numerical surface model 6. Indicates. In FIGS. 5B and 5C, the same elevation is shown with the same hatching. As shown in this figure, in this embodiment, the mesh 30 is set only in the plant region 4. However, as long as the plant region 4 is set, the mesh 30 may be set for the entire calculation target region 2.

  Once the mesh 30 is set, next, a representative value in each cell of the mesh 30 is set. As described above, the numerical terrain model 5 or the like composed of TIN or the like is configured by complementing altitude measurement points with a triangular surface or the like. For this reason, the setting of the representative value is for specifying the altitude coordinate in the above-described unit of the cell 7, and specifically, for example, the altitude coordinate of the center point of the cell is the altitude coordinate of the cell 7, that is, the representative value. Can be adopted as. In addition, the representative value can be appropriately determined in consideration of the accuracy of altitude and the amount of calculation such as the average value of the altitude coordinates in the cell 7 and the median.

  When the representative values are set as described above, the altitude coordinate values of the numerical terrain model 5 and the numerical surface layer model 6 are compared in units of cells 7. As shown in FIG. 3 (b), the numerical terrain model 5 indicates the altitude ht of the ground surface excluding the land cover 31 and the feature 32, and the numerical surface model 6 indicates the surface layer surface including the land cover 31 and the feature 32. Since it has the altitude hs, the height of the desired plant group 1 is given by calculating the difference Δh. If this difference Δh is calculated for all the cells 7 in the plant region 4, the height can be obtained for each plant group 1 constituting each cell 7.

  The above-described height calculation processing of the plant group 1 is performed by the plant height calculation means 17. As shown in FIG. 2, the plant height calculation unit 17 includes a mesh setting unit 33, a representative value setting unit 34, and an altitude difference calculation unit 35. As described above, the mesh setting unit 33 sets the mesh 30 in the numerical terrain model 5 or the like, and is based on, for example, the resolution of the observation image data 3 or the overlay position of the observation image data 3 and the numerical terrain model 5 or the like. Etc. to determine the mesh size and the position of the mesh 30 on the plane coordinates. Further, the representative value setting unit 34 obtains, for example, the plane coordinates of the center point of the cell 7, acquires the altitude value taking the plane coordinates in the numerical terrain model 5, and the like, and obtains the elevation values of the numerical terrain model 5 and the numerical surface model 6. It repeats until an altitude value is acquired in all the cells 7 in the plant region 4 in each. The altitude difference calculation unit 35 subtracts the altitude value of the numerical landform model 5 from the altitude value of the numerical surface layer model 6 for each cell 7 indicating the same point in the numerical landform model 5 and the numerical surface layer model 6, Iterate until it is calculated in all of the cells 7.

After calculating the plant height in this way, finally, the carbon dioxide absorption effect of the plant excluding the forest is calculated (step S4). This calculation is performed by the effect calculation means 18 referring to the plant / carbon dioxide absorption amount table 14 as shown in FIG. As shown in FIG. 1B, the table 14 stores the plant height range section 8 and the annual carbon dioxide absorption amount e per unit area in association with each other. In this embodiment, 1 square meter is set as the unit area. However, it is sufficient to easily calculate the carbon dioxide absorption amount, which will be described later, using the area corresponding to the above-described cell 7 as the unit area.

Further, in this embodiment, the plant height range section 8 includes a grassy-equivalent height range, a farmland-equivalent height range, and a shrub, where carbon dioxide absorption tends to be approximated regardless of the type of plant. Three types of height ranges equivalent to the ground are set, and the forest 23 (or forest land, forest land) is excluded from the calculation target by not setting the height range equivalent to the forest. Each carbon dioxide absorption amount e corresponding to the height range section 8 can be made constant by a predetermined numerical value as shown by the algebra of A, B, and C in FIG. It is desirable to set in accordance with the locality after investigating the types of plants inhabiting the calculation target region 2 in advance. Further, as described above, when the plant height range category 8 is set according to the land use category, for example, when the attribute data of the land use category is set in the GIS data 20 described above, the category is classified using this. Can also be verified. Furthermore, the plant height of the threshold is shown in FIG. 1 (b) for dividing the above height range segment 8 together, on but those determined by experiment, for example, the calculation target region and field survey, actually local It can be set considering the height of the plants that are prosperous.

The calculation procedure will be described with reference to FIG. For the calculation , first, as described above, the plant height of each cell 7 calculated by the plant height calculating means 17 is obtained for each cell 7 (step S4-1), and this exceeds the shrub land equivalent. (Step S4-2) and, if applicable, the forest 7 is not subject to calculation , the cell 7 has a carbon dioxide absorption amount e per unit area. The attribute is registered and set as 0 (step S4-3). On the other hand, if it does not correspond to the above step S4-2, it is similarly determined whether it is the height equivalent to the shrub land (step S4-4) or the height equivalent to the farmland (step S4-6), respectively. If so, the attribute C for the carbon dioxide absorption amount e per unit area corresponding to the shrub land 24 is registered in the cell 7 (step S4-5), or the carbon dioxide absorption amount per unit area corresponding to the farmland 26 Register attribute B for e (step S4-7). If it does not correspond to any of them, since it can be determined that it is grassland 25, the attribute A for the carbon dioxide absorption amount e per unit area corresponding to this is registered (step S4-8). The above processing is repeated until all the cells 7 of the plant area 4 is calculated (step S4-9), multiplied by the carbon dioxide absorption amount e per unit area relative to the area on the ground in all calculation Once you cell 7 Then, after obtaining the carbon dioxide absorption amount of each cell 7, the carbon dioxide absorption amount of the plant region 4 obtained by adding the carbon dioxide absorption amounts of all the cells 7 is calculated (step S4-10).

The calculation result can be displayed on a monitor (not shown) via the output unit 36 shown in FIG. 2, thereby completing the calculation process of the carbon dioxide absorption effect.

FIG. 5D shows the calculation target region 2 by applying different hatching to the cell 7 for each height range section 8 of the plant, in other words, for each difference in the carbon dioxide absorption amount e per unit area. An alternate long and two short dashes line displays the map data 9 in an overlapping manner. As shown in this figure, according to the above processing, the difference in height and distribution of the plant group 1 can be analyzed, and the difference in carbon dioxide absorption amount e per unit area and the distribution can be analyzed. . It is also desirable to output such an image to the output unit 36 in addition to the carbon dioxide absorption amount described above so that the regional distribution of the carbon dioxide absorption capacity can be easily grasped in the calculation result.

In the above, the case where the plant height is determined using the numerical terrain model 5 and the numerical surface model 6 has been described. For example, the existing land cover classification determination method using the normalized vegetation index described above, It is also possible to improve the calculation accuracy by confirming the plant height specified using the numerical terrain model 5 etc., and making corrections as necessary, using field surveys etc. .

The modification of this invention is shown below. In this modification, the carbon dioxide absorption effect obtained by excluding the forest in the calculation target region 2 as described above can be easily calculated in terms of forest. Specifically, the carbon dioxide absorption amount e per unit area per year for the forest is set in advance in the same manner as the shrub land described above, and this is the carbon dioxide absorption amount in the shrub land, farmland, and grassland described above. By dividing each by e, the area corresponding to the forest per unit area such as shrub land is calculated. Using each area calculated in this manner as a conversion factor, the plant / carbon dioxide absorption amount table 14 described above has a one-to-one correspondence for each plant height range section 8 instead of the carbon dioxide absorption amount e per unit area. Then, the area equivalent to the forest in the plant other than the forest in the calculation target region 2 can be obtained by the effect calculation means 18.

Therefore, in this modified example, an area equivalent to the forest can be obtained, and this calculation method was followed when a predetermined calculation method using the area of the forest was established in calculating the carbon dioxide absorption effect. Calculation can be easily obtained.

  In this modification, the case where the area of the shrub land or the like is converted to the area of the forest is shown, but according to the determination factor of the carbon dioxide absorption effect, for example, the volume of the trunk shown in the above-described conventional example, etc. Such other elements may be included in the conversion. In this case, for example, the volume of the trunk of the shrub may be set in advance, and the conversion rate may be set in consideration of this.

DESCRIPTION OF SYMBOLS 1 Vegetation group 2 Calculation object area 3 Observation image data 4 Plant area 5 Numerical landform model 6 Numerical surface layer model 7 Division unit 8 Height range classification 9 Map data 10 Observation image data storage part 11 Numerical landform model storage part 12 Numerical surface layer model storage Section 13 Map data storage section 14 Plant / carbon dioxide absorption table 15 Target area setting means 16 Plant area extraction means 17 Plant height calculation means 18 Effect calculation means ht Ground surface elevation hs Surface elevation altitude Δh Difference e Dioxide per unit area Carbon absorption

Claims (4)

Computer
Extracting the plant region from the calculation target region according to the distribution of pixels in the observed image data obtained by observing the calculation target region from the sky with a predetermined wavelength range capable of distinguishing the plant group;
The difference between the surface elevation of the plant region and the surface elevation based on the numerical terrain model and the numerical surface layer model coordinated to the observed image data is obtained for each division unit obtained by dividing the plant region by a predetermined area unit. Calculating the height of the plant group;
In accordance with the causal relationship between the height per unit area of the plant group and the amount of carbon dioxide absorbed by the plant group, each of the plurality of height range divisions of the plant group set within the height range that is less than forest equivalent Based on the pre-assigned carbon dioxide absorption amount, the carbon dioxide absorption amount of the height range according to the height of the plant group is integrated according to the area over the entire plant region, and the non-forest region in the calculation target region the method of calculating the carbon dioxide absorption effect of performing each step of the step of calculating the carbon dioxide absorption effect in. Instead of assigning a carbon dioxide absorption amount to each of the height range sections of the plant group, a conversion rate for converting the carbon dioxide absorption amount to a forest is assigned, and a forest equivalent area obtained by multiplying the conversion rate by the area is calculated. The method for calculating the carbon dioxide absorption effect according to claim 1. Computer
After setting the calculation target area,
With reference to an observation image data storage unit that stores observation image data obtained by observing the ground surface from the sky with a predetermined wavelength range that is coordinated to map data and distinguishable by a group of plants, the corresponding region on the map of the calculation target region According to the distribution of pixels in the observed image data, extract the plant area from the calculation target area by specifying the position on the map,
Next, with reference to the numerical terrain model storage unit that stores the numerical terrain model coordinated in the map data, and the numerical surface layer model storage unit that stores the numerical surface layer model coordinated in the map data, The difference between the ground surface elevation and the surface elevation is calculated as the height of the plant group for each division unit obtained by dividing the plant region by a predetermined area unit,
Thereafter, the height per unit area of the plant group and the amount of carbon dioxide absorbed by the plant group are set in each of the plurality of height range divisions of the plant group set within the height range not corresponding to the forest. Based on the carbon dioxide absorption amount allocated in advance according to the causal relationship, the carbon dioxide absorption amount in the height range according to the height of the plant group is integrated according to the area in the entire region of the plant region, and within the calculation target region A method for calculating a carbon dioxide absorption effect that executes each process of calculating a carbon dioxide absorption effect in a non-forest area . A map data storage unit for storing map data;
An observation image data storage unit that stores observation image data obtained by observing the ground surface from the sky with a predetermined wavelength range that is coordinated to the map data and that can distinguish the plant group,
A numerical terrain model storage unit for storing a numerical terrain model of a calculation target area coordinated in the map data;
A numerical surface layer model storage unit for storing a numerical surface layer model of a calculation target area coordinated in the map data;
Plant / dioxide that stores the amount of carbon dioxide absorbed per unit area of a plant group in a predetermined height range within a range that is less than the height equivalent to a forest based on the causal relationship with the height of the plant group Carbon absorption table,
Target area setting means for setting the calculation target area;
A plant region extraction unit that refers to the map data storage unit and the observation image data storage unit, and extracts a plant region from the calculation target region according to a distribution of pixels in observation image data for a corresponding region on the map of the calculation target region;
With reference to the numerical terrain model storage unit and the numerical surface model storage unit, the difference between the ground surface elevation of the plant region and the surface elevation is calculated for each division unit obtained by dividing the plant region by a predetermined area unit. Plant height calculating means for calculating the height of the group,
Refer to the plant / carbon dioxide absorption amount table, determine which height range category the height of the plant group corresponds to, and determine the carbon dioxide absorption amount per predetermined unit area of each divided unit. And the carbon dioxide absorption effect calculation apparatus which has the effect calculation means which integrates this carbon dioxide absorption amount according to an area in the whole region of a plant area, and calculates the carbon dioxide absorption effect in the non-forest area | region in a calculation object area | region.

Images