SNAP批量处理Sentinel-1A数据

已于 2022-08-04 17:31:55 修改
阅读量2.5k 收藏 28
点赞数 4
分类专栏： SAR 文章标签： sentinel snap sar
于 2022-06-29 16:21:36 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43262648/article/details/125524468
版权
本文详细介绍了如何使用SNAP的gpt工具，通过预先创建的流程图，批量处理Sentinel-1A的时序SAR数据。步骤包括对参考影像和时序数据进行极化预处理，生成C2矩阵并进行配准。每个步骤都提供了相应的XML配置文件内容，适用于大规模数据处理，提高工作效率。
摘要由CSDN通过智能技术生成

SNAP批量处理Sentinel-1A数据

通常情况，使用SNAP处理单景影像有两种方式，一是在SNAP中按照预处理步骤一个功能一个功能地处理数据，另一种是创建流程图直接运行得到最终的结果。

然而，当遇到时序数据时，上述做法显然不适用，这种情况下，我们可以采用gpt+流程图处理批量数据
.
具体操作流程如下（本实验数据为广州市2018年Sentinel-1A SAR影像，示例步骤为生成每一景影像的C2矩阵并配准）

SNAP批量处理Sentinel-1A数据

打开Windows PowerShell编译器

首先找到snap gpt的位置，选择【文件】->【Windows PowerShell】
在这里插入图片描述在命令行，输入cmd
在这里插入图片描述

运行gpt（总览）

本实验目的是生成一年影像的C2矩阵，并配准
我们可以将此过程分解为三个步骤：
1、对参考影像进行极化预处理，便于后面配准
2、对每一景影像进行极化预处理
3、对极化预处理后的C2矩阵进行配准
总的gpt代码如下：

## 极化预处理——获取C2矩阵
# 打开cmd 
cmd
# 对参考影像进行极化预处理
gpt F:\Data\S1_Guangzhou\C2-graph-main\1Reference.xml
# 对S1数据进行极化预处理
for /r "F:\Data\S1_Guangzhou\Sum" %X in (*.zip) do (gpt F:\Data\S1_Guangzhou\C2-graph-main\1Sum.xml -Pinput=%X  -Poutput="F:\Data\S1_Guangzhou\S2_features\Subset_%~nX_features.dim")
# 对时序S1数据进行配准
for /r "F:\Data\S1_Guangzhou\S2_features" %X in (*.dim) do (gpt F:\Data\S1_Guangzhou\C2-graph-main\1Coregister.xml -Pinput=%X)

！！！注意：
1、循环的语句中的路径需要根据需求修改
2、xml文件可根据需要在SNAP中创建
3、xml文件创建好后，将原有的输入路径改为$ input，输入路径改为$ output，结合for循环即可实现遍历

Step1：对参考影像进行极化预处理

SNAP中创建的流程图以及xml文件内容如下：
在这里插入图片描述

