LU500138B1 - Landsat8-basierte Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses - Google Patents

Landsat8-basierte Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses Download PDF

Info

Publication number
LU500138B1
LU500138B1 LU500138A LU500138A LU500138B1 LU 500138 B1 LU500138 B1 LU 500138B1 LU 500138 A LU500138 A LU 500138A LU 500138 A LU500138 A LU 500138A LU 500138 B1 LU500138 B1 LU 500138B1
Authority
LU
Luxembourg
Prior art keywords
model
soil
soil salinity
salinity
landsat8
Prior art date
Application number
LU500138A
Other languages
English (en)
Other versions
LU500138A1 (de
Inventor
Xiangyun Song
Xinwei Liu
Xiaoguang Zang
Dejie Cui
Original Assignee
Univ Qingdao Agricultural
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Univ Qingdao Agricultural filed Critical Univ Qingdao Agricultural
Publication of LU500138A1 publication Critical patent/LU500138A1/de
Application granted granted Critical
Publication of LU500138B1 publication Critical patent/LU500138B1/de

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • AHUMAN NECESSITIES
    • A01AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01BSOIL WORKING IN AGRICULTURE OR FORESTRY; PARTS, DETAILS, OR ACCESSORIES OF AGRICULTURAL MACHINES OR IMPLEMENTS, IN GENERAL
    • A01B79/00Methods for working soil
    • A01B79/005Precision agriculture
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/194Terrestrial scenes using hyperspectral data, i.e. more or other wavelengths than RGB

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Astronomy & Astrophysics (AREA)
  • Mechanical Engineering (AREA)
  • Soil Sciences (AREA)
  • Environmental Sciences (AREA)
  • Mathematical Physics (AREA)
  • Multimedia (AREA)
  • Pure & Applied Mathematics (AREA)
  • Remote Sensing (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Geometry (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

Dieses Patent erfindet ein Verfahren zur Inversion des Bodensalzgehalts basierend auf Landsat8-Fernerkundungsbildern in den folgenden Schritten. Schritt 1: Sammlung und Erwerbung von Bodensalzgehaltsdaten; Schritt 2: Erwerbung und Verarbeitung von Fernerkundungsbilddaten; Schritt 3: Auslesen und Berechnen der Reflexionsdaten der jeweiligen Frequenzbänder der Landsat8 OLI-Bilder; Schritt 4: Korrelationsanalyse des fernerkundeten Reflexionsgrads und des transformierten fernerkundeten Reflexionsgrads mit dem Bodensalzgehalt, nachdem die entsprechenden lg(R)- und 1/R-Änderungen für jedes Wellenband durchgeführt wurden, und vorläufiges Screening der empfindlichen Wellenbänder; Schritt 5: Durchführen einer multiplen linearen Regressionsanalyse von 1/R, lg(R) des Reflexionsgrads jeder Frequenzbande mit dem in Schritt 1 erhaltenen Bodensalzgehalt, wodurch ein Regressionsmodell erstellt wird; Schritt 6: Anwendung des oben erstellten besten Prognosemodells auf Landsat8-Fernerkundungsbilder eines typischen Untersuchungsgebiets im Delta des Gelben Flusses, wobei Landsat8-basierte Ferninversionskarten des Bodensalzgehalts erstellt werden. Diese Erfindung verwendet das Inversionsmodell, um den Bodensalzgehalt zu invertieren und eine deskriptive statistische Analyse durchzuführen, und vergleicht es mit dem gemessenen Bodensalzgehalt. Fig.1

Description

BESCHREIBUNG Landsat8-basierte Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses Technische Bereiche Diese Erfindung bezieht sich auf ein Inversionsverfahren, insbesondere auf ein Landsat8-basiertes Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses.
Hintergrundtechnik Bei der Bodenversalzung handelt es sich um ein Phänomen der Salzakkumulation in ariden und semiariden Gebieten aufgrund der Anhebung des Grundwasserspiegels und der Verdunstung von der Bodenwasseroberfläche, die größer ist als die Grundwasserneubildung.
Die Versalzung des Bodens ist in China ernsthaft, und die Fläche des salzhaltigen Bodens beträgt etwa 1,0x108ha, von denen der restliche salzhaltige Boden mit bis zu 45% am höchsten, der moderne salzhaltige Boden am zweithäufigsten und der oberflächliche salzhaltige Boden am geringsten ist.
Darunter befinden sich küstennahe salzhaltige Böden hauptsächlich im östlichen Teil Chinas und im Küstenbereich des Südchinesischen Meeres, mit einer Küstenlinie von mehr als 18,000 km Länge, die sich über 24 Grad von Süden nach Norden erstreckt.
Das Delta des Gelben Flusses ist ein typisches Küstenversalzungsgebiet, und der Bodensalzgehalt ist die Hauptursache für die Einschränkung des Getreideertrags und der Getreidequalität in der Region.
Daher ist die rechtzeitige und genaue Ermittlung des
Bodensalzgehalts äußerst wichtig, um die Bodenqualität zu verbessern, den Getreideertrag in der Region zu steigern, das Einkommen der Landwirte zu erhöhen und eine nachhaltige Entwicklung der regionalen Landwirtschaft zu erreichen.
Traditionell wurden Informationen über den Bodensalzgehalt hauptsächlich im Feld durch Erhebungsmethoden gesammelt. Mit der Entwicklung und Anwendung der Fernerkundungstechnologie hat die Extraktion des Salzgehalts mittels Fernerkundungstechnologie eine qualitative Klassifizierung des Bodensalzgehalts auf der Basis von Satellitenbildern erfahren und sich zu einer quantitativen Schätzung des Bodensalzgehalts auf der Basis von hyperspektralen und bodennahen multispektralen Bildern entwickelt. Zhang Tongrui et al. führten Korrelations- und Regressionsanalysen durch, indem sie gemessene Hyperspektraldaten und OLI-Bilddaten in diesem Untersuchungsgebiet miteinander kombinierten, und wählten schließlich das optimale Modell als lineares Modell mit SAVI als abhängiger Variable aus, das den Bodensalzgehalt in Winterweizenanbaugebieten in der Region des Gelben Flussdeltas genauer und zeitnah ermitteln kann; Maimaimaiti - Shawuti benutzte auch die gleiche Methode der Kombination von WorldView-2-Bildern und gemessenen Hyperspektraldaten, um den Bodensalzgehalt des Keliya-Flussgebiets in Xinjiang mit Hilfe der partiellen Kleinstquadrat-Regression bzw. der BP-Neuronalnetz-Methode quantitativ zu invertieren, wobei das R’ des mit dem BP-Neuronalnetz-Modell erstellten Modells bei 0,851 lag, wodurch der Bodensalzgehalt in der Region genauer geschätzt und vorhergesagt werden konnte; Zhao Gengxing et al. erstellten Fernerkundungs-Inversionsmodelle für den Bodensalzgehalt in der Delta-Region des Gelben Flusses, indem sie multispektrale OLI- und hyperspektrale HSI-Bilder mit Hilfe von multipler linearer Regression, partieller Kleinstquadrat-Regression, neuronalem BP-Netzwerk, Support-Vektor-Maschine bzw. Random-Forest-Methoden fusionierten. Durch den Vergleich der Studie wurde festgestellt, dass das durch die BP-Methode des neuronalen Netzwerks konstruierte Modell die höchste Genauigkeit hatte, das Modell R? war so hoch wie 0,966, und das Modell konnte den Bodensalzgehalt in der Region rechtzeitig und genau ermitteln. Da es sich jedoch um ein Dark-Box-Modell handelt, kann die Formel nicht genau geschrieben werden, und es ist für Forscher schwierig, sich auf diese Modellart zu beziehen. Auf der Grundlage der obigen Analyse besteht ein Mangel an Fernerkundungs-Inversionsmodellen, die alle Landnutzungsarten und Anbaupflanzen abdecken und für die Delta-Region des Gelben Flusses einfach zu bedienen sind.
Inhalt der Erfindung Um die oben genannten Probleme zu lösen, bietet die vorliegende Erfindung eine auf Landsat8 basierende Inversionsmethode für den Bodensalzgehalt in größeren Maßstäben im Delta des Gelben Flusses, die Landsat8 OLI-Bilder und den im Labor ermittelten Bodensalzgehalt als Datenquellen verwendet, ein Modell zur Schätzung des Bodensalzgehalts in der Region des Gelben Flussdeltas unter Verwendung mathematischer und statistischer Analysemethoden aufstellt und empfindliche Frequenzbänder und das beste Inversionsmodell für den Bodensalzgehalt untersucht. Dies bietet eine gewisse technische Unterstützung für die rechtzeitige und genaue Erfassung des Bodensalzgehalts in der Region in der Zukunft und stellt eine Grundlage für das Management, die Vorbeugung und Kontrolle des Bodensalzgehalts und die nachhaltige Entwicklung in der Delta-Region des Gelben Flusses dar.
Diese Erfindung wird durch die folgenden technischen Lösungen erreicht, um eine auf Landsat8 basierende Inversionsmethode für den Bodensalzgehalt im Delta des Gelben Flusses bereitzustellen.
Die speziellen Schritte sind wie folgt: Schritt 1: Sammlung und Erwerbung von Bodensalzgehaltsdaten; Schritt 2: Erwerbung und Verarbeitung von Fernerkundungsbilddaten; Schritt 3: Auslese und Berechnung der Reflexionsdaten für die jeweiligen Frequenzbänder aus Landsat8 OLI-Bildern mit der Software ArcGIS10.0; Schritt 4: Vorläufiges Screening empfindlicher Frequenzbänder durch entsprechende 1g(R)- und 1/R-Änderungen für jedes Frequenzband, gefolgt von einer Korrelationsanalyse des fernerkundeten Reflexionsgrads und des transformierten fernerkundeten Reflexionsgrads mit dem jeweiligen Salzgehalt des Bodens mittels SPSS-Software; Schritt 5: Multiple lineare Regressionsanalyse von 1/R, 1g(R) für jede Bandenreflexion mit dem in Schritt 1 erhaltenen Bodensalzgehalt mittels SPSS-Software, um ein Regressionsmodell zu erstellen und zu testen, ob die Reflexionsänderung der einzelnen Bandenwerte mit dem Bodensalzgehalt modelliert werden kann; Schritt 6: Erstellung von Fernerkundungs-Inversionsmodellen des Bodensalzgehalts für log-transformierte bzw. invers-transformierte Fernerkundungs-Bilddaten zum Vergleich des nach dieser einzigen Methode erstellten Fernerkundungs-Bilddatenmodells mit dem Fernerkundungs-Bilddatenmodell, das nach dem Mischen der beiden
Methoden-Transformationen erstellt wurde, und zur Analyse der Auswirkungen der verschiedenen Datentransformationen auf die Modellgenauigkeit; Schritt 7: Das oben erstellte beste Prognosemodell wurde auf Landsat8-Fernerkundungsbilder des Untersuchungsgebiets angewendet, wodurch Landsat8-Fernerkundungs-Inversionskarten des Bodensalzgehalts erstellt wurden.
Insbesondere für Schritt 1 müssen Faktoren wie Bodentyp, Vegetationstyp und die Art der Landnutzung im Untersuchungsgebiet berücksichtigt werden, damit die Probenahmestellen die tatsächlichen Bedingungen im Untersuchungsgebiet genau wiedergeben.
Dies wird speziell durch die folgenden Methoden erreicht: Schritt 1.1: Zunächst wurde ein Rasterlayout erstellt und an die Landnutzungsart sowie den Straßenzugang angepasst.
Insgesamt wurden 86 Probenahmestellen angelegt, und um jede Probenahmestelle herum wurden 3-6 Flecken Boden mit einer Schaufel in der Tiefe von 0-20 cm gesammelt und gründlich gemischt.
Dann wurden 1 kg bis 2 kg in Plastikbeutel gefüllt, beschriftet und die Koordinaten der Probenahmestellen, der Zeitpunkt der Probenahme, die Wetterbedingungen, die Landnutzungsart und der Bepflanzungstyp aufgezeichnet; Schritt 1.2: Die gesammelten Bodenproben wurden in Innenräumen getrocknet.
Nach dem Trocknen wurden die Graswurzeln ausgewählt und die Bodenproben wurden gemahlen und durch ein Imm-Sieb gegeben.
Zur Aufbewahrung für spätere Experimente wurden sie in Weithalsgläsern konserviert; Schritt 1.3: Der in Schritt 1.2 erhaltene Boden wurde mit der Leitfähigkeits-Methode bestimmt. Dies wurde wie folgt durchgeführt: Eine entsprechende Menge an Bodenprobe und destilliertem Wasser wurde gemäß dem Gewichtsverhältnis Boden-Wasser von 1:5 in ein Glasreagenzglas gefüllt und 5 Minuten lang geschüttelt. Anschließend wurde das Filtrat filtriert und die Leitfähigkeit mit einem Leitfähigkeitsmessgerät gemessen. Nachdem die Leitfähigkeit bestimmt wurde, wurden die Daten zur Bodensalinität durch Umrechnung der Formel entsprechend der Beziehung zwischen Leitfähigkeit und Bodensalinität ermittelt. Die spezifische Umrechnungsformel lautete wie folgt: S=3,047EC1:5-0,493 (r=0,981””, p<0,001). In der Formel bezieht sich S auf den Bodensalzgehalt, g/kg; EC1:5 bezieht sich auf das Boden-Wasser-Verhältnis von 1:5, die Sickerwasserleitfähigkeit, ms/c.
Insbesondere wird der genannte Schritt 2 durch das folgende Verfahren erreicht: Es wurden OLI-Daten des Satelliten Landset8, die das Untersuchungsgebiet abdeckten, von der NASA-Website bezogen und anhand der topografischen Karte von Kenli County geometrisch korrigiert, wobei ein quadratisches polynomiales Korrekturmodell und die Methode des „Nearest Neighbor Image Resampling“ verwendet wurden, um den Fehler auf ein Bildelement zu begrenzen. Diese Arbeit wurde mit der Software Envi 5.1 durchgeführt, mit der die Vektorgrenzen des Untersuchungsgebiets von der Abteilung für die Verwaltung natürlicher Ressourcen erhalten wurden. Die korrigierten Fernerkundungsbilder wurden in der ArcGIS 10.0 Plattformsoftware beschnitten, um die Fernerkundungsbilder des Untersuchungsgebiets von Kenli County zu beschneiden.
Insbesondere wird der genannte Schritt 3 durch das folgende Verfahren erreicht: Der
Schlüssel in diesem Schritt ist die Auslese des Spektralwertes des Fernerkundungsbildes, der dem Bodenprobenpunkt entspricht.
In der Software ArcGIS10.0 wird die Funktion „extret value to point“ verwendet, um den dem Bodenprobenpunkt entsprechenden Rasterwert jedes Frequenzbandes des Bildes anhand der Koordinatenposition des Probenpunktes als den dem Bodenprobenpunkt entsprechenden Fernerkundungsreflexionsgrad zu ermitteln.
Insbesondere besteht der Hauptzweck von Schritt 4 darin, zu analysieren und zu filtern, welche Fernerkundungsbänder oder transformierten Bandspektralwerte Form und Bodensalzgehalt eine enge Korrelation aufweisen.
Die Korrelationsanalyse ergab, dass die Korrelation zwischen den Spektralwerten der einzelnen Bänder und dem Salzgehalt gering ist, was darauf hindeutet, dass die Fernerkundungsspektralwerte der unveränderten einzelnen Bänder weniger empfindlich auf den Bodensalzgehalt reagieren.
Daher ist es notwendig, entsprechende mathematische Änderungen an den fernerkundeten Spektralwerten vorzunehmen, einschließlich der inversen Variation(1/R), der logarithmischen Variation(IgR). Die Korrelationsanalyse der transformierten Bandenspektren mit dem Bodensalzgehalt zeigt, dass die Korrelation verbessert wird und ein äußerst signifikanter Zusammenhang zwischen Fernerkundungsdaten und Bodensalzgehalt besteht.
Dies wurde durch die folgenden Verfahren erreicht: Schritt 4.1: Die entsprechende mathematische Transformation der fernerkundeten Spektralwerte der extrahierten Probepunkte, einschließlich inverser Variation, logarithmischer Variation, unter Verwendung von ArcGIS 10.0; Schritt 4.2: Die SPSS-Software wurde verwendet, um Korrelationen zwischen dem
Bodensalzgehalt der Probenahmestellen und den Spektralwerten der Frequenzbänder und den transformierten Spektralformen zu erhalten.
Der Korrelationskoeffizient wurde verwendet, um die Nähe zwischen zwei oder zwei oder mehr Variablengruppen mit Korrelation zu messen, und der Korrelationskoeffizient wurde mit der einfachen Korrelationskoeffizienten-Methode von „Pearson“ verwendet.
Je größer die absolute Zahl des Korrelationskoeffizienten ist, desto stärker ist die Korrelation.
Das heißt, je näher der Korrelationskoeffizient an 1 oder -1 liegt, desto stärker ist die Korrelation; je näher der Korrelationskoeffizient an 0 liegt, desto schwächer ist die Korrelation.
Insbesondere wird der genannte Schritt 5 durch das folgende Verfahren erreicht: Schritt 5.1: Basierend auf den Fernerkundungsvariablen, die in Schritt 4 gescreent wurden, wurde eine multiple lineare Regressionsanalyse von 1/R, 1g(R) der Reflexion in jedem Frequenzband mit dem in Schritt 1 erhaltenen Bodensalzgehalt mittels SPSS-Software durchgeführt, um ein Regressionsmodell zu erstellen.
Schritt 5.2: Um das erstellte Modell zu testen, wurde zuerst der Test auf Homogenität der Varianz (d.h.
F-Test) durchgeführt, und wenn sig<0,01 ist, bedeutet dies, dass der Test bestanden wurde und für die Regressionsanalyse gültig ist.
Dann wird der Präzisionstest durchgeführt, um die Präzision und Stabilität des Modells zu prüfen, um später ein Modell mit besserer Präzision und Stabilität auszuwählen.
Um die Genauigkeit der Gleichung zu testen, werden der Determinationskoeffizient (d. h.
R?) und die mittlere quadratische Abweichung (d. h.
RMSE) verwendet.
Je größer R? und je kleiner RMSE, desto höher ist die Genauigkeit des Modells und desto besser ist die Stabilitét.
Insbesondere wird der genannte Schritt 6 durch das folgende Verfahren erreicht: Schritt 6.1: Die empfindlichen Fernerkundungsvariablen wurden zunächst gemäß Schritt 4 gescreent.
Fine multiple lineare Regressionsanalyse des Logarithmus lg(R) der Reflektion jeder aus den Fernerkundungsdaten erhaltenen Frequenzbänder mit dem in Schritt 1 erhaltenen Bodensalzgehalt wurde mit der SPSS-Software durchgeführt, um eine Koeffizienten-Testtabelle für das Modell zu erhalten; Schritt 6.2: Das erstellte Modell wurde getestet.
F ist 8,452, sig ist 0,005, sig<0,01, was bedeutet, dass das Modell den Test bestanden hat und das Modell gültig ist; der Determinationskoeffizient des Modells ist 0,304 und die mittlere quadratische Abweichung ist 6,23; Schritt 6.3: Die empfindlichen Fernerkundungsvariablen wurden gemäß Schritt 4 gescreent.
Eine multiple lineare Regressionsanalyse des Kehrwerts 1/R der Reflektion jeder aus den Fernerkundungsdaten erhaltenen Frequenzbänder mit dem in Schritt 1 erhaltenen Bodensalzgehalt wurde mit der SPSS-Software durchgeführt, um eine Koeffizienten-Testtabelle für das Modell zu erhalten; Schritt 6.4: Das erstellte Modell wurde getestet.
F ist 7,89, sig ist 0,006, sig<0,01, was bedeutet, dass das Modell den Test bestanden hat und das Modell gültig ist; der Determinationskoeffizient des Modells ist 0,295 und die mittlere quadratische Abweichung ist 6,25; Schritt 6.5: Der Vergleich verschiedener Modellgenauigkeiten ergab, dass das gemischte
Fernerkundungsmodell, das auf dem Kehrwert 1/R und Logarithmus Ig(R) des Reflexionsgrads der einzelnen Frequenzbänder basiert, die höchste Genauigkeit aufweist, und das beste Modell ist Y1= 59,94-43,71* Log(B5) +7105,37* 1/(B11), das für die Überwachung der Bodensalzgehalt-Inversion in der Delta Region des Gelben Flusses verwendet werden kann.
Diese Erfindung stellt eine Landsat8-basierte Inversionsmethode für den Bodensalzgehalt im Delta des Gelben Flusses zur Verfügung, die Landsat8 OLI-Bilder zur Inversion des Bodensalzgehalts in der Region des Gelben Flusses verwendet. Drei spektrale Indikatoren, nämlich die spektrale Reflexion, der Kehrwert der spektralen Reflexion und der Logarithmus der spektralen Reflexion, wurden für die Korrelationsanalyse mit den Bodensalzgehaltsdaten ausgewählt, und die spektralen Indikatoren mit der höchsten Korrelation wurden ausgewählt. Eine multiple lineare Regressionsanalyse wurde mit den Bändern durchgeführt, die das Niveau einer signifikanten Korrelation erreichten, um ein Regressionsmodell zu erstellen. Die Ergebnisse zeigten, dass BS und B11 der OLI-Bilder alle einen signifikanten Korrelationsgrad mit dem gemessenen Bodensalzgehalt erreichten, wobei die inverse Transformation am signifikantesten war; das Regressionsmodell, das den Wert der inversen Transformation des 5. Bandes und den Bodensalzgehalt verwendete, war grundsätzlich in der Lage, die Vorhersage zu erfiillen (R? von 0,081 und RMSE von 6,229). Das Inversmodell wurde verwendet, um den Bodensalzgehalt zu invertieren und eine deskriptive statistische Analyse durchzuführen. Es wurde mit dem gemessenen Bodensalzgehalt verglichen, und es wurde festgestellt, dass der Salzgehalt im Wesentlichen konsistent war, was darauf hindeutet, dass das Inversmodell eine hohe Genauigkeit und gute Stabilität aufweist.
Illustration Abbildung 1 zeigt eine Fernerkundungsinversionskarte des Bodensalzgehalts, die mit einem Verfahren der vorliegenden Erfindung auf der Grundlage von Landsat8-Fernerkundungsbilddaten zur Inversion des Bodensalzgehalts im Delta des Gelben Flusses erhalten wurde.
Dabei stehen 1, 2, 3, 4 und 5 für fünf Salzgehaltsstufen, nämlich nicht salzhaltiger, leicht salzhaltiger, mäBig salzhaltiger, stark salzhaltiger und salzhaltiger Boden.
Spezielle Wege zur Umsetzung In dieser Erfindung werden Landsat8 OLI-Bilder und im Labor gemessene Bodensalzgehalte als Datenquellen verwendet, um ein Modell zur Schätzung des Bodensalzgehalts in der Region des Gelben Flussdeltas zu erstellen.
Dabei werden mathematische und statistische Analysemethoden verwendet, um die empfindlichen Bänder des Bodensalzgehalts und das beste Inversionsmodell zu untersuchen.
Dies bietet eine gewisse technische Unterstützung für die rechtzeitige und genaue Frfassung des Bodensalzgehalts in der Region in der Zukunft und stellt eine Grundlage für das Management, die Vorbeugung und Kontrolle des Bodensalzgehalts und die nachhaltige Entwicklung in der Region des Gelben Flussdeltas dar.
Das Verfahren zur Inversion des Bodensalzgehalts, das durch die vorliegende Erfindung bereitgestellt wird, ist wie folgt spezifiziert: Schritt 1: Sammlung und Erwerbung von Bodensalzgehaltsdaten.
Damit die Probenahmestellen die tatsächliche Situation im Untersuchungsgebiet genau widerspiegeln, müssen der Bodentyp, der Vegetationstyp, die Art der Landnutzung und andere Faktoren im
Untersuchungsgebiet berücksichtigt werden.
Dies wird insbesondere durch die folgenden Methoden erreicht: Schritt 1.1: Zunächst wurde ein Rasterlayout erstellt und an die Landnutzungsart sowie den Straßenzugang angepasst.
Insgesamt wurden 86 Probenahmestellen angelegt, und um jede Probenahmestelle herum wurden 3-6 Flecken Boden mit einer Schaufel in der Tiefe von 0-20 cm gesammelt und gründlich gemischt.
Dann wurden 1 kg bis 2 kg in Plastikbeutel gefüllt, beschriftet und die Koordinaten der Probenahmestellen, der Zeitpunkt der Probenahme, die Wetterbedingungen, die Landnutzungsart und der Bepflanzungstyp aufgezeichnet; Schritt 1.2: Die gesammelten Bodenproben wurden in Innenräumen getrocknet.
Nach dem Trocknen wurden die Graswurzeln ausgewählt und die Bodenproben wurden gemahlen und durch ein Imm-Sieb gegeben.
Zur Aufbewahrung für spätere Experimente wurden sie in Weithalsgläsern konserviert; Schritt 1.3: Der in Schritt 1.2 erhaltene Boden wurde mit der Leitfähigkeits-Methode bestimmt.
Dies wurde wie folgt durchgeführt: Eine entsprechende Menge an Bodenprobe und destilliertem Wasser wurde gemäß dem Gewichtsverhältnis Boden-Wasser von 1:5 in ein Glasreagenzglas gefüllt und 5 Minuten lang geschüttelt.
Anschließend wurde das Filtrat filtriert und die Leitfähigkeit mit einem Leitfähigkeitsmessgerät gemessen.
Nachdem die Leitfähigkeit bestimmt wurde, wurden die Daten zur Bodensalinität durch Umrechnung der Formel entsprechend der Beziehung zwischen Leitfähigkeit und Bodensalinität ermittelt.
Die spezifische Umrechnungsformel lautete wie folgt:
S=3,047EC1:5-0,493 (r=0,981””, p<0,001). In der Formel bezieht sich S auf den Bodensalzgehalt, g/kg; EC1:5 bezieht sich auf das Boden-Wasser-Verhältnis von 1:5, die Sickerwasserleitfähigkeit, ms/c. Die in Schritt 1.3 abgeleiteten Daten werden mit dem K-S-Test auf Übereinstimmung mit der Normalverteilung geprüft. Diejenigen, die nicht der Normalverteilung entsprechen, müssen vorverarbeitet werden, damit sie der Normalverteilung entsprechen, was die Genauigkeit des Modells bis zu einem gewissen Grad verbessern kann. Schritt 2: Erwerbung und Verarbeitung von Fernerkundungsbilddaten. Dies wird insbesondere durch die folgenden Methoden erreicht: Es wurden OLI-Daten des Satelliten Landset8, die das Untersuchungsgebiet abdeckten, von der NASA-Website bezogen und anhand der topografischen Karte von Kenli County geometrisch korrigiert, wobei ein quadratisches polynomiales Korrekturmodell und die Methode des „Nearest Neighbor Image Resampling” verwendet wurden, um den Fehler auf ein Bildelement zu begrenzen. Diese Arbeit wurde mit der Software Envi 5.1 durchgeführt, mit der die Vektorgrenzen des Untersuchungsgebiets von der Abteilung für die Verwaltung natürlicher Ressourcen erhalten wurden. Die korrigierten Fernerkundungsbilder wurden in der ArcGIS 10.0 Plattformsoftware beschnitten, um die Fernerkundungsbilder des Untersuchungsgebiets von Kenli County zu beschneiden.
Schritt 3: Auslese und Berechnung der Reflexionsdaten für die jeweiligen Frequenzbänder aus Landsat8 OLI-Bildern mit der Software ArcGIS10.0. Dies wird insbesondere durch die folgenden Methoden erreicht: Der Schlüssel in diesem Schritt ist die Auslese des Spektralwertes des Fernerkundungsbildes, der dem Bodenprobenpunkt entspricht. In der
Software ArcGIS10.0 wird die Funktion ,extrct value to point“ verwendet, um den dem Bodenprobenpunkt entsprechenden Rasterwert jedes Frequenzbandes des Bildes anhand der Koordinatenposition des Probenpunktes als den dem Bodenprobenpunkt entsprechenden Fernerkundungsreflexionsgrad zu ermitteln, und den lg(R) des Reflexionsgrads jedes Bandes zu erhalten.
Schritt 4: Vorläufiges Screening empfindlicher Frequenzbänder durch entsprechende 1g(R)- und 1/R-Änderungen für jedes Frequenzband, gefolgt von einer Korrelationsanalyse des fernerkundeten Reflexionsgrads und des transformierten fernerkundeten Reflexionsgrads mit dem jeweiligen Salzgehalt des Bodens mittels SPSS-Software. Der Hauptzweck dieses Schritts ist es, zu analysieren und zu filtern, welche Fernerkundungsbänder oder transformierten Bandspektralwerte Form und Bodensalzgehalt eine enge Korrelation aufweisen. Die Korrelationsanalyse ergab, dass die Korrelation zwischen den Spektralwerten der einzelnen Bänder und dem Salzgehalt gering ist, was darauf hindeutet, dass die Fernerkundungsspektralwerte der unveränderten einzelnen Bänder weniger empfindlich auf den Bodensalzgehalt reagieren. Daher ist es notwendig, entsprechende mathematische Änderungen an den fernerkundeten Spektralwerten vorzunehmen, einschließlich der inversen Variation(1/R), der logarithmischen Variation(lgR). Die Korrelationsanalyse der transformierten Bandenspektren mit dem Bodensalzgehalt zeigt, dass die Korrelation verbessert wird und ein äußerst signifikanter Zusammenhang zwischen Fernerkundungsdaten und Bodensalzgehalt besteht. Dies wurde durch die folgenden Methoden erreicht: Schritt 4.1: Die entsprechende mathematische Transformation der fernerkundeten
Spektralwerte der extrahierten Probepunkte, einschließlich inverser Variation, logarithmischer Variation, unter Verwendung von ArcGIS 10.0;
Schritt 4.2: Die SPSS-Software wurde verwendet, um den Bodensalzgehalt mit den Spektralwerten jedes Bandes des Bildes zu korrelieren, wodurch die Korrelation zwischen dem Salzgehalt der Probenahmestellen und den Spektralwerten jedes Bandes und deren transformierten Formen erhalten wurde.
In Tabelle 1 sind die Ergebnisse dargestellt.
Korrelationskoeffizienten werden verwendet, um die Nähe zwischen zwei oder zwei oder mehr Variablengruppen mit Korrelation zu messen, und die wichtigsten Korrelationskoeffizienten sind der einfache Korrelationskoeffizient von Pearson, der konsistente Korrelationskoeffizient von Kendall tua-b und die einfache Korrelationskoeffizientenmethode von Spearman et al.
In der vorliegenden Erfindung wird der einfache Korrelationskoeffizient nach Pearson verwendet.
Je größer der Absolutwert des Korrelationskoeffizienten R ist, desto stärker ist die Korrelation.
Je näher der Korrelationskoeffizient an 1 oder -1 liegt, desto stärker ist die Korrelation, und je näher der Korrelationskoeffizient an 0 liegt, desto schwächer ist die Korrelation.
= a 5 = =
8. £ 8 5 > 2 =
CM > > > ce ba = or S > = + =] a = eo.
& fa = © 7 LR ss = > nn 5 & = = er =] ht) © > 2 > on > = = = S Aus Tabelle 1 geht hervor, dass der Absolutwert des Korrelationskoeffizienten maximal 0,14 und minimal O beträgt. Die Korrelationskoeffizienten liegen alle nahe bei 0, so dass die Korrelation extrem schwach ist, was darauf hindeutet, dass die Spektren der einzelnen Banden gering und unempfindlich gegenüber der Salzreflexion sind. Daher müssen zum Aufbau von Fernerkundungs-Inversionsmodellen aus Bildern entsprechende mathematische Anderungen an der spektralen Reflexion vorgenommen werden, einschließlich der inversen Variante (1/R), der logarithmischen Variante (log(R)). Nachdem die entsprechenden 1g(R), 1/R Varianten für jede Bande vorgenommen wurden, wurde die transformierte spektrale Reflexion dann mit der SPSS-Software mit dem Bodensalzgehalt korreliert, um die Korrelationskoeffizienten zu erhalten Tabelle 2.
—— r= = wd Ya 2 ® Bh = = [= = + El a = — + A La © sd 2 5 © fre) oi = © te 5 Lys 1 = td S = + Lai = 2 S ax + S ta ~ pi = © 3 = & = Es 5 La ç . 2 A ed = be 2 oo a, a © EN ve] = bm om = = S “D + i La So = By ta a
HS = ue I N. = D td iG = < pe 5 x æ = ted Ris oh = = 2 2 Ch La = > el D = © = & = 2 = > to Ed rs De = S “J = ES & m Ri Ë À ps 3 pe
Hinweis: ** steht für eine signifikante Korrelation auf dem 0,01-Niveau, * für eine signifikante Korrelation auf dem 0,05-Niveau.
Wie aus Tabelle 2 hervorgeht, zeigt die Korrelation zwischen 1/R, 1g(R) und dem Bodensalzgehalt in jedem Band der OLI-Bilder, dass der Korrelationskoeffizient zwischen der transformierten spektralen Bodenreflexion und dem Bodensalzgehalt insgesamt zugenommen hat und enger korreliert ist. Bei der inversen Variation hatte der Pearson-Korrelationskoeffizient (R) einen maximalen Absolutwert von 0,304 und einen minimalen Wert von 0. Bei der logarithmischen Variation hatte der Pearson-Korrelationskoeffizient (R) einen maximalen Absolutwert von 0,295 und einen minimalen Wert von 0. Auf dem Niveau P<0,01 waren nach der 1/R, 1g(R)-Variation die 1/R, lg(R) jeder Bande der OLI-Bilder signifikant mit dem Bodensalzgehalt in den B5- und B11-Banden korreliert. Bei der inversen Variation lagen die signifikanten Korrelationskoeffizienten bei -0,304 bzw. -0,238. Bei der logarithmischen Variation lagen die signifikanten Korrelationskoeffizienten bei 0,295 bzw. 0,239.
In der Studie wurde festgestellt, dass durch bestimmte mathematische Anderungen an den Hyperspektraldaten zum Bodensalzgehalt die empfindliche Position der Kurve hervorgehoben werden kann, was die Korrelation zwischen Bodensalzgehalt und spektraler Reflexion erleichtert und die quantitative Inversion des Bodensalzgehalts aus den Bildern weiter verbessert. Daher wurden in dieser Studie die inversen und logarithmischen Anderungen in der spektralen Reflexion der einzelnen Bänder der OLI-Bilder vorgenommen, so dass die spektrale Reflexion der Bänder erhöht wurde und somit die Korrelation signifikant war. Die
Größe und Signifikanz der Korrelationskoeffizienten wurden verwendet, um die logarithmierten, invers-transformierten Werte von Band 5 (B5) und Band 11 (B11) des Bildes zu filtern, die als die fernerkundungsbandsensitiven Variablen für den nächsten Schritt der Erstellung des Fernerkundungsinversionsmodells verwendet werden.
Schritt 5: Multiple lineare Regressionsanalyse von 1/R, 1g(R) für jede Bandenreflexion mit dem in Schritt 1 erhaltenen Bodensalzgehalt mittels SPSS-Software, um ein Regressionsmodell zu erstellen und zu testen, ob die Reflexionsänderung der einzelnen Bandenwerte mit dem Bodensalzgehalt modelliert werden kann.
Dies wurde durch die folgenden Methoden erreicht: Schritt 5.1: Basierend auf den empfindlichen Fernerkundungsvariablen, die zunächst in Schritt 4 gescreent wurden, wurde eine multiple lineare Regressionsanalyse des log 1g(R), invers 1/R des Reflexionsgrads jeder aus den Fernerkundungsdaten erhaltenen Bande mit dem in Schritt 1 erhaltenen Bodensalzgehalt mit der SPSS-Software durchgeführt.
Daraus ergab sich der Koeffiziententest des Modells Tabelle 3.
Tabelle 3 Koeffizient t Sig 59,94 1,75 Konstante 0,08 Log (B5) -43,71 -3,06 0,03 1/ (B11) 7105,37 2,03 0,04 Die t-Statistik ist der Test für die Signifikanz des Koeffizienten; der sig-Wert ist der Wahrscheinlichkeitswert, der der t-Statistik entspricht, so dass t und sig äquivalent sind und es genügt, sig zu betrachten.
Der sig-Wert muss kleiner als das angegebene Signifikanzniveau sein, in der Regel 0,05, 0,01 usw.
Je näher der sig-Wert bei 0 liegt, desto besser.
Aus Tabelle 3 lässt sich ablesen, dass sig nahe bei 0 liegen, so dass das Ergebnis akzeptabel ist.
Schreiben Sie die Gleichung gemäß Tabelle 3, wie in Gleichung (1) unten: Y1= 59,94-43,71* Log(B5) + 7105,37* 1/(B11) (1) In der obigen Gleichung ist Log(B5) der logarithmisch transformierte Wert der Reflektivität des 5. Bandes; 1/(B11) ist der invers transformierte Wert der Reflektivität des 11.
Schritt 5.2: Um das erstellte Modell zu testen, wurde zuerst der Test auf Homogenität der Varianz (d.h.
F-Test) durchgeführt wurde, um Tabelle 4 zu erhalten., und wenn sig<0,05 ist, bedeutet dies, dass der Test bestanden wurde und für die Regressionsanalyse gültig ist.
Dann wird der Präzisionstest durchgeführt, um die Präzision und Stabilität des Modells zu prüfen, um später ein Modell mit besserer Präzision und Stabilität auszuwählen.
Um die Genauigkeit der Gleichung zu testen, werden der Determinationskoeffizient (d. h.
R°) und die mittlere quadratische Abweichung (d. h.
RMSE) verwendet.
Je größer R? und je kleiner RMSE, desto höher ist die Genauigkeit des Modells und desto besser ist die Stabilität.
Tabelle 4 R? RMSE sig F 0,322 6,227 0,011 4,757 In Tabelle 4 ist F 4,757, sig ist 0,011 und sig < 0,05. Das Modell besteht den F-Test und ist für die Regressionsanalyse gültig.
Das R? des Modells ist 0,322 und der RMSE beträgt 6,227. Schritt 6: Erstellung von Fernerkundungs-Inversionsmodellen des Bodensalzgehalts für log-transformierte bzw. invers-transformierte Fernerkundungs-Bilddaten zum Vergleich des nach dieser einzigen Methode erstellten Fernerkundungs-Bilddatenmodells mit dem Fernerkundungs-Bilddatenmodell, das nach dem Mischen der beiden Methoden-Transformationen erstellt wurde, und zur Analyse der Auswirkungen der verschiedenen Datentransformationen auf die Modellgenauigkeit.
Schritt 6.1: Die empfindlichen Fernerkundungsvariablen wurden zunächst gemäß Schritt 4 gescreent. Fine multiple lineare Regressionsanalyse des Logarithmus lg(R) der Reflektion jeder aus den Fernerkundungsdaten erhaltenen Frequenzbänder mit dem in Schritt 1 erhaltenen Bodensalzgehalt wurde mit der SPSS-Software durchgeführt, um den Koeffiziententest des Modells Tabelle 5 zu erhalten.
Tabelle 5 Koeffizient t Sig Konstante 112,19 3,01 0,003 Log (B5) 54.08 -2,91 0,005 Die sig-Werte sowohl für die Konstante als auch für die Variable sind kleiner als die angegebenen 0,01, nahe bei 0. Dies zeigt, dass die Koeffizienten des Modells akzeptabel sind. Die Gleichung wird gemäß Tabelle 5 wie folgt geschrieben (2): Y2= 112,19 - 54,08* Log(B5) (2) In der obigen Gleichung ist Log(B5) der logarithmisch transformierte Wert der Reflektivität des 5. Bandes. Schritt 6.2: Das erstellte Modell wurde getestet. F ist 8,452, sig ist 0,005, sig<0,01, was bedeutet, dass das Modell den Test bestanden hat und das Modell gültig ist. Das R° des Modells ist 0,304 und der RMSE beträgt 6,23.
Schritt 6.3: Die empfindlichen Fernerkundungsvariablen wurden gemäß Schritt 4 gescreent.
Eine multiple lineare Regressionsanalyse des Kehrwerts 1/R der Reflektion jeder aus den Fernerkundungsdaten erhaltenen Frequenzbänder mit dem in Schritt 1 erhaltenen Bodensalzgehalt wurde mit der SPSS-Software durchgeführt, um den Koeffiziententest des Modells Tabelle 6 zu erhalten.
Tabelle 6 Koeffizient t Sig Konstante -17,99 -2,30 0,024 1/ (B11) 2193,07 -2,81 0,006 Die sig-Werte sowohl der Konstanten als auch der Variablen sind kleiner als die angegebenen 0,05, was darauf hinweist, dass die Koeffizienten des Modells akzeptabel sind und die Anforderungen erfüllen.
Die Gleichung wird gemäß Tabelle 5 wie folgt geschrieben (3): Y3= -17,99+2193,07* 1/(B11) (3) In der obigen Gleichung ist 1/(B11) der invers transformierte Wert des Reflexionsgrads des 11. Bandes.
Schritt 6.4: Das erstellte Modell wurde getestet.
F ist 7,89, sig ist 0,006, sig<0,01, was bedeutet, dass das Modell den Test bestanden hat und das Modell gültig ist.
Das R? des
Modells ist 0,295 und der RMSE beträgt 6,25. Schritt 6.5: Der Vergleich verschiedener Modellgenauigkeiten ergab, dass das gemischte Fernerkundungsmodell, das auf dem Kehrwert 1/R und Logarithmus Ig(R) des Reflexionsgrads der einzelnen Frequenzbänder basiert, die höchste Genauigkeit aufweist, und das beste Modell ist Y1= 59,94-43,71* Log(B5) +7105,37* 1/(B11), das für die Überwachung der Bodensalzgehalt-Inversion in der Delta Region des Gelben Flusses verwendet werden kann.
Schritt 7: Das oben erstellte Prognosemodell wurde auf das Landsat8-Fernerkundungsbild des Untersuchungsgebiets angewendet, so dass die Landsat8-Fernerkundungs-Inversionskarte des Bodensalzgehalts entstand, wie in Abbildung 1 dargestellt; eine deskriptive statistische Analyse des Bodensalzgehalts wurde durchgeführt, so dass Tabelle 7 entsteht.
Wie in Abbildung 1 dargestellt, wies der Salzgehalt der Oberflächenböden im Untersuchungsgebiet offensichtliche räumliche Verteilungsunterschiede auf und war deutlich abgestuft.
Die Fläche des Bodens <1g/kg beträgt 37,5 km?, was 1,6 % entspricht und nicht versalzenes Land ist; die Fläche des Bodens 1-2g/kg beträgt 156,14 km?, was 6,7 % entspricht und leicht versalzenes Land ist; die Fläche des Bodens 2-4g/kg beträgt 672,87 km?, was mäßig versalzenes Land ist; die Fläche des Bodens 4-6g/kg beträgt 580,36 km?, das ist stark versalzenes Land; >6g/kg Landfläche ist 884,131 km? das ist noch versalzenes Land.
Je heller die Farbe in der inversen Darstellung ist, desto höher ist der Salzgehalt des dargestellten Bodens.
Aus Abb. 2 ist ersichtlich, dass der Bodensalzgehalt mit zunehmender Entfernung vom Bohai-Meer abnimmt.
Der höchste Bodensalzgehalt findet sich im Nordosten und in den
Küstengebieten, wo der Boden noch salzhaltig ist und eine Fläche von 884,131 km? oder 37,93 % einnimmt; weiter landeinwärts ist der Boden stark salzhaltig und nimmt eine Fläche von 580,36 km? oder 24,9 % ein; und am weitesten landeinwärts ist der Boden mäßig salzhaltig und nimmt 672,87 km? oder 6,7 % ein.
Tabelle 7 Salzgehalt/(g/kg) Messwert Inverswert Maximalwert 6,8 7,8 Mindestwert 0,2 0,6 Durchschnittswert 1,41 1,96 Standardabweichung 1,04 1,33 Schrägungsfaktor 8,02 9,00 Kurtosis-Koeffizient 8,02 8,72 Median 1,1 2,1 Variationskoeffizient 73,76% 74,56% Wie in Tabelle 7 hervorgeht, unterscheiden sich die Messwerte nicht signifikant von den Inverswerten. Daher kann die Inversionsmethode, die durch diese Erfindung bereitgestellt wird, auf die Messung des Bodensalzgehalts angewendet werden, was eine große Arbeits- und Zeitersparnis mit sich bringt. Diese Erfindung verwendet die spektrale Reflexion von Lansat8 zur Korrelation mit dem Bodensalzgehalt und wählt das Band mit signifikanter Korrelation als empfindliches Band aus. Dann wurde eine Regressionsanalyse zwischen den empfindlichen Banden und dem gemessenen Bodensalzgehalt durchgeführt und ein multiples lineares Regressionsmodell erstellt, um die Inversion der räumlichen Verteilung des Bodensalzgehalts in Kenli County zu ermöglichen. Die Ergebnisse zeigten, dass die räumliche Verteilung des Bodensalzgehalts im Untersuchungsgebiet mithilfe von OLI-Bildern recht gut abgeleitet werden konnte.Der Salzgehalt im Untersuchungsgebiet weist eine offensichtliche räumliche Verteilungsvariabilität auf, die zeigt, dass der Salzgehalt von Osten nach Westen allmählich abnimmt; je näher am Gelben Fluss, desto niedriger ist der Bodensalzgehalt; je weiter vom Meer entfernt, desto niedriger ist der Bodensalzgehalt; das östliche Küstengebiet ist vom Meerwasser beeinflusst, mit relativ hoher Grundwassermineralisierung, starker Salzakkumulation und hohem Bodensalzgehalt; das nordöstliche Naturschutzgebiet hat aufgrund der fehlenden anthropogenen offenen Nutzung eine geringe Vegetationsbedeckung sowie hohe Bodensalzgehalte. Die räumliche Variabilität der Bodenversalzung ist groß, und die Hauptgründe für die räumliche Variabilität sind natürliche Faktoren wie Topographie, Klima, bodenbildendes Ausgangsmaterial und Geomorphologie.
Die empfindlichen Spektren des OLI-Bildes für den Salzgehalt des Bodens im Untersuchungsgebiet sind B5 (845-885 nm) und B11 (11500-12510 nm), wobei die B5-Bande die beste Bande darstellt. Die Banden wurden invertiert-transformiert und log-transformiert sowie einzeln der Korrelationsanalyse unterzogen, was die Korrelation mit dem Salzgehalt deutlich verbessern kann, wobei die invertiert-transformierte die deutlichste ist. Mit Hilfe statistischer Methoden wurde ein multiples lineares Regressionsmodell erstellt, und das beste Modell wurde nach dem Prinzip ausgewählt, dass ein größeres R? besser und ein kleinerer RMSE besser ist, und schließlich wurde Y1= 59,94-43,71* Log(B5) + 7105,37* 1/(B11)
gewählt.
Die Koeffizienten des Modells wurden zunächst einem t-Test unterzogen, wobei P=0,001, das Niveau der statistischen Signifikanz auf dem Niveau P<0,01 erreicht wurde; und dann wurde das Modell auf Genauigkeit getestet.
Dieses Modell besitzt eine gute Prognosefähigkeit für den Bodensalzgehalt im Untersuchungsgebiet, wobei der R?-Wert 0,322, der RMSE-Wert 6,229, die Modellgenauigkeit gut und die Stabilität hoch ist.
In Schritt 2 der vorliegenden Erfindung ist Landsat8 der neueste Satellit der Landsat-Reihe, die OLI (mit 9 Bändern) und TIRS (mit 2 separaten thermischen Infrarotbändern) umfasst.
Ein wichtiger Vorteil der OLI-Bilder ist die Anpassung von Band 5, um das Wasserdampf-Absorptionsmerkmal bei 0,825 um auszuschließen, wodurch der Einfluss der atmosphärischen Absorption auf die Daten stark reduziert wird.
Es gibt auch zwei zusätzliche Bänder, das erste ist das blaue Band für die Küstenüberwachung, das andere ist das kurzwellige Infrarotband für die Wolkenüberwachung.

