首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Spatial differences in drought proneness and intensity of drought caused by differences in cropping patterns and crop growing environments within a district indicate the need for agricultural drought assessment at disaggregated level. The objective of this study is to use moderate resolution satellite images for detailed assessment of the agricultural drought situation at different administrative units (blocks) within a district. Monthly time composite NDVI images derived from moderate resolution AWiFS (60 m) and WiFS (180 m) images from Indian Remote Sensing satellites were analysed along with ground data on rainfall and crop sown areas for the kharif seasons (June – November) of 2002 (drought year), 2004 (early season drought) and 2005 (good monsoon year). The impact of the 2002 meteorological drought on crop area in different blocks of the district was assessed. The amplitude of crop condition variability in a severe drought year (2002) and a good year (2005) was used to map the degree of vulnerability of different blocks in the district to agricultural drought. The impact of early season deficit rainfall in 2004 on the agricultural situation and subsequent recovery of the agricultural situation was clearly shown. Agricultural drought assessment at disaggregated level using moderate resolution images is useful for prioritizing the problem areas within a district to undertake, in season drought management plans, such as alternate cropping strategies, as well as for end of the season drought relief management actions. The availability of ground data on rainfall, cropping pattern, crop calendar, irrigation, soil type etc., is very crucial in order to interpret the seasonal NDVI patterns at disaggregated level for drought assessment. The SWIR band of AWiFS sensor is a potential data source for assessing surface drought at the beginning of the season.  相似文献   

2.
Penman–Monteith method adapted to satellite data was used for the estimation of wheat crop evapotranspiration during the entire growth period using satellite data together with ground meteorological measurements. The IRS-1D/IRS-P6 LISS-III sensor data at 23.5 m spatial resolution for path 096 and row 059 covering the study area were used to derive, albedo, normalized difference vegetation index, leaf area index and crop height and then to estimate wheat crop evapotranspiration referred to as actual evapotranspiration (ETact). The ETact varied from 0.86 to 3.41 mm/day during the crop growth period. These values are on an average 16.40 % lower than wheat crop potential evapotranspiration (ETc) estimated as product of reference crop evapotranspiration estimated by Penman–Monteith method and lysimetric crop coefficient (Kc). The deviation of ETact from ETc is significant, when both the values were compared with t test for paired two sample means. Though the observations on ETact were taken from well maintained unstressed experimental plot of 120 × 120 m size, there was significant deviation. This deviation could be attributed to, the satellite images representing the actual crop evapotranspiration as function crop canopy biophysical parameters, condition of the crop stand, climatic and soil conditions and the microclimate variation over area of one hectare. However, Penman–Monteith method represents a flat rate of specific growth stage of the crop.  相似文献   

3.
Monitoring of Agricultural crops using remote sensing data is an emerging tool in recent years. Spatial determination of sowing date is an important input of any crop model. Geostationary satellite has the capability to provide data at high temporal interval to monitor vegetation throughout the entire growth period. A study was conducted to estimate the sowing date of wheat crop in major wheat growing states viz. Punjab, Haryana, Uttar Pradesh (UP), Madhya Pradesh (MP), Rajasthan and Bihar. Data acquired by Charged Couple Detector (CCD) onboard Indian geostationary satellite INSAT 3A have continental (Asia) coverage at 1 km?×?1 km spatial resolution in optical spectral bands with high temporal frequency. Daily operational Normalized Difference Vegetation Index (NDVI) product from INSAT 3A CCD available through Meteorological and Oceanographic Satellite Data Archival Centre (MOSDAC) was used to estimate sowing date of wheat crop in selected six states. Daily NDVI data acquired from September 1, 2010 to December 31, 2010 were used in this study. A composite of 7 days was prepared for further analysis of temporal profile of NDVI. Spatial wheat crop map derived from AWiFS (56 m) were re-sampled at INSAT 3A CCD parent resolution and applied over each 7 day composite. The characteristic temporal profiles of 7 day NDVI composite was used to determine sowing date. NDVI profile showed decreasing trend during maturity of kharif crop, minimum value after harvest and increasing trend after emergence of wheat crop. A mathematical model was made to capture the persistent positive slope of NDVI profile after an inflection point. The change in behavior of NDVI profile was detected on the basis of change in NDVI threshold of 0.3 and sowing date was estimated for wheat crop in six states. Seven days has been deducted after it reached to threshold value with persistent positive slope to get sowing date. The clear distinction between early sowing and late sowing regions was observed in study area. Variation of sowing date was observed ranging from November 1 to December 20. The estimated sowing date was validated with the reported sowing date for the known wheat crop regions. The RMSD of 3.2 (n?=?45) has been observed for wheat sowing date. This methodology can also be applied over different crops with the availability of crop maps.  相似文献   