<graph id="Graph">
  <version>1.0</version>
  <node id="Read">
    <operator>Read</operator>
    <sources/>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <file>F:\Data\S1_Guangzhou\S1A_IW_SLC__1SDV_20200108T103351_20200108T103419_030708_03852A_7855.zip</file>
    </parameters>
  </node>
  <node id="Apply-Orbit-File">
    <operator>Apply-Orbit-File</operator>
    <sources>
      <sourceProduct refid="Read"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <orbitType>Sentinel Precise (Auto Download)</orbitType>
      <polyDegree>3</polyDegree>
      <continueOnFail>true</continueOnFail>
    </parameters>
  </node>
  <node id="TOPSAR-Split">
    <operator>TOPSAR-Split</operator>
    <sources>
      <sourceProduct refid="Apply-Orbit-File"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <subswath>IW2</subswath>
      <selectedPolarisations/>
      <firstBurstIndex>3</firstBurstIndex>
      <lastBurstIndex>6</lastBurstIndex>
      <wktAoi/>
    </parameters>
  </node>
  <node id="Calibration">
    <operator>Calibration</operator>
    <sources>
      <sourceProduct refid="TOPSAR-Split"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <sourceBands/>
      <auxFile>Latest Auxiliary File</auxFile>
      <externalAuxFile/>
      <outputImageInComplex>true</outputImageInComplex>
      <outputImageScaleInDb>false</outputImageScaleInDb>
      <createGammaBand>false</createGammaBand>
      <createBetaBand>false</createBetaBand>
      <selectedPolarisations/>
      <outputSigmaBand>true</outputSigmaBand>
      <outputGammaBand>false</outputGammaBand>
      <outputBetaBand>false</outputBetaBand>
    </parameters>
  </node>
  <node id="TOPSAR-Deburst">
    <operator>TOPSAR-Deburst</operator>
    <sources>
      <sourceProduct refid="Calibration"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <selectedPolarisations/>
    </parameters>
  </node>
  <node id="Subset">
    <operator>Subset</operator>
    <sources>
      <sourceProduct refid="TOPSAR-Deburst"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <sourceBands/>
      <region>0,0,0,0</region>
      <referenceBand/>
      <geoRegion>POLYGON ((113.17552185058594 23.3626651763916, 113.56804656982422 23.3626651763916, 113.56804656982422 22.93994140625, 113.17552185058594 22.93994140625, 113.17552185058594 23.3626651763916, 113.17552185058594 23.3626651763916))</geoRegion>
      <subSamplingX>1</subSamplingX>
      <subSamplingY>1</subSamplingY>
      <fullSwath>false</fullSwath>
      <tiePointGrids/>
      <copyMetadata>true</copyMetadata>
    </parameters>
  </node>
  <node id="Polarimetric-Matrices">
    <operator>Polarimetric-Matrices</operator>
    <sources>
      <sourceProduct refid="Subset"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <matrix>C2</matrix>
    </parameters>
  </node>
  <node id="Multilook">
    <operator>Multilook</operator>
    <sources>
      <sourceProduct refid="Polarimetric-Matrices"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <sourceBands/>
      <nRgLooks>4</nRgLooks>
      <nAzLooks>1</nAzLooks>
      <outputIntensity>false</outputIntensity>
      <grSquarePixel>true</grSquarePixel>
    </parameters>
  </node>
  <node id="Terrain-Correction">
    <operator>Terrain-Correction</operator>
    <sources>
      <sourceProduct refid="Polarimetric-Speckle-Filter"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <sourceBands/>
      <demName>SRTM 1Sec HGT</demName>
      <externalDEMFile/>
      <externalDEMNoDataValue>0.0</externalDEMNoDataValue>
      <externalDEMApplyEGM>true</externalDEMApplyEGM>
      <demResamplingMethod>BILINEAR_INTERPOLATION</demResamplingMethod>
      <imgResamplingMethod>BILINEAR_INTERPOLATION</imgResamplingMethod>
      <pixelSpacingInMeter>14.73</pixelSpacingInMeter>
      <pixelSpacingInDegree>1.3232184135080554E-4</pixelSpacingInDegree>
      <mapProjection>GEOGCS[&quot;WGS84(DD)&quot;, &#xd;
  DATUM[&quot;WGS84&quot;, &#xd;
    SPHEROID[&quot;WGS84&quot;, 6378137.0, 298.257223563]], &#xd;
  PRIMEM[&quot;Greenwich&quot;, 0.0], &#xd;
  UNIT[&quot;degree&quot;, 0.017453292519943295], &#xd;
  AXIS[&quot;Geodetic longitude&quot;, EAST], &#xd;
  AXIS[&quot;Geodetic latitude&quot;, NORTH]]</mapProjection>
      <alignToStandardGrid>false</alignToStandardGrid>
      <standardGridOriginX>0.0</standardGridOriginX>
      <standardGridOriginY>0.0</standardGridOriginY>
      <nodataValueAtSea>true</nodataValueAtSea>
      <saveDEM>false</saveDEM>
      <saveLatLon>false</saveLatLon>
      <saveIncidenceAngleFromEllipsoid>false</saveIncidenceAngleFromEllipsoid>
      <saveLocalIncidenceAngle>false</saveLocalIncidenceAngle>
      <saveProjectedLocalIncidenceAngle>false</saveProjectedLocalIncidenceAngle>
      <saveSelectedSourceBand>true</saveSelectedSourceBand>
      <saveLayoverShadowMask>false</saveLayoverShadowMask>
      <applyRadiometricNormalization>false</applyRadiometricNormalization>
      <saveSigmaNought>false</saveSigmaNought>
      <saveGammaNought>false</saveGammaNought>
      <saveBetaNought>false</saveBetaNought>
      <incidenceAngleForSigma0>Use projected local incidence angle from DEM</incidenceAngleForSigma0>
      <incidenceAngleForGamma0>Use projected local incidence angle from DEM</incidenceAngleForGamma0>
      <auxFile>Latest Auxiliary File</auxFile>
      <externalAuxFile/>
    </parameters>
  </node>
  <node id="Polarimetric-Speckle-Filter">
    <operator>Polarimetric-Speckle-Filter</operator>
    <sources>
      <sourceProduct refid="Multilook"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <filter>Refined Lee Filter</filter>
      <filterSize>5</filterSize>
      <numLooksStr>1</numLooksStr>
      <windowSize>7x7</windowSize>
      <targetWindowSizeStr>3x3</targetWindowSizeStr>
      <anSize>50</anSize>
      <sigmaStr>0.9</sigmaStr>
      <searchWindowSizeStr>15</searchWindowSizeStr>
      <patchSizeStr>5</patchSizeStr>
      <scaleSizeStr>1</scaleSizeStr>
    </parameters>
  </node>
  <node id="Write">
    <operator>Write</operator>
    <sources>
      <sourceProduct refid="Terrain-Correction"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <file>F:\Data\S1_Guangzhou\Reference.dim</file>
      <formatName>BEAM-DIMAP</formatName>
    </parameters>
  </node>
  <applicationData id="Presentation">
    <Description/>
    <node id="Read">
            <displayPosition x="14.0" y="52.0"/>
    </node>
    <node id="Apply-Orbit-File">
      <displayPosition x="91.0" y="53.0"/>
    </node>
    <node id="TOPSAR-Split">
      <displayPosition x="201.0" y="54.0"/>
    </node>
    <node id="Calibration">
      <displayPosition x="301.0" y="55.0"/>
    </node>
    <node id="TOPSAR-Deburst">
      <displayPosition x="409.0" y="56.0"/>
    </node>
    <node id="Subset">
      <displayPosition x="436.0" y="106.0"/>
    </node>
    <node id="Polarimetric-Matrices">
      <displayPosition x="403.0" y="150.0"/>
    </node>
    <node id="Multilook">
      <displayPosition x="437.0" y="195.0"/>
    </node>
    <node id="Terrain-Correction">
      <displayPosition x="568.0" y="249.0"/>
    </node>
    <node id="Polarimetric-Speckle-Filter">
      <displayPosition x="390.0" y="248.0"/>
    </node>
    <node id="Write">
      <displayPosition x="704.0" y="250.0"/>
    </node>
  </applicationData>