Claims (7)

PATENTANSPRÜCHE
1. Ein Verfahren zur Inversion des Bodensalzgehalts im Delta des Gelben Flusses auf der Basis von Landsat8, dadurch gekennzeichnet, dass der Schritt umfasst: Schritt 1: Sammlung und Erwerbung von Bodensalzgehaltsdaten; Schritt 2: Erwerbung und Verarbeitung von Fernerkundungsbilddaten; Schritt 3: Auslese und Berechnung der Reflexionsdaten für die jeweiligen Frequenzbänder aus Landsat8 OLI-Bildern mit der Software ArcGIS10.0; Schritt 4: Vorläufiges Screening empfindlicher Frequenzbänder durch entsprechende 1g(R)- und 1/R-Änderungen für jedes Frequenzband, gefolgt von einer Korrelationsanalyse des fernerkundeten Reflexionsgrads und des transformierten fernerkundeten Reflexionsgrads mit dem jeweiligen Salzgehalt des Bodens mittels SPSS-Software; Schritt 5: Multiple lineare Regressionsanalyse von 1/R, 1g(R) für jede Bandenreflexion mit dem in Schritt 1 erhaltenen Bodensalzgehalt mittels SPSS-Software, um ein Regressionsmodell zu erstellen und zu testen, ob die Reflexionsänderung der einzelnen Bandenwerte mit dem Bodensalzgehalt modelliert werden kann; Schritt 6: Erstellung von Fernerkundungs-Inversionsmodellen des Bodensalzgehalts für log-transformierte bzw. invers-transformierte Fernerkundungs-Bilddaten zum Vergleich des nach dieser einzigen Methode erstellten Fernerkundungs-Bilddatenmodells mit dem Fernerkundungs-Bilddatenmodell, das nach dem Mischen der beiden Methoden-Transformationen erstellt wurde, und zur Analyse der Auswirkungen der verschiedenen Datentransformationen auf die Modellgenauigkeit;
Schritt 7: Das oben erstellte beste Prognosemodell wurde auf Landsat8-Fernerkundungsbilder des Untersuchungsgebiets angewendet, wodurch Landsat8-Fernerkundungs-Inversionskarten des Bodensalzgehalts erstellt wurden.
2. Das Landsat8-basierende Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses gemäß Anspruch 1, dadurch gekennzeichnet, dass der Schritt 1 speziell durch folgende Methoden erreicht wird: Schritt 1.1: Zunächst wurde ein Rasterlayout erstellt und an die Landnutzungsart sowie den Straßenzugang angepasst; insgesamt wurden 86 Probenahmestellen angelegt, und um jede Probenahmestelle herum wurden 3-6 Flecken Boden mit einer Schaufel in der Tiefe von 0-20 cm gesammelt und gründlich gemischt; dann wurden 1 kg bis 2 kg in Plastikbeutel gefüllt, beschriftet und die Koordinaten der Probenahmestellen, der Zeitpunkt der Probenahme, die Wetterbedingungen, die Landnutzungsart und der Bepflanzungstyp aufgezeichnet; Schritt 1.2: Die gesammelten Bodenproben wurden in Innenräumen getrocknet; nach dem Trocknen wurden die Graswurzeln ausgewählt und die Bodenproben wurden gemahlen und durch ein Imm-Sieb gegeben; zur Aufbewahrung für spätere Experimente wurden sie in Weithalsgläsern konserviert; Schritt 1.3: Der in Schritt 1.2 erhaltene Boden wurde mit der Leitfähigkeits-Methode bestimmt; Dies wurde wie folgt durchgeführt: Eine entsprechende Menge an Bodenprobe und destilliertem Wasser wurde gemäß dem Gewichtsverhältnis Boden-Wasser von 1:5 in ein Glasreagenzglas gefüllt und 5 Minuten lang geschüttelt; anschließend wurde das Filtrat filtriert und die Leitfähigkeit mit einem Leitfähigkeitsmessgerät gemessen; nachdem die Leitfähigkeit bestimmt wurde, wurden die Daten zur Bodensalinität durch Umrechnung der Formel entsprechend der Beziehung zwischen Leitfähigkeit und Bodensalinität ermittelt; die spezifische Umrechnungsformel lautete wie folgt: S=3,047EC1:5-0,493 (r=0,981"", p<0,001). In der Formel bezieht sich S auf den Bodensalzgehalt, g/kg; EC1:5 bezieht sich auf das Boden-Wasser-Verhältnis von 1:5, die Sickerwasserleitfähigkeit, ms/c.
3. Das Landsat8-basierende Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses gemäß Anspruch 1, dadurch gekennzeichnet, dass der Schritt 2 speziell durch folgende Methoden erreicht wird: Es wurden OLI-Daten des Satelliten Landset8, die das Untersuchungsgebiet abdeckten, von der NASA-Website bezogen und anhand der topografischen Karte von Kenli County geometrisch korrigiert, wobei ein quadratisches polynomiales Korrekturmodell und die Methode des „Nearest Neighbor Image Resampling” verwendet wurden, um den Fehler auf ein Bildelement zu begrenzen; diese Arbeit wurde mit der Software Envi 5.1 durchgeführt, mit der die Vektorgrenzen des Untersuchungsgebiets von der Abteilung für die Verwaltung natürlicher Ressourcen erhalten wurden; die korrigierten Fernerkundungsbilder wurden in der ArcGIS 10.0 Plattformsoftware beschnitten, um die Fernerkundungsbilder des Untersuchungsgebiets von Kenli County zu beschneiden.
4. Das Landsat8-basierende Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses gemäß Anspruch 3, dadurch gekennzeichnet, dass der Schritt 3 speziell durch folgende Methoden erreicht wird: der Schlüssel in diesem Schritt ist die Auslese des Spektralwertes des Fernerkundungsbildes, der dem Bodenprobenpunkt entspricht; in der Software ArcGIS10.0 wird die Funktion ,extrct value to point“ verwendet, um den dem Bodenprobenpunkt entsprechenden Rasterwert jedes Frequenzbandes des Bildes anhand der Koordinatenposition des Probenpunktes als den dem Bodenprobenpunkt entsprechenden Fernerkundungsreflexionsgrad zu ermitteln.
5. Das Landsat8-basierende Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses gemäß Anspruch 4, dadurch gekennzeichnet, dass der Schritt 4 speziell durch folgende Methoden erreicht wird: Schritt 4.1: Die entsprechende mathematische Transformation der fernerkundeten Spektralwerte der extrahierten Probepunkte, einschließlich inverser Variation, logarithmischer Variation, unter Verwendung von ArcGIS 10.0; Schritt 4.2: Die SPSS-Software wurde verwendet, um Korrelationen zwischen dem Bodensalzgehalt der Probenahmestellen und den Spektralwerten der Frequenzbänder und den transformierten Spektralformen zu erhalten; der Korrelationskoeffizient wurde verwendet, um die Nahe zwischen zwei oder zwei oder mehr Variablengruppen mit Korrelation zu messen, und der Korrelationskoeffizient wurde mit der einfachen Korrelationskoeffizienten-Methode von Pearson verwendet; je größer die absolute Zahl des Korrelationskoeffizienten ist, desto stärker ist die Korrelation; das heißt, je näher der Korrelationskoeffizient an 1 oder -1 liegt, desto stärker ist die Korrelation; je näher der Korrelationskoeffizient an O liegt, desto schwächer ist die Korrelation.
6. Das Landsat8-basierende Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses gemäß Anspruch 5, dadurch gekennzeichnet, dass der Schritt 5 speziell durch folgende Methoden erreicht wird: Schritt 5.1: Basierend auf den Fernerkundungsvariablen, die in Schritt 4 gescreent wurden, wurde eine multiple lineare Regressionsanalyse von 1/R, 1g(R) der Reflexion in jedem Frequenzband mit dem in Schritt 1 erhaltenen Bodensalzgehalt mittels SPSS-Software durchgeführt, um ein Regressionsmodell zu erstellen; Schritt 5.2: Um das erstellte Modell zu testen, wurde zuerst der F-Test durchgeführt, und wenn sig<0,01 ist, bedeutet dies, dass der Test bestanden wurde und für die Regressionsanalyse gültig ist; dann wird der Präzisionstest durchgeführt, um die Präzision und Stabilität des Modells zu prüfen, um später ein Modell mit besserer Präzision und Stabilität auszuwählen; um die Genauigkeit der Gleichung zu testen, werden der Determinationskoeffizient und die mittlere quadratische Abweichung verwendet; je größer der Determinationskoeffizient und je kleiner die mittlere quadratische Abweichung, desto höher ist die Genauigkeit des Modells und desto besser ist die Stabilität.
7. Das Landsat8-basierende Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses gemäß Anspruch 6, dadurch gekennzeichnet, dass der Schritt 6 speziell durch folgende Methoden erreicht wird: Schritt 6.1: Die empfindlichen Fernerkundungsvariablen wurden zunächst gemäß Schritt 4 gescreent; Eine multiple lineare Regressionsanalyse des Logarithmus lg(R) der Reflektion jeder aus den Fernerkundungsdaten erhaltenen Frequenzbänder mit dem in Schritt 1 erhaltenen Bodensalzgehalt wurde mit der SPSS-Software durchgeführt, um eine Koeffizienten-Testtabelle für das Modell zu erhalten;
Schritt 6.2: Das erstellte Modell wurde getestet; F ist 8,452, sig ist 0,005, sig<0,01, was bedeutet, dass das Modell den Test bestanden hat und das Modell gültig ist; der Determinationskoeffizient des Modells ist 0,304 und die mittlere quadratische Abweichung ist 6,23;
Schritt 6.3: Die empfindlichen Fernerkundungsvariablen wurden gemäß Schritt 4 gescreent; Eine multiple lineare Regressionsanalyse des Kehrwerts 1/R der Reflektion jeder aus den Fernerkundungsdaten erhaltenen Frequenzbänder mit dem in Schritt 1 erhaltenen Bodensalzgehalt wurde mit der SPSS-Software durchgeführt, um eine Koeffizienten-Testtabelle für das Modell zu erhalten;
Schritt 6.4: Das erstellte Modell wurde getestet; F ist 7,89, sig ist 0,006, sig<0,01, was bedeutet, dass das Modell den Test bestanden hat und das Modell gültig ist; der Determinationskoeffizient des Modells ist 0,295 und die mittlere quadratische Abweichung ist 6,25;
Schritt 6.5: Der Vergleich verschiedener Modellgenauigkeiten ergab, dass das gemischte Fernerkundungsmodell, das auf dem Kehrwert 1/R und Logarithmus Ig(R) des Reflexionsgrads der einzelnen Frequenzbänder basiert, die höchste Genauigkeit aufweist, und das beste Modell ist Y1= 59,94-43,71* Log(B5) +7105,37* 1/(B11), das für die Überwachung der Bodensalzgehalt-Inversion in der Delta Region des Gelben Flusses verwendet werden kann.
LU500138A 2020-06-19 2021-05-07 Landsat8-basierte Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses LU500138B1 (de)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010567384.0A CN111783288B (zh) 2020-06-19 2020-06-19 基于Landsat8对黄河三角洲土壤盐分的反演方法

