Identiﬁcation of Paddy Varieties from Landsat 8 Satellite Image Data Using Spectral Unmixing Method in Indramayu Regency, Indonesia

: Indramayu Regency is the highest rice producer in West Java province, Indonesia. Accord-ing to the Central Statistics Agency (BPS), in 2021, rice production in 2020 reached 1,365,435.39 tons of GKG (milled dry grain). Technological developments in the food sector produce various kinds of premium quality rice and rice varieties resistant to climate change, such as Ciherang, Inpari 32 HDB and IR 64. The regular monitoring of speciﬁc rice varieties over large areas effectively maintains the quality and quantity of rice production. This study used remote sensing data to monitor rice conditions and distribution based on the spectral unmixing method. The spectral unmixing method was used to identify the percentage of the presence of a pure object in a pixel. The results obtained in this study were images of the endmember fractions of rice varieties and areas of dominant rice varieties used in the Indramayu district. The dominant variety detected with the processing results was the Inpari 32 HDB variety, with an area of 30,738.64 hectares. In comparison, varieties other than Inpari 32 HDB were also detected in several areas in the Indramayu district, with an area of 12,192.68 hectares.


Introduction
Southeast Asia is located on the equator, where most of the population relies on the agricultural sector as a livelihood source and development support [1]. Rice (Oryza sativa L.) is a species of the genus Oryza, which belongs to the family Poaceae (Gramineae). In Indonesia, 90% of the people consume rice as a staple food [2]. The high level of rice consumption must be balanced with appropriate production so as to not cause a shortage of supplies, resulting in food insecurity [3].
Indramayu Regency is a reasonably fertile area. Of the 204,011 ha area, 41.90% is paddy fields [4]. Based on data from the Central Statistics Agency (BPS), in 2021, rice production in Indramayu Regency in 2020 reached 1,365,435.39 tons of GKG (milled dry grain), with a productivity level of 6.07 tons/ha. These statistical data place Indramayu Regency in the first rank as the highest rice producer in West Java province. Other data show that there may be a bias in the statistics on area and rice productivity ranging from 12.61% to 20.72% for the rice harvested area and 0.23% to 0.25% for rice productivity, as reported by Hidayat (2018). This bias is caused by the data source, which are the harvested area (from monthly reports of agricultural statistics; Statistik Pertanian (SP)) and productivity from the Ubinan Survey [5]. Rice plants generally take 3-4 months to grow, from planting to harvesting, depending on the type of variety and the environmental conditions in which the rice grows.
Technological developments in the food sector produce various kinds of premium-quality rice and rice varieties resistant to global climate change caused by global warming [6]. Rice plants can be distinguished based on their varieties. There are wide varieties of rice plants, and almost every year they appear with better genetic traits [7]. Based on the annual data issued by the Directorate of Food Crops Germination in 2020, the dominant superior rice varieties grown in Indonesia include Ciherang, Inpari 32 HDB and IR 64.
Information on the distribution of dominant rice varieties using remote sensing techniques can be used as a reference for efforts to control and guarantee the quality of rice products to be marketed based on the type of variety. Many of them are quantitative measurements based on specific methods. For example, spectral classification uses the linear spectral unmixing (LSU) method [8]. The "spectral unmixing" of satellite images is one of the most widely used methods for obtaining mixed-pixel information. The spectral unmixing method is used to identify the percentage of the presence of a pure object (endmember) in a pixel. Therefore, the spectral unmixing method requires the endmember of each object to be identified.
The spectral mixture analysis was originally developed to interpret hyperspectral data using a "High Spectral Resolution Advanced Visible/Infrared Image Spectrometer" (HSR AVIRIS). The spectral unmixing method in its development has also been applied to analyze multispectral images [9]. For example, Liming Bai (2012) used Landsat 7 ETM+ imagery data to estimate forest cover in the Pingnan area. Gandharum (2015) utilized Landsat 5 TM to categorize rice plants infected with bacterial leaf blight. Tompolidi (2020) mapped a hydrothermal field using spectral unmixing from ASTER, Landsat 8/OLI and Sentinel-2 MSI data. Finally, Wu (2021) utilized spectral unmixing from stacked multitemporal remote sensing images with a variability of endmembers for change detection [9][10][11][12]. The objective of this study is to identify varieties of rice plants using Landsat 8 image data in Indramayu Regency using the spectral unmixing method. There remain some issues in the existing method for estimating paddy areas and rice production. Hopefully, this research can be used to estimate rice production accurately by estimating paddy areas of certain varieties.