4.
Assessment of area under agroforestry in Tehri district of North Western Himalaya, Uttarakhand, India has been done using GIS and remote sensing technology. The study district characterized by hilly terrain with varying elevations from 288 m to more than 2800 m and generally gentle slopes, valleys, flat land covers and agricultural terraces. High-resolution satellite imageries (spatial resolution 5.8 m) were used in this study for land uses and land covers classification. According to unsupervised classification, highest area was found under forest class (65.22%) followed by cropland (20.41%). Considerable area was also found under snow cover (9.45%) in the district. Area under agroforestry was estimated to be 5572.26 ha (1.53%) by this method, whereas it was estimated to be 7029.06 ha (1.93%) by supervised classification. Estimated cropland area comes out to be about 20.0%. An accuracy of 86.5% was found in this classification for agroforestry class. Highest area under agroforestry of 3707.36 ha was obtained in 1200–2000 m elevations followed by 2231.26 ha in 288–1200 m elevations. Negligible area was found on high elevation zones of more than 2800 m. The major agroforestry systems of dominated by Grewia oppositifolia (Bhimal), Celtis australis (Kharik) and Quercus leucotrichophora (Banj) were identified and mapped and remaining systems were grouped as others class. Estimated area under G. oppositifolia, C. australis and Q. leucotrichophora based systems come out to be 2330.82, 1456.80 and 1129.10 ha, respectively. These systems are multiple usufructs are food, fuelwood, fodder, fiber and small timber. It has been observed from the accuracy assessment that the estimates of area under agroforestry obtained under this study are reliable.  相似文献   

5.
6.
Crop yield estimation has an important role on economy development and its accuracy and speed influence yield price and helps in deciding the excess or deficit production conditions. The water productivity evaluates the irrigation command through water use efficiency (WUE). Remote sensing (RS) and geographical information system (GIS) techniques were used for crop yield and water productivity estimation of wheat crop (Triticum aestivum) grown in Tarafeni South Main Canal (TSMC) irrigation command of West Bengal State in India. One IRS P6 image and four wide field sensor (WiFS) images for different months of winter season were used to determine the Normalized Difference Vegetation Index (NDVI) and Soil Adjusted Vegetation Index (SAVI) for area under wheat crop. The temporally and spatially distributed spectral growth profile and AREASUM of NDVI (ANDVI) and SAVI (ASAVI) with time after sowing of wheat crop were developed and correlated with actual crop yield of wheat (Yact). The developed relationships between ASAVI and Yact resulted high correlation in comparison to that of ANDVI. Using the developed model the RS based wheat yield (YRS) predicted from ASAVI varied on entire TSMC irrigation command from 22.67 to 33.13 q ha−1 respectively, which gave an average yield of 26.50 q ha−1. The RS generated yield based water use efficiency (WUEYRS) for water supplied from canal of TSMC irrigation command was found to be 6.69 kg ha−1 mm−1.  相似文献   

7.
Field experiment was conducted in a sandy loam soil of Indian Agricultural Research Institute, New Delhi during the year 2011–13 to see the effect of irrigation, mulch and nitrogen on canopy spectral reflectance indices and their use in predicting the grain and biomass yield of wheat. The canopy reflectances were measured using a hand held ASD FieldSpec Spectroradiometer at booting stage of wheat. Four spectral reflectance indices (SRIs) viz. RNDVI (Red Normalized Difference Vegetation Index), GNDVI (Green Normalized Difference Vegetation Index), SR (Simple Ratio) and WI (Water Index) were computed using the spectral reflectance data. Out of these four indices, RNDVI, GNDVI and SR were significantly and positively related with the grain and biomass yield of wheat whereas WI was significantly and negatively related with the grain and biomass yield of wheat. Calibration with the second year data showed that among the SRIs, WI could account for respectively, 85 % and 86 % variation in grain and biomass yield of wheat with least RMSE (395 kg ha?1 (15 %) for grain yield and 1609 kg ha?1 (20 %) for biomass yield) and highest d index (0.95 for grain yield and 0.91 for biomass yield). Therefore it can be concluded that WI measured at booting stage can be successfully used for prediction of grain and biomass yield of wheat.  相似文献   

