MintPy icon indicating copy to clipboard operation
MintPy copied to clipboard

Geocode 4D dataset

Open bjmarfito opened this issue 1 year ago • 0 comments

Description of the desired feature

Request to geocode 4D dataset such as timeseriesDecorCov.h5 since geocode.py can only geocode 3D dataset.

Is your feature request related to a problem? Please describe

The error is generated by geocode.py when I try to geocode timeseriesDecorCov.h5.

reading timeseries in block (2032, 6116, 6465, 6605) from timeseriesDecorCov.h5 ...
nearest resampling with pyresample.kd_tree using 1 CPU cores in 1 segments ...
Traceback (most recent call last):
  File "/home/bryanjim001/tools/MintPy/src/mintpy/cli/geocode.py", line 269, in <module>
    main(sys.argv[1:])
  File "/home/bryanjim001/tools/MintPy/src/mintpy/cli/geocode.py", line 264, in main
    run_geocode(inps)
  File "/home/bryanjim001/tools/MintPy/src/mintpy/geocode.py", line 121, in run_geocode
    data = res_obj.run_resample(src_data=data, box_ind=i)
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/home/bryanjim001/tools/MintPy/src/mintpy/objects/resample.py", line 214, in run_resample
    dest_data = self.run_pyresample(
                ^^^^^^^^^^^^^^^^^^^^
  File "/home/bryanjim001/tools/MintPy/src/mintpy/objects/resample.py", line 724, in run_pyresample
    dest_data = pr.kd_tree.resample_nearest(
                ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/home/bryanjim001/.conda/envs/isce2expe/lib/python3.11/site-packages/pyresample/kd_tree.py", line 107, in resample_nearest
    return _resample(source_geo_def, data, target_geo_def, 'nn',
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/home/bryanjim001/.conda/envs/isce2expe/lib/python3.11/site-packages/pyresample/kd_tree.py", line 270, in _resample
    return get_sample_from_neighbour_info(resample_type,
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/home/bryanjim001/.conda/envs/isce2expe/lib/python3.11/site-packages/pyresample/kd_tree.py", line 613, in get_sample_from_neighbour_info
    raise ValueError('Mismatch between geometry and dataset')
ValueError: Mismatch between geometry and dataset

Describe the solution you'd like

Describe alternatives you have considered

Additional context

Are you willing to help implement and maintain this feature?

  • [ ] Yes
  • [ ] No

bjmarfito avatar Nov 22 '24 05:11 bjmarfito