Study Area
Indramayu Regency is one of the regencies in West Java province, Indonesia, as presented in Figure 1. It has geographic coordinates between 107 • 52 E and 108 • 36 E and 6 • 15 S and 6 • 40 S. Administratively, it is divided into 31 subdistricts with 317 villages. The location of Indramayu Regency stretches along the northern coast of Java Island with its topographic appearance. Most of it is plains or sloping areas with an average slope of 0-2%.
planting to harvesting, depending on the type of variety and the environmental conditions in which the rice grows. Technological developments in the food sector produce various kinds of premium-quality rice and rice varieties resistant to global climate change caused by global warming [6]. Rice plants can be distinguished based on their varieties. There are wide varieties of rice plants, and almost every year they appear with better genetic traits [7]. Based on the annual data issued by the Directorate of Food Crops Germination in 2020, the dominant superior rice varieties grown in Indonesia include Ciherang, Inpari 32 HDB and IR 64.
Information on the distribution of dominant rice varieties using remote sensing techniques can be used as a reference for efforts to control and guarantee the quality of rice products to be marketed based on the type of variety. Many of them are quantitative measurements based on specific methods. For example, spectral classification uses the linear spectral unmixing (LSU) method [8]. The "spectral unmixing" of satellite images is one of the most widely used methods for obtaining mixed-pixel information. The spectral unmixing method is used to identify the percentage of the presence of a pure object (endmember) in a pixel. Therefore, the spectral unmixing method requires the endmember of each object to be identified.
The spectral mixture analysis was originally developed to interpret hyperspectral data using a "High Spectral Resolution Advanced Visible/Infrared Image Spectrometer" (HSR AVIRIS). The spectral unmixing method in its development has also been applied to analyze multispectral images [9]. For example, Liming Bai (2012) used Landsat 7 ETM+ imagery data to estimate forest cover in the Pingnan area. Gandharum (2015) utilized Landsat 5 TM to categorize rice plants infected with bacterial leaf blight. Tompolidi (2020) mapped a hydrothermal field using spectral unmixing from ASTER, Landsat 8/OLI and Sentinel-2 MSI data. Finally, Wu (2021) utilized spectral unmixing from stacked multitemporal remote sensing images with a variability of endmembers for change detection [9][10][11][12]. The objective of this study is to identify varieties of rice plants using Landsat 8 image data in Indramayu Regency using the spectral unmixing method. There remain some issues in the existing method for estimating paddy areas and rice production. Hopefully, this research can be used to estimate rice production accurately by estimating paddy areas of certain varieties.

Study Area
Indramayu Regency is one of the regencies in West Java province, Indonesia, as presented in Figure 1. It has geographic coordinates between 107°52′ E and 108°36′ E and 6°15′ S and 6°40′ S. Administratively, it is divided into 31 subdistricts with 317 villages. The location of Indramayu Regency stretches along the northern coast of Java Island with its topographic appearance. Most of it is plains or sloping areas with an average slope of 0-2%.