8.
Cotton aphid (Aphis gossypii) is considered as one of the most important agriculture pest for the cotton production. However, it is generally labor-intensive and time-consuming to obtain some information of Cotton aphid with conventional methods through direct measurement by sampling in the field. This study explores the potential of using a new method to obtain information of the Cotton aphid rapidly. In our study, the cotton canopy spectral indices (NDVI, VI_2, REDrefc, NIRrefc) and chlorophyll concentration, obtained from hand-held high spectrometer GreenSeeker and chlorophyll meter SPAD-502 and Cotton aphid amount derived from the artificial field-based survey were used to uncover the relationship between Cotton aphid amount and canopy spectral index and SPAD value of the cotton in city of Shihezi, China. The results showed that NDVI and NIRrefc were negatively related to Cotton aphid amount. VI_2 content had a significant and positive relationship with its amount. The non-linear three cubic models with alate Aphid amount as independent variables have been established between VI_2 value and alatae Aphid amount, which could explain 92.37 % of VI_2 value variance. SPAD values were also significantly and negatively correlated to the Aphid amount. The non-linear logarithm model with wingless Aphid amount as independent variables was the best for uncovering the relationship between SPAD value and wingless Aphid amount, which could explain 85.48 % of SPAD value variance. The results demonstrate the establishment of the function model provides a theoretical basis and techniques for indirect and rapid monitoring and management of Cotton aphid.  相似文献   

9.
Chlorophyll fluorescence is an indicator of plant photosynthetic activity and has been used to monitor the health status of vegetation. Several studies have exploited the application of red/far-red chlorophyll fluorescence ratio in detecting the impact of various types of stresses in plants. Recently, sunlight-induced chlorophyll fluorescence imaging has been used to detect and discriminate different stages of mosaic virus infection in potted cassava plants with a multi-spectral imaging system (MSIS). In this study, the MSIS is used to investigate the impact of drought and herbicide stress in field grown crop plants. Towards this control and treatment groups of colocasia and sweet potato plants were grown in laterite soil beds and the reflectance images of these crop plants were recorded up to 14-days of treatment at the Fraunhofer lines of O2 B at 687 nm and O2 A at 759.5 nm and the off-lines at 684 and 757.5 nm. The recorded images were analyzed using the Fraunhofer Line Discrimination technique to extract the sunlight-induced chlorophyll fluorescence component from the reflectance images of the plant leaves. As compared to the control group, the chlorophyll fluorescence image ratio (F 687/F 760) in the treatment groups of both the plant varieties shows an increasing trend with increase in the extent of stress. Further, the F 687/F 760 ratio was found to correlate with the net photosynthetic rate (Pn) and stomatal conductance (gs) of leaves. The correlation coefficient (R 2) for the relationship of F 687/F 760 ratio with Pn were found to be 0.78, 0.79 and 0.78, respectively for the control, herbicide treated and drought treated colocasia plants, while these were 0.77, 0.86 and 0.88, respectively for sweet potato plants. The results presented show the potential of proximal remote sensing and the application F 687/F 760 fluorescence image ratio for effective monitoring of stress-induced changes in field grown plants.  相似文献   

10.
Considering the contribution of the hardware biases to the estimated clock errors, an improved method for estimating the satellite inter-frequency clock bias (IFCB) is presented, i.e., the difference in the satellite clock error as computed from ionospheric-free pseudorange and carrier phase observations using L1/L2 and P1/P2 versus L1/L5 and P1/P5. The IFCB is composed of a constant and a variable part. The constant part is the inter-frequency hardware bias (IFHB). It contains the satellite and receiver hardware delays and can be expressed as a function of the DCBs [DCB (P1 ? P2) and DCB (P1 ? P5)]. When a reference satellite is selected, the satellite IFHB can be computed but is biased by a reference satellite IFHB. This bias will not affect the utilization of IFCB in positioning since it can be absorbed by the receiver clock error. Triple-frequency observations of 30 IGS stations between June 1, 2013, and May 31, 2014, were processed to show the variations of the IFHB. The IFHB values show a long-term variation with time. When a linear and a fourth-order harmonic function are used to model the estimated IFCB, which contains contributions of the hardware delays and clock errors, the results show that 89 % of the IFCB can be corrected given the current five triple-frequency GPS satellites with the averaged fitting RMS of 1.35 cm. Five days of data are processed to test the estimated satellite clock errors using the strategy presented. The residuals of P1/P5 and L1/L5 have a STD of <0.27 m and 0.97 cm, respectively. In addition, most predicted satellite IFCBs reach an accuracy of centimeter level and its mean accuracy of 5 days is better than 7 cm.  相似文献   

