Ingesting wrapped interferogram from topsStack causes an error

[bjmarfito]

Description of the problem

I loaded wrapped interferograms generated from topsStack using smallbaselineApp.py and it causes a loading error

Full script that generated the error

smallbaselineApp.py --dostep load_data

Full error message

create dataset /wrapPhase        of <class 'numpy.float32'>   in size of (6, 13374, 3607) with compression = None
[=======>                17%                       ] 20250403_20250415    0s /     0sTraceback (most recent call last):
  File "/home/bryanjim001/tools/MintPy/src/mintpy/cli/smallbaselineApp.py", line 214, in <module>
    main(sys.argv[1:])
  File "/home/bryanjim001/tools/MintPy/src/mintpy/cli/smallbaselineApp.py", line 209, in main
    run_smallbaselineApp(inps)
  File "/home/bryanjim001/tools/MintPy/src/mintpy/smallbaselineApp.py", line 1155, in run_smallbaselineApp
    app.run(steps=inps.runSteps)
  File "/home/bryanjim001/tools/MintPy/src/mintpy/smallbaselineApp.py", line 908, in run
    self.run_load_data(sname)
  File "/home/bryanjim001/tools/MintPy/src/mintpy/smallbaselineApp.py", line 176, in run_load_data
    mintpy.cli.load_data.main(iargs)
  File "/home/bryanjim001/tools/MintPy/src/mintpy/cli/load_data.py", line 121, in main
    load_data(inps)
  File "/home/bryanjim001/tools/MintPy/src/mintpy/load_data.py", line 881, in load_data
    stack_obj.write2hdf5(
  File "/home/bryanjim001/tools/MintPy/src/mintpy/objects/stackDict.py", line 215, in write2hdf5
    ds[i, :, :] = data
    ~~^^^^^^^^^
  File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper
  File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper
  File "/home/bryanjim001/.conda/envs/isce2expe/lib/python3.11/site-packages/h5py/_hl/dataset.py", line 1043, in __setitem__
    self.id.write(mspace, fspace, val, mtype, dxpl=self._dxpl)
  File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper
  File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper
  File "h5py/h5d.pyx", line 281, in h5py.h5d.DatasetID.write
  File "h5py/_proxy.pyx", line 106, in h5py._proxy.dset_rw
  File "h5py/_proxy.pyx", line 187, in h5py._proxy.make_reduced_type
ValueError: Cannot return member number (operation not supported for type class)

System information

• Operating system: Linux
• Python environment: Python 3.11
• MintPy version: MintPy version 1.6.1.post40, date 2025-06-10
• InSAR processor/product: isce2 v2.6.4
• Your custom / default template file (if the bug is related to a specific dataset): smallbaselineApp.txt

[yunjunz]

Hi @bjmarfito, I have not tested the wrapped phase loading for a long time, since we do not need it: we could re-wrap the unwrapped phase to get the wrapped phase.

As for your error, it looks like related to the 20250403_20250415 interferogram, could you try to plot its .int file with view.py to ensure everything with the file itself is normal?

[bjmarfito]

Hi @yunjunz, I see so the wrapped phase is not needed anymore in MintPy. Please correct me if I'm wrong.

Here is the wrapped interferogram (.int) for your reference, please.

[ditafaith]

Indeed, no need wrapped since SBAS method use the unwrapped one.

[yunjunz]

Right. I will remove that option in the future to avoid potential confusion.

[yunjunz]

And your .int file looks normal to me.

[bjmarfito]

Thank you @yunjunz and @ditafaith for your answers. Closing this discussion.