Data
The data used in this study were Landsat 8 surface reflectance tier 1 satellite images with two scenes to cover the study area obtained through the Google Earth Engine platform (https://earthengine.google.com/, accessed on 25 January 2022) with the image recording time adjusted to the reproductive phase of the rice plant, namely, on 17 March 2021, in the Indramayu Regency area for identifying rice varieties. In addition, this study also used MODIS satellite imagery with a resolution of 500 m in a total of 730 scene images obtained through the Google Earth Engine platform (https://earthengine.google.com/, accessed on 25 January 2022) with an acquisition time from January to December 2021. MODIS satellite image data were used to determine the reproductive phase of rice plants and the threshold value of the reproductive phase of rice plants. This study also used vector data for raw rice fields in Indramayu Regency from the Geospatial Information Agency (BIG) to delineate the study area, and used spectral data for each rice variety obtained through direct field measurements in the Indramayu Regency area for the spectral library development. The software used in this research was Google Earth Engine and Quantum GIS.

Methods
The overall workflow of our study was structured into three sections as presented in Figure 2: First, we determined the reproductive phase of rice plants and determined the threshold value of the reproductive phase of rice plants from MODIS satellite imagery data. Second, the threshold value of the reproductive phase of rice was used in Landsat 8 satellite imagery for the phase separation of rice plants. Third, we determined the endmember used in the spectral unmixing algorithm. Finally, we implemented the spectral unmixing algorithm to the Landsat 8 data.

Data
The data used in this study were Landsat 8 surface reflectance tier 1 satellite images with two scenes to cover the study area obtained through the Google Earth Engine platform (https://earthengine.google.com/, accessed on 25 January 2022) with the image recording time adjusted to the reproductive phase of the rice plant, namely, on 17 March 2021, in the Indramayu Regency area for identifying rice varieties. In addition, this study also used MODIS satellite imagery with a resolution of 500 m in a total of 730 scene images obtained through the Google Earth Engine platform (https://earthengine.google.com/, accessed on 25 January 2022) with an acquisition time from January to December 2021. MODIS satellite image data were used to determine the reproductive phase of rice plants and the threshold value of the reproductive phase of rice plants. This study also used vector data for raw rice fields in Indramayu Regency from the Geospatial Information Agency (BIG) to delineate the study area, and used spectral data for each rice variety obtained through direct field measurements in the Indramayu Regency area for the spectral library development. The software used in this research was Google Earth Engine and Quantum GIS.

Methods
The overall workflow of our study was structured into three sections as presented in Figure 2: First, we determined the reproductive phase of rice plants and determined the threshold value of the reproductive phase of rice plants from MODIS satellite imagery data. Second, the threshold value of the reproductive phase of rice was used in Landsat 8 satellite imagery for the phase separation of rice plants. Third, we determined the endmember used in the spectral unmixing algorithm. Finally, we implemented the spectral unmixing algorithm to the Landsat 8 data.

Determination of the Reproductive Phase and Threshold Value of the Reproductive Phase of Rice Plants
The reproductive phase of rice plants in this study was determined using the NDVI and MNDWI spectral index. NDVI and MNDWI spectral index calculations were performed using BRDF MODIS/Terra+Aqua Nadir image data. The fashionable image data from January 2020 to December 2021 comprised of 730 image scenes.
NDVI is one of the vegetation indices commonly used for vegetation monitoring in remote sensing. This index was introduced by Rouse [13]. NDVI compares the subtraction value of the NIR surface reflectance band and the red band surface reflectance with the total value of the NIR surface reflectance band and the red band surface reflectance.
NDVI generally has a range of values from −1 to 1. A low NDVI value indicates that the vegetation on the pixel has low health [14]. Another index called the NDWI was proposed by McFeeters in 1996, and is often used for detecting and monitoring minor changes in the water content of bodies of water. The NDWI can enhance water bodies in satellite images by using NIR (near-infrared) and green (visible green) spectral bands. The index's disadvantage is that it is sensitive to constructed buildings, which can lead to an overestimation of water bodies [15,16].
To improve open water features while effectively reducing and even eliminating builtup land noise, as well as vegetation and soil noise, a modified NDWI called the MNDWI was proposed by Xu in 2006 [17]. The method compares the subtraction value of the green band surface reflectance and the SWIR band surface reflectance with the total value of the green band surface reflectance and SWIR band surface reflectance.
The MNDWI method has advantages over NDWI. When compared to NDWI, it can be utilized to better define water features combined with vegetation, since it can reduce soil and vegetation disruptions to identify water bodies [17][18][19].
The determination of the reproductive phase of rice plants in this study was obtained through the extraction of the spectral index values of NDVI and MNDWI using 5 sample points spread over several rice fields. The determination of these 5 sample points was carried out randomly in rice fields in Indramayu Regency. The results of the extraction of the spectral index values were determined in a time series graph from January 2020 to December 2021 to identify the reproductive phase of rice plants and determine the threshold value of the reproductive phase of rice plants. The threshold value was determined based on observations in the time series graph of the NDVI and MNDWI spectral index values of the MODIS Terra/Aqua data.

Calculation of NDVI and MNDWI Spectral Index with Landsat 8 Satellite
The calculation of the NDVI and MNDWI spectral indices was carried out; then, the spectral index value was masked with the threshold value obtained from calculations on MODIS satellite imagery. The results produced were satellite images of the reproductive phase of rice plants. The images of the reproductive phase were calculated using the spectral unmixing method.

Endmember Determination
The endmember is the pure reflectance spectral value of an object [20]. Endmember values usually correspond to macroscopic objects in a study area, such as water, soil, metal or any natural materials [21]. To perform the spectral unmixing process, endmember data had to be observed. The number of endmember data used was better than the number of available spectral channels [10]. Endmember data could be obtained through (1) direct measurements of the field using a spectrometer, (2) spectral libraries or (3) through image transformation [22].
The determination of the endmembers in this study was based on several types of rice varieties, namely, Ciherang, IR64 and Inpari 32 HDB. The endmember value was obtained through direct measurements in Indramayu Regency rice fields using an Ocean Optics Spectrometer. Based on field observations, rice leaves started to become dense in the reproductive phase; thus, the selection of endmembers was narrowed down to the reproductive phase, as shown in Figure 3.
The determination of the endmembers in this study was based on several types rice varieties, namely, Ciherang, IR64 and Inpari 32 HDB. The endmember value was o tained through direct measurements in Indramayu Regency rice fields using an Oce Optics Spectrometer. Based on field observations, rice leaves started to become dense the reproductive phase; thus, the selection of endmembers was narrowed down to t reproductive phase, as shown in Figure 3.

Spectral Unmixing Algorithm Calculation
Spectral unmixing is a procedure in which the spectral values of mixed pixels a decomposed into a collection of constituent spectra or endmembers and a fraction map gray map, showing the proportion of each endmember contained in the pixel [22].
In the results of high-resolution satellite imagery, one pixel can represent one obje However, in the results of medium-resolution satellite images, one pixel can contain v ious objects. Therefore, it was necessary to simplify the mixed pixel using the endmemb value.
The results of this method were highly dependent on the endmember data inp Unmixing spectral classification produced a pair of gray-level images. Each consisted an abundance map and a quantitative pixel error image for all endmembers. The resulti endmember fraction map was in the form of the percentage of the presence of an obj (endmember) in a pixel, indicated by a bright hue (white) if the object's presence was hi and a dark hue (black) if the object's presence was low [9].
We calculated the spectral unmixing using reflectance band 1, band 2, band 3, ba 4 and band 5 on Landsat 8 surface reflectance collection 2 tier 1 images. The value of t resulting fraction image ranged from 0 to 1. Data processing was carried out on the Goog Earth Engine platform.

Rice Plant Phase Identification
The generative phase of rice plants could be identified by decreasing the NDVI val from the peak value. Based on Figure 4, it could be identified that the generative phase rice was in April 2020, August 2020, March 2021 and August 2021. Based on the range values owned by the two indices in the graph formed in Figures 4 and 5, the thresho value for the reproductive phase could be determined as values of NDVI > 0.7 a MNDWI < −0.4. The threshold value was determined based on observations of the tim series graph of the formed NDVI and MNDWI spectral index values. The threshold val

Spectral Unmixing Algorithm Calculation
Spectral unmixing is a procedure in which the spectral values of mixed pixels are decomposed into a collection of constituent spectra or endmembers and a fraction map or gray map, showing the proportion of each endmember contained in the pixel [21].
In the results of high-resolution satellite imagery, one pixel can represent one object. However, in the results of medium-resolution satellite images, one pixel can contain various objects. Therefore, it was necessary to simplify the mixed pixel using the endmember value.
The results of this method were highly dependent on the endmember data input. Unmixing spectral classification produced a pair of gray-level images. Each consisted of an abundance map and a quantitative pixel error image for all endmembers. The resulting endmember fraction map was in the form of the percentage of the presence of an object (endmember) in a pixel, indicated by a bright hue (white) if the object's presence was high and a dark hue (black) if the object's presence was low [9].
We calculated the spectral unmixing using reflectance band 1, band 2, band 3, band 4 and band 5 on Landsat 8 surface reflectance collection 2 tier 1 images. The value of the resulting fraction image ranged from 0 to 1. Data processing was carried out on the Google Earth Engine platform.

Rice Plant Phase Identification
The generative phase of rice plants could be identified by decreasing the NDVI value from the peak value. Based on Figure 4, it could be identified that the generative phase of rice was in April 2020, August 2020, March 2021 and August 2021. Based on the range of values owned by the two indices in the graph formed in Figures 4 and 5, the threshold value for the reproductive phase could be determined as values of NDVI > 0.7 and MNDWI < −0.4. The threshold value was determined based on observations of the time series graph of the formed NDVI and MNDWI spectral index values. The threshold value indicated that the rice fields were in the reproductive phase where the NDVI index value was high while the MNDWI index value was low. indicated that the rice fields were in the reproductive phase where the NDVI index value was high while the MNDWI index value was low.

Landsat 8 Satellite Imagery Reproductive Phase of Rice Plants
The reproductive phase of rice plants in Indramayu Regency was identified using MODIS data. It occurred in April 2020, August 2020, March 2021 and August 2021. We utilized the reproductive phase in March 2021 for this research. The Landsat 8 satellite image used was taken on 17 March 2021. The Landsat 8 image was subsetted with vector data for raw rice fields in Indramayu Regency. The results of the subsetted image can be seen in Figure 6. indicated that the rice fields were in the reproductive phase where the NDVI index value was high while the MNDWI index value was low.

Landsat 8 Satellite Imagery Reproductive Phase of Rice Plants
The reproductive phase of rice plants in Indramayu Regency was identified using MODIS data. It occurred in April 2020, August 2020, March 2021 and August 2021. We utilized the reproductive phase in March 2021 for this research. The Landsat 8 satellite image used was taken on 17 March 2021. The Landsat 8 image was subsetted with vecto data for raw rice fields in Indramayu Regency. The results of the subsetted image can be seen in Figure 6.

Landsat 8 Satellite Imagery Reproductive Phase of Rice Plants
The reproductive phase of rice plants in Indramayu Regency was identified using MODIS data. It occurred in April 2020, August 2020, March 2021 and August 2021. We utilized the reproductive phase in March 2021 for this research. The Landsat 8 satellite image used was taken on 17 March 2021. The Landsat 8 image was subsetted with vector data for raw rice fields in Indramayu Regency. The results of the subsetted image can be seen in Figure 6.
The NDVI and MNDWI spectral indices were calculated and then masked based on the threshold value obtained from the calculation on the MODIS image. The resulting results were in the form of satellite images of the reproductive phase of rice plants. The results of satellite images in the reproductive phase of rice plants can be seen in Figure 7. The reproductive phase image was calculated using the spectral unmixing method. The NDVI and MNDWI spectral indices were calculated and then masked based on the threshold value obtained from the calculation on the MODIS image. The resulting results were in the form of satellite images of the reproductive phase of rice plants. The results of satellite images in the reproductive phase of rice plants can be seen in Figure 7. The reproductive phase image was calculated using the spectral unmixing method.

Distribution of Rice Plant Varieties in Indramayu Regency
The fraction map generated from the spectral unmixing algorithm was the percentage of endmembers for each pixel. The resulting endmember percentage values range was in the 0-1 or 0-100% range. The higher the endmember percentage, the brighter the pixels. The results of the processing with the spectral unmixing algorithm can be seen in Figure 8.

Distribution of Rice Plant Varieties in Indramayu Regency
The fraction map generated from the spectral unmixing algorithm was the percentage of endmembers for each pixel. The resulting endmember percentage values range was in the 0-1 or 0-100% range. The higher the endmember percentage, the brighter the pixels. The results of the processing with the spectral unmixing algorithm can be seen in Figure  8.  Figure 8 shows the results of the spectral unmixing algorithm. Based on Figure 8, the dominant variety used in Indramayu Regency was the Inpari 32 HDB variety, which was indicated with a bright hue (red) if the presence of the object was high and a dark hue (black) if the presence of the object was low. Figure 8 shows that there were no Ciherang varieties and IR 64 varieties in the data processing results, which were indicated by the absence of other colors on the map. Figures 9-11 show images of the endmember fractions of each rice variety. Based on the statistical values generated from the endmember fraction map, it was obtained that the Inpari 32 HDB variety had a value range of 0 to 0.978. This indicated that there was rice of the Inpari 32 HDB variety in some pixels. This variety had dominant appearance with 50% larger endmembers indicated in blue and light blue. Because the endmember was less than 50%, other variants (IR 64 and Ciherang) had no pixels with these colors. The Ciherang variety had a value of 0 in each pixel. This indicated that the Ciherang variety did not exist in any pixels. The IR 64 variety had a value range of 0 to 0.440, which indicated that in some pixels there was rice of the IR 64 variety, but because the value range was still below 0.5 or 50%, this indicated that there were other dominant varieties in that pixel.  Figure 8 shows the results of the spectral unmixing algorithm. Based on Figure 8, the dominant variety used in Indramayu Regency was the Inpari 32 HDB variety, which was indicated with a bright hue (red) if the presence of the object was high and a dark hue (black) if the presence of the object was low. Figure 8 shows that there were no Ciherang varieties and IR 64 varieties in the data processing results, which were indicated by the absence of other colors on the map. Figures 9-11 show images of the endmember fractions of each rice variety. Based on the statistical values generated from the endmember fraction map, it was obtained that the Inpari 32 HDB variety had a value range of 0 to 0.978. This indicated that there was rice of the Inpari 32 HDB variety in some pixels. This variety had dominant appearance with 50% larger endmembers indicated in blue and light blue. Because the endmember was less than 50%, other variants (IR 64 and Ciherang) had no pixels with these colors. The Ciherang variety had a value of 0 in each pixel. This indicated that the Ciherang variety did not exist in any pixels. The IR 64 variety had a value range of 0 to 0.440, which indicated that in some pixels there was rice of the IR 64 variety, but because the value range was still below 0.5 or 50%, this indicated that there were other dominant varieties in that pixel.     Based on Table 1 and Figure 12, the dominant variety used in Indramayu Regency was the Inpari 32 HDB variety. However, the range of values indicated that there were still pixels whose values were below 50%. The low values of the pixels indicated that there were other varieties in the pixel. Therefore, this study determined a threshold value to distinguish the Inpari 32 HDB variety and other varieties (besides Inpari 32 HDB, Ciherang and IR 64). The threshold value used to detect the presence of the Inpari 32 HDB variety in the fraction image was ≥50%. Meanwhile, the threshold value used to detect the presence of other varieties in the fractional image was < 50%. The threshold value was based on the presence of 32 HDB Inpari varieties in one pixel. If the image fraction value was below 50%, it indicated other varieties in that pixel. Using a threshold value above 50% indicated that in that pixel the dominant variety was the Inpari 32 HDB variety.  Based on Table 1 and Figure 12, the dominant variety used in Indramayu Regency was the Inpari 32 HDB variety. However, the range of values indicated that there were still pixels whose values were below 50%. The low values of the pixels indicated that there were other varieties in the pixel. Therefore, this study determined a threshold value to distinguish the Inpari 32 HDB variety and other varieties (besides Inpari 32 HDB, Ciherang and IR 64). The threshold value used to detect the presence of the Inpari 32 HDB variety in the fraction image was ≥50%. Meanwhile, the threshold value used to detect the presence of other varieties in the fractional image was <50%. The threshold value was based on the presence of 32 HDB Inpari varieties in one pixel. If the image fraction value was below 50%, it indicated other varieties in that pixel. Using a threshold value above 50% indicated that in that pixel the dominant variety was the Inpari 32 HDB variety. Based on Table 2, the area of paddy fields planted with the Inpari 32 HDB variety with a fraction value ≥ 50% was 30,738.65 hectares. Meanwhile, the area of paddy fields planted with other varieties (besides Inpari 32 HDB, Ciherang and IR 64) with a fraction value < 50% was 12,192.68 hectares. Based on the results of data processing, it could be seen that the dominant variety used by farmers in Indramayu Regency in March 2021 was the Inpari 32 HDB variety. Based on Table 2, the area of paddy fields planted with the Inpari 32 HDB variety with a fraction value ≥ 50% was 30,738.65 hectares. Meanwhile, the area of paddy fields planted with other varieties (besides Inpari 32 HDB, Ciherang and IR 64) with a fraction value < 50% was 12,192.68 hectares. Based on the results of data processing, it could be seen that the dominant variety used by farmers in Indramayu Regency in March 2021 was the Inpari 32 HDB variety.

Validation of Rice Fields
Information on the distribution of dominant rice varieties from the results of Landsat 8 satellite image processing needed to be validated with data on the use of rice varieties obtained from direct field surveys to ensure the suitability of the information on the time and location of the use of rice varieties. Field validation data in the form of rice field coordinate points and information regarding the use of rice plant varieties were obtained. Validation points were taken randomly in rice fields in Indramayu Regency. The validation results are presented in Table 3 and Figure 13.

Validation of Rice Fields
Information on the distribution of dominant rice varieties from the results of Landsat 8 satellite image processing needed to be validated with data on the use of rice varieties obtained from direct field surveys to ensure the suitability of the information on the time and location of the use of rice varieties. Field validation data in the form of rice field coordinate points and information regarding the use of rice plant varieties were obtained. Validation points were taken randomly in rice fields in Indramayu Regency. The validation results are presented in Table 3 and Figure 13.  Based on the comparison of the results between the image of the endmember fraction and the data on the use of rice varieties in Table 3, it can be seen that there were several differences in the data. Based on sample point number three in Table 3 in the image of the endmember fraction, it was identified that the rice field used the Inpari 32 HDB rice plant  Based on the comparison of the results between the image of the endmember fraction and the data on the use of rice varieties in Table 3, it can be seen that there were several differences in the data. Based on sample point number three in Table 3 in the image of the endmember fraction, it was identified that the rice field used the Inpari 32 HDB rice plant variety, but the data on the use of rice plant varieties stated that the rice field used another rice plant variety, namely, the Borang variety (Inpari 32). Based on the results of the processing in this study, there were differences in the results produced by the field validation data.
The difference in the results between the endmember fraction image and data on rice varieties could be influenced by several factors, including the image data used. The spectral unmixing method was intended for hyperspectral images with a high spatial resolution, but this study used images with a medium resolution. Another factor could be the use of spectral libraries, given that there are wide varieties of rice plants in Indonesia. In contrast, in this study, only three types of spectral libraries were used, which could have affected the results obtained.

Conclusions
After applying the threshold value to distinguish between the Inpari 32 HDB variety and other varieties, rice fields planted with the Inpari 32 HDB variety comprised 30,738.65 hectares. Meanwhile, the area of paddy fields planted with other varieties was 12,192.68. Based on the data processing results, it was found that the dominant variety used by farmers in Indramayu Regency in March 2021 was the Inpari 32 HDB variety. Based on the comparison of the results of data processing with data on the use of rice varieties, it was found that there were several differences in the results produced. This showed that the spectral unmixing method using Landsat 8 satellite images still struggled to identify rice varieties. This could have been caused by several factors, including the image data used having a low spatial resolution. Other factors could be in the form of the spectral libraries used, where this study only used three types of spectral libraries. These factors could have affected the obtained results.
Author Contributions: Iqbal Maulana Cipta, Lalu Muhamad Jaelani and Hartanto Sanjaya designed the paper. Iqbal Maulana Cipta conducted the data processing, statistical analysis and wrote the paper. Lalu Muhamad Jaelani and Hartanto Sanjaya contributed great efforts to the research methodology, substantial editing, final writing and figure optimization. All the authors contributed to the paper revision. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement:
The Landsat 8 surface reflectance tier 1 images were obtained with the Google Earth Engine (https://earthengine.google.com/, accessed on 25 January 2022). The MODIS satellite imagery was obtained with the Google Earth Engine (https://earthengine.google.com/, accessed on 25 January 2022). Spectral data for each rice variety were obtained through direct field measurements.