11.
Soil moisture estimation from satellite earth observation has emerged effectively advantageous due to the high temporal resolution, spatial resolution, coverage, and processing convenience it affords. In this paper, we present a study carried out to estimate soil moisture level at every location within Enugu State Nigeria from satellite earth observation. Comparative analysis of multiple indices for soil moisture estimation was carried out with a view to evaluating the robustness, correlation, appropriateness and accuracy of the indices in estimating the spatial distribution of soil moisture level in Enugu State. Results were correlated and validated with In-Situ soil moisture observations from multi-sample points. To achieve this, the Topographic Wetness Index (TWI), based on digital elevation data, the Temperature Vegetation Dryness Index (TVDI) and an improved TVDI (iTVDI) incorporating air temperature and a Digital Elevation Model (DEM) were calculated from ASTER global DEM and Landsat images. Possible dependencies of the indices on land cover type, topography, and precipitation were explored. In-Situ soil moisture data were used to validate the derived indices. The results showed that there was a positive significant relationship between iTVDI versus TVDI (R = 0.53, P value < 0.05), while in iTVDI versus TWI (R = 0.00, P value > 0.05) and TVDI versus TWI (R = ?0.01, P value > 0.05) no significant relationship existed. There was a strong relationship between iTVDI and topography, land cover type, and precipitation than other indices (TVDI, TWI). In situ measured soil moisture values showed negative significant relationship with TVDI (R = ?0.52, P value < 0.05) and iTVDI (R = ?0.63, P value < 0.05) but not with TWI (R = ?0.10, P value > 0.05). The iTVDI outperformed the other two index; having a stronger relationship with topography, precipitation, land cover classes and soil moisture. It concludes that although iTVDI outperformed other indices (TVDI, TWI) in soil moisture estimation, the decision of which index to apply is dependent on available data, the intent of usage and spatial scale.  相似文献   

12.
Present study was designed to determine the effect of various growing environments on sucking pest population dynamics in cotton and to work out their relation with spectral indices. Crop spectral reflectance in four IRS bands was measured with ground truth radiometer during 1000–1200 h in all the treatment combinations. Incidence of sucking pest in cotton was found out to be highly influenced by growing environments. The leafhopper and whitefly population was highest in 15 May sown cotton crop and was lowest in 15 April sown crop. Cultivar HS-6 was highly affected by both the sucking pest than the other cultivar H-1226. The spectral indices (SR, NDVI and TVI) were highest in 15 April sown crop at all the phenophases and were lowest in 15 May sown crop. The cultivar H-1226 showed higher values of spectral indices as compared to HS-6. The relationship of pests’ population with various spectral indices was established. Multiple regression models based on spectral indices can be used for prediction of sucking pest population more than 69 and 74 % accuracy in leafhopper and whitefly, respectively in cotton crop.  相似文献   

13.
The existence of mixed pixels in the satellite images has always been an area of concern. Adding to the challenge is an occurrence of non-linearity between the classes, which is generally overlooked. The study makes an attempt to solve the two frequently occurring problems by kernel based fuzzy approach. This research work deals with Possibilistic c-Means (PCM) classifier with local, global, spectral angle and hyper tangent kernels for wheat crop (Triticum aestivum) identification in Haridwar, Uttarakhand, India. The multi-temporal vegetation index data of Formosat-2 have been used which covers the whole phenology of wheat crop. The additional sensor Landsat-8 OLI imagery has been filled the crucial gap of Formosat-2 temporal datasets. Nine types of proposed kernels based PCM classifier have been applied on three temporal datasets (four, five and six date combinations) to classify two classes early sown and late sown wheat crop. These test results have been concluded that at optimized weighted constant KMOD and polynomial kernel was found effective to separate wheat crop. The five and six date combination were sufficient to discriminate early sown and late sown wheat crop.  相似文献   

