Error by reading Sentinel-1 bands #209
Answered
by
remi-braun
bcerripromethee
asked this question in
Q&A
|
Hey, I'm trying to read a Sentinel-1 product with this code : reader = Reader()
prod = reader.open(os.path.abspath('S1A_IW_GRDH_1SDV_20250331T051040_20250331T051105_058544_073E9F_F12E.SAFE.zip'), constellation= Constellation.S1, remove_tmp=True)
bands = [VV, VV_DSPK]
prod.load(bands=bands, pixel_size=30.)But i get this error message : 2025-04-11 11:57:49,628 - [DEBUG] - Loading bands ['VV', 'VV_DSPK']
-> [1214](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/geoseries.py:1214) return self.values.estimate_utm_crs(datum_name)
File /usr/local/lib/python3.11/dist-packages/geopandas/array.py:942, in GeometryArray.estimate_utm_crs(self, datum_name)
[937](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:937) else:
[938](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:938) x_center = np.mean([minx, maxx])
[940](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:940) utm_crs_list = query_utm_crs_info(
[941](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:941) datum_name=datum_name,
--> [942](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:942) area_of_interest=AreaOfInterest(
[943](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:943) west_lon_degree=x_center,
[944](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:944) south_lat_degree=y_center,
[945](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:945) east_lon_degree=x_center,
[946](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:946) north_lat_degree=y_center,
[947](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:947) ),
[948](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:948) )
[949](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:949) try:
[950](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/geopandas/array.py:950) return CRS.from_epsg(utm_crs_list[0].code)
File <string>:7, in __init__(self, west_lon_degree, south_lat_degree, east_lon_degree, north_lat_degree)
File /usr/local/lib/python3.11/dist-packages/pyproj/aoi.py:38, in AreaOfInterest.__post_init__(self)
[31](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/pyproj/aoi.py:31) def __post_init__(self):
[32](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/pyproj/aoi.py:32) if (
[33](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/pyproj/aoi.py:33) is_null(self.west_lon_degree)
[34](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/pyproj/aoi.py:34) or is_null(self.south_lat_degree)
[35](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/pyproj/aoi.py:35) or is_null(self.east_lon_degree)
[36](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/pyproj/aoi.py:36) or is_null(self.north_lat_degree)
[37](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/pyproj/aoi.py:37) ):
---> [38](https://vscode-remote+attached-002dcontainer-002b7b22636f6e7461696e65724e616d65223a222f736572766963652d707269736d65732d707269736d65732d31222c2273657474696e6773223a7b22636f6e74657874223a226465736b746f702d6c696e7578227d7d.vscode-resource.vscode-cdn.net/usr/local/lib/python3.11/dist-packages/pyproj/aoi.py:38) raise ValueError("NaN or None values are not allowed.")
ValueError: NaN or None values are not allowed.Do you have an idea ? PS : My eoreader version is : '0.22.1' |
Answered by
remi-braun
Apr 17, 2025
Replies: 2 comments 24 replies
|
Hello @bcerripromethee, I will test this tomorrow, but can I have access to your data? |
18 replies
|
Do you manage to open properly |
|
Actually, i'm not interested in reading map-overlay. I just want to preprocess sentinel-1 images. But no worries, i'll script it with |
|
The thing is that I am trying to get the extent of the S1 product when pre-processing the data. If it fails when reading |
|
After reflection, I think that it's related with gdal (as you were saying). Even if i can do things like gdalinfo, executing a Traceback (most recent call last):
File "/usr/lib/python3/dist-packages/osgeo/__init__.py", line 30, in swig_import_helper
return importlib.import_module(mname)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/usr/lib/python3.11/importlib/__init__.py", line 126, in import_module
return _bootstrap._gcd_import(name[level:], package, level)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "<frozen importlib._bootstrap>", line 1206, in _gcd_import
File "<frozen importlib._bootstrap>", line 1178, in _find_and_load
File "<frozen importlib._bootstrap>", line 1142, in _find_and_load_unlocked
ModuleNotFoundError: No module named 'osgeo._gdal'So it's not eoreader's fault ! |
|
Good luck debugging all that :'( |
|
Ok i finally did it ! Eoreader works ! I would like to use a personal graph : <graph id="Graph">
<version>1.0</version>
<node id="Read">
<operator>Read</operator>
<sources/>
<parameters class="com.bc.ceres.binding.dom.XppDomElement">
<useAdvancedOptions>true</useAdvancedOptions>
<file>${file}</file>
<geometryRegion>${geo_region}</geometryRegion>
<copyMetadata>true</copyMetadata>
<bandNames>Amplitude_VH,Intensity_VH,Amplitude_VV,Intensity_VV</bandNames>
<maskNames/>
</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>false</continueOnFail>
</parameters>
</node>
<node id="Remove-GRD-Border-Noise">
<operator>Remove-GRD-Border-Noise</operator>
<sources>
<sourceProduct refid="Apply-Orbit-File"/>
</sources>
<parameters class="com.bc.ceres.binding.dom.XppDomElement">
<selectedPolarisations/>
<borderLimit>500</borderLimit>
<trimThreshold>0.5</trimThreshold>
</parameters>
</node>
<node id="ThermalNoiseRemoval">
<operator>ThermalNoiseRemoval</operator>
<sources>
<sourceProduct refid="Remove-GRD-Border-Noise"/>
</sources>
<parameters class="com.bc.ceres.binding.dom.XppDomElement">
<selectedPolarisations/>
<removeThermalNoise>true</removeThermalNoise>
<outputNoise>false</outputNoise>
<reIntroduceThermalNoise>false</reIntroduceThermalNoise>
</parameters>
</node>
<node id="Calibration">
<operator>Calibration</operator>
<sources>
<sourceProduct refid="ThermalNoiseRemoval"/>
</sources>
<parameters class="com.bc.ceres.binding.dom.XppDomElement">
<sourceBands/>
<auxFile>Product Auxiliary File</auxFile>
<externalAuxFile/>
<outputImageInComplex>false</outputImageInComplex>
<outputImageScaleInDb>false</outputImageScaleInDb>
<createGammaBand>false</createGammaBand>
<createBetaBand>false</createBetaBand>
<selectedPolarisations/>
<outputSigmaBand>false</outputSigmaBand>
<outputGammaBand>false</outputGammaBand>
<outputBetaBand>true</outputBetaBand>
</parameters>
</node>
<node id="Terrain-Flattening">
<operator>Terrain-Flattening</operator>
<sources>
<sourceProduct refid="Calibration"/>
</sources>
<parameters class="com.bc.ceres.binding.dom.XppDomElement">
<sourceBands/>
<demName>${dem_name}</demName>
<demResamplingMethod>BILINEAR_INTERPOLATION</demResamplingMethod>
<externalDEMFile>${dem_path}</externalDEMFile>
<externalDEMNoDataValue>0.0</externalDEMNoDataValue>
<externalDEMApplyEGM>false</externalDEMApplyEGM>
<outputSimulatedImage>false</outputSimulatedImage>
<outputSigma0>false</outputSigma0>
<nodataValueAtSea>true</nodataValueAtSea>
<additionalOverlap>0.1</additionalOverlap>
<oversamplingMultiple>1.0</oversamplingMultiple>
</parameters>
</node>
<node id="Terrain-Correction">
<operator>Terrain-Correction</operator>
<sources>
<sourceProduct refid="LinearToFromdB"/>
</sources>
<parameters class="com.bc.ceres.binding.dom.XppDomElement">
<sourceBands/>
<demName>${dem_name}</demName>
<externalDEMFile>${dem_path}</externalDEMFile>
<externalDEMNoDataValue>0.0</externalDEMNoDataValue>
<externalDEMApplyEGM>true</externalDEMApplyEGM>
<demResamplingMethod>BILINEAR_INTERPOLATION</demResamplingMethod>
<imgResamplingMethod>BILINEAR_INTERPOLATION</imgResamplingMethod>
<pixelSpacingInMeter>${res_m}</pixelSpacingInMeter>
<pixelSpacingInDegree>${res_deg}</pixelSpacingInDegree>
<mapProjection>${crs}</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>
<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="LinearToFromdB">
<operator>LinearToFromdB</operator>
<sources>
<sourceProduct refid="Terrain-Flattening"/>
</sources>
<parameters class="com.bc.ceres.binding.dom.XppDomElement">
<sourceBands/>
</parameters>
</node>
<node id="Write">
<operator>Write</operator>
<sources>
<sourceProduct refid="Terrain-Correction"/>
</sources>
<parameters class="com.bc.ceres.binding.dom.XppDomElement">
<file>${out}</file>
<formatName>BEAM-DIMAP</formatName>
</parameters>
</node>
<applicationData id="Presentation">
<Description/>
<node id="Read">
<displayPosition x="16.0" y="16.0"/>
</node>
<node id="Apply-Orbit-File">
<displayPosition x="122.0" y="16.0"/>
</node>
<node id="Remove-GRD-Border-Noise">
<displayPosition x="91.0" y="78.0"/>
</node>
<node id="ThermalNoiseRemoval">
<displayPosition x="292.0" y="76.0"/>
</node>
<node id="Calibration">
<displayPosition x="326.0" y="132.0"/>
</node>
<node id="Terrain-Flattening">
<displayPosition x="440.0" y="132.0"/>
</node>
<node id="Terrain-Correction">
<displayPosition x="595.0" y="186.0"/>
</node>
<node id="LinearToFromdB">
<displayPosition x="452.0" y="185.0"/>
</node>
<node id="Write">
<displayPosition x="626.0" y="248.0"/>
</node>
</applicationData>
</graph>I just have to replace the environment variable like that ? : os.environ['EOREADER_PP_GRAPH'] = PATH_TO_MY_GRAPHand it's supposed to work ? |
6 replies
|
Ok ! I have the feeling that you subset is lonely in your graph, no ? <node id="Subset">
<operator>Subset</operator>
<sources>
<sourceProduct refid="Subset"/>
</sources>
<parameters class="com.bc.ceres.binding.dom.XppDomElement">
<sourceBands/>
<tiePointGrids/>
<region/>
<referenceBand/>
<region>${region}</region>
<geoRegion>${geo_region}</geoRegion>
<subSamplingX>1</subSamplingX>
<subSamplingY>1</subSamplingY>
<fullSwath>false</fullSwath>
<copyMetadata>true</copyMetadata>
</parameters>
</node> |
|
Hahaha well this is a copy pasting error 😂 |
|
Ok ! No worries 😂 |
|
I've just pushed the corrected versions ;) |
Answer selected by
remi-braun
|
Okay thanks !!! I'll try to adapt the graph to my needs, it's quite difficult but i'll do it, thanks again ! |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
I've just pushed the corrected versions ;)