</graph>

Step2：对时序S1数据进行极化预处理

SNAP中创建的流程图以及xml文件内容如下：
在这里插入图片描述

在这里插入代码片<graph id="Graph">
  <version>1.0</version>
  <node id="Read">
    <operator>Read</operator>
    <sources/>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <file>$input</file>
    </parameters>
  </node>
  <node id="Apply-Orbit-File">
    <operator>Apply-Orbit-File</operator>
    <sources>
      <sourceProduct refid="Read"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <orbitType>Sentinel Precise (Auto Download)</orbitType>
      <polyDegree>3</polyDegree>
      <continueOnFail>true</continueOnFail>
    </parameters>
  </node>
  <node id="TOPSAR-Split">
    <operator>TOPSAR-Split</operator>
    <sources>
      <sourceProduct refid="Apply-Orbit-File"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <subswath>IW2</subswath>
      <selectedPolarisations/>
      <firstBurstIndex>3</firstBurstIndex>
      <lastBurstIndex>6</lastBurstIndex>
      <wktAoi/>
    </parameters>
  </node>
  <node id="Calibration">
    <operator>Calibration</operator>
    <sources>
      <sourceProduct refid="TOPSAR-Split"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <sourceBands/>
      <auxFile>Latest Auxiliary File</auxFile>
      <externalAuxFile/>
      <outputImageInComplex>true</outputImageInComplex>
      <outputImageScaleInDb>false</outputImageScaleInDb>
      <createGammaBand>false</createGammaBand>
      <createBetaBand>false</createBetaBand>
      <selectedPolarisations/>
      <outputSigmaBand>true</outputSigmaBand>
      <outputGammaBand>false</outputGammaBand>
      <outputBetaBand>false</outputBetaBand>
    </parameters>
  </node>
  <node id="TOPSAR-Deburst">
    <operator>TOPSAR-Deburst</operator>
    <sources>
      <sourceProduct refid="Calibration"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <selectedPolarisations/>
    </parameters>
  </node>
  <node id="Subset">
    <operator>Subset</operator>
    <sources>
      <sourceProduct refid="TOPSAR-Deburst"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <sourceBands/>
      <region>0,0,0,0</region>
      <referenceBand/>
      <geoRegion>POLYGON ((113.17552185058594 23.3626651763916, 113.56804656982422 23.3626651763916, 113.56804656982422 22.93994140625, 113.17552185058594 22.93994140625, 113.17552185058594 23.3626651763916, 113.17552185058594 23.3626651763916))</geoRegion>
      <subSamplingX>1</subSamplingX>
      <subSamplingY>1</subSamplingY>
      <fullSwath>false</fullSwath>
      <tiePointGridNames/>
      <copyMetadata>true</copyMetadata>
    </parameters>
  </node>
  <node id="Polarimetric-Matrices">
    <operator>Polarimetric-Matrices</operator>
    <sources>
      <sourceProduct refid="Subset"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <matrix>C2</matrix>
    </parameters>
  </node>
  <node id="Multilook">
    <operator>Multilook</operator>
    <sources>
      <sourceProduct refid="Polarimetric-Matrices"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <sourceBands/>
      <nRgLooks>4</nRgLooks>
      <nAzLooks>1</nAzLooks>
      <outputIntensity>false</outputIntensity>
      <grSquarePixel>true</grSquarePixel>
    </parameters>
  </node>
  <node id="Terrain-Correction">
    <operator>Terrain-Correction</operator>
    <sources>
      <sourceProduct refid="Polarimetric-Speckle-Filter"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <sourceBands/>
      <demName>SRTM 1Sec HGT</demName>
      <externalDEMFile/>
      <externalDEMNoDataValue>0.0</externalDEMNoDataValue>
      <externalDEMApplyEGM>true</externalDEMApplyEGM>
      <demResamplingMethod>BILINEAR_INTERPOLATION</demResamplingMethod>
      <imgResamplingMethod>BILINEAR_INTERPOLATION</imgResamplingMethod>
      <pixelSpacingInMeter>14.73</pixelSpacingInMeter>
      <pixelSpacingInDegree>1.3232184135080554E-4</pixelSpacingInDegree>
      <mapProjection>GEOGCS[&quot;WGS84(DD)&quot;, &#xd;
  DATUM[&quot;WGS84&quot;, &#xd;
    SPHEROID[&quot;WGS84&quot;, 6378137.0, 298.257223563]], &#xd;
  PRIMEM[&quot;Greenwich&quot;, 0.0], &#xd;
  UNIT[&quot;degree&quot;, 0.017453292519943295], &#xd;
  AXIS[&quot;Geodetic longitude&quot;, EAST], &#xd;
  AXIS[&quot;Geodetic latitude&quot;, NORTH]]</mapProjection>
      <alignToStandardGrid>false</alignToStandardGrid>
      <standardGridOriginX>0.0</standardGridOriginX>
      <standardGridOriginY>0.0</standardGridOriginY>
      <nodataValueAtSea>true</nodataValueAtSea>
      <saveDEM>false</saveDEM>
      <saveLatLon>false</saveLatLon>
      <saveIncidenceAngleFromEllipsoid>false</saveIncidenceAngleFromEllipsoid>
      <saveLocalIncidenceAngle>false</saveLocalIncidenceAngle>
      <saveProjectedLocalIncidenceAngle>false</saveProjectedLocalIncidenceAngle>
      <saveSelectedSourceBand>true</saveSelectedSourceBand>
      <outputComplex>false</outputComplex>
      <applyRadiometricNormalization>false</applyRadiometricNormalization>
      <saveSigmaNought>false</saveSigmaNought>
      <saveGammaNought>false</saveGammaNought>
      <saveBetaNought>false</saveBetaNought>
      <incidenceAngleForSigma0>Use projected local incidence angle from DEM</incidenceAngleForSigma0>
      <incidenceAngleForGamma0>Use projected local incidence angle from DEM</incidenceAngleForGamma0>
      <auxFile>Latest Auxiliary File</auxFile>
      <externalAuxFile/>
    </parameters>
  </node>
  <node id="Polarimetric-Speckle-Filter">
    <operator>Polarimetric-Speckle-Filter</operator>
    <sources>
      <sourceProduct refid="Multilook"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <filter>Refined Lee Filter</filter>
      <filterSize>5</filterSize>
      <numLooksStr>1</numLooksStr>
      <windowSize>7x7</windowSize>
      <targetWindowSizeStr>3x3</targetWindowSizeStr>
      <anSize>50</anSize>
      <sigmaStr>0.9</sigmaStr>
      <searchWindowSizeStr>15</searchWindowSizeStr>
      <patchSizeStr>5</patchSizeStr>
      <scaleSizeStr>1</scaleSizeStr>
    </parameters>
  </node>
  <node id="Write">
    <operator>Write</operator>
    <sources>
      <sourceProduct refid="Terrain-Correction"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <file>$output</file>
      <formatName>BEAM-DIMAP</formatName>
    </parameters>
  </node>
  <applicationData id="Presentation">
    <Description/>
    <node id="Read">
            <displayPosition x="14.0" y="52.0"/>
    </node>
    <node id="Apply-Orbit-File">
      <displayPosition x="91.0" y="53.0"/>
    </node>
    <node id="TOPSAR-Split">
      <displayPosition x="201.0" y="54.0"/>
    </node>
    <node id="Calibration">
      <displayPosition x="301.0" y="55.0"/>
    </node>
    <node id="TOPSAR-Deburst">
      <displayPosition x="409.0" y="56.0"/>
    </node>
    <node id="Subset">
      <displayPosition x="436.0" y="106.0"/>
    </node>
    <node id="Polarimetric-Matrices">
      <displayPosition x="403.0" y="150.0"/>
    </node>
    <node id="Multilook">
      <displayPosition x="437.0" y="195.0"/>
    </node>
    <node id="Terrain-Correction">
      <displayPosition x="568.0" y="249.0"/>
    </node>
    <node id="Polarimetric-Speckle-Filter">
      <displayPosition x="390.0" y="248.0"/>
    </node>
    <node id="Write">
      <displayPosition x="704.0" y="250.0"/>
    </node>
  </applicationData>