14.
The aim of this study was to monitor changes in leaf spectral reflectance due to phytoaccumulation of trace elements (Cd, Pb, and As) in sunflower mutant (M5 mutant line 38/R4-R6/15-35-190-04-M5) grown in spiked and in situ metal-contaminated potted soils. Reflectance spectra (350–2500 nm) of leaves were collected using portable ASD spectroradiometer, and respective leaves sample were analyzed for total metal contents. The spectral changes were quite noticeable and showed increased visible and decreased NIR reflectance for sunflower grown in soil spiked with 900 mg As kg?1, and in in situ metal-contaminated soils. These changes also involved a blue-shift feature of red-edge position in the first derivatives spectra, studied vegetation indices and continuum removed absorption features at 495, 680, 970, 1165, 1435, 1780, and 1925 nm wavelength. Correlograms of leaf-metal concentration and reflectance values show highest degrees of overall correlation for visible, near-infrared, and water-sensitive wavelengths. Partial least square and multiple linear regression statistical models (cross-validated), respectively, based on Savitzky–Golay filter first-order derivative spectra and combination of spectral feature such as vegetation indices and band depths yielded good prediction of leaf-metal concentrations.  相似文献   

15.
Remote sensing technology becomes an effective and inexpensive technique for detecting disease in vegetation. In this study, an attempt has been done to discriminate healthy and late blight affected crop using remote sensing based indices such as NDVI and LSWI. NDVI and LSWI spectral profiles between healthy and late blight affected crop shows large difference. Mean difference in reflectance between two acquired dates Jan. 10 and 29, 2009 crop clusters varied from 31.28 % in red band, 7.7 % in NIR band and 6.23 % in SWIR bands in healthy crops while in late blight affected crops it is ?15.5 % in red, 44.4 % in NIR and ?14.61 % in SWIR bands. Negative percentage differences in reflectance indicate reflectance increases from Jan. 10, 2009 to Jan. 29, 2009, while positive difference indicate decrease in reflectance between the two dates. Since potato is an irrigated crop, these differences in reflectance are attributed to prevalent disease at that time. It is found that severely affected areas are Bardhman, Arambag, Bishnupur, Ghatal and Hugli taluka with crop damage areas are 4036.66, 1138.68, 2025.23, 469.15, and 380.08 ha, respectively.  相似文献   

16.
Detection of crop water stress is crucial for efficient irrigation water management. Potential of Satellite data to provide spatial and temporal dynamics of crop growth conditions makes it possible to monitor crop water stress at regional level. This study was conducted in parts of western Uttar Pradesh and Haryana. Multi-temporal Landsat data were used for detecting wheat crop water stress using vegetation indices (VIs), viz. vegetation water stress index (VWSI) and land surface wetness index water stress factor (Ws_LSWI). The estimated water stress from satellite data-based VIs was validated by water stress factor (Ws) derived from flux-tower data. The study observed Ws_LSWI to be better index for water stress detection. The results indicated that Ws_LSWI was superior over other index showing RMSE = 0.12, R2 = 0.65, whereas VWSI showed overestimated values with mean RD 4%.  相似文献   

17.
Field experiment was conducted during 2009–10 and 2010–11 rabi season at research farm of IARI, New Delhi for assessing the aphid infestation in mustard. In aphid infested plant the LAI was 67 to 94% lower than healthy plant. Chlorophyll concentration decreased to 50% in infested plant as compared to healthy plant. Infestation was more severe in late sown crop and due to aphid infestation the percentage oil content and yield was reduced significantly. The spectral reflectance of aphid infested canopy and healthy canopy taken in the laboratory had significant difference in NIR region. In the visible region, the reflectance peak occurred in healthy canopy at around 550–560 nm while this peak was lower by 31% in the aphid infested canopy. The reflectance for healthy crop was found to be more in visible as well as NIR region as compared to aphid infested canopy. The most significant spectral bands for the aphid infestation in mustard are in visible (550–560 nm) and near infrared regions (700–1250 nm and 1950–2450 nm). The different level of aphid infestation can be identified in 1950–2450 nm spectral regions. Spectral indices viz NDVI, RVI, AI and SIPI had significant correlation with aphid infestation. Hence these indices could be used for identifying aphid infestation in mustard.  相似文献   