Publications (2)

Publication Number Publication Date
LU500138A1 LU500138A1 (de) 2021-12-20
LU500138B1 true LU500138B1 (de) 2022-01-06

Family

ID=72757435

Family Applications (1)

Application Number Title Priority Date Filing Date
LU500138A LU500138B1 (de) 2020-06-19 2021-05-07 Landsat8-basierte Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses

Country Status (2)

Country Link
CN (1) CN111783288B (de)
LU (1) LU500138B1 (de)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113011372B (zh) * 2021-04-01 2022-08-09 清华大学 一种盐碱化土地自动监测和识别方法
CN113435640B (zh) * 2021-06-24 2022-07-26 中国科学院东北地理与农业生态研究所 预测苏打盐碱地水稻主要生育期不同耕层土壤原位ec的方法
CN114062439B (zh) * 2021-11-18 2022-08-23 浙江大学 一种利用时间序列遥感影像联合估算土壤剖面盐分的方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10107770B2 (en) * 2014-06-18 2018-10-23 Texas Tech University System Portable apparatus for soil chemical characterization
CN106918566A (zh) * 2017-03-09 2017-07-04 长江勘测规划设计研究有限责任公司 基于高光谱数据的土壤水分、盐分信息联合提取方法
CN109342337A (zh) * 2018-12-19 2019-02-15 山东农业大学 一种重度盐渍土盐分获取方法、系统及装置
CN109738380B (zh) * 2019-01-25 2022-09-30 西北农林科技大学 一种土壤盐渍化程度的高光谱遥感判断方法