</graph>

Step3：对时序C2数据进行配准

SNAP中创建的流程图以及xml文件内容如下：
在这里插入图片描述

<graph id="Graph">
  <version>1.0</version>
  <node id="CreateStack">
    <operator>CreateStack</operator>
    <sources>
      <sourceProduct.2 refid="ProductSet-Reader"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <masterBands/>
      <sourceBands/>
      <resamplingType>NONE</resamplingType>
      <extent>Master</extent>
      <initialOffsetMethod>Orbit</initialOffsetMethod>
    </parameters>
  </node>
  <node id="Cross-Correlation">
    <operator>Cross-Correlation</operator>
    <sources>
      <sourceProduct refid="CreateStack"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <numGCPtoGenerate>200</numGCPtoGenerate>
      <coarseRegistrationWindowWidth>64</coarseRegistrationWindowWidth>
      <coarseRegistrationWindowHeight>64</coarseRegistrationWindowHeight>
      <rowInterpFactor>2</rowInterpFactor>
      <columnInterpFactor>2</columnInterpFactor>
      <maxIteration>2</maxIteration>
      <gcpTolerance>0.5</gcpTolerance>
      <applyFineRegistration>false</applyFineRegistration>
      <inSAROptimized>false</inSAROptimized>
      <fineRegistrationWindowWidth>64</fineRegistrationWindowWidth>
      <fineRegistrationWindowHeight>64</fineRegistrationWindowHeight>
      <fineRegistrationWindowAccAzimuth/>
      <fineRegistrationWindowAccRange/>
      <fineRegistrationOversampling/>
      <coherenceWindowSize>3</coherenceWindowSize>
      <coherenceThreshold>0.6</coherenceThreshold>
      <useSlidingWindow>false</useSlidingWindow>
      <computeOffset>false</computeOffset>
      <onlyGCPsOnLand>false</onlyGCPsOnLand>
    </parameters>
  </node>
  <node id="Warp">
    <operator>Warp</operator>
    <sources>
      <sourceProduct refid="Cross-Correlation"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <rmsThreshold>1.0</rmsThreshold>
      <warpPolynomialOrder>1</warpPolynomialOrder>
      <interpolationMethod>Bilinear interpolation</interpolationMethod>
      <demRefinement/>
      <demName/>
      <excludeMaster>false</excludeMaster>
      <openResidualsFile>false</openResidualsFile>
    </parameters>
  </node>
  <node id="Write">
    <operator>Write</operator>
    <sources>
      <sourceProduct refid="Warp"/>
    </sources>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <file>F:\Data\S1_Guangzhou\C2_features_Stack.dim</file>
      <formatName>BEAM-DIMAP</formatName>
    </parameters>
  </node>
  <node id="ProductSet-Reader">
    <operator>ProductSet-Reader</operator>
    <sources/>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <fileList>F:\Data\S1_Guangzhou\Reference.dim,$input</fileList>
    </parameters>
  </node>
  <applicationData id="Presentation">
    <Description/>
    <node id="CreateStack">
      <displayPosition x="154.0" y="101.0"/>
    </node>
    <node id="Cross-Correlation">
      <displayPosition x="248.0" y="100.0"/>
    </node>
    <node id="Warp">
      <displayPosition x="373.0" y="100.0"/>
    </node>
    <node id="Write">
      <displayPosition x="449.0" y="101.0"/>
    </node>
    <node id="ProductSet-Reader">
      <displayPosition x="24.0" y="99.0"/>
    </node>
  </applicationData>