18.
An Effective Model to Retrieve Soil Moisture from L- and C-Band SAR Data   总被引:1,自引:0,他引:1  
This study investigated an appropriate method for soil moisture retrieval from radar images and coincident ground measurements acquired over bare soil and sparsely vegetated regions. The adopted approach based on a single scattering integral equation method (IEM) was developed to establish the relationship between backscatter coefficient and surface soil parameters including volumetric soil moisture content and surface roughness. The performance of IEM in 0–7.6 cm is better than that in 0–20 cm. Moreover, IEM can simulate correctly the backscatter coefficients only for the root mean square (RMS) height s < 1.5 cm at C-band and s < 2.5 cm at L-band by using an exponential correlation function and for s > 1.5 cm at C-band and s > 2.5 cm at L-band by using Gaussian function. However, due to the difficulties involved in the parameterization of soil surface roughness, the estimated accuracy is not satisfactory for the inversion of IEM. This paper used a combined roughness parameter and Fresnel reflection coefficient to develop an empirical model. Simulations were performed to support experimental results and to highlight soil moisture content and surface roughness effects in different polarizations. Results showed that a good agreement was found between the IEM simulations and the SAR measurements over a wide range of soil moisture and surface roughness characteristics. The model had a significant operational advantage in soil moisture retrieval. The correlation coefficients were 77.03 % at L-band and 81.45 % at C-band with the RMSEs of 0.515 and 0.4996 dB, respectively. Additionally, this work offered insight into the required application accuracy of soil moisture retrieval at a large area of arid regions.  相似文献   

19.
Real time, accurate and reliable estimation of maize yield is valuable to policy makers in decision making. The current study was planned for yield estimation of spring maize using remote sensing and crop modeling. In crop modeling, the CERES-Maize model was calibrated and evaluated with the field experiment data and after calibration and evaluation, this model was used to forecast maize yield. A Field survey of 64 farm was also conducted in Faisalabad to collect data on initial field conditions and crop management data. These data were used to forecast maize yield using crop model at farmers’ field. While in remote sensing, peak season Landsat 8 images were classified for landcover classification using machine learning algorithm. After classification, time series normalized difference vegetation index (NDVI) and land surface temperature (LST) of the surveyed 64 farms were calculated. Principle component analysis were run to correlate the indicators with maize yield. The selected LSTs and NDVIs were used to develop yield forecasting equations using least absolute shrinkage and selection operator (LASSO) regression. Calibrated and evaluated results of CERES-Maize showed the mean absolute % error (MAPE) of 0.35–6.71% for all recorded variables. In remote sensing all machine learning algorithms showed the accuracy greater the 90%, however support vector machine (SVM-radial basis) showed the higher accuracy of 97%, that was used for classification of maize area. The accuracy of area estimated through SVM-radial basis was 91%, when validated with crop reporting service. Yield forecasting results of crop model were precise with RMSE of 255 kg ha?1, while remote sensing showed the RMSE of 397 kg ha?1. Overall strength of relationship between estimated and actual grain yields were good with R2 of 0.94 in both techniques. For regional yield forecasting remote sensing could be used due greater advantages of less input dataset and if focus is to assess specific stress, and interaction of plant genetics to soil and environmental conditions than crop model is very useful tool.  相似文献   

20.
Remote sensing and FAO 56 crop water model are used for estimating crop water requirement for paddy crop located in the main branch canal of Bhadra Command Area in Karnataka, India. The estimation of crop-water requirement depends on the meteorological factors, soil type and crop coefficients. The result obtained showed that water requirements of rabi crops higher than those of the kariff crops. The total irrigated area estimated from the IRS image is 29,353 ha. It is found that the total paddy crop acreage is 18,257 ha covering 62 % in the total irrigated area of the command area, Arecanut 20 %, coconut 15 % and sugarcane with other crops 3 %. The water requirement for paddy is 1180.4 mm for its entire growth period. The total water requirement for irrigation supply for crops in the entire command area is 5,790 at a demand of 0.10501 cusecs per ha.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号