Also Published As

Publication number Publication date
LU500138A1 (de) 2021-12-20
CN111783288A (zh) 2020-10-16
CN111783288B (zh) 2022-12-02

Similar Documents

Publication Publication Date Title
LU500138B1 (de) Landsat8-basierte Inversionsverfahren für den Bodensalzgehalt im Delta des Gelben Flusses
Teshaev et al. The soil-adjusted vegetation index for soil salinity assessment in Uzbekistan
DE102020122010B4 (de) Verfahren und system zur bestimmung einer beschaffenheit einer geografischen linie
DE102020122052B3 (de) Verfahren und system zur kodierung von bilddaten, die eine geografische linie repräsentieren
CN105809148A (zh) 基于遥感时空谱融合的作物干旱识别及风险评估方法
Zarei et al. Monitoring meteorological drought in Iran using remote sensing and drought indices
CN109342697A (zh) 基于随机森林—普通克里格法的土壤有机碳含量预测方法
CH716789B1 (de) Verfahren zur Extraktion verborgener geologischer Strukturinformationen aus einem Fernerkundung-Bildpaar.
CN114819737B (zh) 公路路域植被的碳储量估算方法、系统及存储介质
Belete The impact of sedimentation and climate variability on the hydrological status of Lake Hawassa, South Ethiopia
DE102011109503B4 (de) Qualitätsabhängige Modellierung von Radar-Rückstreuung
Ludwig et al. Mesoscale water balance modelling in the Upper Danube watershed using sub-scale land cover information derived from NOAA-AVHRR imagery and GIS-techniques
EP1695029A1 (de) Verfahren und vorrichtung zur mindestens teilweise automatisierten auswertung von fernerkundungsdaten
AU2021100533A4 (en) Soil salinity at Yellow River Delta Inversion Method based on Landsat 8
Bannari et al. Potential of WorldView-3 for soil salinity modeling and mapping in an arid environment
Mandlburger et al. Topographische Daten aus Laserscanning als Grundlage für Hydrologie und Wasserwirtschaft
DE102020008065A1 (de) Verfahren und System zur Bestimmung einer Beschaffenheit einer geografischen Linie
CH717361A2 (de) Verfahren und System zum Extrahieren von Winkelanomalie-Informationen aus Fernerkundungsdaten.
Asam Potential of high resolution remote sensing data for Leaf Area Index derivation using statistical and physical models
Janalipour et al. Satellite-based monitoring of lake Urmia basin health with focus on its ecosystem regulating services by applying total ecosystem health (TEH) model
Wu et al. Dynamic Change Monitoring and Assessment for Sandy Land Based on Quantitative Remote Sensing
DE202022100123U1 (de) Intelligentes System zur Überwachung landwirtschaftlicher Felder mit IOT-Sensoren und maschinellem Lernen
WO2017050548A1 (de) Verfahren zur bestimmung bodenartspezifischer parameter durch erfassung der bodenfeuchte von oberflächennahen bodenschichten
EP1986127B1 (de) Erzeugung von Endmember-Spektren
Mozhaeva et al. RELATIONSHIPS BETWEEN VEGETATION INDICES AND RAINFALL AND PET AT DIFFERENT TIME-LAGS: A STUDY AT A MEDITERRANEAN TO ARID GRADIENT

Legal Events

Date Code Title Description
FG Patent granted

Effective date: 20220106