</graph>

所有的文件可在C2_features_Stack文件夹中找到，后续处理可在matlab中进行

暮归山泉
关注
  • 4
    点赞
  • 踩
  • 28
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
snap如何打开哨兵二号数据_SNAP批量处理Sentinel-2数据
weixin_39637920的博客
11-23 2975
Sentinel-2数据以其高时空分辨率获得了越来越多的应用，然而由于Sentinel-2数据空间分辨率高，光谱分辨率高，数据量也较Landsat数据大幅上升。好在ESA SNAP软件批处理功能也很强大，可以支撑起大数据量的应用。下面就介绍一下如何使用SNAP批量处理Sentinel-2数据。Sentinel-2数据的下载和由L1C数据处理为L2A数据过程请见往期文章：Sentinel-2数据简介...
SNAP处理Sentinel-1 SAR影像流程
weixin_51610638的博客
03-19 3841
本文的处理流程适用于Sentinel-1数据关于数据和数据的区别和介绍，很容易找到，本文不再赘述。本文的处理步骤包含了辐射定标、噪声去除、轨道校正、地形校正（地理编码），这些步骤完成后SAR影像便可正常投入使用。
定制ARM 版本的ubuntu core (16.04)
热门推荐
Qing的专栏
09-25 1万+
1.  先从snap 包讲起 ubuntu core 基本上是诸多snap 包堆积起来的一个系统，正如传统的ubuntu 是debian 堆积起来的一样。但不同的是ubuntu core 也是一个snap 包， 甚至 kernel， uboot 也分别是一个snap 包。因此，就让我们先 了解下snap 包 是个甚么东东吧。 1.1 什么是snap ？ A snap : is a s...
08-SNAP的命令行处理工具gpt及其批处理（Sentinel-1和Sentinel-2为例）-附件资源
03-05
08-SNAP的命令行处理工具gpt及其批处理（Sentinel-1和Sentinel-2为例）-附件资源
08-SNAP的命令行处理工具gpt及其批处理（Sentinel-1和Sentinel-2为例）
lidahuilidahui的博客
07-02 1万+
08-SNAP的命令行处理工具gpt及其批处理前言gpt工具简介gpt实现单一操作单个数据集单一操作处理多个数据集批量单一操作处理gpt搭配流程图实现流程处理单个数据集流程图处理多个数据集批量流程图处理后语 前言 gpt工具简介 gpt实现单一操作 单个数据集单一操作处理 多个数据集批量单一操作处理 gpt搭配流程图实现流程处理 单个数据集流程图处理 多个数据集批量流程图处理 后语 ...
使用 SNAP 处理 SENTINEL-1
m0_60730723的博客
12-24 5953
使用 SNAP 处理 SENTINEL-1
SNAP处理Sentinel-1数据的详细指导
01-04
Sentinel-1是一颗双卫星系统（Sentinel-1A和Sentinel-1B），运行在太阳同步轨道上，提供全天候、全天时的C波段合成孔径雷达（SAR）数据。SAR技术能穿透云层和暗夜，使其在任何天气条件下都能进行观测。IW模式是其...
SNAP处理Sentinel-1 IW SLC数据
08-08
SNAP处理Sentinel-1 IW SLC数据 Sentinel-1是欧空局开发的雷达卫星系统，主要用于地表监测和环境监测。IW模式是Sentinel-1卫星的一种工作模式，SLC（Single Look Complex）是 Sentinel-1 IW模式下的数据产品之一。...
Sentinel3AltimetryL2_sentinel-3数据处理_
10-03
本文将详细介绍如何读取和处理Sentinel-3的测高数据，以Sentinel3AltimetryL2为例。 Sentinel-3 Altimetry Level-2数据代表了经过初步处理的高度测量结果，包含了原始雷达回波信号处理后的海面高度、冰雪覆盖区域的...
PolSARpro v5.1.3 处理Sentinel-1A SLC数据-附件资源
03-05
PolSARpro v5.1.3 处理Sentinel-1A SLC数据-附件资源
sentinel-1_SNAP数据处理教程.pdf
04-21
简单易懂，初学者可快速学会SNAP的数据预处理操作。从数据的读取到数据的预处理步骤详细有说明。对于哨兵数据的处理SNAP是一个很好的处理工具，大致内容介绍有Calibration（校准，标度） Speckle Filtering（几何配准，相干斑滤波） Terrain Correction（地形辐射校正） Mosaicking（拼接、镶嵌） Polarimetry（偏振测定） Interferometry（干涉量度分析法） Classification(分类)
ESA 哨兵数据处理软件 SNAP
11-25
ESA（欧空局）提供的哨兵系列数据处理软件，可对哨兵系列数据、landsat系列数据、spot系列数据等光学遥感数据和SAR数据进行诸如几何校正、辐射定标、滤波、裁剪、镶嵌等一系列处理。
哨兵一号数据Snap预处理指导书
12-03
我们对于欧洲空间局下载的哨兵一号的数据进行预处理的步骤写在了知道书中，希望可以帮助到同学。
遥感【SAR数据处理】【SNAP使用指导】【哨兵1号2号】如何预处理sentinel-1和2数据
06-10
本篇文章将重点介绍如何使用SNAP（Sentinels Application Platform）工具对Sentinel-1和Sentinel-2数据进行预处理。 **一、Sentinel-1 SAR数据处理** Sentinel-1数据处理主要包括辐射定标、几何校正、滤波和多视...
SNAP：如何批量预处理Sentinel2 L2A数据集并输出为TIFF文件？
最新发布
炒茄子的博客
03-03 1508
我目前就是希望下载哨兵2号数据，然后在SNAP中进行批量提取真彩色波段并输出为TIFF文件。
SNAP 处理Sentinel-1A数据（极化预处理，生成C2矩阵）
~Up~
04-07 5999
SNAP 处理Sentinel-1A数据（极化预处理）
C# 如何获取checkedListBox 的选中项
weixin_45002375的博客
09-30 6069
如图已经选择 2 5 8项 怎么获取其选中项 执行效果
Sentinel-1A 数据预处理
09-11
您好！对于Sentinel-1A数据的预处理，一般的步骤包括数据下载、校正、滤波和几何校正等。以下是一个常见的Sentinel-1A数据预处理流程： 1. 数据下载：从ESA的Sentinel数据存档中下载需要的Sentinel-1A数据。 2. 校正：进行辐射定标和表面反射率校正，以消除大气和地表特性的影响。这可以通过使用标准的Sentinel-1A SAR处理工具，如SNAP（Sentinel Application Platform）或其他开源软件来完成。 3. 滤波：在校正后，可以应用滤波技术来减少噪声和杂散信号，提高图像质量。常用的滤波方法包括Lee滤波器、Frost滤波器等。 4. 几何校正：进行几何校正以纠正因地球曲率、拓扑效应和传感器位置误差引起的图像畸变。这可以通过使用附带的DEM（数字高程模型）数据进行配准来实现。 5. 拼接和合成：如果需要，可以将多个Sentinel-1A图像拼接在一起，并生成复合图像，以获得更广阔的覆盖区域或更长的时间序列。 这里只是一个基本的概述，具体的预处理步骤可能会因应用需求而有所不同。希望以上信息对您有所帮助！如果您有任何其他问题，请随时提问。
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值